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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKY


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	-708.0220	1581.2292	0.4777	0.4777
HKY	59	-710.5837	1582.6031	0.2403	0.7181
HKYI	60	-708.9056	1582.9964	0.1974	0.9155
HKYIG	61	-708.0213	1585.0240	0.0716	0.9871
GTRG	64	-704.6925	1590.0433	0.0058	0.9930
GTR	63	-707.1619	1591.0408	0.0035	0.9965
GTRI	64	-705.4561	1591.5705	0.0027	0.9992
GTRIG	65	-704.6919	1594.0334	0.0008	1.0000
K2PG	57	-726.8675	1607.8077	0.0000	1.0000
K2PI	57	-727.6660	1609.4046	0.0000	1.0000
K2PIG	58	-726.8680	1611.4676	0.0000	1.0000
K2P	56	-730.5178	1611.4934	0.0000	1.0000
F81G	59	-727.7434	1616.9224	0.0000	1.0000
F81I	59	-727.9586	1617.3528	0.0000	1.0000
F81	58	-729.8686	1617.4690	0.0000	1.0000
SYMG	61	-725.3467	1619.6747	0.0000	1.0000
F81IG	60	-727.7597	1620.7046	0.0000	1.0000
SYMI	61	-726.0177	1621.0167	0.0000	1.0000
SYM	60	-728.9499	1623.0849	0.0000	1.0000
SYMIG	62	-725.3503	1623.5256	0.0000	1.0000
JC69G	56	-743.1870	1636.8319	0.0000	1.0000
JC69I	56	-743.6324	1637.7227	0.0000	1.0000
JC69	55	-746.0258	1638.9378	0.0000	1.0000
JC69IG	57	-743.1944	1640.4616	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-708.0220	1536.0440	0.3282	0.3282
GTRG	64	-704.6925	1537.3850	0.1679	0.4961
HKYI	60	-708.9056	1537.8112	0.1356	0.6317
HKYIG	61	-708.0213	1538.0426	0.1208	0.7525
GTRI	64	-705.4561	1538.9122	0.0782	0.8307
HKY	59	-710.5837	1539.1675	0.0688	0.8996
GTRIG	65	-704.6919	1539.3837	0.0618	0.9614
GTR	63	-707.1619	1540.3238	0.0386	1.0000
K2PG	57	-726.8675	1567.7350	0.0000	1.0000
K2PI	57	-727.6660	1569.3319	0.0000	1.0000
K2PIG	58	-726.8680	1569.7359	0.0000	1.0000
SYMG	61	-725.3467	1572.6933	0.0000	1.0000
K2P	56	-730.5178	1573.0355	0.0000	1.0000
F81G	59	-727.7434	1573.4868	0.0000	1.0000
F81I	59	-727.9586	1573.9172	0.0000	1.0000
SYMI	61	-726.0177	1574.0354	0.0000	1.0000
SYMIG	62	-725.3503	1574.7006	0.0000	1.0000
F81IG	60	-727.7597	1575.5194	0.0000	1.0000
F81	58	-729.8686	1575.7373	0.0000	1.0000
SYM	60	-728.9499	1577.8997	0.0000	1.0000
JC69G	56	-743.1870	1598.3741	0.0000	1.0000
JC69I	56	-743.6324	1599.2648	0.0000	1.0000
JC69IG	57	-743.1944	1600.3888	0.0000	1.0000
JC69	55	-746.0258	1602.0516	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	59	-710.5837	1740.1906	0.4375	0.4375
HKYG	60	-708.0220	1740.4743	0.3796	0.8171
HKYI	60	-708.9056	1742.2415	0.1569	0.9740
HKYIG	61	-708.0213	1745.8801	0.0254	0.9994
GTR	63	-707.1619	1754.9756	0.0003	0.9997
GTRG	64	-704.6925	1755.4440	0.0002	0.9999
GTRI	64	-705.4561	1756.9712	0.0001	1.0000
GTRIG	65	-704.6919	1760.8499	0.0000	1.0000
K2PG	57	-726.8675	1761.9438	0.0000	1.0000
K2PI	57	-727.6660	1763.5407	0.0000	1.0000
K2P	56	-730.5178	1763.8372	0.0000	1.0000
K2PIG	58	-726.8680	1767.3519	0.0000	1.0000
F81	58	-729.8686	1773.3532	0.0000	1.0000
F81G	59	-727.7434	1774.5099	0.0000	1.0000
F81I	59	-727.9586	1774.9404	0.0000	1.0000
F81IG	60	-727.7597	1779.9497	0.0000	1.0000
SYMG	61	-725.3467	1780.5308	0.0000	1.0000
SYMI	61	-726.0177	1781.8729	0.0000	1.0000
SYM	60	-728.9499	1782.3300	0.0000	1.0000
SYMIG	62	-725.3503	1785.9453	0.0000	1.0000
JC69G	56	-743.1870	1789.1757	0.0000	1.0000
JC69	55	-746.0258	1789.4460	0.0000	1.0000
JC69I	56	-743.6324	1790.0664	0.0000	1.0000
JC69IG	57	-743.1944	1794.5976	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/3551_hg18.phylip = 1 - 223;
 Partition Dummy = 1:261-aggressive-phylip/3551_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/3551_hg18.phylip = 1 - 223;
 Partition Dummy = 1:261-aggressive-phylip/3551_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 (HKY)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/3551_hg18.phylip = 1 - 223;
 Partition Dummy = 1:261-aggressive-phylip/3551_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

