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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: F81G


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
GTRG	64	-652.7346	1465.7173	0.4308	0.4308
GTRI	64	-652.7988	1465.8457	0.4040	0.8349
GTRIG	65	-652.6673	1468.7199	0.0960	0.9309
GTR	63	-656.4682	1470.0716	0.0488	0.9797
F81G	59	-664.7941	1474.5083	0.0053	0.9850
F81I	59	-664.8448	1474.6098	0.0051	0.9901
HKYG	60	-663.6959	1475.3307	0.0035	0.9936
HKYI	60	-663.7560	1475.4509	0.0033	0.9969
F81IG	60	-664.7130	1477.3650	0.0013	0.9982
HKYIG	61	-663.6049	1478.1906	0.0008	0.9990
F81	58	-668.5218	1478.9678	0.0006	0.9996
HKY	59	-667.3810	1479.6821	0.0004	1.0000
JC69G	56	-683.6360	1503.2721	0.0000	1.0000
JC69I	56	-683.6750	1503.3500	0.0000	1.0000
SYMG	61	-676.3342	1503.6492	0.0000	1.0000
SYMI	61	-676.3713	1503.7234	0.0000	1.0000
K2PG	57	-683.0029	1504.9567	0.0000	1.0000
K2PI	57	-683.0445	1505.0400	0.0000	1.0000
JC69IG	57	-683.5598	1506.0706	0.0000	1.0000
SYMIG	62	-676.3036	1506.6533	0.0000	1.0000
SYM	60	-679.8233	1507.5855	0.0000	1.0000
K2PIG	58	-682.9196	1507.7635	0.0000	1.0000
JC69	55	-687.4841	1508.0393	0.0000	1.0000
K2P	56	-686.8506	1509.7013	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-652.7346	1433.4692	0.4166	0.4166
GTRI	64	-652.7988	1433.5977	0.3906	0.8072
GTRIG	65	-652.6673	1435.3347	0.1639	0.9711
GTR	63	-656.4682	1438.9365	0.0271	0.9982
HKYG	60	-663.6959	1447.3917	0.0004	0.9986
HKYI	60	-663.7560	1447.5120	0.0004	0.9989
F81G	59	-664.7941	1447.5882	0.0004	0.9993
F81I	59	-664.8448	1447.6896	0.0003	0.9996
HKYIG	61	-663.6049	1449.2098	0.0002	0.9998
F81IG	60	-664.7130	1449.4261	0.0001	0.9999
HKY	59	-667.3810	1452.7619	0.0000	1.0000
F81	58	-668.5218	1453.0436	0.0000	1.0000
SYMG	61	-676.3342	1474.6684	0.0000	1.0000
SYMI	61	-676.3713	1474.7425	0.0000	1.0000
SYMIG	62	-676.3036	1476.6071	0.0000	1.0000
JC69G	56	-683.6360	1479.2721	0.0000	1.0000
JC69I	56	-683.6750	1479.3500	0.0000	1.0000
SYM	60	-679.8233	1479.6466	0.0000	1.0000
K2PG	57	-683.0029	1480.0057	0.0000	1.0000
K2PI	57	-683.0445	1480.0891	0.0000	1.0000
JC69IG	57	-683.5598	1481.1197	0.0000	1.0000
K2PIG	58	-682.9196	1481.8393	0.0000	1.0000
JC69	55	-687.4841	1484.9681	0.0000	1.0000
K2P	56	-686.8506	1485.7013	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
F81G	59	-664.7941	1670.4697	0.3251	0.3251
F81I	59	-664.8448	1670.5711	0.3090	0.6341
F81	58	-668.5218	1672.1474	0.1405	0.7746
HKYG	60	-663.6959	1674.0509	0.0542	0.8288
HKYI	60	-663.7560	1674.1711	0.0511	0.8799
GTRG	64	-652.7346	1675.2390	0.0299	0.9098
GTRI	64	-652.7988	1675.3674	0.0281	0.9379
HKY	59	-667.3810	1675.6434	0.0245	0.9624
F81IG	60	-664.7130	1676.0852	0.0196	0.9820
GTR	63	-656.4682	1676.9286	0.0129	0.9949
HKYIG	61	-663.6049	1679.6466	0.0033	0.9982
GTRIG	65	-652.6673	1680.8821	0.0018	1.0000
JC69G	56	-683.6360	1690.8206	0.0000	1.0000
JC69I	56	-683.6750	1690.8986	0.0000	1.0000
JC69	55	-687.4841	1692.7390	0.0000	1.0000
K2PG	57	-683.0029	1695.3319	0.0000	1.0000
K2PI	57	-683.0445	1695.4153	0.0000	1.0000
JC69IG	57	-683.5598	1696.4458	0.0000	1.0000
K2P	56	-686.8506	1697.2498	0.0000	1.0000
K2PIG	58	-682.9196	1700.9431	0.0000	1.0000
SYMG	61	-676.3342	1705.1052	0.0000	1.0000
SYMI	61	-676.3713	1705.1793	0.0000	1.0000
SYM	60	-679.8233	1706.3057	0.0000	1.0000
SYMIG	62	-676.3036	1710.8216	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/7141_hg18.phylip = 1 - 323;
 Partition Dummy = 1:261-aggressive-phylip/7141_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 (GTRG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7141_hg18.phylip = 1 - 323;
 Partition Dummy = 1:261-aggressive-phylip/7141_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 BIC model (F81G)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7141_hg18.phylip = 1 - 323;
 Partition Dummy = 1:261-aggressive-phylip/7141_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=1 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

