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

Input data from file "591-tomodel-phylip/chr2_5526.phylip" (7 taxa, 281 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	-833.3595	1700.7795	0.4007	0.4007
HKYG	16	-834.1726	1702.4057	0.1777	0.5784
GTRI	20	-829.6850	1702.6008	0.1612	0.7395
HKYIG	17	-833.3789	1703.0849	0.1265	0.8661
GTRG	20	-830.3421	1703.9150	0.0835	0.9496
GTRIG	21	-829.6998	1704.9672	0.0494	0.9990
HKY	15	-841.0170	1713.8454	0.0006	0.9996
GTR	19	-836.7749	1714.4617	0.0004	1.0000
K2PI	13	-849.2027	1725.7686	0.0000	1.0000
K2PG	13	-850.1743	1727.7120	0.0000	1.0000
K2PIG	14	-849.2284	1728.0357	0.0000	1.0000
F81I	15	-849.0106	1729.8325	0.0000	1.0000
SYMI	17	-846.7822	1729.8915	0.0000	1.0000
F81G	15	-849.8003	1731.4119	0.0000	1.0000
SYMG	17	-847.7343	1731.7955	0.0000	1.0000
F81IG	16	-849.0313	1732.1232	0.0000	1.0000
SYMIG	18	-846.8076	1732.2259	0.0000	1.0000
K2P	12	-858.2585	1741.6811	0.0000	1.0000
F81	14	-856.3677	1742.3144	0.0000	1.0000
SYM	16	-855.7746	1745.6097	0.0000	1.0000
JC69I	12	-862.8777	1750.9197	0.0000	1.0000
JC69G	12	-863.7167	1752.5976	0.0000	1.0000
JC69IG	13	-862.9007	1753.1647	0.0000	1.0000
JC69	11	-871.0870	1765.1554	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-833.3595	1698.7189	0.3159	0.3159
GTRI	20	-829.6850	1699.3700	0.2281	0.5440
HKYG	16	-834.1726	1700.3451	0.1401	0.6842
GTRG	20	-830.3421	1700.6842	0.1183	0.8024
HKYIG	17	-833.3789	1700.7579	0.1140	0.9164
GTRIG	21	-829.6998	1701.3996	0.0827	0.9991
GTR	19	-836.7749	1711.5498	0.0005	0.9996
HKY	15	-841.0170	1712.0340	0.0004	1.0000
K2PI	13	-849.2027	1724.4053	0.0000	1.0000
K2PG	13	-850.1743	1726.3487	0.0000	1.0000
K2PIG	14	-849.2284	1726.4568	0.0000	1.0000
SYMI	17	-846.7822	1727.5645	0.0000	1.0000
F81I	15	-849.0106	1728.0212	0.0000	1.0000
SYMG	17	-847.7343	1729.4685	0.0000	1.0000
F81G	15	-849.8003	1729.6006	0.0000	1.0000
SYMIG	18	-846.8076	1729.6152	0.0000	1.0000
F81IG	16	-849.0313	1730.0626	0.0000	1.0000
K2P	12	-858.2585	1740.5169	0.0000	1.0000
F81	14	-856.3677	1740.7355	0.0000	1.0000
SYM	16	-855.7746	1743.5491	0.0000	1.0000
JC69I	12	-862.8777	1749.7555	0.0000	1.0000
JC69G	12	-863.7167	1751.4334	0.0000	1.0000
JC69IG	13	-862.9007	1751.8014	0.0000	1.0000
JC69	11	-871.0870	1764.1740	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-833.3595	1756.9326	0.6616	0.6616
HKYG	16	-834.1726	1758.5588	0.2934	0.9549
HKYIG	17	-833.3789	1762.6099	0.0387	0.9936
HKY	15	-841.0170	1766.6094	0.0052	0.9989
K2PI	13	-849.2027	1771.7039	0.0004	0.9993
GTRI	20	-829.6850	1772.1371	0.0003	0.9996
GTRG	20	-830.3421	1773.4513	0.0002	0.9998
K2PG	13	-850.1743	1773.6473	0.0002	0.9999
K2PIG	14	-849.2284	1777.3937	0.0000	1.0000
GTRIG	21	-829.6998	1777.8050	0.0000	1.0000
GTR	19	-836.7749	1780.6786	0.0000	1.0000
F81I	15	-849.0106	1782.5965	0.0000	1.0000
F81G	15	-849.8003	1784.1759	0.0000	1.0000
K2P	12	-858.2585	1784.1772	0.0000	1.0000
F81IG	16	-849.0313	1788.2762	0.0000	1.0000
SYMI	17	-846.7822	1789.4165	0.0000	1.0000
SYMG	17	-847.7343	1791.3206	0.0000	1.0000
F81	14	-856.3677	1791.6724	0.0000	1.0000
JC69I	12	-862.8777	1793.4158	0.0000	1.0000
JC69G	12	-863.7167	1795.0937	0.0000	1.0000
SYMIG	18	-846.8076	1795.1056	0.0000	1.0000
JC69IG	13	-862.9007	1799.1001	0.0000	1.0000
SYM	16	-855.7746	1801.7628	0.0000	1.0000
JC69	11	-871.0870	1804.1959	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/chr2_5526.phylip = 1 - 281;
 Partition Dummy = 1:591-tomodel-phylip/chr2_5526.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/chr2_5526.phylip = 1 - 281;
 Partition Dummy = 1:591-tomodel-phylip/chr2_5526.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/chr2_5526.phylip = 1 - 281;
 Partition Dummy = 1:591-tomodel-phylip/chr2_5526.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

