
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:30 08/01/2011

Input data from file "261-aggressive-phylip/3831_hg18.phylip" (29 taxa, 209 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-731.5276	1632.5148	0.5556	0.5556
K2PG	57	-738.6691	1635.1264	0.1505	0.7061
K2PI	57	-739.0195	1635.8271	0.1060	0.8122
HKYI	60	-733.4138	1636.2871	0.0843	0.8964
HKYIG	61	-731.4815	1636.4189	0.0789	0.9753
K2PIG	58	-738.6672	1638.9611	0.0221	0.9974
GTRG	64	-729.6087	1644.9951	0.0011	0.9985
GTRI	64	-730.0277	1645.8331	0.0007	0.9992
SYMG	61	-736.8922	1647.2401	0.0004	0.9996
SYMI	61	-737.2274	1647.9106	0.0003	0.9998
GTRIG	65	-729.5782	1649.1564	0.0001	1.0000
SYMIG	62	-736.8893	1651.2854	0.0000	1.0000
F81G	59	-747.7855	1661.0879	0.0000	1.0000
F81I	59	-748.1751	1661.8669	0.0000	1.0000
K2P	56	-753.9497	1661.8994	0.0000	1.0000
HKY	59	-748.7470	1663.0108	0.0000	1.0000
F81IG	60	-747.7860	1665.0314	0.0000	1.0000
JC69G	56	-755.7716	1665.5433	0.0000	1.0000
JC69I	56	-756.1367	1666.2734	0.0000	1.0000
JC69IG	57	-755.7703	1669.3286	0.0000	1.0000
SYM	60	-751.8018	1673.0631	0.0000	1.0000
GTR	63	-746.9588	1675.5313	0.0000	1.0000
JC69	55	-770.7807	1691.8228	0.0000	1.0000
F81	58	-765.1431	1691.9129	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-731.5276	1583.0553	0.5471	0.5471
HKYIG	61	-731.4815	1584.9631	0.2108	0.7579
HKYI	60	-733.4138	1586.8277	0.0830	0.8409
GTRG	64	-729.6087	1587.2173	0.0683	0.9092
GTRI	64	-730.0277	1588.0554	0.0449	0.9541
GTRIG	65	-729.5782	1589.1564	0.0259	0.9800
K2PG	57	-738.6691	1591.3383	0.0087	0.9887
K2PI	57	-739.0195	1592.0390	0.0061	0.9948
K2PIG	58	-738.6672	1593.3344	0.0032	0.9980
SYMG	61	-736.8922	1595.7843	0.0009	0.9990
SYMI	61	-737.2274	1596.4548	0.0007	0.9997
SYMIG	62	-736.8893	1597.7786	0.0003	1.0000
F81G	59	-747.7855	1613.5711	0.0000	1.0000
F81I	59	-748.1751	1614.3501	0.0000	1.0000
HKY	59	-748.7470	1615.4940	0.0000	1.0000
F81IG	60	-747.7860	1615.5720	0.0000	1.0000
K2P	56	-753.9497	1619.8994	0.0000	1.0000
GTR	63	-746.9588	1619.9175	0.0000	1.0000
JC69G	56	-755.7716	1623.5433	0.0000	1.0000
SYM	60	-751.8018	1623.6036	0.0000	1.0000
JC69I	56	-756.1367	1624.2734	0.0000	1.0000
JC69IG	57	-755.7703	1625.5405	0.0000	1.0000
F81	58	-765.1431	1646.2863	0.0000	1.0000
JC69	55	-770.7807	1651.5614	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-738.6691	1781.8514	0.4375	0.4375
K2PI	57	-739.0195	1782.5520	0.3082	0.7456
HKYG	60	-731.5276	1783.5954	0.1829	0.9285
K2PIG	58	-738.6672	1787.1898	0.0303	0.9589
HKYI	60	-733.4138	1787.3677	0.0277	0.9866
HKYIG	61	-731.4815	1788.8455	0.0132	0.9998
SYMG	61	-736.8922	1799.6667	0.0001	0.9999
SYMI	61	-737.2274	1800.3372	0.0000	0.9999
GTRG	64	-729.6087	1801.1267	0.0000	1.0000
GTRI	64	-730.0277	1801.9648	0.0000	1.0000
SYMIG	62	-736.8893	1805.0033	0.0000	1.0000
GTRIG	65	-729.5782	1806.4082	0.0000	1.0000
K2P	56	-753.9497	1807.0701	0.0000	1.0000
JC69G	56	-755.7716	1810.7140	0.0000	1.0000
F81G	59	-747.7855	1810.7688	0.0000	1.0000
JC69I	56	-756.1367	1811.4441	0.0000	1.0000
F81I	59	-748.1751	1811.5478	0.0000	1.0000
HKY	59	-748.7470	1812.6917	0.0000	1.0000
JC69IG	57	-755.7703	1816.0536	0.0000	1.0000
F81IG	60	-747.7860	1816.1120	0.0000	1.0000
SYM	60	-751.8018	1824.1437	0.0000	1.0000
GTR	63	-746.9588	1830.4846	0.0000	1.0000
JC69	55	-770.7807	1835.3897	0.0000	1.0000
F81	58	-765.1431	1840.1416	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 261-aggressive-phylip/3831_hg18.phylip = 1 - 209;
 Partition Dummy = 1:261-aggressive-phylip/3831_hg18.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 261-aggressive-phylip/3831_hg18.phylip = 1 - 209;
 Partition Dummy = 1:261-aggressive-phylip/3831_hg18.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 (K2PG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/3831_hg18.phylip = 1 - 209;
 Partition Dummy = 1:261-aggressive-phylip/3831_hg18.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

