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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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	-3496.6682	7129.1121	0.7432	0.7432
GTRIG	65	-3496.6766	7131.3794	0.2392	0.9824
HKYG	60	-3505.1903	7137.1963	0.0131	0.9955
HKYIG	61	-3505.1989	7139.4472	0.0042	0.9997
K2PG	57	-3513.1671	7146.4735	0.0001	0.9999
SYMG	61	-3509.2343	7147.5180	0.0001	0.9999
K2PIG	58	-3513.1757	7148.7120	0.0000	1.0000
SYMIG	62	-3509.2428	7149.7730	0.0000	1.0000
F81G	59	-3558.3633	7241.3126	0.0000	1.0000
F81IG	60	-3558.3718	7243.5593	0.0000	1.0000
JC69G	56	-3563.6905	7245.3030	0.0000	1.0000
JC69IG	57	-3563.6991	7247.5375	0.0000	1.0000
GTRI	64	-3578.3287	7292.4331	0.0000	1.0000
HKYI	60	-3586.6004	7300.0164	0.0000	1.0000
K2PI	57	-3594.3428	7308.8249	0.0000	1.0000
SYMI	61	-3590.1417	7309.3327	0.0000	1.0000
GTR	63	-3589.9716	7313.4727	0.0000	1.0000
HKY	59	-3597.4866	7319.5592	0.0000	1.0000
SYM	60	-3600.2416	7327.2988	0.0000	1.0000
K2P	56	-3604.6936	7327.3093	0.0000	1.0000
F81I	59	-3638.9371	7402.4603	0.0000	1.0000
JC69I	56	-3644.3158	7406.5537	0.0000	1.0000
F81	58	-3649.5450	7421.4507	0.0000	1.0000
JC69	55	-3654.7899	7425.2888	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-3496.6682	7121.3364	0.7247	0.7247
GTRIG	65	-3496.6766	7123.3532	0.2644	0.9891
HKYG	60	-3505.1903	7130.3807	0.0079	0.9970
HKYIG	61	-3505.1989	7132.3978	0.0029	0.9999
K2PG	57	-3513.1671	7140.3343	0.0001	0.9999
SYMG	61	-3509.2343	7140.4686	0.0001	1.0000
K2PIG	58	-3513.1757	7142.3514	0.0000	1.0000
SYMIG	62	-3509.2428	7142.4857	0.0000	1.0000
F81G	59	-3558.3633	7234.7266	0.0000	1.0000
F81IG	60	-3558.3718	7236.7437	0.0000	1.0000
JC69G	56	-3563.6905	7239.3809	0.0000	1.0000
JC69IG	57	-3563.6991	7241.3982	0.0000	1.0000
GTRI	64	-3578.3287	7284.6574	0.0000	1.0000
HKYI	60	-3586.6004	7293.2007	0.0000	1.0000
SYMI	61	-3590.1417	7302.2833	0.0000	1.0000
K2PI	57	-3594.3428	7302.6856	0.0000	1.0000
GTR	63	-3589.9716	7305.9433	0.0000	1.0000
HKY	59	-3597.4866	7312.9731	0.0000	1.0000
SYM	60	-3600.2416	7320.4831	0.0000	1.0000
K2P	56	-3604.6936	7321.3872	0.0000	1.0000
F81I	59	-3638.9371	7395.8742	0.0000	1.0000
JC69I	56	-3644.3158	7400.6316	0.0000	1.0000
F81	58	-3649.5450	7415.0901	0.0000	1.0000
JC69	55	-3654.7899	7419.5798	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-3513.1671	7427.2944	0.9024	0.9024
HKYG	60	-3505.1903	7432.4440	0.0687	0.9711
K2PIG	58	-3513.1757	7434.3459	0.0266	0.9977
HKYIG	61	-3505.1989	7439.4954	0.0020	0.9997
GTRG	64	-3496.6682	7443.5372	0.0003	1.0000
SYMG	61	-3509.2343	7447.5663	0.0000	1.0000
GTRIG	65	-3496.6766	7450.5885	0.0000	1.0000
SYMIG	62	-3509.2428	7454.6177	0.0000	1.0000
JC69G	56	-3563.6905	7521.3067	0.0000	1.0000
JC69IG	57	-3563.6991	7528.3584	0.0000	1.0000
F81G	59	-3558.3633	7531.7554	0.0000	1.0000
F81IG	60	-3558.3718	7538.8070	0.0000	1.0000
K2PI	57	-3594.3428	7589.6457	0.0000	1.0000
HKYI	60	-3586.6004	7595.2640	0.0000	1.0000
K2P	56	-3604.6936	7603.3130	0.0000	1.0000
GTRI	64	-3578.3287	7606.8583	0.0000	1.0000
SYMI	61	-3590.1417	7609.3810	0.0000	1.0000
HKY	59	-3597.4866	7610.0020	0.0000	1.0000
SYM	60	-3600.2416	7622.5464	0.0000	1.0000
GTR	63	-3589.9716	7623.1097	0.0000	1.0000
JC69I	56	-3644.3158	7682.5573	0.0000	1.0000
F81I	59	-3638.9371	7692.9031	0.0000	1.0000
JC69	55	-3654.7899	7696.4712	0.0000	1.0000
F81	58	-3649.5450	7707.0846	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/7140_hg18.phylip = 1 - 1135;
 Partition Dummy = 1:261-aggressive-phylip/7140_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/7140_hg18.phylip = 1 - 1135;
 Partition Dummy = 1:261-aggressive-phylip/7140_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 (K2PG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7140_hg18.phylip = 1 - 1135;
 Partition Dummy = 1:261-aggressive-phylip/7140_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

