
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:40 02/17/2011

Input data from file "183-phylip/chr13_714.phylip" (29 taxa, 226 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	60	-725.6751	1615.7138	0.5005	0.5005
HKYI	60	-726.0547	1616.4730	0.3424	0.8430
HKYIG	61	-725.6704	1619.4628	0.0768	0.9198
HKY	59	-729.7139	1620.0784	0.0565	0.9762
GTRG	64	-722.2332	1624.1435	0.0074	0.9836
GTRI	64	-722.5256	1624.7283	0.0055	0.9892
K2PG	57	-735.8355	1625.0282	0.0048	0.9939
K2PI	57	-736.3603	1626.0778	0.0028	0.9967
GTR	63	-726.1136	1628.0049	0.0011	0.9978
GTRIG	65	-722.2278	1628.0806	0.0010	0.9988
K2PIG	58	-735.8277	1628.6374	0.0008	0.9996
K2P	56	-740.1609	1630.0970	0.0004	1.0000
SYMG	61	-734.6723	1637.4665	0.0000	1.0000
SYMI	61	-735.1557	1638.4334	0.0000	1.0000
SYMIG	62	-734.6494	1641.2252	0.0000	1.0000
SYM	60	-738.4921	1641.3479	0.0000	1.0000
F81G	59	-741.6465	1643.9436	0.0000	1.0000
F81I	59	-742.1315	1644.9136	0.0000	1.0000
F81IG	60	-741.6427	1647.6490	0.0000	1.0000
F81	58	-745.6226	1648.2272	0.0000	1.0000
JC69G	56	-749.7750	1649.3251	0.0000	1.0000
JC69I	56	-750.3226	1650.4203	0.0000	1.0000
JC69IG	57	-749.7701	1652.8973	0.0000	1.0000
JC69	55	-753.9697	1654.1748	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-725.6751	1571.3502	0.2988	0.2988
HKYI	60	-726.0547	1572.1094	0.2044	0.5033
GTRG	64	-722.2332	1572.4665	0.1710	0.6743
GTRI	64	-722.5256	1573.0513	0.1276	0.8019
HKYIG	61	-725.6704	1573.3408	0.1104	0.9123
GTRIG	65	-722.2278	1574.4556	0.0633	0.9756
HKY	59	-729.7139	1577.4278	0.0143	0.9899
GTR	63	-726.1136	1578.2271	0.0096	0.9995
K2PG	57	-735.8355	1585.6710	0.0002	0.9997
K2PI	57	-736.3603	1586.7207	0.0001	0.9999
K2PIG	58	-735.8277	1587.6553	0.0001	1.0000
SYMG	61	-734.6723	1591.3445	0.0000	1.0000
SYMI	61	-735.1557	1592.3114	0.0000	1.0000
K2P	56	-740.1609	1592.3219	0.0000	1.0000
SYMIG	62	-734.6494	1593.2988	0.0000	1.0000
SYM	60	-738.4921	1596.9843	0.0000	1.0000
F81G	59	-741.6465	1601.2930	0.0000	1.0000
F81I	59	-742.1315	1602.2630	0.0000	1.0000
F81IG	60	-741.6427	1603.2854	0.0000	1.0000
F81	58	-745.6226	1607.2451	0.0000	1.0000
JC69G	56	-749.7750	1611.5499	0.0000	1.0000
JC69I	56	-750.3226	1612.6452	0.0000	1.0000
JC69IG	57	-749.7701	1613.5402	0.0000	1.0000
JC69	55	-753.9697	1617.9395	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-725.6751	1776.5823	0.4422	0.4422
HKYI	60	-726.0547	1777.3415	0.3026	0.7448
HKY	59	-729.7139	1779.2394	0.1171	0.8619
K2PG	57	-735.8355	1780.6415	0.0581	0.9200
K2PI	57	-736.3603	1781.6912	0.0344	0.9544
HKYIG	61	-725.6704	1781.9934	0.0296	0.9840
K2P	56	-740.1609	1783.8718	0.0116	0.9955
K2PIG	58	-735.8277	1786.0463	0.0039	0.9994
GTRG	64	-722.2332	1791.3807	0.0003	0.9997
GTRI	64	-722.5256	1791.9655	0.0002	0.9999
GTR	63	-726.1136	1793.7208	0.0001	1.0000
GTRIG	65	-722.2278	1796.7904	0.0000	1.0000
SYMG	61	-734.6723	1799.9972	0.0000	1.0000
SYMI	61	-735.1557	1800.9641	0.0000	1.0000
SYM	60	-738.4921	1802.2164	0.0000	1.0000
JC69G	56	-749.7750	1803.0999	0.0000	1.0000
F81G	59	-741.6465	1803.1046	0.0000	1.0000
F81I	59	-742.1315	1804.0746	0.0000	1.0000
JC69I	56	-750.3226	1804.1951	0.0000	1.0000
SYMIG	62	-734.6494	1805.3720	0.0000	1.0000
F81	58	-745.6226	1805.6362	0.0000	1.0000
JC69	55	-753.9697	1806.0689	0.0000	1.0000
JC69IG	57	-749.7701	1808.5107	0.0000	1.0000
F81IG	60	-741.6427	1808.5175	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 183-phylip/chr13_714.phylip = 1 - 226;
 Partition Dummy = 1:183-phylip/chr13_714.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 183-phylip/chr13_714.phylip = 1 - 226;
 Partition Dummy = 1:183-phylip/chr13_714.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 183-phylip/chr13_714.phylip = 1 - 226;
 Partition Dummy = 1:183-phylip/chr13_714.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

