
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:46 04/04/2011

Input data from file "591-tomodel-phylip/chr5_1674.phylip" (7 taxa, 1156 characters)

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


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-3509.3630	7050.7260	0.2594	0.2594
HKYI	16	-3509.4509	7050.9017	0.2375	0.4969
GTRG	20	-3505.6936	7051.3872	0.1863	0.6832
GTRI	20	-3505.8052	7051.6104	0.1667	0.8499
HKYIG	17	-3509.4499	7052.8998	0.0875	0.9374
GTRIG	21	-3505.7840	7053.5680	0.0626	1.0000
K2PI	13	-3542.1986	7110.3972	0.0000	1.0000
K2PG	13	-3542.2909	7110.5817	0.0000	1.0000
K2PIG	14	-3542.3825	7112.7650	0.0000	1.0000
SYMI	17	-3539.8193	7113.6385	0.0000	1.0000
SYMG	17	-3540.0492	7114.0984	0.0000	1.0000
SYMIG	18	-3540.1421	7116.2842	0.0000	1.0000
HKY	15	-3544.9056	7119.8113	0.0000	1.0000
GTR	19	-3541.8262	7121.6523	0.0000	1.0000
K2P	12	-3577.0753	7178.1506	0.0000	1.0000
SYM	16	-3574.9061	7181.8122	0.0000	1.0000
F81I	15	-3587.6367	7205.2733	0.0000	1.0000
F81G	15	-3587.8791	7205.7582	0.0000	1.0000
F81IG	16	-3587.9462	7207.8924	0.0000	1.0000
JC69I	12	-3614.0569	7252.1138	0.0000	1.0000
JC69G	12	-3614.3585	7252.7170	0.0000	1.0000
JC69IG	13	-3614.4268	7254.8536	0.0000	1.0000
F81	14	-3619.0182	7266.0365	0.0000	1.0000
JC69	11	-3645.5392	7313.0785	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3509.3630	7131.5695	0.5147	0.5147
HKYI	16	-3509.4509	7131.7453	0.4714	0.9861
HKYIG	17	-3509.4499	7138.7960	0.0139	1.0000
GTRG	20	-3505.6936	7152.4416	0.0000	1.0000
GTRI	20	-3505.8052	7152.6648	0.0000	1.0000
GTRIG	21	-3505.7840	7159.6751	0.0000	1.0000
K2PI	13	-3542.1986	7176.0826	0.0000	1.0000
K2PG	13	-3542.2909	7176.2671	0.0000	1.0000
K2PIG	14	-3542.3825	7183.5031	0.0000	1.0000
HKY	15	-3544.9056	7195.6021	0.0000	1.0000
SYMI	17	-3539.8193	7199.5348	0.0000	1.0000
SYMG	17	-3540.0492	7199.9947	0.0000	1.0000
SYMIG	18	-3540.1421	7207.2332	0.0000	1.0000
GTR	19	-3541.8262	7217.6540	0.0000	1.0000
K2P	12	-3577.0753	7238.7832	0.0000	1.0000
SYM	16	-3574.9061	7262.6558	0.0000	1.0000
F81I	15	-3587.6367	7281.0641	0.0000	1.0000
F81G	15	-3587.8791	7281.5490	0.0000	1.0000
F81IG	16	-3587.9462	7288.7360	0.0000	1.0000
JC69I	12	-3614.0569	7312.7464	0.0000	1.0000
JC69G	12	-3614.3585	7313.3497	0.0000	1.0000
JC69IG	13	-3614.4268	7320.5390	0.0000	1.0000
F81	14	-3619.0182	7336.7746	0.0000	1.0000
JC69	11	-3645.5392	7368.6584	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_1674.phylip = 1 - 1156;
 Partition Dummy = 1:591-tomodel-phylip/chr5_1674.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 AICc model (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_1674.phylip = 1 - 1156;
 Partition Dummy = 1:591-tomodel-phylip/chr5_1674.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 591-tomodel-phylip/chr5_1674.phylip = 1 - 1156;
 Partition Dummy = 1:591-tomodel-phylip/chr5_1674.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

