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

Input data from file "261-aggressive-phylip/12814_hg18.phylip" (29 taxa, 434 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	-1295.2402	2741.0278	0.4501	0.4501
GTRI	64	-1296.0593	2742.6659	0.1984	0.6486
GTRIG	65	-1295.0652	2743.4457	0.1344	0.7829
HKYG	60	-1302.0030	2743.6307	0.1225	0.9054
HKYI	60	-1302.8315	2745.2876	0.0535	0.9589
HKYIG	61	-1301.8117	2745.9567	0.0383	0.9972
K2PG	57	-1310.6580	2752.9011	0.0012	0.9984
K2PI	57	-1311.4561	2754.4972	0.0005	0.9990
K2PIG	58	-1310.4274	2755.1054	0.0004	0.9994
SYMG	61	-1306.5011	2755.3356	0.0004	0.9997
SYMI	61	-1307.1800	2756.6934	0.0002	0.9999
SYMIG	62	-1306.2943	2757.6452	0.0001	1.0000
GTR	63	-1307.9683	2763.7311	0.0000	1.0000
HKY	59	-1315.2251	2767.3808	0.0000	1.0000
K2P	56	-1324.9854	2778.9046	0.0000	1.0000
SYM	60	-1320.2021	2780.0288	0.0000	1.0000
F81G	59	-1326.2504	2789.4312	0.0000	1.0000
F81I	59	-1327.0533	2791.0371	0.0000	1.0000
F81IG	60	-1326.0387	2791.7020	0.0000	1.0000
JC69G	56	-1333.9939	2796.9215	0.0000	1.0000
JC69I	56	-1334.7903	2798.5142	0.0000	1.0000
JC69IG	57	-1333.7597	2799.1046	0.0000	1.0000
F81	58	-1339.7169	2813.6844	0.0000	1.0000
JC69	55	-1348.0771	2822.4505	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1295.2402	2718.4804	0.5003	0.5003
GTRI	64	-1296.0593	2720.1185	0.2206	0.7209
GTRIG	65	-1295.0652	2720.1304	0.2192	0.9401
HKYG	60	-1302.0030	2724.0060	0.0316	0.9717
HKYIG	61	-1301.8117	2725.6233	0.0141	0.9857
HKYI	60	-1302.8315	2725.6629	0.0138	0.9995
SYMG	61	-1306.5011	2735.0022	0.0001	0.9997
K2PG	57	-1310.6580	2735.3160	0.0001	0.9998
SYMI	61	-1307.1800	2736.3601	0.0001	0.9998
SYMIG	62	-1306.2943	2736.5886	0.0001	0.9999
K2PIG	58	-1310.4274	2736.8548	0.0001	0.9999
K2PI	57	-1311.4561	2736.9121	0.0000	1.0000
GTR	63	-1307.9683	2741.9365	0.0000	1.0000
HKY	59	-1315.2251	2748.4503	0.0000	1.0000
SYM	60	-1320.2021	2760.4041	0.0000	1.0000
K2P	56	-1324.9854	2761.9709	0.0000	1.0000
F81G	59	-1326.2504	2770.5007	0.0000	1.0000
F81IG	60	-1326.0387	2772.0774	0.0000	1.0000
F81I	59	-1327.0533	2772.1066	0.0000	1.0000
JC69G	56	-1333.9939	2779.9878	0.0000	1.0000
JC69IG	57	-1333.7597	2781.5195	0.0000	1.0000
JC69I	56	-1334.7903	2781.5805	0.0000	1.0000
F81	58	-1339.7169	2795.4337	0.0000	1.0000
JC69	55	-1348.0771	2806.1542	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1310.6580	2967.4796	0.4058	0.4058
HKYG	60	-1302.0030	2968.3887	0.2576	0.6633
K2PI	57	-1311.4561	2969.0757	0.1827	0.8460
HKYI	60	-1302.8315	2970.0456	0.1125	0.9585
K2PIG	58	-1310.4274	2973.0914	0.0245	0.9830
HKYIG	61	-1301.8117	2974.0790	0.0150	0.9980
GTRG	64	-1295.2402	2979.1552	0.0012	0.9992
GTRI	64	-1296.0593	2980.7934	0.0005	0.9997
SYMG	61	-1306.5011	2983.4579	0.0001	0.9998
SYMI	61	-1307.1800	2984.8158	0.0001	0.9999
GTRIG	65	-1295.0652	2984.8783	0.0001	1.0000
HKY	59	-1315.2251	2988.7599	0.0000	1.0000
SYMIG	62	-1306.2943	2989.1173	0.0000	1.0000
K2P	56	-1324.9854	2990.0614	0.0000	1.0000
GTR	63	-1307.9683	2998.5383	0.0000	1.0000
SYM	60	-1320.2021	3004.7868	0.0000	1.0000
JC69G	56	-1333.9939	3008.0783	0.0000	1.0000
JC69I	56	-1334.7903	3009.6710	0.0000	1.0000
F81G	59	-1326.2504	3010.8104	0.0000	1.0000
F81I	59	-1327.0533	3012.4163	0.0000	1.0000
JC69IG	57	-1333.7597	3013.6830	0.0000	1.0000
F81IG	60	-1326.0387	3016.4601	0.0000	1.0000
JC69	55	-1348.0771	3030.1716	0.0000	1.0000
F81	58	-1339.7169	3031.6703	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/12814_hg18.phylip = 1 - 434;
 Partition Dummy = 1:261-aggressive-phylip/12814_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/12814_hg18.phylip = 1 - 434;
 Partition Dummy = 1:261-aggressive-phylip/12814_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/12814_hg18.phylip = 1 - 434;
 Partition Dummy = 1:261-aggressive-phylip/12814_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

