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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYG


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


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

Model	df	lnL		AICc		wAICc	cwAICc
GTRG	20	-2144.7727	4330.7025	0.2730	0.2730
HKYG	16	-2149.1130	4330.9712	0.2386	0.5116
HKYI	16	-2149.5374	4331.8200	0.1561	0.6677
GTRI	20	-2145.3548	4331.8666	0.1525	0.8202
GTRIG	21	-2144.7692	4332.8128	0.0950	0.9153
HKYIG	17	-2149.1017	4333.0429	0.0847	1.0000
GTR	19	-2155.5323	4350.1100	0.0000	1.0000
HKY	15	-2160.1738	4351.0042	0.0000	1.0000
K2PG	13	-2164.6575	4355.8116	0.0000	1.0000
K2PI	13	-2164.7829	4356.0624	0.0000	1.0000
K2PIG	14	-2165.4171	4359.4079	0.0000	1.0000
SYMG	17	-2162.7234	4360.2864	0.0000	1.0000
SYMI	17	-2162.9717	4360.7829	0.0000	1.0000
SYMIG	18	-2163.4252	4363.7899	0.0000	1.0000
K2P	12	-2177.6467	4379.7185	0.0000	1.0000
SYM	16	-2175.1790	4383.1031	0.0000	1.0000
F81I	15	-2215.7282	4462.1131	0.0000	1.0000
F81G	15	-2215.8046	4462.2658	0.0000	1.0000
F81IG	16	-2215.7205	4464.1861	0.0000	1.0000
F81	14	-2226.7655	4482.1047	0.0000	1.0000
JC69I	12	-2229.6764	4483.7779	0.0000	1.0000
JC69G	12	-2229.7607	4483.9465	0.0000	1.0000
JC69IG	13	-2229.7631	4486.0227	0.0000	1.0000
JC69	11	-2241.4806	4505.3204	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2144.7727	4329.5455	0.2967	0.2967
HKYG	16	-2149.1130	4330.2260	0.2112	0.5079
GTRI	20	-2145.3548	4330.7096	0.1658	0.6737
HKYI	16	-2149.5374	4331.0748	0.1381	0.8118
GTRIG	21	-2144.7692	4331.5383	0.1096	0.9214
HKYIG	17	-2149.1017	4332.2034	0.0786	1.0000
GTR	19	-2155.5323	4349.0646	0.0000	1.0000
HKY	15	-2160.1738	4350.3476	0.0000	1.0000
K2PG	13	-2164.6575	4355.3150	0.0000	1.0000
K2PI	13	-2164.7829	4355.5658	0.0000	1.0000
K2PIG	14	-2165.4171	4358.8341	0.0000	1.0000
SYMG	17	-2162.7234	4359.4469	0.0000	1.0000
SYMI	17	-2162.9717	4359.9434	0.0000	1.0000
SYMIG	18	-2163.4252	4362.8503	0.0000	1.0000
K2P	12	-2177.6467	4379.2934	0.0000	1.0000
SYM	16	-2175.1790	4382.3579	0.0000	1.0000
F81I	15	-2215.7282	4461.4564	0.0000	1.0000
F81G	15	-2215.8046	4461.6092	0.0000	1.0000
F81IG	16	-2215.7205	4463.4409	0.0000	1.0000
F81	14	-2226.7655	4481.5310	0.0000	1.0000
JC69I	12	-2229.6764	4483.3529	0.0000	1.0000
JC69G	12	-2229.7607	4483.5215	0.0000	1.0000
JC69IG	13	-2229.7631	4485.5261	0.0000	1.0000
JC69	11	-2241.4806	4504.9612	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2149.1130	4404.0830	0.5887	0.5887
HKYI	16	-2149.5374	4404.9318	0.3851	0.9738
HKYIG	17	-2149.1017	4410.6765	0.0218	0.9956
K2PG	13	-2164.6575	4415.3239	0.0021	0.9977
K2PI	13	-2164.7829	4415.5747	0.0019	0.9996
HKY	15	-2160.1738	4419.5886	0.0003	0.9998
GTRG	20	-2144.7727	4421.8668	0.0001	0.9999
GTRI	20	-2145.3548	4423.0309	0.0000	1.0000
K2PIG	14	-2165.4171	4423.4591	0.0000	1.0000
GTRIG	21	-2144.7692	4428.4757	0.0000	1.0000
K2P	12	-2177.6467	4434.6862	0.0000	1.0000
GTR	19	-2155.5323	4436.7698	0.0000	1.0000
SYMG	17	-2162.7234	4437.9200	0.0000	1.0000
SYMI	17	-2162.9717	4438.4165	0.0000	1.0000
SYMIG	18	-2163.4252	4445.9395	0.0000	1.0000
SYM	16	-2175.1790	4456.2150	0.0000	1.0000
F81I	15	-2215.7282	4530.6974	0.0000	1.0000
F81G	15	-2215.8046	4530.8502	0.0000	1.0000
F81IG	16	-2215.7205	4537.2980	0.0000	1.0000
JC69I	12	-2229.6764	4538.7457	0.0000	1.0000
JC69G	12	-2229.7607	4538.9142	0.0000	1.0000
JC69IG	13	-2229.7631	4545.5349	0.0000	1.0000
F81	14	-2226.7655	4546.1559	0.0000	1.0000
JC69	11	-2241.4806	4555.7380	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/chr8_4947.phylip = 1 - 747;
 Partition Dummy = 1:591-tomodel-phylip/chr8_4947.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/chr8_4947.phylip = 1 - 747;
 Partition Dummy = 1:591-tomodel-phylip/chr8_4947.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/chr8_4947.phylip = 1 - 747;
 Partition Dummy = 1:591-tomodel-phylip/chr8_4947.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

