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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-844.3973	1722.7168	0.5417	0.5417
HKYIG	17	-844.0849	1724.3399	0.2406	0.7822
HKYI	16	-845.7054	1725.3330	0.1464	0.9287
GTRG	20	-842.4129	1727.8366	0.0419	0.9705
GTRIG	21	-842.1120	1729.5478	0.0178	0.9883
GTRI	20	-843.7314	1730.4735	0.0112	0.9996
K2PG	13	-855.4046	1738.0820	0.0002	0.9998
K2PIG	14	-855.0872	1739.6481	0.0001	0.9999
K2PI	13	-856.7392	1740.7511	0.0001	1.0000
SYMG	17	-854.4363	1745.0427	0.0000	1.0000
HKY	15	-857.0391	1745.7684	0.0000	1.0000
SYMIG	18	-854.1048	1746.6437	0.0000	1.0000
SYMI	17	-855.7309	1747.6320	0.0000	1.0000
GTR	19	-855.3855	1751.4852	0.0000	1.0000
K2P	12	-867.8712	1760.8294	0.0000	1.0000
F81G	15	-868.0997	1767.8896	0.0000	1.0000
SYM	16	-867.0229	1767.9681	0.0000	1.0000
F81IG	16	-867.9375	1769.7972	0.0000	1.0000
F81I	15	-869.1835	1770.0571	0.0000	1.0000
JC69G	12	-876.7288	1778.5448	0.0000	1.0000
JC69IG	13	-876.5510	1780.3748	0.0000	1.0000
JC69I	12	-877.8492	1780.7855	0.0000	1.0000
F81	14	-879.6843	1788.8422	0.0000	1.0000
JC69	11	-888.3593	1799.6353	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-844.3973	1720.7945	0.4978	0.4978
HKYIG	17	-844.0849	1722.1697	0.2503	0.7481
HKYI	16	-845.7054	1723.4108	0.1346	0.8827
GTRG	20	-842.4129	1724.8259	0.0663	0.9490
GTRIG	21	-842.1120	1726.2241	0.0330	0.9819
GTRI	20	-843.7314	1727.4628	0.0177	0.9997
K2PG	13	-855.4046	1736.8093	0.0002	0.9999
K2PIG	14	-855.0872	1738.1744	0.0001	0.9999
K2PI	13	-856.7392	1739.4784	0.0000	1.0000
SYMG	17	-854.4363	1742.8725	0.0000	1.0000
HKY	15	-857.0391	1744.0782	0.0000	1.0000
SYMIG	18	-854.1048	1744.2095	0.0000	1.0000
SYMI	17	-855.7309	1745.4618	0.0000	1.0000
GTR	19	-855.3855	1748.7710	0.0000	1.0000
K2P	12	-867.8712	1759.7423	0.0000	1.0000
SYM	16	-867.0229	1766.0459	0.0000	1.0000
F81G	15	-868.0997	1766.1994	0.0000	1.0000
F81IG	16	-867.9375	1767.8749	0.0000	1.0000
F81I	15	-869.1835	1768.3670	0.0000	1.0000
JC69G	12	-876.7288	1777.4577	0.0000	1.0000
JC69IG	13	-876.5510	1779.1021	0.0000	1.0000
JC69I	12	-877.8492	1779.6984	0.0000	1.0000
F81	14	-879.6843	1787.3686	0.0000	1.0000
JC69	11	-888.3593	1798.7186	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-844.3973	1780.0550	0.6826	0.6826
HKYI	16	-845.7054	1782.6713	0.1845	0.8671
K2PG	13	-855.4046	1784.9585	0.0588	0.9259
HKYIG	17	-844.0849	1785.1340	0.0539	0.9797
K2PI	13	-856.7392	1787.6276	0.0155	0.9952
K2PIG	14	-855.0872	1790.0273	0.0047	0.9999
GTRG	20	-842.4129	1798.9015	0.0001	0.9999
HKY	15	-857.0391	1799.6350	0.0000	1.0000
GTRI	20	-843.7314	1801.5384	0.0000	1.0000
GTRIG	21	-842.1120	1804.0035	0.0000	1.0000
K2P	12	-867.8712	1804.1877	0.0000	1.0000
SYMG	17	-854.4363	1805.8368	0.0000	1.0000
SYMI	17	-855.7309	1808.4261	0.0000	1.0000
SYMIG	18	-854.1048	1810.8776	0.0000	1.0000
GTR	19	-855.3855	1819.1428	0.0000	1.0000
F81G	15	-868.0997	1821.7562	0.0000	1.0000
JC69G	12	-876.7288	1821.9030	0.0000	1.0000
F81I	15	-869.1835	1823.9237	0.0000	1.0000
JC69I	12	-877.8492	1824.1437	0.0000	1.0000
SYM	16	-867.0229	1825.3064	0.0000	1.0000
F81IG	16	-867.9375	1827.1355	0.0000	1.0000
JC69IG	13	-876.5510	1827.2513	0.0000	1.0000
F81	14	-879.6843	1839.2215	0.0000	1.0000
JC69	11	-888.3593	1839.4602	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/chr5_5183.phylip = 1 - 300;
 Partition Dummy = 1:591-tomodel-phylip/chr5_5183.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/chr5_5183.phylip = 1 - 300;
 Partition Dummy = 1:591-tomodel-phylip/chr5_5183.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 (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_5183.phylip = 1 - 300;
 Partition Dummy = 1:591-tomodel-phylip/chr5_5183.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

