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

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

Minimum AIC  model: HKYIG
Minimum AICc model: HKYIG
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
HKYIG	61	-1532.2175	3210.4477	0.5279	0.5279
HKYG	60	-1533.8086	3210.7818	0.4467	0.9746
GTRIG	65	-1530.2683	3218.1250	0.0114	0.9860
HKYI	60	-1537.8619	3218.8883	0.0078	0.9937
GTRG	64	-1532.3230	3219.3127	0.0063	1.0000
GTRI	64	-1542.2265	3239.1196	0.0000	1.0000
K2PG	57	-1559.4017	3253.5307	0.0000	1.0000
K2PIG	58	-1558.6400	3254.8020	0.0000	1.0000
K2PI	57	-1563.0912	3260.9096	0.0000	1.0000
SYMG	61	-1559.8144	3265.6415	0.0000	1.0000
SYMIG	62	-1558.6973	3266.2736	0.0000	1.0000
F81G	59	-1565.6430	3271.6205	0.0000	1.0000
F81IG	60	-1564.3784	3271.9213	0.0000	1.0000
SYMI	61	-1564.2539	3274.5205	0.0000	1.0000
F81I	59	-1570.1518	3280.6379	0.0000	1.0000
JC69G	56	-1585.3847	3302.7194	0.0000	1.0000
JC69IG	57	-1584.4353	3303.5979	0.0000	1.0000
JC69I	56	-1589.5225	3310.9950	0.0000	1.0000
HKY	59	-1596.7882	3333.9108	0.0000	1.0000
GTR	63	-1598.0134	3347.7904	0.0000	1.0000
K2P	56	-1620.9992	3373.9484	0.0000	1.0000
SYM	60	-1622.7941	3388.7528	0.0000	1.0000
F81	58	-1627.7813	3393.0845	0.0000	1.0000
JC69	55	-1647.2063	3423.6026	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	61	-1532.2175	3186.4350	0.5758	0.5758
HKYG	60	-1533.8086	3187.6172	0.3188	0.8946
GTRIG	65	-1530.2683	3190.5366	0.0741	0.9687
GTRG	64	-1532.3230	3192.6460	0.0258	0.9945
HKYI	60	-1537.8619	3195.7237	0.0055	1.0000
GTRI	64	-1542.2265	3212.4529	0.0000	1.0000
K2PG	57	-1559.4017	3232.8034	0.0000	1.0000
K2PIG	58	-1558.6400	3233.2800	0.0000	1.0000
K2PI	57	-1563.0912	3240.1824	0.0000	1.0000
SYMIG	62	-1558.6973	3241.3947	0.0000	1.0000
SYMG	61	-1559.8144	3241.6288	0.0000	1.0000
F81IG	60	-1564.3784	3248.7567	0.0000	1.0000
F81G	59	-1565.6430	3249.2861	0.0000	1.0000
SYMI	61	-1564.2539	3250.5078	0.0000	1.0000
F81I	59	-1570.1518	3258.3035	0.0000	1.0000
JC69G	56	-1585.3847	3282.7694	0.0000	1.0000
JC69IG	57	-1584.4353	3282.8706	0.0000	1.0000
JC69I	56	-1589.5225	3291.0450	0.0000	1.0000
HKY	59	-1596.7882	3311.5764	0.0000	1.0000
GTR	63	-1598.0134	3322.0268	0.0000	1.0000
K2P	56	-1620.9992	3353.9984	0.0000	1.0000
SYM	60	-1622.7941	3365.5882	0.0000	1.0000
F81	58	-1627.7813	3371.5625	0.0000	1.0000
JC69	55	-1647.2063	3404.4125	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1533.8086	3423.5520	0.7872	0.7872
HKYIG	61	-1532.2175	3426.3019	0.1990	0.9863
HKYI	60	-1537.8619	3431.6584	0.0137	1.0000
GTRG	64	-1532.3230	3444.3097	0.0000	1.0000
GTRIG	65	-1530.2683	3446.1325	0.0000	1.0000
K2PG	57	-1559.4017	3456.9414	0.0000	1.0000
K2PIG	58	-1558.6400	3461.3502	0.0000	1.0000
GTRI	64	-1542.2265	3464.1166	0.0000	1.0000
K2PI	57	-1563.0912	3464.3203	0.0000	1.0000
F81G	59	-1565.6430	3481.2885	0.0000	1.0000
SYMG	61	-1559.8144	3481.4958	0.0000	1.0000
F81IG	60	-1564.3784	3484.6915	0.0000	1.0000
SYMIG	62	-1558.6973	3485.1939	0.0000	1.0000
F81I	59	-1570.1518	3490.3060	0.0000	1.0000
SYMI	61	-1564.2539	3490.3748	0.0000	1.0000
JC69G	56	-1585.3847	3502.9752	0.0000	1.0000
JC69IG	57	-1584.4353	3507.0086	0.0000	1.0000
JC69I	56	-1589.5225	3511.2507	0.0000	1.0000
HKY	59	-1596.7882	3543.5789	0.0000	1.0000
GTR	63	-1598.0134	3569.7582	0.0000	1.0000
K2P	56	-1620.9992	3574.2042	0.0000	1.0000
F81	58	-1627.7813	3599.6327	0.0000	1.0000
SYM	60	-1622.7941	3601.5229	0.0000	1.0000
JC69	55	-1647.2063	3620.6860	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/6562_hg18.phylip = 1 - 377;
 Partition Dummy = 1:261-aggressive-phylip/6562_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/6562_hg18.phylip = 1 - 377;
 Partition Dummy = 1:261-aggressive-phylip/6562_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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


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

