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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-2674.8922	5495.1902	0.7608	0.7608
GTRIG	65	-2674.8277	5497.6427	0.2232	0.9840
GTRI	64	-2679.0644	5503.5347	0.0117	0.9957
HKYG	60	-2685.4513	5506.0893	0.0033	0.9990
HKYIG	61	-2685.4177	5508.5610	0.0010	1.0000
HKYI	60	-2689.7872	5514.7611	0.0000	1.0000
SYMG	61	-2695.5756	5528.8767	0.0000	1.0000
SYMIG	62	-2695.5580	5531.3911	0.0000	1.0000
SYMI	61	-2699.1367	5535.9990	0.0000	1.0000
K2PG	57	-2707.5419	5542.7168	0.0000	1.0000
K2PIG	58	-2707.5525	5545.2454	0.0000	1.0000
K2PI	57	-2710.9497	5549.5323	0.0000	1.0000
GTR	63	-2727.1252	5597.0856	0.0000	1.0000
HKY	59	-2736.2772	5605.2127	0.0000	1.0000
SYM	60	-2744.6145	5624.4156	0.0000	1.0000
F81G	59	-2748.8124	5630.2831	0.0000	1.0000
F81IG	60	-2749.7033	5634.5933	0.0000	1.0000
F81I	59	-2752.3631	5637.3846	0.0000	1.0000
K2P	56	-2759.0061	5643.1480	0.0000	1.0000
JC69G	56	-2777.1171	5679.3700	0.0000	1.0000
JC69IG	57	-2777.1361	5681.9051	0.0000	1.0000
JC69I	56	-2780.0774	5685.2905	0.0000	1.0000
F81	58	-2797.9125	5725.9655	0.0000	1.0000
JC69	55	-2823.9399	5770.5286	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2674.8922	5477.7844	0.7093	0.7093
GTRIG	65	-2674.8277	5479.6553	0.2783	0.9877
GTRI	64	-2679.0644	5486.1288	0.0109	0.9986
HKYG	60	-2685.4513	5490.9026	0.0010	0.9996
HKYIG	61	-2685.4177	5492.8354	0.0004	1.0000
HKYI	60	-2689.7872	5499.5743	0.0000	1.0000
SYMG	61	-2695.5756	5513.1512	0.0000	1.0000
SYMIG	62	-2695.5580	5515.1161	0.0000	1.0000
SYMI	61	-2699.1367	5520.2735	0.0000	1.0000
K2PG	57	-2707.5419	5529.0838	0.0000	1.0000
K2PIG	58	-2707.5525	5531.1049	0.0000	1.0000
K2PI	57	-2710.9497	5535.8994	0.0000	1.0000
GTR	63	-2727.1252	5580.2505	0.0000	1.0000
HKY	59	-2736.2772	5590.5543	0.0000	1.0000
SYM	60	-2744.6145	5609.2289	0.0000	1.0000
F81G	59	-2748.8124	5615.6248	0.0000	1.0000
F81IG	60	-2749.7033	5619.4066	0.0000	1.0000
F81I	59	-2752.3631	5622.7262	0.0000	1.0000
K2P	56	-2759.0061	5630.0122	0.0000	1.0000
JC69G	56	-2777.1171	5666.2342	0.0000	1.0000
JC69IG	57	-2777.1361	5668.2721	0.0000	1.0000
JC69I	56	-2780.0774	5672.1547	0.0000	1.0000
F81	58	-2797.9125	5711.8250	0.0000	1.0000
JC69	55	-2823.9399	5757.8797	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2685.4513	5748.7291	0.8360	0.8360
GTRG	64	-2674.8922	5752.7994	0.1092	0.9453
HKYIG	61	-2685.4177	5754.9591	0.0371	0.9824
HKYI	60	-2689.7872	5757.4009	0.0109	0.9933
GTRIG	65	-2674.8277	5758.9674	0.0050	0.9983
GTRI	64	-2679.0644	5761.1438	0.0017	1.0000
K2PG	57	-2707.5419	5774.0191	0.0000	1.0000
SYMG	61	-2695.5756	5775.2748	0.0000	1.0000
K2PIG	58	-2707.5525	5780.3373	0.0000	1.0000
K2PI	57	-2710.9497	5780.8346	0.0000	1.0000
SYMIG	62	-2695.5580	5781.5368	0.0000	1.0000
SYMI	61	-2699.1367	5782.3971	0.0000	1.0000
HKY	59	-2736.2772	5844.0838	0.0000	1.0000
GTR	63	-2727.1252	5850.9684	0.0000	1.0000
SYM	60	-2744.6145	5867.0555	0.0000	1.0000
F81G	59	-2748.8124	5869.1542	0.0000	1.0000
K2P	56	-2759.0061	5870.6503	0.0000	1.0000
F81I	59	-2752.3631	5876.2557	0.0000	1.0000
F81IG	60	-2749.7033	5877.2332	0.0000	1.0000
JC69G	56	-2777.1171	5906.8723	0.0000	1.0000
JC69I	56	-2780.0774	5912.7929	0.0000	1.0000
JC69IG	57	-2777.1361	5913.2074	0.0000	1.0000
F81	58	-2797.9125	5961.0573	0.0000	1.0000
JC69	55	-2823.9399	5994.2207	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/8668_hg18.phylip = 1 - 543;
 Partition Dummy = 1:261-aggressive-phylip/8668_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/8668_hg18.phylip = 1 - 543;
 Partition Dummy = 1:261-aggressive-phylip/8668_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 (HKYG)]

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

