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

Input data from file "683-tomodel-phylip/chr2_18959.phylip" (7 taxa, 892 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	-2078.5679	4197.1358	0.4249	0.4249
GTRI	20	-2078.6383	4197.2765	0.3961	0.8210
GTRIG	21	-2078.5536	4199.1073	0.1586	0.9796
HKYI	16	-2086.4410	4204.8821	0.0088	0.9884
HKYG	16	-2086.4811	4204.9622	0.0085	0.9969
HKYIG	17	-2086.4811	4206.9623	0.0031	1.0000
GTR	19	-2114.4025	4266.8050	0.0000	1.0000
F81G	15	-2119.3445	4268.6889	0.0000	1.0000
F81I	15	-2119.3603	4268.7206	0.0000	1.0000
F81IG	16	-2119.3439	4270.6879	0.0000	1.0000
HKY	15	-2121.5725	4273.1450	0.0000	1.0000
F81	14	-2152.4956	4332.9912	0.0000	1.0000
K2PG	13	-2167.5862	4361.1724	0.0000	1.0000
K2PI	13	-2167.5965	4361.1929	0.0000	1.0000
K2PIG	14	-2167.5841	4363.1682	0.0000	1.0000
SYMG	17	-2166.4963	4366.9926	0.0000	1.0000
SYMI	17	-2166.5145	4367.0290	0.0000	1.0000
SYMIG	18	-2166.4727	4368.9454	0.0000	1.0000
JC69G	12	-2192.4541	4408.9081	0.0000	1.0000
JC69I	12	-2192.4989	4408.9978	0.0000	1.0000
JC69IG	13	-2192.4046	4410.8091	0.0000	1.0000
K2P	12	-2204.0601	4432.1201	0.0000	1.0000
SYM	16	-2202.7951	4437.5902	0.0000	1.0000
JC69	11	-2227.7180	4477.4360	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2086.4410	4281.5775	0.5002	0.5002
HKYG	16	-2086.4811	4281.6577	0.4805	0.9807
HKYIG	17	-2086.4811	4288.4512	0.0161	0.9968
GTRG	20	-2078.5679	4293.0051	0.0017	0.9984
GTRI	20	-2078.6383	4293.1458	0.0015	0.9999
GTRIG	21	-2078.5536	4299.7701	0.0001	1.0000
F81G	15	-2119.3445	4340.5909	0.0000	1.0000
F81I	15	-2119.3603	4340.6226	0.0000	1.0000
HKY	15	-2121.5725	4345.0470	0.0000	1.0000
F81IG	16	-2119.3439	4347.3834	0.0000	1.0000
GTR	19	-2114.4025	4357.8808	0.0000	1.0000
F81	14	-2152.4956	4400.0997	0.0000	1.0000
K2PG	13	-2167.5862	4423.4875	0.0000	1.0000
K2PI	13	-2167.5965	4423.5080	0.0000	1.0000
K2PIG	14	-2167.5841	4430.2767	0.0000	1.0000
SYMG	17	-2166.4963	4448.4816	0.0000	1.0000
SYMI	17	-2166.5145	4448.5179	0.0000	1.0000
SYMIG	18	-2166.4727	4455.2278	0.0000	1.0000
JC69G	12	-2192.4541	4466.4297	0.0000	1.0000
JC69I	12	-2192.4989	4466.5194	0.0000	1.0000
JC69IG	13	-2192.4046	4473.1242	0.0000	1.0000
K2P	12	-2204.0601	4489.6417	0.0000	1.0000
SYM	16	-2202.7951	4514.2857	0.0000	1.0000
JC69	11	-2227.7180	4530.1641	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 683-tomodel-phylip/chr2_18959.phylip = 1 - 892;
 Partition Dummy = 1:683-tomodel-phylip/chr2_18959.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 683-tomodel-phylip/chr2_18959.phylip = 1 - 892;
 Partition Dummy = 1:683-tomodel-phylip/chr2_18959.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 683-tomodel-phylip/chr2_18959.phylip = 1 - 892;
 Partition Dummy = 1:683-tomodel-phylip/chr2_18959.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

