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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYI


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2627.5118	5295.0235	0.4756	0.4756
GTRI	20	-2627.8155	5295.6310	0.3510	0.8267
GTRIG	21	-2627.5221	5297.0442	0.1732	0.9998
HKYI	16	-2640.3057	5312.6113	0.0001	0.9999
HKYG	16	-2640.5489	5313.0978	0.0001	1.0000
HKYIG	17	-2640.3089	5314.6179	0.0000	1.0000
GTR	19	-2644.4729	5326.9459	0.0000	1.0000
SYMG	17	-2647.9837	5329.9674	0.0000	1.0000
SYMI	17	-2648.1739	5330.3477	0.0000	1.0000
SYMIG	18	-2648.0191	5332.0382	0.0000	1.0000
K2PI	13	-2655.3371	5336.6742	0.0000	1.0000
K2PG	13	-2655.4759	5336.9517	0.0000	1.0000
K2PIG	14	-2655.4965	5338.9930	0.0000	1.0000
HKY	15	-2657.7292	5345.4585	0.0000	1.0000
SYM	16	-2665.3520	5362.7039	0.0000	1.0000
K2P	12	-2673.3674	5370.7347	0.0000	1.0000
F81G	15	-2687.5004	5405.0007	0.0000	1.0000
F81I	15	-2687.6029	5405.2058	0.0000	1.0000
F81IG	16	-2687.5107	5407.0214	0.0000	1.0000
JC69G	12	-2700.6211	5425.2423	0.0000	1.0000
JC69I	12	-2700.7373	5425.4745	0.0000	1.0000
JC69IG	13	-2700.6286	5427.2572	0.0000	1.0000
F81	14	-2704.2109	5436.4217	0.0000	1.0000
JC69	11	-2718.4706	5458.9411	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2640.3057	5388.9989	0.3752	0.3752
HKYG	16	-2640.5489	5389.4854	0.2942	0.6694
GTRG	20	-2627.5118	5390.5080	0.1764	0.8458
GTRI	20	-2627.8155	5391.1155	0.1302	0.9760
HKYIG	17	-2640.3089	5395.7797	0.0126	0.9886
GTRIG	21	-2627.5221	5397.3029	0.0059	0.9945
K2PI	13	-2655.3371	5398.7391	0.0029	0.9974
K2PG	13	-2655.4759	5399.0167	0.0025	0.9999
K2PIG	14	-2655.4965	5405.8321	0.0001	1.0000
SYMG	17	-2647.9837	5411.1292	0.0000	1.0000
SYMI	17	-2648.1739	5411.5095	0.0000	1.0000
HKY	15	-2657.7292	5417.0718	0.0000	1.0000
GTR	19	-2644.4729	5417.6561	0.0000	1.0000
SYMIG	18	-2648.0191	5417.9742	0.0000	1.0000
K2P	12	-2673.3674	5428.0254	0.0000	1.0000
SYM	16	-2665.3520	5439.0915	0.0000	1.0000
F81G	15	-2687.5004	5476.6141	0.0000	1.0000
F81I	15	-2687.6029	5476.8191	0.0000	1.0000
JC69G	12	-2700.6211	5482.5329	0.0000	1.0000
JC69I	12	-2700.7373	5482.7652	0.0000	1.0000
F81IG	16	-2687.5107	5483.4090	0.0000	1.0000
JC69IG	13	-2700.6286	5489.3221	0.0000	1.0000
F81	14	-2704.2109	5503.2609	0.0000	1.0000
JC69	11	-2718.4706	5511.4576	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_371.phylip = 1 - 875;
 Partition Dummy = 1:591-tomodel-phylip/chr3_371.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_371.phylip = 1 - 875;
 Partition Dummy = 1:591-tomodel-phylip/chr3_371.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 (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_371.phylip = 1 - 875;
 Partition Dummy = 1:591-tomodel-phylip/chr3_371.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

