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

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

Minimum AIC  model: HKY
Minimum AICc model: HKY
Minimum BIC  model: HKY


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
HKY	15	-757.7795	1547.5262	0.4124	0.4124
HKYI	16	-757.0202	1548.2791	0.2830	0.6955
HKYG	16	-757.3536	1548.9458	0.2028	0.8983
HKYIG	17	-757.0326	1550.5942	0.0889	0.9872
GTR	19	-757.4871	1556.1409	0.0056	0.9928
GTRI	20	-756.7588	1557.0322	0.0036	0.9963
GTRG	20	-757.0810	1557.6767	0.0026	0.9989
GTRIG	21	-756.7709	1559.4242	0.0011	1.0000
K2P	12	-774.1440	1573.5512	0.0000	1.0000
K2PI	13	-773.4473	1574.3743	0.0000	1.0000
K2PG	13	-773.7322	1574.9440	0.0000	1.0000
K2PIG	14	-773.4584	1576.6310	0.0000	1.0000
SYM	16	-772.8173	1579.8732	0.0000	1.0000
SYMI	17	-772.0745	1580.6778	0.0000	1.0000
F81	14	-775.6432	1581.0007	0.0000	1.0000
SYMG	17	-772.3261	1581.1810	0.0000	1.0000
F81I	15	-774.9373	1581.8417	0.0000	1.0000
F81G	15	-775.2569	1582.4811	0.0000	1.0000
SYMIG	18	-772.0841	1583.0064	0.0000	1.0000
F81IG	16	-774.9495	1584.1377	0.0000	1.0000
JC69	11	-790.4167	1603.8980	0.0000	1.0000
JC69I	12	-789.7874	1604.8379	0.0000	1.0000
JC69G	12	-790.0695	1605.4021	0.0000	1.0000
JC69IG	13	-789.7989	1607.0775	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKY	15	-757.7795	1545.5590	0.3706	0.3706
HKYI	16	-757.0202	1546.0404	0.2913	0.6619
HKYG	16	-757.3536	1546.7071	0.2087	0.8706
HKYIG	17	-757.0326	1548.0653	0.1058	0.9764
GTR	19	-757.4871	1552.9743	0.0091	0.9855
GTRI	20	-756.7588	1553.5176	0.0069	0.9925
GTRG	20	-757.0810	1554.1621	0.0050	0.9975
GTRIG	21	-756.7709	1555.5418	0.0025	1.0000
K2P	12	-774.1440	1572.2880	0.0000	1.0000
K2PI	13	-773.4473	1572.8946	0.0000	1.0000
K2PG	13	-773.7322	1573.4643	0.0000	1.0000
K2PIG	14	-773.4584	1574.9167	0.0000	1.0000
SYM	16	-772.8173	1577.6346	0.0000	1.0000
SYMI	17	-772.0745	1578.1489	0.0000	1.0000
SYMG	17	-772.3261	1578.6521	0.0000	1.0000
F81	14	-775.6432	1579.2864	0.0000	1.0000
F81I	15	-774.9373	1579.8745	0.0000	1.0000
SYMIG	18	-772.0841	1580.1682	0.0000	1.0000
F81G	15	-775.2569	1580.5138	0.0000	1.0000
F81IG	16	-774.9495	1581.8991	0.0000	1.0000
JC69	11	-790.4167	1602.8335	0.0000	1.0000
JC69I	12	-789.7874	1603.5748	0.0000	1.0000
JC69G	12	-790.0695	1604.1390	0.0000	1.0000
JC69IG	13	-789.7989	1605.5978	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	15	-757.7795	1598.9692	0.8091	0.8091
HKYI	16	-757.0202	1603.0113	0.1072	0.9163
HKYG	16	-757.3536	1603.6780	0.0768	0.9931
HKYIG	17	-757.0326	1608.5968	0.0066	0.9997
K2P	12	-774.1440	1615.0162	0.0003	0.9999
K2PI	13	-773.4473	1619.1835	0.0000	1.0000
K2PG	13	-773.7322	1619.7532	0.0000	1.0000
GTR	19	-757.4871	1620.6272	0.0000	1.0000
GTRI	20	-756.7588	1624.7312	0.0000	1.0000
K2PIG	14	-773.4584	1624.7663	0.0000	1.0000
GTRG	20	-757.0810	1625.3757	0.0000	1.0000
F81	14	-775.6432	1629.1360	0.0000	1.0000
GTRIG	21	-756.7709	1630.3161	0.0000	1.0000
F81I	15	-774.9373	1633.2847	0.0000	1.0000
F81G	15	-775.2569	1633.9241	0.0000	1.0000
SYM	16	-772.8173	1634.6055	0.0000	1.0000
SYMI	17	-772.0745	1638.6805	0.0000	1.0000
F81IG	16	-774.9495	1638.8700	0.0000	1.0000
SYMG	17	-772.3261	1639.1837	0.0000	1.0000
JC69	11	-790.4167	1642.0010	0.0000	1.0000
SYMIG	18	-772.0841	1644.2605	0.0000	1.0000
JC69I	12	-789.7874	1646.3030	0.0000	1.0000
JC69G	12	-790.0695	1646.8671	0.0000	1.0000
JC69IG	13	-789.7989	1651.8866	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_17991.phylip = 1 - 260;
 Partition Dummy = 1:591-tomodel-phylip/chr3_17991.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_17991.phylip = 1 - 260;
 Partition Dummy = 1:591-tomodel-phylip/chr3_17991.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_17991.phylip = 1 - 260;
 Partition Dummy = 1:591-tomodel-phylip/chr3_17991.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


-------------------------------------------------------------
End of Output

