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

Input data from file "683-tomodel-phylip/chr2_7180.phylip" (7 taxa, 975 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	-3070.8020	6173.6039	0.3433	0.3433
GTRI	20	-3067.2931	6174.5861	0.2101	0.5533
HKYG	16	-3071.6092	6175.2183	0.1531	0.7065
HKYIG	17	-3070.8172	6175.6344	0.1244	0.8309
GTRG	20	-3068.1089	6176.2178	0.0929	0.9238
GTRIG	21	-3067.3065	6176.6131	0.0762	1.0000
HKY	15	-3085.7701	6201.5401	0.0000	1.0000
GTR	19	-3083.1305	6204.2611	0.0000	1.0000
K2PI	13	-3116.3324	6258.6647	0.0000	1.0000
K2PG	13	-3116.9836	6259.9671	0.0000	1.0000
K2PIG	14	-3116.3449	6260.6898	0.0000	1.0000
SYMI	17	-3115.6815	6265.3630	0.0000	1.0000
SYMG	17	-3116.2770	6266.5540	0.0000	1.0000
SYMIG	18	-3115.6926	6267.3852	0.0000	1.0000
F81I	15	-3123.7389	6277.4777	0.0000	1.0000
F81G	15	-3123.7560	6277.5119	0.0000	1.0000
F81IG	16	-3123.6739	6279.3478	0.0000	1.0000
K2P	12	-3130.7205	6285.4409	0.0000	1.0000
SYM	16	-3129.7362	6291.4723	0.0000	1.0000
F81	14	-3136.0360	6300.0721	0.0000	1.0000
JC69G	12	-3161.2990	6346.5980	0.0000	1.0000
JC69I	12	-3161.3916	6346.7831	0.0000	1.0000
JC69IG	13	-3161.2845	6348.5691	0.0000	1.0000
JC69	11	-3173.5567	6369.1133	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3070.8020	6251.7229	0.6767	0.6767
HKYG	16	-3071.6092	6253.3373	0.3019	0.9786
HKYIG	17	-3070.8172	6258.6359	0.0213	1.0000
GTRI	20	-3067.2931	6272.2349	0.0000	1.0000
GTRG	20	-3068.1089	6273.8666	0.0000	1.0000
HKY	15	-3085.7701	6274.7767	0.0000	1.0000
GTRIG	21	-3067.3065	6279.1443	0.0000	1.0000
GTR	19	-3083.1305	6297.0274	0.0000	1.0000
K2PI	13	-3116.3324	6322.1364	0.0000	1.0000
K2PG	13	-3116.9836	6323.4388	0.0000	1.0000
K2PIG	14	-3116.3449	6329.0439	0.0000	1.0000
K2P	12	-3130.7205	6344.0302	0.0000	1.0000
SYMI	17	-3115.6815	6348.3644	0.0000	1.0000
SYMG	17	-3116.2770	6349.5555	0.0000	1.0000
F81I	15	-3123.7389	6350.7143	0.0000	1.0000
F81G	15	-3123.7560	6350.7485	0.0000	1.0000
SYMIG	18	-3115.6926	6355.2691	0.0000	1.0000
F81IG	16	-3123.6739	6357.4668	0.0000	1.0000
F81	14	-3136.0360	6368.4262	0.0000	1.0000
SYM	16	-3129.7362	6369.5913	0.0000	1.0000
JC69G	12	-3161.2990	6405.1873	0.0000	1.0000
JC69I	12	-3161.3916	6405.3724	0.0000	1.0000
JC69IG	13	-3161.2845	6412.0408	0.0000	1.0000
JC69	11	-3173.5567	6422.8202	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/chr2_7180.phylip = 1 - 975;
 Partition Dummy = 1:683-tomodel-phylip/chr2_7180.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/chr2_7180.phylip = 1 - 975;
 Partition Dummy = 1:683-tomodel-phylip/chr2_7180.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/chr2_7180.phylip = 1 - 975;
 Partition Dummy = 1:683-tomodel-phylip/chr2_7180.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

