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

Input data from file "591-tomodel-phylip/chr1_28747.phylip" (7 taxa, 927 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	-3045.7939	6131.5877	0.4243	0.4243
GTRI	20	-3045.8091	6131.6183	0.4179	0.8421
GTRIG	21	-3045.7941	6133.5882	0.1560	0.9982
HKYG	16	-3056.0544	6144.1087	0.0008	0.9990
HKYI	16	-3056.2104	6144.4207	0.0007	0.9997
HKYIG	17	-3056.0545	6146.1091	0.0003	1.0000
K2PG	13	-3096.3985	6218.7970	0.0000	1.0000
SYMG	17	-3092.5688	6219.1376	0.0000	1.0000
K2PI	13	-3097.2884	6220.5768	0.0000	1.0000
SYMI	17	-3093.3962	6220.7924	0.0000	1.0000
K2PIG	14	-3096.3985	6220.7970	0.0000	1.0000
SYMIG	18	-3092.5688	6221.1376	0.0000	1.0000
F81G	15	-3135.3217	6300.6434	0.0000	1.0000
F81I	15	-3135.7276	6301.4553	0.0000	1.0000
GTR	19	-3131.9389	6301.8778	0.0000	1.0000
F81IG	16	-3135.3219	6302.6437	0.0000	1.0000
HKY	15	-3141.5690	6313.1380	0.0000	1.0000
JC69G	12	-3166.2765	6356.5529	0.0000	1.0000
JC69I	12	-3166.7796	6357.5592	0.0000	1.0000
JC69IG	13	-3166.2766	6358.5531	0.0000	1.0000
SYM	16	-3175.7541	6383.5081	0.0000	1.0000
K2P	12	-3180.2267	6384.4534	0.0000	1.0000
F81	14	-3214.8826	6457.7652	0.0000	1.0000
JC69	11	-3245.4032	6512.8064	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3056.0544	6221.4200	0.5114	0.5114
HKYI	16	-3056.2104	6221.7320	0.4375	0.9489
GTRG	20	-3045.7939	6228.2268	0.0170	0.9659
HKYIG	17	-3056.0545	6228.2523	0.0168	0.9827
GTRI	20	-3045.8091	6228.2574	0.0168	0.9994
GTRIG	21	-3045.7941	6235.0592	0.0006	1.0000
K2PG	13	-3096.3985	6281.6124	0.0000	1.0000
K2PI	13	-3097.2884	6283.3922	0.0000	1.0000
K2PIG	14	-3096.3985	6288.4444	0.0000	1.0000
SYMG	17	-3092.5688	6301.2808	0.0000	1.0000
SYMI	17	-3093.3962	6302.9356	0.0000	1.0000
SYMIG	18	-3092.5688	6308.1128	0.0000	1.0000
F81G	15	-3135.3217	6373.1227	0.0000	1.0000
F81I	15	-3135.7276	6373.9346	0.0000	1.0000
F81IG	16	-3135.3219	6379.9550	0.0000	1.0000
HKY	15	-3141.5690	6385.6173	0.0000	1.0000
GTR	19	-3131.9389	6393.6849	0.0000	1.0000
JC69G	12	-3166.2765	6414.5364	0.0000	1.0000
JC69I	12	-3166.7796	6415.5426	0.0000	1.0000
JC69IG	13	-3166.2766	6421.3685	0.0000	1.0000
K2P	12	-3180.2267	6442.4369	0.0000	1.0000
SYM	16	-3175.7541	6460.8194	0.0000	1.0000
F81	14	-3214.8826	6525.4125	0.0000	1.0000
JC69	11	-3245.4032	6565.9578	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/chr1_28747.phylip = 1 - 927;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28747.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/chr1_28747.phylip = 1 - 927;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28747.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/chr1_28747.phylip = 1 - 927;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28747.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

