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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-757.3553	1557.1524	0.3948	0.3948
GTRI	20	-757.4717	1557.3853	0.3514	0.7463
GTRIG	21	-757.0857	1558.8652	0.1677	0.9140
HKYG	16	-764.7625	1563.0882	0.0203	0.9343
HKYI	16	-764.8448	1563.2529	0.0187	0.9530
SYMG	17	-764.1663	1564.0963	0.0123	0.9652
SYMI	17	-764.2593	1564.2823	0.0112	0.9764
HKYIG	17	-764.5709	1564.9055	0.0082	0.9846
SYMIG	18	-763.9836	1565.9441	0.0049	0.9895
K2PG	13	-769.6610	1566.3590	0.0040	0.9934
K2PI	13	-769.7376	1566.5122	0.0037	0.9971
K2PIG	14	-769.4926	1568.1853	0.0016	0.9987
GTR	19	-764.3141	1568.8312	0.0011	0.9998
HKY	15	-771.2268	1573.8290	0.0001	0.9999
SYM	16	-770.5735	1574.7101	0.0001	1.0000
K2P	12	-775.9879	1576.8621	0.0000	1.0000
F81G	15	-774.6659	1580.7071	0.0000	1.0000
F81I	15	-774.7738	1580.9229	0.0000	1.0000
JC69G	12	-778.8207	1582.5277	0.0000	1.0000
F81IG	16	-774.5430	1582.6492	0.0000	1.0000
JC69I	12	-778.9279	1582.7422	0.0000	1.0000
JC69IG	13	-778.7021	1584.4413	0.0000	1.0000
F81	14	-780.4217	1590.0434	0.0000	1.0000
JC69	11	-784.5461	1591.8401	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-757.3553	1554.7106	0.3973	0.3973
GTRI	20	-757.4717	1554.9434	0.3536	0.7510
GTRIG	21	-757.0857	1556.1714	0.1914	0.9424
HKYG	16	-764.7625	1561.5250	0.0132	0.9555
HKYI	16	-764.8448	1561.6896	0.0121	0.9677
SYMG	17	-764.1663	1562.3326	0.0088	0.9764
SYMI	17	-764.2593	1562.5186	0.0080	0.9845
HKYIG	17	-764.5709	1563.1418	0.0059	0.9903
SYMIG	18	-763.9836	1563.9672	0.0039	0.9942
K2PG	13	-769.6610	1565.3220	0.0020	0.9962
K2PI	13	-769.7376	1565.4752	0.0018	0.9980
GTR	19	-764.3141	1566.6283	0.0010	0.9990
K2PIG	14	-769.4926	1566.9853	0.0009	0.9999
HKY	15	-771.2268	1572.4536	0.0001	0.9999
SYM	16	-770.5735	1573.1469	0.0000	1.0000
K2P	12	-775.9879	1575.9757	0.0000	1.0000
F81G	15	-774.6659	1579.3318	0.0000	1.0000
F81I	15	-774.7738	1579.5476	0.0000	1.0000
F81IG	16	-774.5430	1581.0860	0.0000	1.0000
JC69G	12	-778.8207	1581.6413	0.0000	1.0000
JC69I	12	-778.9279	1581.8558	0.0000	1.0000
JC69IG	13	-778.7021	1583.4042	0.0000	1.0000
F81	14	-780.4217	1588.8434	0.0000	1.0000
JC69	11	-784.5461	1591.0922	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-769.6610	1616.0206	0.4833	0.4833
K2PI	13	-769.7376	1616.1738	0.4476	0.9309
K2PIG	14	-769.4926	1621.5838	0.0299	0.9608
K2P	12	-775.9879	1622.7745	0.0165	0.9773
HKYG	16	-764.7625	1623.9234	0.0093	0.9866
HKYI	16	-764.8448	1624.0880	0.0086	0.9952
JC69G	12	-778.8207	1628.4401	0.0010	0.9961
SYMG	17	-764.1663	1628.6308	0.0009	0.9970
JC69I	12	-778.9279	1628.6546	0.0009	0.9979
SYMI	17	-764.2593	1628.8168	0.0008	0.9987
HKYIG	17	-764.5709	1629.4401	0.0006	0.9993
HKY	15	-771.2268	1630.9521	0.0003	0.9996
GTRG	20	-757.3553	1632.7085	0.0001	0.9997
GTRI	20	-757.4717	1632.9414	0.0001	0.9998
JC69	11	-784.5461	1633.9911	0.0001	0.9998
JC69IG	13	-778.7021	1634.1029	0.0001	0.9999
SYMIG	18	-763.9836	1634.1654	0.0001	0.9999
SYM	16	-770.5735	1635.5453	0.0000	1.0000
F81G	15	-774.6659	1637.8302	0.0000	1.0000
F81I	15	-774.7738	1638.0460	0.0000	1.0000
GTRIG	21	-757.0857	1638.0692	0.0000	1.0000
GTR	19	-764.3141	1640.7263	0.0000	1.0000
F81	14	-780.4217	1643.4420	0.0000	1.0000
F81IG	16	-774.5430	1643.4843	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_18608.phylip = 1 - 365;
 Partition Dummy = 1:591-tomodel-phylip/chr2_18608.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best AICc model (GTRG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_18608.phylip = 1 - 365;
 Partition Dummy = 1:591-tomodel-phylip/chr2_18608.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,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/chr2_18608.phylip = 1 - 365;
 Partition Dummy = 1:591-tomodel-phylip/chr2_18608.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

