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

Input data from file "591-tomodel-phylip/chr8_7536.phylip" (7 taxa, 327 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	-772.7484	1579.2516	0.3658	0.3658
HKYI	16	-772.7760	1579.3068	0.3558	0.7215
HKYIG	17	-772.7500	1581.4807	0.1200	0.8415
GTRG	20	-769.8945	1582.5340	0.0709	0.9124
GTRI	20	-769.9979	1582.7409	0.0639	0.9763
GTRIG	21	-769.8780	1584.7854	0.0230	0.9993
HKY	15	-780.5044	1592.5521	0.0005	0.9998
GTR	19	-776.7406	1593.9569	0.0002	1.0000
F81G	15	-785.9142	1603.3717	0.0000	1.0000
F81I	15	-785.9910	1603.5255	0.0000	1.0000
F81IG	16	-785.8917	1605.5382	0.0000	1.0000
K2PG	13	-794.3261	1615.8152	0.0000	1.0000
K2PI	13	-794.3500	1615.8630	0.0000	1.0000
F81	14	-793.3210	1615.9882	0.0000	1.0000
K2PIG	14	-794.3381	1618.0223	0.0000	1.0000
SYMG	17	-793.0222	1622.0250	0.0000	1.0000
SYMI	17	-793.1227	1622.2260	0.0000	1.0000
SYMIG	18	-792.9837	1624.1882	0.0000	1.0000
K2P	12	-802.4821	1629.9579	0.0000	1.0000
SYM	16	-800.5810	1634.9168	0.0000	1.0000
JC69G	12	-805.3268	1635.6473	0.0000	1.0000
JC69I	12	-805.4161	1635.8259	0.0000	1.0000
JC69IG	13	-805.2887	1637.7404	0.0000	1.0000
JC69	11	-812.9480	1648.7341	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-772.7484	1577.4967	0.3262	0.3262
HKYI	16	-772.7760	1577.5519	0.3173	0.6435
HKYIG	17	-772.7500	1579.5001	0.1198	0.7633
GTRG	20	-769.8945	1579.7889	0.1037	0.8670
GTRI	20	-769.9979	1579.9958	0.0935	0.9605
GTRIG	21	-769.8780	1581.7559	0.0388	0.9993
HKY	15	-780.5044	1591.0087	0.0004	0.9997
GTR	19	-776.7406	1591.4813	0.0003	1.0000
F81G	15	-785.9142	1601.8283	0.0000	1.0000
F81I	15	-785.9910	1601.9821	0.0000	1.0000
F81IG	16	-785.8917	1603.7834	0.0000	1.0000
F81	14	-793.3210	1614.6420	0.0000	1.0000
K2PG	13	-794.3261	1614.6523	0.0000	1.0000
K2PI	13	-794.3500	1614.7000	0.0000	1.0000
K2PIG	14	-794.3381	1616.6761	0.0000	1.0000
SYMG	17	-793.0222	1620.0444	0.0000	1.0000
SYMI	17	-793.1227	1620.2454	0.0000	1.0000
SYMIG	18	-792.9837	1621.9675	0.0000	1.0000
K2P	12	-802.4821	1628.9643	0.0000	1.0000
SYM	16	-800.5810	1633.1620	0.0000	1.0000
JC69G	12	-805.3268	1634.6537	0.0000	1.0000
JC69I	12	-805.4161	1634.8322	0.0000	1.0000
JC69IG	13	-805.2887	1636.5775	0.0000	1.0000
JC69	11	-812.9480	1647.8960	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-772.7484	1638.1361	0.4911	0.4911
HKYI	16	-772.7760	1638.1913	0.4778	0.9689
HKYIG	17	-772.7500	1643.9294	0.0271	0.9960
HKY	15	-780.5044	1647.8581	0.0038	0.9998
GTRG	20	-769.8945	1655.5881	0.0001	0.9999
GTRI	20	-769.9979	1655.7950	0.0001	1.0000
F81G	15	-785.9142	1658.6777	0.0000	1.0000
F81I	15	-785.9910	1658.8315	0.0000	1.0000
GTRIG	21	-769.8780	1661.3451	0.0000	1.0000
GTR	19	-776.7406	1663.4905	0.0000	1.0000
K2PG	13	-794.3261	1663.9217	0.0000	1.0000
K2PI	13	-794.3500	1663.9695	0.0000	1.0000
F81IG	16	-785.8917	1664.4228	0.0000	1.0000
F81	14	-793.3210	1667.7015	0.0000	1.0000
K2PIG	14	-794.3381	1669.7356	0.0000	1.0000
K2P	12	-802.4821	1674.4438	0.0000	1.0000
JC69G	12	-805.3268	1680.1332	0.0000	1.0000
JC69I	12	-805.4161	1680.3117	0.0000	1.0000
SYMG	17	-793.0222	1684.4737	0.0000	1.0000
SYMI	17	-793.1227	1684.6747	0.0000	1.0000
JC69IG	13	-805.2887	1685.8470	0.0000	1.0000
JC69	11	-812.9480	1689.5856	0.0000	1.0000
SYMIG	18	-792.9837	1690.1867	0.0000	1.0000
SYM	16	-800.5810	1693.8013	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/chr8_7536.phylip = 1 - 327;
 Partition Dummy = 1:591-tomodel-phylip/chr8_7536.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/chr8_7536.phylip = 1 - 327;
 Partition Dummy = 1:591-tomodel-phylip/chr8_7536.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/chr8_7536.phylip = 1 - 327;
 Partition Dummy = 1:591-tomodel-phylip/chr8_7536.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

