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

Input data from file "591-tomodel-phylip/chr7_10298.phylip" (7 taxa, 878 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	-2493.6774	5027.3549	0.5697	0.5697
GTRI	20	-2494.6398	5029.2795	0.2176	0.7873
GTRIG	21	-2493.6711	5029.3422	0.2109	0.9982
HKYG	16	-2504.1236	5040.2472	0.0009	0.9991
HKYI	16	-2504.6375	5041.2750	0.0005	0.9997
HKYIG	17	-2504.1307	5042.2614	0.0003	1.0000
GTR	19	-2512.1339	5062.2678	0.0000	1.0000
HKY	15	-2520.7149	5071.4298	0.0000	1.0000
F81G	15	-2543.9236	5117.8471	0.0000	1.0000
F81I	15	-2544.4567	5118.9135	0.0000	1.0000
F81IG	16	-2543.9361	5119.8722	0.0000	1.0000
K2PG	13	-2552.7547	5131.5095	0.0000	1.0000
K2PI	13	-2553.2232	5132.4464	0.0000	1.0000
K2PIG	14	-2552.7754	5133.5507	0.0000	1.0000
SYMG	17	-2550.5182	5135.0364	0.0000	1.0000
SYMI	17	-2550.9076	5135.8152	0.0000	1.0000
SYMIG	18	-2550.5372	5137.0744	0.0000	1.0000
F81	14	-2559.5714	5147.1429	0.0000	1.0000
K2P	12	-2568.8868	5161.7736	0.0000	1.0000
SYM	16	-2566.2714	5164.5428	0.0000	1.0000
JC69G	12	-2584.7517	5193.5035	0.0000	1.0000
JC69I	12	-2585.2367	5194.4734	0.0000	1.0000
JC69IG	13	-2584.7723	5195.5446	0.0000	1.0000
JC69	11	-2599.7827	5221.5654	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2504.1236	5116.6895	0.5900	0.5900
HKYI	16	-2504.6375	5117.7173	0.3529	0.9429
GTRG	20	-2493.6774	5122.9078	0.0263	0.9693
HKYIG	17	-2504.1307	5123.4814	0.0198	0.9890
GTRI	20	-2494.6398	5124.8325	0.0101	0.9991
GTRIG	21	-2493.6711	5129.6728	0.0009	1.0000
HKY	15	-2520.7149	5143.0945	0.0000	1.0000
GTR	19	-2512.1339	5153.0431	0.0000	1.0000
F81G	15	-2543.9236	5189.5118	0.0000	1.0000
F81I	15	-2544.4567	5190.5782	0.0000	1.0000
K2PG	13	-2552.7547	5193.6189	0.0000	1.0000
K2PI	13	-2553.2232	5194.5558	0.0000	1.0000
F81IG	16	-2543.9361	5196.3145	0.0000	1.0000
K2PIG	14	-2552.7754	5200.4378	0.0000	1.0000
F81	14	-2559.5714	5214.0300	0.0000	1.0000
SYMG	17	-2550.5182	5216.2564	0.0000	1.0000
SYMI	17	-2550.9076	5217.0352	0.0000	1.0000
K2P	12	-2568.8868	5219.1054	0.0000	1.0000
SYMIG	18	-2550.5372	5223.0721	0.0000	1.0000
SYM	16	-2566.2714	5240.9852	0.0000	1.0000
JC69G	12	-2584.7517	5250.8352	0.0000	1.0000
JC69I	12	-2585.2367	5251.8052	0.0000	1.0000
JC69IG	13	-2584.7723	5257.6540	0.0000	1.0000
JC69	11	-2599.7827	5274.1195	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/chr7_10298.phylip = 1 - 878;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10298.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/chr7_10298.phylip = 1 - 878;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10298.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/chr7_10298.phylip = 1 - 878;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10298.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

