
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:11 04/04/2011

Input data from file "683-tomodel-phylip/chr1_32597.phylip" (7 taxa, 1145 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKY


(Output sorted by AIC. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2473.8304	4979.6609	0.2650	0.2650
HKYI	16	-2473.8501	4979.7002	0.2598	0.5248
GTRG	20	-2470.5524	4981.1048	0.1287	0.6535
GTRI	20	-2470.5664	4981.1328	0.1269	0.7805
HKYIG	17	-2473.8313	4981.6625	0.0974	0.8779
HKY	15	-2476.4595	4982.9190	0.0520	0.9298
GTRIG	21	-2470.5524	4983.1048	0.0474	0.9772
GTR	19	-2473.2825	4984.5649	0.0228	1.0000
K2P	12	-2525.8856	5075.7712	0.0000	1.0000
K2PG	13	-2524.9796	5075.9592	0.0000	1.0000
K2PI	13	-2524.9876	5075.9751	0.0000	1.0000
K2PIG	14	-2524.9797	5077.9593	0.0000	1.0000
SYM	16	-2524.2103	5080.4205	0.0000	1.0000
SYMG	17	-2523.4560	5080.9120	0.0000	1.0000
SYMI	17	-2523.4667	5080.9333	0.0000	1.0000
F81I	15	-2525.8748	5081.7495	0.0000	1.0000
F81G	15	-2526.0710	5082.1419	0.0000	1.0000
SYMIG	18	-2523.4554	5082.9108	0.0000	1.0000
F81IG	16	-2525.8810	5083.7619	0.0000	1.0000
F81	14	-2527.9529	5083.9058	0.0000	1.0000
JC69I	12	-2572.4873	5168.9747	0.0000	1.0000
JC69G	12	-2572.6765	5169.3531	0.0000	1.0000
JC69	11	-2574.3911	5170.7821	0.0000	1.0000
JC69IG	13	-2572.4939	5170.9877	0.0000	1.0000



(Output sorted by BIC. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		BIC		wBIC	cwBIC
HKY	15	-2476.4595	5058.5664	0.5484	0.5484
HKYG	16	-2473.8304	5060.3514	0.2246	0.7731
HKYI	16	-2473.8501	5060.3907	0.2203	0.9933
HKYIG	17	-2473.8313	5067.3963	0.0066	1.0000
GTR	19	-2473.2825	5080.3849	0.0000	1.0000
GTRG	20	-2470.5524	5081.9680	0.0000	1.0000
GTRI	20	-2470.5664	5081.9960	0.0000	1.0000
GTRIG	21	-2470.5524	5089.0111	0.0000	1.0000
K2P	12	-2525.8856	5136.2891	0.0000	1.0000
K2PG	13	-2524.9796	5141.5202	0.0000	1.0000
K2PI	13	-2524.9876	5141.5362	0.0000	1.0000
K2PIG	14	-2524.9797	5148.5636	0.0000	1.0000
F81	14	-2527.9529	5154.5100	0.0000	1.0000
F81I	15	-2525.8748	5157.3969	0.0000	1.0000
F81G	15	-2526.0710	5157.7893	0.0000	1.0000
SYM	16	-2524.2103	5161.1111	0.0000	1.0000
F81IG	16	-2525.8810	5164.4525	0.0000	1.0000
SYMG	17	-2523.4560	5166.6457	0.0000	1.0000
SYMI	17	-2523.4667	5166.6671	0.0000	1.0000
SYMIG	18	-2523.4554	5173.6877	0.0000	1.0000
JC69	11	-2574.3911	5226.2569	0.0000	1.0000
JC69I	12	-2572.4873	5229.4926	0.0000	1.0000
JC69G	12	-2572.6765	5229.8710	0.0000	1.0000
JC69IG	13	-2572.4939	5236.5488	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr1_32597.phylip = 1 - 1145;
 Partition Dummy = 1:683-tomodel-phylip/chr1_32597.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best AICc model (HKYG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr1_32597.phylip = 1 - 1145;
 Partition Dummy = 1:683-tomodel-phylip/chr1_32597.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best BIC model (HKY)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr1_32597.phylip = 1 - 1145;
 Partition Dummy = 1:683-tomodel-phylip/chr1_32597.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


-------------------------------------------------------------
End of Output

