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

Input data from file "261-aggressive-phylip/10084_hg18.phylip" (29 taxa, 586 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	-2674.4515	5492.8723	0.4381	0.4381
GTRIG	65	-2673.7508	5494.0016	0.2491	0.6872
HKYIG	61	-2678.9709	5494.3769	0.2065	0.8937
HKYG	60	-2681.1257	5496.1943	0.0832	0.9769
SYMG	61	-2681.6033	5499.6418	0.0148	0.9918
SYMIG	62	-2681.4277	5501.7924	0.0051	0.9969
K2PIG	58	-2687.2158	5503.4182	0.0022	0.9991
K2PG	57	-2689.3632	5505.2492	0.0009	1.0000
GTRI	64	-2687.1769	5518.3231	0.0000	1.0000
HKYI	60	-2693.3182	5520.5792	0.0000	1.0000
SYMI	61	-2694.5441	5525.5234	0.0000	1.0000
K2PI	57	-2701.2971	5529.1170	0.0000	1.0000
GTR	63	-2781.3490	5704.1462	0.0000	1.0000
HKY	59	-2788.3832	5708.2265	0.0000	1.0000
F81IG	60	-2787.2777	5708.4983	0.0000	1.0000
SYM	60	-2788.0732	5710.0892	0.0000	1.0000
F81G	59	-2789.4187	5710.2976	0.0000	1.0000
K2P	56	-2796.1902	5716.4484	0.0000	1.0000
F81I	59	-2799.1062	5729.6724	0.0000	1.0000
JC69G	56	-2803.8106	5731.6893	0.0000	1.0000
JC69IG	57	-2804.0105	5734.5438	0.0000	1.0000
JC69I	56	-2812.0770	5748.2220	0.0000	1.0000
F81	58	-2890.6364	5910.2594	0.0000	1.0000
JC69	55	-2905.6188	5932.8602	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2674.4515	5476.9030	0.4870	0.4870
GTRIG	65	-2673.7508	5477.5016	0.3610	0.8481
HKYIG	61	-2678.9709	5479.9418	0.1066	0.9546
HKYG	60	-2681.1257	5482.2515	0.0336	0.9882
SYMG	61	-2681.6033	5485.2066	0.0077	0.9959
SYMIG	62	-2681.4277	5486.8555	0.0034	0.9993
K2PIG	58	-2687.2158	5490.4315	0.0006	0.9998
K2PG	57	-2689.3632	5492.7265	0.0002	1.0000
GTRI	64	-2687.1769	5502.3538	0.0000	1.0000
HKYI	60	-2693.3182	5506.6363	0.0000	1.0000
SYMI	61	-2694.5441	5511.0883	0.0000	1.0000
K2PI	57	-2701.2971	5516.5943	0.0000	1.0000
GTR	63	-2781.3490	5688.6979	0.0000	1.0000
F81IG	60	-2787.2777	5694.5554	0.0000	1.0000
HKY	59	-2788.3832	5694.7664	0.0000	1.0000
SYM	60	-2788.0732	5696.1463	0.0000	1.0000
F81G	59	-2789.4187	5696.8375	0.0000	1.0000
K2P	56	-2796.1902	5704.3803	0.0000	1.0000
F81I	59	-2799.1062	5716.2123	0.0000	1.0000
JC69G	56	-2803.8106	5719.6213	0.0000	1.0000
JC69IG	57	-2804.0105	5722.0211	0.0000	1.0000
JC69I	56	-2812.0770	5736.1539	0.0000	1.0000
F81	58	-2890.6364	5897.2727	0.0000	1.0000
JC69	55	-2905.6188	5921.2375	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-2689.3632	5742.0057	0.5804	0.5804
K2PIG	58	-2687.2158	5744.0841	0.2053	0.7857
HKYG	60	-2681.1257	5744.6507	0.1547	0.9403
HKYIG	61	-2678.9709	5746.7143	0.0551	0.9955
SYMG	61	-2681.6033	5751.9791	0.0040	0.9994
GTRG	64	-2674.4515	5756.7955	0.0004	0.9998
SYMIG	62	-2681.4277	5758.0013	0.0002	1.0000
GTRIG	65	-2673.7508	5761.7674	0.0000	1.0000
K2PI	57	-2701.2971	5765.8735	0.0000	1.0000
HKYI	60	-2693.3182	5769.0355	0.0000	1.0000
SYMI	61	-2694.5441	5777.8608	0.0000	1.0000
GTRI	64	-2687.1769	5782.2463	0.0000	1.0000
K2P	56	-2796.1902	5949.2862	0.0000	1.0000
HKY	59	-2788.3832	5952.7923	0.0000	1.0000
F81G	59	-2789.4187	5954.8633	0.0000	1.0000
F81IG	60	-2787.2777	5956.9546	0.0000	1.0000
SYM	60	-2788.0732	5958.5455	0.0000	1.0000
GTR	63	-2781.3490	5964.2170	0.0000	1.0000
JC69G	56	-2803.8106	5964.5272	0.0000	1.0000
JC69IG	57	-2804.0105	5971.3003	0.0000	1.0000
F81I	59	-2799.1062	5974.2382	0.0000	1.0000
JC69I	56	-2812.0770	5981.0598	0.0000	1.0000
F81	58	-2890.6364	6150.9252	0.0000	1.0000
JC69	55	-2905.6188	6161.7701	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/10084_hg18.phylip = 1 - 586;
 Partition Dummy = 1:261-aggressive-phylip/10084_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/10084_hg18.phylip = 1 - 586;
 Partition Dummy = 1:261-aggressive-phylip/10084_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/10084_hg18.phylip = 1 - 586;
 Partition Dummy = 1:261-aggressive-phylip/10084_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

