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

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

Minimum AIC  model: GTRG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3477.6945	6995.3890	0.3552	0.3552
HKYG	16	-3481.7371	6995.4742	0.3403	0.6955
GTRIG	21	-3477.6848	6997.3696	0.1319	0.8274
HKYIG	17	-3481.7272	6997.4543	0.1265	0.9539
HKYI	16	-3484.4240	7000.8479	0.0232	0.9770
GTRI	20	-3480.4332	7000.8664	0.0230	1.0000
HKY	15	-3507.2958	7044.5917	0.0000	1.0000
GTR	19	-3503.6334	7045.2669	0.0000	1.0000
K2PG	13	-3541.7052	7109.4104	0.0000	1.0000
K2PIG	14	-3541.7048	7111.4096	0.0000	1.0000
K2PI	13	-3543.6613	7113.3225	0.0000	1.0000
SYMG	17	-3539.8590	7113.7180	0.0000	1.0000
SYMIG	18	-3539.8500	7115.7000	0.0000	1.0000
SYMI	17	-3541.7033	7117.4066	0.0000	1.0000
F81G	15	-3551.4337	7132.8673	0.0000	1.0000
F81IG	16	-3551.4338	7134.8675	0.0000	1.0000
F81I	15	-3553.8339	7137.6678	0.0000	1.0000
K2P	12	-3567.6162	7159.2323	0.0000	1.0000
SYM	16	-3566.3459	7164.6919	0.0000	1.0000
F81	14	-3575.1853	7178.3707	0.0000	1.0000
JC69G	12	-3601.8083	7227.6166	0.0000	1.0000
JC69IG	13	-3601.8034	7229.6069	0.0000	1.0000
JC69I	12	-3604.0591	7232.1182	0.0000	1.0000
JC69	11	-3625.9199	7273.8397	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3481.7371	7075.6824	0.9104	0.9104
HKYI	16	-3484.4240	7081.0562	0.0620	0.9724
HKYIG	17	-3481.7272	7082.6756	0.0276	1.0000
GTRG	20	-3477.6945	7095.6493	0.0000	1.0000
GTRI	20	-3480.4332	7101.1267	0.0000	1.0000
GTRIG	21	-3477.6848	7102.6429	0.0000	1.0000
HKY	15	-3507.2958	7119.7869	0.0000	1.0000
GTR	19	-3503.6334	7140.5142	0.0000	1.0000
K2PG	13	-3541.7052	7174.5796	0.0000	1.0000
K2PI	13	-3543.6613	7178.4917	0.0000	1.0000
K2PIG	14	-3541.7048	7181.5918	0.0000	1.0000
SYMG	17	-3539.8590	7198.9392	0.0000	1.0000
SYMI	17	-3541.7033	7202.6278	0.0000	1.0000
SYMIG	18	-3539.8500	7205.9343	0.0000	1.0000
F81G	15	-3551.4337	7208.0626	0.0000	1.0000
F81I	15	-3553.8339	7212.8630	0.0000	1.0000
F81IG	16	-3551.4338	7215.0758	0.0000	1.0000
K2P	12	-3567.6162	7219.3885	0.0000	1.0000
SYM	16	-3566.3459	7244.9002	0.0000	1.0000
F81	14	-3575.1853	7248.5529	0.0000	1.0000
JC69G	12	-3601.8083	7287.7728	0.0000	1.0000
JC69I	12	-3604.0591	7292.2744	0.0000	1.0000
JC69IG	13	-3601.8034	7294.7761	0.0000	1.0000
JC69	11	-3625.9199	7328.9829	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_27321.phylip = 1 - 1111;
 Partition Dummy = 1:683-tomodel-phylip/chr2_27321.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_27321.phylip = 1 - 1111;
 Partition Dummy = 1:683-tomodel-phylip/chr2_27321.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 (HKYG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_27321.phylip = 1 - 1111;
 Partition Dummy = 1:683-tomodel-phylip/chr2_27321.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;


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

