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

Input data from file "683-tomodel-phylip/chr1_32446.phylip" (7 taxa, 957 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	-2113.4969	4266.9938	0.4672	0.4672
GTRI	20	-2114.1148	4268.2297	0.2518	0.7191
GTRIG	21	-2113.4601	4268.9202	0.1783	0.8974
HKYG	16	-2119.7844	4271.5688	0.0474	0.9448
HKYI	16	-2120.3358	4272.6716	0.0273	0.9721
HKYIG	17	-2119.7638	4273.5276	0.0178	0.9899
GTR	19	-2118.4675	4274.9349	0.0088	0.9987
HKY	15	-2124.4180	4278.8359	0.0013	1.0000
F81G	15	-2150.6912	4331.3825	0.0000	1.0000
F81I	15	-2151.0851	4332.1702	0.0000	1.0000
F81IG	16	-2150.6821	4333.3642	0.0000	1.0000
F81	14	-2154.9474	4337.8948	0.0000	1.0000
K2PG	13	-2181.9184	4389.8368	0.0000	1.0000
K2PI	13	-2182.3753	4390.7506	0.0000	1.0000
K2PIG	14	-2181.9074	4391.8149	0.0000	1.0000
SYMG	17	-2179.0154	4392.0307	0.0000	1.0000
SYMI	17	-2179.4702	4392.9404	0.0000	1.0000
SYMIG	18	-2179.0036	4394.0072	0.0000	1.0000
K2P	12	-2186.3523	4396.7046	0.0000	1.0000
SYM	16	-2183.2519	4398.5039	0.0000	1.0000
JC69G	12	-2207.5941	4439.1882	0.0000	1.0000
JC69I	12	-2207.9498	4439.8997	0.0000	1.0000
JC69IG	13	-2207.5874	4441.1747	0.0000	1.0000
JC69	11	-2211.7591	4445.5182	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2119.7844	4349.3897	0.5233	0.5233
HKYI	16	-2120.3358	4350.4925	0.3015	0.8248
HKY	15	-2124.4180	4351.7930	0.1574	0.9822
HKYIG	17	-2119.7638	4356.2123	0.0173	0.9995
GTRG	20	-2113.4969	4364.2698	0.0003	0.9998
GTRI	20	-2114.1148	4365.5057	0.0002	0.9999
GTR	19	-2118.4675	4367.3472	0.0001	1.0000
GTRIG	21	-2113.4601	4371.0601	0.0000	1.0000
F81G	15	-2150.6912	4404.3396	0.0000	1.0000
F81I	15	-2151.0851	4405.1273	0.0000	1.0000
F81	14	-2154.9474	4405.9880	0.0000	1.0000
F81IG	16	-2150.6821	4411.1851	0.0000	1.0000
K2PG	13	-2181.9184	4453.0662	0.0000	1.0000
K2PI	13	-2182.3753	4453.9800	0.0000	1.0000
K2P	12	-2186.3523	4455.0702	0.0000	1.0000
K2PIG	14	-2181.9074	4459.9081	0.0000	1.0000
SYMG	17	-2179.0154	4474.7154	0.0000	1.0000
SYMI	17	-2179.4702	4475.6251	0.0000	1.0000
SYM	16	-2183.2519	4476.3247	0.0000	1.0000
SYMIG	18	-2179.0036	4481.5557	0.0000	1.0000
JC69G	12	-2207.5941	4497.5538	0.0000	1.0000
JC69I	12	-2207.9498	4498.2653	0.0000	1.0000
JC69	11	-2211.7591	4499.0200	0.0000	1.0000
JC69IG	13	-2207.5874	4504.4042	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/chr1_32446.phylip = 1 - 957;
 Partition Dummy = 1:683-tomodel-phylip/chr1_32446.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/chr1_32446.phylip = 1 - 957;
 Partition Dummy = 1:683-tomodel-phylip/chr1_32446.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 683-tomodel-phylip/chr1_32446.phylip = 1 - 957;
 Partition Dummy = 1:683-tomodel-phylip/chr1_32446.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

