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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	16	-754.9029	1543.3780	0.5526	0.5526
HKYIG	17	-754.7457	1545.2653	0.2151	0.7677
HKYG	16	-756.1295	1545.8312	0.1621	0.9297
GTRI	20	-752.9509	1548.3580	0.0458	0.9756
GTRIG	21	-752.9509	1550.6114	0.0148	0.9904
GTRG	20	-754.7039	1551.8639	0.0079	0.9984
K2PI	13	-764.7627	1556.5684	0.0008	0.9991
K2PIG	14	-764.4937	1558.1943	0.0003	0.9994
K2PG	13	-765.6712	1558.3854	0.0003	0.9997
F81I	15	-764.8177	1561.0186	0.0001	0.9998
SYMI	17	-762.9430	1561.6598	0.0001	0.9999
F81G	15	-765.6307	1562.6447	0.0000	0.9999
F81IG	16	-764.5792	1562.7307	0.0000	1.0000
SYMIG	18	-762.8165	1563.6214	0.0000	1.0000
SYMG	17	-764.1907	1564.1553	0.0000	1.0000
HKY	15	-768.5350	1568.4534	0.0000	1.0000
JC69I	12	-773.4420	1571.7754	0.0000	1.0000
JC69G	12	-774.1541	1573.1996	0.0000	1.0000
JC69IG	13	-773.1972	1573.4374	0.0000	1.0000
GTR	19	-767.5947	1575.4052	0.0000	1.0000
K2P	12	-777.4215	1579.7344	0.0000	1.0000
F81	14	-777.0101	1583.2271	0.0000	1.0000
SYM	16	-776.2937	1586.1597	0.0000	1.0000
JC69	11	-785.3109	1593.3739	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-754.9029	1541.8058	0.5195	0.5195
HKYIG	17	-754.7457	1543.4913	0.2236	0.7431
HKYG	16	-756.1295	1544.2589	0.1524	0.8955
GTRI	20	-752.9509	1545.9018	0.0670	0.9625
GTRIG	21	-752.9509	1547.9018	0.0247	0.9871
GTRG	20	-754.7039	1549.4078	0.0116	0.9987
K2PI	13	-764.7627	1555.5254	0.0005	0.9993
K2PIG	14	-764.4937	1556.9874	0.0003	0.9995
K2PG	13	-765.6712	1557.3425	0.0002	0.9998
F81I	15	-764.8177	1559.6353	0.0001	0.9998
SYMI	17	-762.9430	1559.8859	0.0001	0.9999
F81IG	16	-764.5792	1561.1584	0.0000	0.9999
F81G	15	-765.6307	1561.2615	0.0000	1.0000
SYMIG	18	-762.8165	1561.6330	0.0000	1.0000
SYMG	17	-764.1907	1562.3814	0.0000	1.0000
HKY	15	-768.5350	1567.0701	0.0000	1.0000
JC69I	12	-773.4420	1570.8840	0.0000	1.0000
JC69G	12	-774.1541	1572.3082	0.0000	1.0000
JC69IG	13	-773.1972	1572.3945	0.0000	1.0000
GTR	19	-767.5947	1573.1894	0.0000	1.0000
K2P	12	-777.4215	1578.8430	0.0000	1.0000
F81	14	-777.0101	1582.0202	0.0000	1.0000
SYM	16	-776.2937	1584.5874	0.0000	1.0000
JC69	11	-785.3109	1592.6218	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-754.9029	1604.1162	0.5292	0.5292
K2PI	13	-764.7627	1606.1526	0.1912	0.7203
HKYG	16	-756.1295	1606.5694	0.1552	0.8755
K2PG	13	-765.6712	1607.9697	0.0771	0.9526
HKYIG	17	-754.7457	1609.6962	0.0325	0.9851
K2PIG	14	-764.4937	1611.5090	0.0131	0.9982
JC69I	12	-773.4420	1617.6169	0.0006	0.9988
F81I	15	-764.8177	1618.0514	0.0005	0.9993
JC69G	12	-774.1541	1619.0411	0.0003	0.9996
F81G	15	-765.6307	1619.6775	0.0002	0.9999
JC69IG	13	-773.1972	1623.0217	0.0000	0.9999
F81IG	16	-764.5792	1623.4688	0.0000	0.9999
GTRI	20	-752.9509	1623.7899	0.0000	1.0000
HKY	15	-768.5350	1625.4861	0.0000	1.0000
K2P	12	-777.4215	1625.5758	0.0000	1.0000
SYMI	17	-762.9430	1626.0908	0.0000	1.0000
GTRG	20	-754.7039	1627.2959	0.0000	1.0000
SYMG	17	-764.1907	1628.5863	0.0000	1.0000
GTRIG	21	-752.9509	1629.6842	0.0000	1.0000
SYMIG	18	-762.8165	1631.7323	0.0000	1.0000
JC69	11	-785.3109	1635.4602	0.0000	1.0000
F81	14	-777.0101	1636.5418	0.0000	1.0000
SYM	16	-776.2937	1646.8979	0.0000	1.0000
GTR	19	-767.5947	1647.1831	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_3180.phylip = 1 - 363;
 Partition Dummy = 1:591-tomodel-phylip/chr3_3180.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_3180.phylip = 1 - 363;
 Partition Dummy = 1:591-tomodel-phylip/chr3_3180.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_3180.phylip = 1 - 363;
 Partition Dummy = 1:591-tomodel-phylip/chr3_3180.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

