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

Input data from file "683-tomodel-phylip/chr2_7985.phylip" (7 taxa, 930 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	-2447.9514	4927.9029	0.4795	0.4795
HKYI	16	-2448.5650	4929.1301	0.2596	0.7391
HKYIG	17	-2447.9003	4929.8005	0.1857	0.9247
GTRG	20	-2446.4706	4932.9413	0.0386	0.9633
GTRI	20	-2447.0330	4934.0661	0.0220	0.9853
GTRIG	21	-2446.4389	4934.8779	0.0147	1.0000
HKY	15	-2472.2486	4974.4973	0.0000	1.0000
GTR	19	-2470.7370	4979.4741	0.0000	1.0000
K2PG	13	-2493.0698	5012.1397	0.0000	1.0000
K2PI	13	-2493.5287	5013.0573	0.0000	1.0000
K2PIG	14	-2493.0251	5014.0502	0.0000	1.0000
SYMG	17	-2491.1294	5016.2588	0.0000	1.0000
SYMI	17	-2491.5549	5017.1097	0.0000	1.0000
F81G	15	-2493.7160	5017.4319	0.0000	1.0000
F81I	15	-2494.0032	5018.0063	0.0000	1.0000
SYMIG	18	-2491.0904	5018.1808	0.0000	1.0000
F81IG	16	-2493.7236	5019.4472	0.0000	1.0000
F81	14	-2515.3114	5058.6228	0.0000	1.0000
K2P	12	-2518.4717	5060.9435	0.0000	1.0000
SYM	16	-2516.6465	5065.2931	0.0000	1.0000
JC69G	12	-2533.8467	5091.6934	0.0000	1.0000
JC69I	12	-2534.0924	5092.1848	0.0000	1.0000
JC69IG	13	-2533.8559	5093.7118	0.0000	1.0000
JC69	11	-2555.8575	5133.7150	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2447.9514	5005.2658	0.6346	0.6346
HKYI	16	-2448.5650	5006.4930	0.3435	0.9781
HKYIG	17	-2447.9003	5011.9987	0.0219	1.0000
GTRG	20	-2446.4706	5029.6450	0.0000	1.0000
GTRI	20	-2447.0330	5030.7698	0.0000	1.0000
GTRIG	21	-2446.4389	5036.4168	0.0000	1.0000
HKY	15	-2472.2486	5047.0250	0.0000	1.0000
GTR	19	-2470.7370	5071.3426	0.0000	1.0000
K2PG	13	-2493.0698	5074.9971	0.0000	1.0000
K2PI	13	-2493.5287	5075.9147	0.0000	1.0000
K2PIG	14	-2493.0251	5081.7428	0.0000	1.0000
F81G	15	-2493.7160	5089.9597	0.0000	1.0000
F81I	15	-2494.0032	5090.5341	0.0000	1.0000
F81IG	16	-2493.7236	5096.8101	0.0000	1.0000
SYMG	17	-2491.1294	5098.4569	0.0000	1.0000
SYMI	17	-2491.5549	5099.3079	0.0000	1.0000
SYMIG	18	-2491.0904	5105.2141	0.0000	1.0000
K2P	12	-2518.4717	5118.9657	0.0000	1.0000
F81	14	-2515.3114	5126.3153	0.0000	1.0000
SYM	16	-2516.6465	5142.6560	0.0000	1.0000
JC69G	12	-2533.8467	5149.7156	0.0000	1.0000
JC69I	12	-2534.0924	5150.2070	0.0000	1.0000
JC69IG	13	-2533.8559	5156.5692	0.0000	1.0000
JC69	11	-2555.8575	5186.9021	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 683-tomodel-phylip/chr2_7985.phylip = 1 - 930;
 Partition Dummy = 1:683-tomodel-phylip/chr2_7985.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 683-tomodel-phylip/chr2_7985.phylip = 1 - 930;
 Partition Dummy = 1:683-tomodel-phylip/chr2_7985.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_7985.phylip = 1 - 930;
 Partition Dummy = 1:683-tomodel-phylip/chr2_7985.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

