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

Input data from file "591-tomodel-phylip/chr12_5837.phylip" (7 taxa, 916 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	-2090.2746	4212.5493	0.3683	0.3683
HKYI	16	-2090.8786	4213.7573	0.2013	0.5696
HKYIG	17	-2089.9525	4213.9049	0.1870	0.7566
GTRG	20	-2087.4137	4214.8274	0.1179	0.8745
GTRIG	21	-2087.0078	4216.0157	0.0651	0.9396
GTRI	20	-2088.0818	4216.1636	0.0604	1.0000
HKY	15	-2104.8515	4239.7031	0.0000	1.0000
GTR	19	-2102.5144	4243.0288	0.0000	1.0000
K2PG	13	-2108.7560	4243.5119	0.0000	1.0000
K2PI	13	-2109.3699	4244.7398	0.0000	1.0000
K2PIG	14	-2108.4067	4244.8135	0.0000	1.0000
SYMG	17	-2106.3590	4246.7181	0.0000	1.0000
SYMIG	18	-2105.9334	4247.8668	0.0000	1.0000
SYMI	17	-2107.1032	4248.2064	0.0000	1.0000
K2P	12	-2123.1487	4270.2974	0.0000	1.0000
SYM	16	-2121.0517	4274.1033	0.0000	1.0000
F81G	15	-2128.5563	4287.1126	0.0000	1.0000
F81I	15	-2129.0607	4288.1214	0.0000	1.0000
F81IG	16	-2128.2956	4288.5912	0.0000	1.0000
F81	14	-2142.7367	4313.4734	0.0000	1.0000
JC69G	12	-2147.8263	4319.6526	0.0000	1.0000
JC69I	12	-2148.3882	4320.7763	0.0000	1.0000
JC69IG	13	-2147.5491	4321.0982	0.0000	1.0000
JC69	11	-2161.6789	4345.3578	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2090.2746	4289.6695	0.6279	0.6279
HKYI	16	-2090.8786	4290.8775	0.3432	0.9711
HKYIG	17	-2089.9525	4295.8452	0.0286	0.9997
K2PG	13	-2108.7560	4306.1721	0.0002	0.9999
K2PI	13	-2109.3699	4307.4001	0.0001	1.0000
GTRG	20	-2087.4137	4311.2277	0.0000	1.0000
HKY	15	-2104.8515	4312.0033	0.0000	1.0000
K2PIG	14	-2108.4067	4312.2937	0.0000	1.0000
GTRI	20	-2088.0818	4312.5639	0.0000	1.0000
GTRIG	21	-2087.0078	4317.2360	0.0000	1.0000
K2P	12	-2123.1487	4328.1376	0.0000	1.0000
SYMG	17	-2106.3590	4328.6584	0.0000	1.0000
SYMI	17	-2107.1032	4330.1467	0.0000	1.0000
GTR	19	-2102.5144	4334.6091	0.0000	1.0000
SYMIG	18	-2105.9334	4334.6271	0.0000	1.0000
SYM	16	-2121.0517	4351.2236	0.0000	1.0000
F81G	15	-2128.5563	4359.4129	0.0000	1.0000
F81I	15	-2129.0607	4360.4217	0.0000	1.0000
F81IG	16	-2128.2956	4365.7114	0.0000	1.0000
JC69G	12	-2147.8263	4377.4928	0.0000	1.0000
JC69I	12	-2148.3882	4378.6165	0.0000	1.0000
F81	14	-2142.7367	4380.9537	0.0000	1.0000
JC69IG	13	-2147.5491	4383.7584	0.0000	1.0000
JC69	11	-2161.6789	4398.3780	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/chr12_5837.phylip = 1 - 916;
 Partition Dummy = 1:591-tomodel-phylip/chr12_5837.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/chr12_5837.phylip = 1 - 916;
 Partition Dummy = 1:591-tomodel-phylip/chr12_5837.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/chr12_5837.phylip = 1 - 916;
 Partition Dummy = 1:591-tomodel-phylip/chr12_5837.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

