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

Input data from file "683-tomodel-phylip/chr8_7425.phylip" (7 taxa, 1197 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	-3110.4693	6260.9387	0.2974	0.2974
HKYI	16	-3114.7279	6261.4559	0.2296	0.5271
GTRG	20	-3111.0395	6262.0789	0.1682	0.6952
HKYG	16	-3115.4322	6262.8643	0.1136	0.8088
GTRIG	21	-3110.4810	6262.9620	0.1081	0.9169
HKYIG	17	-3114.7447	6263.4895	0.0831	1.0000
GTR	19	-3129.4602	6296.9203	0.0000	1.0000
HKY	15	-3134.2473	6298.4946	0.0000	1.0000
K2PI	13	-3141.2547	6308.5094	0.0000	1.0000
K2PG	13	-3141.8741	6309.7483	0.0000	1.0000
SYMI	17	-3138.8470	6311.6940	0.0000	1.0000
K2PIG	14	-3141.9932	6311.9864	0.0000	1.0000
SYMG	17	-3139.4342	6312.8683	0.0000	1.0000
SYMIG	18	-3139.5588	6315.1176	0.0000	1.0000
K2P	12	-3163.1530	6350.3061	0.0000	1.0000
SYM	16	-3160.6985	6353.3969	0.0000	1.0000
F81I	15	-3231.2733	6492.5466	0.0000	1.0000
F81G	15	-3232.2136	6494.4271	0.0000	1.0000
F81IG	16	-3231.2990	6494.5981	0.0000	1.0000
F81	14	-3249.4040	6526.8081	0.0000	1.0000
JC69I	12	-3254.7213	6533.4425	0.0000	1.0000
JC69G	12	-3255.6059	6535.2119	0.0000	1.0000
JC69IG	13	-3255.6967	6537.3934	0.0000	1.0000
JC69	11	-3273.8527	6569.7055	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3114.7279	6342.8570	0.6566	0.6566
HKYG	16	-3115.4322	6344.2655	0.3247	0.9813
HKYIG	17	-3114.7447	6349.9782	0.0187	0.9999
GTRI	20	-3110.4693	6362.6902	0.0000	1.0000
GTRG	20	-3111.0395	6363.8304	0.0000	1.0000
GTRIG	21	-3110.4810	6369.8010	0.0000	1.0000
K2PI	13	-3141.2547	6374.6478	0.0000	1.0000
HKY	15	-3134.2473	6374.8082	0.0000	1.0000
K2PG	13	-3141.8741	6375.8867	0.0000	1.0000
K2PIG	14	-3141.9932	6383.2125	0.0000	1.0000
GTR	19	-3129.4602	6393.5842	0.0000	1.0000
SYMI	17	-3138.8470	6398.1827	0.0000	1.0000
SYMG	17	-3139.4342	6399.3571	0.0000	1.0000
SYMIG	18	-3139.5588	6406.6939	0.0000	1.0000
K2P	12	-3163.1530	6411.3570	0.0000	1.0000
SYM	16	-3160.6985	6434.7981	0.0000	1.0000
F81I	15	-3231.2733	6568.8602	0.0000	1.0000
F81G	15	-3232.2136	6570.7407	0.0000	1.0000
F81IG	16	-3231.2990	6575.9993	0.0000	1.0000
JC69I	12	-3254.7213	6594.4934	0.0000	1.0000
JC69G	12	-3255.6059	6596.2628	0.0000	1.0000
F81	14	-3249.4040	6598.0341	0.0000	1.0000
JC69IG	13	-3255.6967	6603.5319	0.0000	1.0000
JC69	11	-3273.8527	6625.6688	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/chr8_7425.phylip = 1 - 1197;
 Partition Dummy = 1:683-tomodel-phylip/chr8_7425.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/chr8_7425.phylip = 1 - 1197;
 Partition Dummy = 1:683-tomodel-phylip/chr8_7425.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/chr8_7425.phylip = 1 - 1197;
 Partition Dummy = 1:683-tomodel-phylip/chr8_7425.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

