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

Input data from file "183-phylip/chr11_3432.phylip" (29 taxa, 266 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	-741.5528	1638.8128	0.5762	0.5762
HKYI	60	-742.5410	1640.7893	0.2145	0.7906
HKYIG	61	-741.1306	1641.3397	0.1629	0.9535
GTRG	64	-738.2205	1645.8340	0.0172	0.9707
F81G	59	-747.2922	1646.9534	0.0098	0.9806
GTRI	64	-739.0715	1647.5360	0.0074	0.9879
GTRIG	65	-737.7170	1648.3340	0.0049	0.9929
F81I	59	-748.1015	1648.5719	0.0044	0.9972
F81IG	60	-746.8969	1649.5012	0.0028	1.0000
HKY	59	-755.1575	1662.6839	0.0000	1.0000
K2PG	57	-760.8886	1667.5656	0.0000	1.0000
K2PI	57	-761.4487	1668.6859	0.0000	1.0000
K2PIG	58	-760.6053	1670.2735	0.0000	1.0000
GTR	63	-752.2170	1670.3548	0.0000	1.0000
F81	58	-761.0716	1671.2059	0.0000	1.0000
JC69G	56	-767.5026	1677.5507	0.0000	1.0000
SYMG	61	-759.2618	1677.6020	0.0000	1.0000
JC69I	56	-768.0237	1678.5929	0.0000	1.0000
SYMI	61	-759.8196	1678.7176	0.0000	1.0000
JC69IG	57	-767.2455	1680.2796	0.0000	1.0000
SYMIG	62	-758.9952	1680.4732	0.0000	1.0000
K2P	56	-774.3666	1691.2787	0.0000	1.0000
JC69	55	-780.7892	1700.9117	0.0000	1.0000
SYM	60	-772.6316	1700.9705	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-741.5528	1603.1055	0.3339	0.3339
HKYIG	61	-741.1306	1604.2612	0.1874	0.5213
GTRG	64	-738.2205	1604.4410	0.1713	0.6925
HKYI	60	-742.5410	1605.0820	0.1243	0.8168
GTRIG	65	-737.7170	1605.4340	0.1042	0.9211
GTRI	64	-739.0715	1606.1430	0.0731	0.9942
F81G	59	-747.2922	1612.5845	0.0029	0.9971
F81IG	60	-746.8969	1613.7938	0.0016	0.9987
F81I	59	-748.1015	1614.2030	0.0013	1.0000
HKY	59	-755.1575	1628.3150	0.0000	1.0000
GTR	63	-752.2170	1630.4340	0.0000	1.0000
K2PG	57	-760.8886	1635.7771	0.0000	1.0000
K2PI	57	-761.4487	1636.8974	0.0000	1.0000
K2PIG	58	-760.6053	1637.2107	0.0000	1.0000
F81	58	-761.0716	1638.1431	0.0000	1.0000
SYMG	61	-759.2618	1640.5236	0.0000	1.0000
SYMI	61	-759.8196	1641.6391	0.0000	1.0000
SYMIG	62	-758.9952	1641.9904	0.0000	1.0000
JC69G	56	-767.5026	1647.0052	0.0000	1.0000
JC69I	56	-768.0237	1648.0474	0.0000	1.0000
JC69IG	57	-767.2455	1648.4911	0.0000	1.0000
K2P	56	-774.3666	1660.7332	0.0000	1.0000
SYM	60	-772.6316	1665.2632	0.0000	1.0000
JC69	55	-780.7892	1671.5783	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-741.5528	1818.1153	0.6464	0.6464
HKYI	60	-742.5410	1820.0918	0.2406	0.8870
HKYIG	61	-741.1306	1822.8545	0.0605	0.9475
F81G	59	-747.2922	1824.0107	0.0339	0.9814
F81I	59	-748.1015	1825.6293	0.0151	0.9965
F81IG	60	-746.8969	1828.8036	0.0031	0.9996
GTRG	64	-738.2205	1833.7847	0.0003	0.9998
GTRI	64	-739.0715	1835.4868	0.0001	0.9999
GTRIG	65	-737.7170	1838.3613	0.0000	1.0000
HKY	59	-755.1575	1839.7412	0.0000	1.0000
K2PG	57	-760.8886	1840.0364	0.0000	1.0000
K2PI	57	-761.4487	1841.1567	0.0000	1.0000
K2PIG	58	-760.6053	1845.0535	0.0000	1.0000
F81	58	-761.0716	1845.9859	0.0000	1.0000
JC69G	56	-767.5026	1847.6810	0.0000	1.0000
JC69I	56	-768.0237	1848.7232	0.0000	1.0000
JC69IG	57	-767.2455	1852.7504	0.0000	1.0000
GTR	63	-752.2170	1856.1942	0.0000	1.0000
SYMG	61	-759.2618	1859.1169	0.0000	1.0000
SYMI	61	-759.8196	1860.2324	0.0000	1.0000
K2P	56	-774.3666	1861.4090	0.0000	1.0000
SYMIG	62	-758.9952	1864.1672	0.0000	1.0000
JC69	55	-780.7892	1868.6706	0.0000	1.0000
SYM	60	-772.6316	1880.2730	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/chr11_3432.phylip = 1 - 266;
 Partition Dummy = 1:183-phylip/chr11_3432.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/chr11_3432.phylip = 1 - 266;
 Partition Dummy = 1:183-phylip/chr11_3432.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/chr11_3432.phylip = 1 - 266;
 Partition Dummy = 1:183-phylip/chr11_3432.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

