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

Input data from file "683-tomodel-phylip/chr2_24734.phylip" (7 taxa, 529 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-1631.4262	3295.9149	0.6532	0.6532
HKYIG	17	-1631.4076	3298.0129	0.2288	0.8820
HKYI	16	-1633.5424	3300.1473	0.0787	0.9607
GTRG	20	-1630.3198	3302.2932	0.0269	0.9876
GTRIG	21	-1630.3023	3304.4271	0.0093	0.9969
GTRI	20	-1632.4717	3306.5969	0.0031	1.0000
HKY	15	-1643.4551	3317.8459	0.0000	1.0000
GTR	19	-1642.6446	3324.7824	0.0000	1.0000
K2PG	13	-1652.8468	3332.4004	0.0000	1.0000
K2PIG	14	-1652.8268	3334.4707	0.0000	1.0000
K2PI	13	-1654.8041	3336.3150	0.0000	1.0000
SYMG	17	-1651.7864	3338.7705	0.0000	1.0000
SYMIG	18	-1651.7598	3340.8609	0.0000	1.0000
SYMI	17	-1653.6348	3342.4674	0.0000	1.0000
K2P	12	-1664.6338	3353.8723	0.0000	1.0000
F81G	15	-1663.9474	3358.8305	0.0000	1.0000
SYM	16	-1663.8761	3360.8147	0.0000	1.0000
F81IG	16	-1663.9412	3360.9449	0.0000	1.0000
F81I	15	-1665.4368	3361.8093	0.0000	1.0000
F81	14	-1674.6001	3378.0173	0.0000	1.0000
JC69G	12	-1681.4570	3387.5187	0.0000	1.0000
JC69IG	13	-1681.4514	3389.6095	0.0000	1.0000
JC69I	12	-1682.8977	3390.4000	0.0000	1.0000
JC69	11	-1692.0827	3406.6760	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1631.4262	3294.8524	0.6338	0.6338
HKYIG	17	-1631.4076	3296.8153	0.2375	0.8713
HKYI	16	-1633.5424	3299.0848	0.0764	0.9477
GTRG	20	-1630.3198	3300.6396	0.0351	0.9828
GTRIG	21	-1630.3023	3302.6046	0.0131	0.9959
GTRI	20	-1632.4717	3304.9434	0.0041	1.0000
HKY	15	-1643.4551	3316.9102	0.0000	1.0000
GTR	19	-1642.6446	3323.2893	0.0000	1.0000
K2PG	13	-1652.8468	3331.6936	0.0000	1.0000
K2PIG	14	-1652.8268	3333.6536	0.0000	1.0000
K2PI	13	-1654.8041	3335.6082	0.0000	1.0000
SYMG	17	-1651.7864	3337.5728	0.0000	1.0000
SYMIG	18	-1651.7598	3339.5197	0.0000	1.0000
SYMI	17	-1653.6348	3341.2697	0.0000	1.0000
K2P	12	-1664.6338	3353.2677	0.0000	1.0000
F81G	15	-1663.9474	3357.8948	0.0000	1.0000
SYM	16	-1663.8761	3359.7522	0.0000	1.0000
F81IG	16	-1663.9412	3359.8824	0.0000	1.0000
F81I	15	-1665.4368	3360.8736	0.0000	1.0000
F81	14	-1674.6001	3377.2001	0.0000	1.0000
JC69G	12	-1681.4570	3386.9141	0.0000	1.0000
JC69IG	13	-1681.4514	3388.9027	0.0000	1.0000
JC69I	12	-1682.8977	3389.7954	0.0000	1.0000
JC69	11	-1692.0827	3406.1653	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1631.4262	3363.1882	0.8584	0.8584
HKYI	16	-1633.5424	3367.4206	0.1034	0.9618
HKYIG	17	-1631.4076	3369.4221	0.0380	0.9999
HKY	15	-1643.4551	3380.9751	0.0001	1.0000
GTRG	20	-1630.3198	3386.0594	0.0000	1.0000
K2PG	13	-1652.8468	3387.2164	0.0000	1.0000
GTRI	20	-1632.4717	3390.3632	0.0000	1.0000
K2PI	13	-1654.8041	3391.1311	0.0000	1.0000
GTRIG	21	-1630.3023	3392.2953	0.0000	1.0000
K2PIG	14	-1652.8268	3393.4474	0.0000	1.0000
GTR	19	-1642.6446	3404.4381	0.0000	1.0000
K2P	12	-1664.6338	3404.5195	0.0000	1.0000
SYMG	17	-1651.7864	3410.1796	0.0000	1.0000
SYMI	17	-1653.6348	3413.8765	0.0000	1.0000
SYMIG	18	-1651.7598	3416.3975	0.0000	1.0000
F81G	15	-1663.9474	3421.9596	0.0000	1.0000
F81I	15	-1665.4368	3424.9384	0.0000	1.0000
SYM	16	-1663.8761	3428.0880	0.0000	1.0000
F81IG	16	-1663.9412	3428.2182	0.0000	1.0000
F81	14	-1674.6001	3436.9940	0.0000	1.0000
JC69G	12	-1681.4570	3438.1659	0.0000	1.0000
JC69I	12	-1682.8977	3441.0472	0.0000	1.0000
JC69IG	13	-1681.4514	3444.4256	0.0000	1.0000
JC69	11	-1692.0827	3453.1462	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_24734.phylip = 1 - 529;
 Partition Dummy = 1:683-tomodel-phylip/chr2_24734.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;


[Mrbayes block for the best AICc model (HKYG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_24734.phylip = 1 - 529;
 Partition Dummy = 1:683-tomodel-phylip/chr2_24734.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;


[Mrbayes block for the best BIC model (HKYG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_24734.phylip = 1 - 529;
 Partition Dummy = 1:683-tomodel-phylip/chr2_24734.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

