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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-1393.5667	2920.2752	0.7591	0.7591
HKYIG	61	-1393.5768	2922.7580	0.2194	0.9784
GTRG	64	-1392.4259	2927.8970	0.0168	0.9952
GTRIG	65	-1392.4361	2930.4156	0.0048	1.0000
K2PG	57	-1409.7985	2945.4041	0.0000	1.0000
K2PIG	58	-1409.8087	2947.8607	0.0000	1.0000
SYMG	61	-1408.7359	2953.0761	0.0000	1.0000
SYMIG	62	-1408.7461	2955.5679	0.0000	1.0000
HKY	59	-1414.5379	2959.7641	0.0000	1.0000
F81G	59	-1415.6259	2961.9400	0.0000	1.0000
HKYI	60	-1414.5449	2962.2316	0.0000	1.0000
F81IG	60	-1415.6361	2964.4140	0.0000	1.0000
GTR	63	-1413.4857	2967.5273	0.0000	1.0000
GTRI	64	-1413.4928	2970.0307	0.0000	1.0000
JC69G	56	-1429.8403	2983.0602	0.0000	1.0000
JC69IG	57	-1429.8505	2985.5081	0.0000	1.0000
K2P	56	-1431.1926	2985.7649	0.0000	1.0000
K2PI	57	-1431.1997	2988.2066	0.0000	1.0000
SYM	60	-1430.7200	2994.5818	0.0000	1.0000
SYMI	61	-1430.7271	2997.0585	0.0000	1.0000
F81	58	-1436.2374	3000.7180	0.0000	1.0000
F81I	59	-1436.2444	3003.1770	0.0000	1.0000
JC69	55	-1451.1316	3023.2241	0.0000	1.0000
JC69I	56	-1451.1387	3025.6570	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1393.5667	2907.1334	0.6933	0.6933
HKYIG	61	-1393.5768	2909.1537	0.2525	0.9458
GTRG	64	-1392.4259	2912.8518	0.0397	0.9855
GTRIG	65	-1392.4361	2914.8721	0.0145	1.0000
K2PG	57	-1409.7985	2933.5969	0.0000	1.0000
K2PIG	58	-1409.8087	2935.6174	0.0000	1.0000
SYMG	61	-1408.7359	2939.4718	0.0000	1.0000
SYMIG	62	-1408.7461	2941.4922	0.0000	1.0000
HKY	59	-1414.5379	2947.0759	0.0000	1.0000
HKYI	60	-1414.5449	2949.0898	0.0000	1.0000
F81G	59	-1415.6259	2949.2518	0.0000	1.0000
F81IG	60	-1415.6361	2951.2722	0.0000	1.0000
GTR	63	-1413.4857	2952.9713	0.0000	1.0000
GTRI	64	-1413.4928	2954.9855	0.0000	1.0000
JC69G	56	-1429.8403	2971.6806	0.0000	1.0000
JC69IG	57	-1429.8505	2973.7010	0.0000	1.0000
K2P	56	-1431.1926	2974.3852	0.0000	1.0000
K2PI	57	-1431.1997	2976.3995	0.0000	1.0000
SYM	60	-1430.7200	2981.4400	0.0000	1.0000
SYMI	61	-1430.7271	2983.4542	0.0000	1.0000
F81	58	-1436.2374	2988.4747	0.0000	1.0000
F81I	59	-1436.2444	2990.4889	0.0000	1.0000
JC69	55	-1451.1316	3012.2633	0.0000	1.0000
JC69I	56	-1451.1387	3014.2773	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1393.5667	3172.7227	0.9604	0.9604
HKYIG	61	-1393.5768	3179.1695	0.0382	0.9986
K2PG	57	-1409.7985	3185.9068	0.0013	0.9999
K2PIG	58	-1409.8087	3192.3537	0.0001	1.0000
GTRG	64	-1392.4259	3196.1470	0.0000	1.0000
GTRIG	65	-1392.4361	3202.5939	0.0000	1.0000
HKY	59	-1414.5379	3208.2387	0.0000	1.0000
SYMG	61	-1408.7359	3209.4876	0.0000	1.0000
F81G	59	-1415.6259	3210.4146	0.0000	1.0000
HKYI	60	-1414.5449	3214.6791	0.0000	1.0000
SYMIG	62	-1408.7461	3215.9345	0.0000	1.0000
F81IG	60	-1415.6361	3216.8615	0.0000	1.0000
JC69G	56	-1429.8403	3219.5639	0.0000	1.0000
K2P	56	-1431.1926	3222.2686	0.0000	1.0000
JC69IG	57	-1429.8505	3226.0108	0.0000	1.0000
K2PI	57	-1431.1997	3228.7093	0.0000	1.0000
GTR	63	-1413.4857	3231.8401	0.0000	1.0000
GTRI	64	-1413.4928	3238.2808	0.0000	1.0000
F81	58	-1436.2374	3245.2111	0.0000	1.0000
SYM	60	-1430.7200	3247.0293	0.0000	1.0000
F81I	59	-1436.2444	3251.6517	0.0000	1.0000
SYMI	61	-1430.7271	3253.4700	0.0000	1.0000
JC69	55	-1451.1316	3255.7201	0.0000	1.0000
JC69I	56	-1451.1387	3262.1607	0.0000	1.0000



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

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

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


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

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


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

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

