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

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

Minimum AIC  model: GTRG
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	60	-665.2172	1484.6400	0.4062	0.4062
HKYG	60	-665.3255	1484.8567	0.3645	0.7706
HKYIG	61	-664.7807	1487.0732	0.1203	0.8910
GTRG	64	-660.5664	1488.7519	0.0520	0.9429
GTRI	64	-660.8500	1489.3190	0.0391	0.9821
GTRIG	65	-659.9156	1490.8839	0.0179	1.0000
HKY	59	-678.0133	1506.9569	0.0000	1.0000
GTR	63	-672.0610	1508.3400	0.0000	1.0000
F81I	59	-678.8519	1508.6340	0.0000	1.0000
F81G	59	-678.9484	1508.8270	0.0000	1.0000
F81IG	60	-678.5497	1511.3050	0.0000	1.0000
K2PI	57	-684.7812	1514.0324	0.0000	1.0000
K2PG	57	-684.9205	1514.3111	0.0000	1.0000
K2PIG	58	-684.5058	1516.6967	0.0000	1.0000
SYMG	61	-682.6814	1522.8746	0.0000	1.0000
SYMI	61	-682.6865	1522.8848	0.0000	1.0000
SYMIG	62	-682.3447	1525.5385	0.0000	1.0000
F81	58	-691.2090	1530.1032	0.0000	1.0000
JC69I	56	-695.9727	1533.2299	0.0000	1.0000
JC69G	56	-696.0827	1533.4497	0.0000	1.0000
JC69IG	57	-695.7150	1535.8999	0.0000	1.0000
K2P	56	-697.4614	1536.2073	0.0000	1.0000
SYM	60	-694.1280	1542.4617	0.0000	1.0000
JC69	55	-708.4052	1554.9382	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-660.5664	1449.1329	0.2670	0.2670
GTRI	64	-660.8500	1449.6999	0.2011	0.4681
GTRIG	65	-659.9156	1449.8312	0.1883	0.6564
HKYI	60	-665.2172	1450.4344	0.1393	0.7957
HKYG	60	-665.3255	1450.6511	0.1250	0.9207
HKYIG	61	-664.7807	1451.5614	0.0793	1.0000
GTR	63	-672.0610	1470.1220	0.0000	1.0000
HKY	59	-678.0133	1474.0267	0.0000	1.0000
F81I	59	-678.8519	1475.7038	0.0000	1.0000
F81G	59	-678.9484	1475.8968	0.0000	1.0000
F81IG	60	-678.5497	1477.0994	0.0000	1.0000
K2PI	57	-684.7812	1483.5623	0.0000	1.0000
K2PG	57	-684.9205	1483.8410	0.0000	1.0000
K2PIG	58	-684.5058	1485.0116	0.0000	1.0000
SYMG	61	-682.6814	1487.3629	0.0000	1.0000
SYMI	61	-682.6865	1487.3730	0.0000	1.0000
SYMIG	62	-682.3447	1488.6894	0.0000	1.0000
F81	58	-691.2090	1498.4180	0.0000	1.0000
JC69I	56	-695.9727	1503.9455	0.0000	1.0000
JC69G	56	-696.0827	1504.1653	0.0000	1.0000
JC69IG	57	-695.7150	1505.4299	0.0000	1.0000
K2P	56	-697.4614	1506.9229	0.0000	1.0000
SYM	60	-694.1280	1508.2561	0.0000	1.0000
JC69	55	-708.4052	1526.8104	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	60	-665.2172	1667.4407	0.5017	0.5017
HKYG	60	-665.3255	1667.6573	0.4502	0.9518
HKYIG	61	-664.7807	1672.1845	0.0468	0.9986
GTRG	64	-660.5664	1680.6062	0.0007	0.9993
GTRI	64	-660.8500	1681.1733	0.0005	0.9999
GTRIG	65	-659.9156	1684.9214	0.0001	0.9999
HKY	59	-678.0133	1687.4162	0.0000	1.0000
F81I	59	-678.8519	1689.0933	0.0000	1.0000
F81G	59	-678.9484	1689.2863	0.0000	1.0000
K2PI	57	-684.7812	1689.7183	0.0000	1.0000
K2PG	57	-684.9205	1689.9970	0.0000	1.0000
F81IG	60	-678.5497	1694.1057	0.0000	1.0000
K2PIG	58	-684.5058	1694.7843	0.0000	1.0000
GTR	63	-672.0610	1697.9786	0.0000	1.0000
JC69I	56	-695.9727	1706.4847	0.0000	1.0000
JC69G	56	-696.0827	1706.7045	0.0000	1.0000
SYMG	61	-682.6814	1707.9859	0.0000	1.0000
SYMI	61	-682.6865	1707.9961	0.0000	1.0000
F81	58	-691.2090	1708.1907	0.0000	1.0000
K2P	56	-697.4614	1709.4620	0.0000	1.0000
JC69IG	57	-695.7150	1711.5859	0.0000	1.0000
SYMIG	62	-682.3447	1712.9292	0.0000	1.0000
SYM	60	-694.1280	1725.2623	0.0000	1.0000
JC69	55	-708.4052	1725.7328	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 261-aggressive-phylip/12098_hg18.phylip = 1 - 275;
 Partition Dummy = 1:261-aggressive-phylip/12098_hg18.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 (HKYI)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/12098_hg18.phylip = 1 - 275;
 Partition Dummy = 1:261-aggressive-phylip/12098_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/12098_hg18.phylip = 1 - 275;
 Partition Dummy = 1:261-aggressive-phylip/12098_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

