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

Input data from file "591-tomodel-phylip/chr2_12995.phylip" (7 taxa, 534 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	-1547.7749	3137.1872	0.5623	0.5623
GTRI	20	-1548.6108	3138.8590	0.2438	0.8061
GTRIG	21	-1547.7585	3139.3216	0.1934	0.9995
HKYG	16	-1559.9843	3153.0208	0.0002	0.9998
HKYI	16	-1560.2943	3153.6408	0.0002	0.9999
HKYIG	17	-1559.9152	3155.0164	0.0001	1.0000
SYMG	17	-1561.9034	3158.9929	0.0000	1.0000
SYMI	17	-1562.4738	3160.1337	0.0000	1.0000
SYMIG	18	-1561.8614	3161.0510	0.0000	1.0000
K2PG	13	-1569.6345	3165.9690	0.0000	1.0000
K2PI	13	-1570.0828	3166.8655	0.0000	1.0000
K2PIG	14	-1569.5847	3167.9787	0.0000	1.0000
GTR	19	-1564.8593	3169.1973	0.0000	1.0000
HKY	15	-1576.4452	3183.8170	0.0000	1.0000
SYM	16	-1577.2847	3187.6215	0.0000	1.0000
K2P	12	-1585.1805	3194.9598	0.0000	1.0000
F81G	15	-1598.6602	3228.2471	0.0000	1.0000
F81I	15	-1599.0706	3229.0678	0.0000	1.0000
F81IG	16	-1598.6365	3230.3252	0.0000	1.0000
JC69G	12	-1607.0093	3238.6175	0.0000	1.0000
JC69I	12	-1607.4492	3239.4973	0.0000	1.0000
JC69IG	13	-1606.9838	3240.6676	0.0000	1.0000
F81	14	-1612.2413	3253.2918	0.0000	1.0000
JC69	11	-1620.4603	3263.4264	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1547.7749	3135.5498	0.5531	0.5531
GTRI	20	-1548.6108	3137.2216	0.2398	0.7928
GTRIG	21	-1547.7585	3137.5169	0.2068	0.9997
HKYG	16	-1559.9843	3151.9686	0.0002	0.9998
HKYI	16	-1560.2943	3152.5886	0.0001	0.9999
HKYIG	17	-1559.9152	3153.8304	0.0001	1.0000
SYMG	17	-1561.9034	3157.8069	0.0000	1.0000
SYMI	17	-1562.4738	3158.9477	0.0000	1.0000
SYMIG	18	-1561.8614	3159.7228	0.0000	1.0000
K2PG	13	-1569.6345	3165.2690	0.0000	1.0000
K2PI	13	-1570.0828	3166.1655	0.0000	1.0000
K2PIG	14	-1569.5847	3167.1695	0.0000	1.0000
GTR	19	-1564.8593	3167.7187	0.0000	1.0000
HKY	15	-1576.4452	3182.8903	0.0000	1.0000
SYM	16	-1577.2847	3186.5693	0.0000	1.0000
K2P	12	-1585.1805	3194.3609	0.0000	1.0000
F81G	15	-1598.6602	3227.3205	0.0000	1.0000
F81I	15	-1599.0706	3228.1411	0.0000	1.0000
F81IG	16	-1598.6365	3229.2730	0.0000	1.0000
JC69G	12	-1607.0093	3238.0186	0.0000	1.0000
JC69I	12	-1607.4492	3238.8985	0.0000	1.0000
JC69IG	13	-1606.9838	3239.9676	0.0000	1.0000
F81	14	-1612.2413	3252.4825	0.0000	1.0000
JC69	11	-1620.4603	3262.9206	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1559.9843	3220.4549	0.2399	0.2399
K2PG	13	-1569.6345	3220.9141	0.1907	0.4306
HKYI	16	-1560.2943	3221.0750	0.1759	0.6065
GTRG	20	-1547.7749	3221.1577	0.1688	0.7753
K2PI	13	-1570.0828	3221.8107	0.1218	0.8971
GTRI	20	-1548.6108	3222.8295	0.0732	0.9703
HKYIG	17	-1559.9152	3226.5971	0.0111	0.9814
K2PIG	14	-1569.5847	3227.0950	0.0087	0.9901
GTRIG	21	-1547.7585	3227.4053	0.0074	0.9975
SYMG	17	-1561.9034	3230.5736	0.0015	0.9991
SYMI	17	-1562.4738	3231.7144	0.0009	0.9999
SYMIG	18	-1561.8614	3236.7700	0.0001	1.0000
K2P	12	-1585.1805	3245.7257	0.0000	1.0000
HKY	15	-1576.4452	3247.0963	0.0000	1.0000
GTR	19	-1564.8593	3249.0462	0.0000	1.0000
SYM	16	-1577.2847	3255.0556	0.0000	1.0000
JC69G	12	-1607.0093	3289.3834	0.0000	1.0000
JC69I	12	-1607.4492	3290.2632	0.0000	1.0000
F81G	15	-1598.6602	3291.5264	0.0000	1.0000
F81I	15	-1599.0706	3292.3471	0.0000	1.0000
JC69IG	13	-1606.9838	3295.6128	0.0000	1.0000
F81IG	16	-1598.6365	3297.7594	0.0000	1.0000
JC69	11	-1620.4603	3310.0050	0.0000	1.0000
F81	14	-1612.2413	3312.4081	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/chr2_12995.phylip = 1 - 534;
 Partition Dummy = 1:591-tomodel-phylip/chr2_12995.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/chr2_12995.phylip = 1 - 534;
 Partition Dummy = 1:591-tomodel-phylip/chr2_12995.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/chr2_12995.phylip = 1 - 534;
 Partition Dummy = 1:591-tomodel-phylip/chr2_12995.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

