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

Input data from file "591-tomodel-phylip/chr9_1164.phylip" (7 taxa, 977 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	-2782.2266	5596.4531	0.3956	0.3956
HKYIG	17	-2782.0441	5598.0881	0.1747	0.5703
GTRG	20	-2779.0829	5598.1658	0.1680	0.7383
HKYI	16	-2783.2709	5598.5417	0.1392	0.8775
GTRIG	21	-2778.8687	5599.7374	0.0766	0.9541
GTRI	20	-2780.3802	5600.7604	0.0459	1.0000
K2PG	13	-2807.1328	5640.2656	0.0000	1.0000
K2PIG	14	-2806.9434	5641.8868	0.0000	1.0000
K2PI	13	-2808.0051	5642.0101	0.0000	1.0000
SYMG	17	-2805.7833	5645.5667	0.0000	1.0000
SYMIG	18	-2805.6313	5647.2627	0.0000	1.0000
SYMI	17	-2806.6339	5647.2679	0.0000	1.0000
GTR	19	-2819.4780	5676.9559	0.0000	1.0000
HKY	15	-2823.4862	5676.9724	0.0000	1.0000
K2P	12	-2849.9100	5723.8200	0.0000	1.0000
SYM	16	-2848.4341	5728.8682	0.0000	1.0000
F81G	15	-2851.8373	5733.6747	0.0000	1.0000
F81I	15	-2852.7287	5735.4575	0.0000	1.0000
F81IG	16	-2851.7520	5735.5040	0.0000	1.0000
JC69G	12	-2874.5593	5773.1186	0.0000	1.0000
JC69I	12	-2875.3692	5774.7384	0.0000	1.0000
JC69IG	13	-2874.4760	5774.9521	0.0000	1.0000
F81	14	-2890.4980	5808.9961	0.0000	1.0000
JC69	11	-2914.2635	5850.5269	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2782.2266	5674.6049	0.7192	0.7192
HKYI	16	-2783.2709	5676.6935	0.2531	0.9724
HKYIG	17	-2782.0441	5681.1244	0.0276	1.0000
GTRG	20	-2779.0829	5695.8556	0.0000	1.0000
GTRI	20	-2780.3802	5698.4501	0.0000	1.0000
GTRIG	21	-2778.8687	5702.3116	0.0000	1.0000
K2PG	13	-2807.1328	5703.7639	0.0000	1.0000
K2PI	13	-2808.0051	5705.5084	0.0000	1.0000
K2PIG	14	-2806.9434	5710.2696	0.0000	1.0000
SYMG	17	-2805.7833	5728.6030	0.0000	1.0000
SYMI	17	-2806.6339	5730.3042	0.0000	1.0000
SYMIG	18	-2805.6313	5735.1834	0.0000	1.0000
HKY	15	-2823.4862	5750.2397	0.0000	1.0000
GTR	19	-2819.4780	5769.7612	0.0000	1.0000
K2P	12	-2849.9100	5782.4339	0.0000	1.0000
F81G	15	-2851.8373	5806.9420	0.0000	1.0000
SYM	16	-2848.4341	5807.0199	0.0000	1.0000
F81I	15	-2852.7287	5808.7248	0.0000	1.0000
F81IG	16	-2851.7520	5813.6558	0.0000	1.0000
JC69G	12	-2874.5593	5831.7324	0.0000	1.0000
JC69I	12	-2875.3692	5833.3522	0.0000	1.0000
JC69IG	13	-2874.4760	5838.4504	0.0000	1.0000
F81	14	-2890.4980	5877.3789	0.0000	1.0000
JC69	11	-2914.2635	5904.2563	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/chr9_1164.phylip = 1 - 977;
 Partition Dummy = 1:591-tomodel-phylip/chr9_1164.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/chr9_1164.phylip = 1 - 977;
 Partition Dummy = 1:591-tomodel-phylip/chr9_1164.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/chr9_1164.phylip = 1 - 977;
 Partition Dummy = 1:591-tomodel-phylip/chr9_1164.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

