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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2120.2156	4280.4312	0.4241	0.4241
GTRI	20	-2120.2242	4280.4483	0.4205	0.8446
GTRIG	21	-2120.2255	4282.4510	0.1545	0.9991
HKYG	16	-2131.2108	4294.4216	0.0004	0.9995
HKYI	16	-2131.3642	4294.7283	0.0003	0.9998
HKYIG	17	-2131.0364	4296.0728	0.0002	1.0000
SYMG	17	-2142.0353	4318.0706	0.0000	1.0000
SYMI	17	-2142.1192	4318.2383	0.0000	1.0000
GTR	19	-2140.1919	4318.3838	0.0000	1.0000
SYMIG	18	-2141.9364	4319.8728	0.0000	1.0000
K2PI	13	-2149.7132	4325.4264	0.0000	1.0000
K2PG	13	-2151.1158	4328.2316	0.0000	1.0000
HKY	15	-2149.8669	4329.7337	0.0000	1.0000
K2PIG	14	-2150.9038	4329.8075	0.0000	1.0000
F81G	15	-2156.6829	4343.3658	0.0000	1.0000
F81I	15	-2156.7455	4343.4909	0.0000	1.0000
F81IG	16	-2156.6526	4345.3052	0.0000	1.0000
SYM	16	-2160.5006	4353.0012	0.0000	1.0000
K2P	12	-2167.1791	4358.3581	0.0000	1.0000
JC69G	12	-2173.1009	4370.2019	0.0000	1.0000
JC69I	12	-2173.1819	4370.3638	0.0000	1.0000
F81	14	-2171.3352	4370.6705	0.0000	1.0000
JC69IG	13	-2173.0581	4372.1162	0.0000	1.0000
JC69	11	-2187.6750	4397.3500	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2131.2108	4372.2759	0.4932	0.4932
HKYI	16	-2131.3642	4372.5826	0.4231	0.9162
GTRG	20	-2120.2156	4377.7490	0.0320	0.9482
GTRI	20	-2120.2242	4377.7661	0.0317	0.9798
HKYIG	17	-2131.0364	4378.7929	0.0190	0.9988
GTRIG	21	-2120.2255	4384.6348	0.0010	0.9998
K2PI	13	-2149.7132	4388.6830	0.0001	1.0000
K2PG	13	-2151.1158	4391.4881	0.0000	1.0000
K2PIG	14	-2150.9038	4397.9300	0.0000	1.0000
SYMG	17	-2142.0353	4400.7907	0.0000	1.0000
SYMI	17	-2142.1192	4400.9585	0.0000	1.0000
HKY	15	-2149.8669	4402.7221	0.0000	1.0000
SYMIG	18	-2141.9364	4407.4589	0.0000	1.0000
GTR	19	-2140.1919	4410.8357	0.0000	1.0000
F81G	15	-2156.6829	4416.3542	0.0000	1.0000
F81I	15	-2156.7455	4416.4793	0.0000	1.0000
K2P	12	-2167.1791	4416.7488	0.0000	1.0000
F81IG	16	-2156.6526	4423.1595	0.0000	1.0000
JC69G	12	-2173.1009	4428.5926	0.0000	1.0000
JC69I	12	-2173.1819	4428.7545	0.0000	1.0000
SYM	16	-2160.5006	4430.8554	0.0000	1.0000
JC69IG	13	-2173.0581	4435.3728	0.0000	1.0000
F81	14	-2171.3352	4438.7930	0.0000	1.0000
JC69	11	-2187.6750	4450.8748	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 591-tomodel-phylip/chr3_17883.phylip = 1 - 959;
 Partition Dummy = 1:591-tomodel-phylip/chr3_17883.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 (GTRG)]

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

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

