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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	65	-3500.6594	7142.9448	0.8027	0.8027
GTRG	64	-3503.2479	7145.7543	0.1970	0.9997
HKYIG	61	-3513.3113	7158.8166	0.0003	1.0000
HKYG	60	-3516.2417	7162.3354	0.0000	1.0000
GTRI	64	-3520.9194	7181.0973	0.0000	1.0000
SYMIG	62	-3527.5085	7189.5595	0.0000	1.0000
K2PIG	58	-3533.3297	7191.8460	0.0000	1.0000
SYMG	61	-3530.2588	7192.7117	0.0000	1.0000
HKYI	60	-3531.9883	7193.8286	0.0000	1.0000
K2PG	57	-3537.0925	7197.0484	0.0000	1.0000
SYMI	61	-3544.7329	7221.6598	0.0000	1.0000
K2PI	57	-3551.7542	7226.3717	0.0000	1.0000
F81IG	60	-3585.5595	7300.9709	0.0000	1.0000
F81G	59	-3588.1835	7303.8832	0.0000	1.0000
JC69IG	57	-3603.7706	7330.4045	0.0000	1.0000
JC69G	56	-3606.4250	7333.3962	0.0000	1.0000
F81I	59	-3603.0080	7333.5321	0.0000	1.0000
JC69I	56	-3618.3506	7357.2473	0.0000	1.0000
GTR	63	-3616.6757	7370.2488	0.0000	1.0000
HKY	59	-3629.2083	7385.9327	0.0000	1.0000
K2P	56	-3646.9194	7414.3850	0.0000	1.0000
SYM	60	-3644.3421	7418.5361	0.0000	1.0000
F81	58	-3700.0827	7525.3520	0.0000	1.0000
JC69	55	-3718.4910	7555.2173	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-3500.6594	7131.3188	0.8303	0.8303
GTRG	64	-3503.2479	7134.4958	0.1696	0.9998
HKYIG	61	-3513.3113	7148.6225	0.0001	1.0000
HKYG	60	-3516.2417	7152.4834	0.0000	1.0000
GTRI	64	-3520.9194	7169.8388	0.0000	1.0000
SYMIG	62	-3527.5085	7179.0169	0.0000	1.0000
SYMG	61	-3530.2588	7182.5176	0.0000	1.0000
K2PIG	58	-3533.3297	7182.6595	0.0000	1.0000
HKYI	60	-3531.9883	7183.9767	0.0000	1.0000
K2PG	57	-3537.0925	7188.1851	0.0000	1.0000
SYMI	61	-3544.7329	7211.4657	0.0000	1.0000
K2PI	57	-3551.7542	7217.5084	0.0000	1.0000
F81IG	60	-3585.5595	7291.1189	0.0000	1.0000
F81G	59	-3588.1835	7294.3670	0.0000	1.0000
JC69IG	57	-3603.7706	7321.5412	0.0000	1.0000
F81I	59	-3603.0080	7324.0160	0.0000	1.0000
JC69G	56	-3606.4250	7324.8500	0.0000	1.0000
JC69I	56	-3618.3506	7348.7012	0.0000	1.0000
GTR	63	-3616.6757	7359.3515	0.0000	1.0000
HKY	59	-3629.2083	7376.4166	0.0000	1.0000
K2P	56	-3646.9194	7405.8389	0.0000	1.0000
SYM	60	-3644.3421	7408.6841	0.0000	1.0000
F81	58	-3700.0827	7516.1654	0.0000	1.0000
JC69	55	-3718.4910	7546.9820	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-3516.2417	7433.8594	0.3759	0.3759
GTRG	64	-3503.2479	7434.6302	0.2557	0.6316
HKYIG	61	-3513.3113	7434.6881	0.2484	0.8800
GTRIG	65	-3500.6594	7436.1428	0.1200	1.0000
K2PIG	58	-3533.3297	7454.6562	0.0000	1.0000
K2PG	57	-3537.0925	7455.4922	0.0000	1.0000
HKYI	60	-3531.9883	7465.3526	0.0000	1.0000
SYMG	61	-3530.2588	7468.5832	0.0000	1.0000
SYMIG	62	-3527.5085	7469.7721	0.0000	1.0000
GTRI	64	-3520.9194	7469.9732	0.0000	1.0000
K2PI	57	-3551.7542	7484.8156	0.0000	1.0000
SYMI	61	-3544.7329	7497.5313	0.0000	1.0000
F81G	59	-3588.1835	7571.0534	0.0000	1.0000
F81IG	60	-3585.5595	7572.4949	0.0000	1.0000
JC69G	56	-3606.4250	7587.4675	0.0000	1.0000
JC69IG	57	-3603.7706	7588.8484	0.0000	1.0000
F81I	59	-3603.0080	7600.7024	0.0000	1.0000
JC69I	56	-3618.3506	7611.3187	0.0000	1.0000
HKY	59	-3629.2083	7653.1029	0.0000	1.0000
GTR	63	-3616.6757	7654.7962	0.0000	1.0000
K2P	56	-3646.9194	7668.4564	0.0000	1.0000
SYM	60	-3644.3421	7690.0601	0.0000	1.0000
F81	58	-3700.0827	7788.1621	0.0000	1.0000
JC69	55	-3718.4910	7804.9100	0.0000	1.0000



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

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

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


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

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


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

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

