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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-4025.3931	8090.7863	0.6647	0.6647
GTRIG	21	-4025.4081	8092.8161	0.2409	0.9056
GTRG	20	-4027.3470	8094.6940	0.0942	0.9998
HKYI	16	-4037.8698	8107.7396	0.0001	0.9999
HKYIG	17	-4037.8959	8109.7918	0.0000	1.0000
HKYG	16	-4039.8391	8111.6783	0.0000	1.0000
GTR	19	-4056.1028	8150.2055	0.0000	1.0000
HKY	15	-4065.4830	8160.9660	0.0000	1.0000
SYMI	17	-4073.9949	8181.9898	0.0000	1.0000
K2PI	13	-4078.0221	8182.0441	0.0000	1.0000
SYMIG	18	-4074.0172	8184.0345	0.0000	1.0000
K2PIG	14	-4078.0440	8184.0880	0.0000	1.0000
K2PG	13	-4079.6640	8185.3280	0.0000	1.0000
SYMG	17	-4075.6871	8185.3741	0.0000	1.0000
K2P	12	-4103.5891	8231.1781	0.0000	1.0000
SYM	16	-4099.9723	8231.9446	0.0000	1.0000
F81I	15	-4135.2621	8300.5241	0.0000	1.0000
F81IG	16	-4135.3168	8302.6336	0.0000	1.0000
F81G	15	-4137.9380	8305.8759	0.0000	1.0000
F81	14	-4158.3621	8344.7242	0.0000	1.0000
JC69I	12	-4164.8636	8353.7272	0.0000	1.0000
JC69IG	13	-4164.9138	8355.8276	0.0000	1.0000
JC69G	12	-4167.3391	8358.6783	0.0000	1.0000
JC69	11	-4187.4263	8396.8526	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-4037.8698	8188.4302	0.7130	0.7130
GTRI	20	-4025.3931	8191.6495	0.1426	0.8556
HKYG	16	-4039.8391	8192.3688	0.0995	0.9551
HKYIG	17	-4037.8959	8195.5255	0.0205	0.9756
GTRG	20	-4027.3470	8195.5572	0.0202	0.9958
GTRIG	21	-4025.4081	8198.7225	0.0042	1.0000
HKY	15	-4065.4830	8236.6134	0.0000	1.0000
GTR	19	-4056.1028	8246.0256	0.0000	1.0000
K2PI	13	-4078.0221	8247.6052	0.0000	1.0000
K2PG	13	-4079.6640	8250.8890	0.0000	1.0000
K2PIG	14	-4078.0440	8254.6922	0.0000	1.0000
SYMI	17	-4073.9949	8267.7235	0.0000	1.0000
SYMG	17	-4075.6871	8271.1078	0.0000	1.0000
SYMIG	18	-4074.0172	8274.8114	0.0000	1.0000
K2P	12	-4103.5891	8291.6961	0.0000	1.0000
SYM	16	-4099.9723	8312.6352	0.0000	1.0000
F81I	15	-4135.2621	8376.1715	0.0000	1.0000
F81G	15	-4137.9380	8381.5233	0.0000	1.0000
F81IG	16	-4135.3168	8383.3242	0.0000	1.0000
JC69I	12	-4164.8636	8414.2452	0.0000	1.0000
F81	14	-4158.3621	8415.3285	0.0000	1.0000
JC69G	12	-4167.3391	8419.1962	0.0000	1.0000
JC69IG	13	-4164.9138	8421.3887	0.0000	1.0000
JC69	11	-4187.4263	8452.3274	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_6216.phylip = 1 - 1145;
 Partition Dummy = 1:683-tomodel-phylip/chr2_6216.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_6216.phylip = 1 - 1145;
 Partition Dummy = 1:683-tomodel-phylip/chr2_6216.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,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_6216.phylip = 1 - 1145;
 Partition Dummy = 1:683-tomodel-phylip/chr2_6216.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

