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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: K2PG


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	-854.1635	1741.6192	0.2554	0.2554
HKYI	16	-854.2150	1741.7221	0.2426	0.4981
GTRG	20	-850.3572	1742.7288	0.1467	0.6447
GTRI	20	-850.4137	1742.8418	0.1386	0.7833
HKYIG	17	-854.0362	1743.5295	0.0983	0.8816
GTRIG	21	-850.3000	1744.8212	0.0515	0.9331
HKY	15	-857.5093	1746.1560	0.0264	0.9596
GTR	19	-853.4898	1746.7978	0.0192	0.9787
K2PG	13	-860.9238	1748.7061	0.0074	0.9861
K2PI	13	-860.9805	1748.8196	0.0070	0.9931
K2PIG	14	-860.7948	1750.5825	0.0029	0.9960
SYMG	17	-858.3529	1752.1629	0.0013	0.9973
SYMI	17	-858.4128	1752.2827	0.0012	0.9985
K2P	12	-864.2078	1753.1498	0.0008	0.9993
SYMIG	18	-858.2707	1754.1738	0.0005	0.9998
SYM	16	-861.4754	1756.2429	0.0002	1.0000
F81G	15	-865.9581	1763.0537	0.0000	1.0000
F81I	15	-866.0272	1763.1919	0.0000	1.0000
F81IG	16	-865.8594	1765.0110	0.0000	1.0000
F81	14	-869.0560	1767.1048	0.0000	1.0000
JC69G	12	-872.2577	1769.2496	0.0000	1.0000
JC69I	12	-872.3257	1769.3856	0.0000	1.0000
JC69IG	13	-872.1550	1771.1684	0.0000	1.0000
JC69	11	-875.3826	1773.3849	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-854.1635	1740.3271	0.2195	0.2195
HKYI	16	-854.2150	1740.4299	0.2085	0.4280
GTRG	20	-850.3572	1740.7145	0.1808	0.6088
GTRI	20	-850.4137	1740.8274	0.1709	0.7798
HKYIG	17	-854.0362	1742.0724	0.0917	0.8715
GTRIG	21	-850.3000	1742.6000	0.0704	0.9419
GTR	19	-853.4898	1744.9796	0.0214	0.9634
HKY	15	-857.5093	1745.0185	0.0210	0.9844
K2PG	13	-860.9238	1747.8476	0.0051	0.9895
K2PI	13	-860.9805	1747.9611	0.0048	0.9943
K2PIG	14	-860.7948	1749.5896	0.0021	0.9965
SYMG	17	-858.3529	1750.7058	0.0012	0.9977
SYMI	17	-858.4128	1750.8256	0.0012	0.9988
K2P	12	-864.2078	1752.4157	0.0005	0.9994
SYMIG	18	-858.2707	1752.5413	0.0005	0.9998
SYM	16	-861.4754	1754.9508	0.0001	1.0000
F81G	15	-865.9581	1761.9162	0.0000	1.0000
F81I	15	-866.0272	1762.0545	0.0000	1.0000
F81IG	16	-865.8594	1763.7188	0.0000	1.0000
F81	14	-869.0560	1766.1119	0.0000	1.0000
JC69G	12	-872.2577	1768.5155	0.0000	1.0000
JC69I	12	-872.3257	1768.6515	0.0000	1.0000
JC69IG	13	-872.1550	1770.3099	0.0000	1.0000
JC69	11	-875.3826	1772.7652	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-860.9238	1800.9164	0.3287	0.3287
K2PI	13	-860.9805	1801.0299	0.3105	0.6392
K2P	12	-864.2078	1801.4023	0.2578	0.8970
HKYG	16	-854.1635	1805.6426	0.0309	0.9279
HKYI	16	-854.2150	1805.7454	0.0294	0.9573
HKY	15	-857.5093	1806.2518	0.0228	0.9801
K2PIG	14	-860.7948	1806.7406	0.0179	0.9980
HKYIG	17	-854.0362	1811.4701	0.0017	0.9997
JC69G	12	-872.2577	1817.5021	0.0001	0.9997
JC69I	12	-872.3257	1817.6381	0.0001	0.9998
JC69	11	-875.3826	1817.6696	0.0001	0.9999
SYMG	17	-858.3529	1820.1035	0.0000	0.9999
SYMI	17	-858.4128	1820.2233	0.0000	0.9999
SYM	16	-861.4754	1820.2663	0.0000	1.0000
GTRG	20	-850.3572	1822.3588	0.0000	1.0000
GTRI	20	-850.4137	1822.4718	0.0000	1.0000
GTR	19	-853.4898	1822.5418	0.0000	1.0000
F81G	15	-865.9581	1823.1495	0.0000	1.0000
F81	14	-869.0560	1823.2630	0.0000	1.0000
F81I	15	-866.0272	1823.2877	0.0000	1.0000
JC69IG	13	-872.1550	1823.3788	0.0000	1.0000
SYMIG	18	-858.2707	1826.0213	0.0000	1.0000
GTRIG	21	-850.3000	1828.3266	0.0000	1.0000
F81IG	16	-865.8594	1829.0343	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 591-tomodel-phylip/chr3_6231.phylip = 1 - 438;
 Partition Dummy = 1:591-tomodel-phylip/chr3_6231.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 591-tomodel-phylip/chr3_6231.phylip = 1 - 438;
 Partition Dummy = 1:591-tomodel-phylip/chr3_6231.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 (K2PG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_6231.phylip = 1 - 438;
 Partition Dummy = 1:591-tomodel-phylip/chr3_6231.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

