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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-3478.6508	6989.3016	0.2398	0.2398
GTRI	20	-3474.7404	6989.4808	0.2192	0.4590
HKYG	16	-3478.8010	6989.6020	0.2064	0.6654
GTRG	20	-3475.1121	6990.2243	0.1512	0.8166
HKYIG	17	-3478.5443	6991.0885	0.0981	0.9147
GTRIG	21	-3474.6845	6991.3691	0.0853	1.0000
HKY	15	-3511.6798	7053.3597	0.0000	1.0000
GTR	19	-3509.2337	7056.4673	0.0000	1.0000
K2PG	13	-3568.9625	7163.9251	0.0000	1.0000
K2PIG	14	-3569.0642	7166.1285	0.0000	1.0000
K2PI	13	-3571.9369	7169.8739	0.0000	1.0000
SYMIG	18	-3568.3086	7172.6172	0.0000	1.0000
SYMG	17	-3570.4312	7174.8625	0.0000	1.0000
SYMI	17	-3571.3880	7176.7761	0.0000	1.0000
F81G	15	-3578.1947	7186.3894	0.0000	1.0000
F81I	15	-3578.9112	7187.8223	0.0000	1.0000
F81IG	16	-3578.2289	7188.4577	0.0000	1.0000
K2P	12	-3607.8589	7239.7178	0.0000	1.0000
F81	14	-3606.8686	7241.7371	0.0000	1.0000
SYM	16	-3607.3582	7246.7163	0.0000	1.0000
JC69G	12	-3668.2726	7360.5453	0.0000	1.0000
JC69I	12	-3669.1569	7362.3138	0.0000	1.0000
JC69IG	13	-3668.3064	7362.6128	0.0000	1.0000
JC69	11	-3699.0735	7420.1471	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3478.6508	7070.3650	0.5282	0.5282
HKYG	16	-3478.8010	7070.6655	0.4546	0.9828
HKYIG	17	-3478.5443	7077.2185	0.0172	1.0000
GTRI	20	-3474.7404	7090.8102	0.0000	1.0000
GTRG	20	-3475.1121	7091.5536	0.0000	1.0000
GTRIG	21	-3474.6845	7097.7649	0.0000	1.0000
HKY	15	-3511.6798	7129.3567	0.0000	1.0000
GTR	19	-3509.2337	7152.7302	0.0000	1.0000
K2PG	13	-3568.9625	7229.7891	0.0000	1.0000
K2PI	13	-3571.9369	7235.7380	0.0000	1.0000
K2PIG	14	-3569.0642	7237.0590	0.0000	1.0000
SYMG	17	-3570.4312	7260.9924	0.0000	1.0000
F81G	15	-3578.1947	7262.3864	0.0000	1.0000
SYMI	17	-3571.3880	7262.9060	0.0000	1.0000
SYMIG	18	-3568.3086	7263.8136	0.0000	1.0000
F81I	15	-3578.9112	7263.8193	0.0000	1.0000
F81IG	16	-3578.2289	7269.5212	0.0000	1.0000
K2P	12	-3607.8589	7300.5154	0.0000	1.0000
F81	14	-3606.8686	7312.6676	0.0000	1.0000
SYM	16	-3607.3582	7327.7798	0.0000	1.0000
JC69G	12	-3668.2726	7421.3429	0.0000	1.0000
JC69I	12	-3669.1569	7423.1114	0.0000	1.0000
JC69IG	13	-3668.3064	7428.4769	0.0000	1.0000
JC69	11	-3699.0735	7475.8782	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr1_39151.phylip = 1 - 1172;
 Partition Dummy = 1:683-tomodel-phylip/chr1_39151.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr1_39151.phylip = 1 - 1172;
 Partition Dummy = 1:683-tomodel-phylip/chr1_39151.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr1_39151.phylip = 1 - 1172;
 Partition Dummy = 1:683-tomodel-phylip/chr1_39151.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

