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

Input data from file "261-aggressive-phylip/9099_hg18.phylip" (29 taxa, 684 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	-1239.3632	2620.1675	0.6277	0.6277
GTRIG	65	-1239.3723	2622.6280	0.1834	0.8111
HKYG	60	-1245.6700	2623.0895	0.1456	0.9568
HKYIG	61	-1245.6789	2625.5186	0.0432	1.0000
GTR	63	-1255.7983	2650.6031	0.0000	1.0000
GTRI	64	-1255.8060	2653.0530	0.0000	1.0000
HKY	59	-1262.4118	2654.1697	0.0000	1.0000
HKYI	60	-1262.4196	2656.5888	0.0000	1.0000
F81G	59	-1266.3519	2662.0499	0.0000	1.0000
F81IG	60	-1266.3606	2664.4707	0.0000	1.0000
K2PG	57	-1276.8209	2678.2040	0.0000	1.0000
SYMG	61	-1273.3893	2680.9393	0.0000	1.0000
K2PIG	58	-1278.6081	2684.1667	0.0000	1.0000
SYMIG	62	-1275.1503	2686.8803	0.0000	1.0000
F81	58	-1283.1488	2693.2480	0.0000	1.0000
F81I	59	-1283.1565	2695.6591	0.0000	1.0000
K2P	56	-1294.5212	2711.2243	0.0000	1.0000
SYM	60	-1290.4227	2712.5949	0.0000	1.0000
JC69G	56	-1295.8288	2713.8394	0.0000	1.0000
K2PI	57	-1296.4074	2717.3772	0.0000	1.0000
SYMI	61	-1292.3035	2718.7678	0.0000	1.0000
JC69IG	57	-1297.6836	2719.9294	0.0000	1.0000
JC69	55	-1313.0843	2745.9776	0.0000	1.0000
JC69I	56	-1314.9859	2752.1537	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1239.3632	2606.7265	0.6665	0.6665
GTRIG	65	-1239.3723	2608.7445	0.2430	0.9094
HKYG	60	-1245.6700	2611.3399	0.0664	0.9758
HKYIG	61	-1245.6789	2613.3578	0.0242	1.0000
GTR	63	-1255.7983	2637.5966	0.0000	1.0000
GTRI	64	-1255.8060	2639.6119	0.0000	1.0000
HKY	59	-1262.4118	2642.8235	0.0000	1.0000
HKYI	60	-1262.4196	2644.8392	0.0000	1.0000
F81G	59	-1266.3519	2650.7037	0.0000	1.0000
F81IG	60	-1266.3606	2652.7211	0.0000	1.0000
K2PG	57	-1276.8209	2667.6417	0.0000	1.0000
SYMG	61	-1273.3893	2668.7785	0.0000	1.0000
K2PIG	58	-1278.6081	2673.2163	0.0000	1.0000
SYMIG	62	-1275.1503	2674.3006	0.0000	1.0000
F81	58	-1283.1488	2682.2976	0.0000	1.0000
F81I	59	-1283.1565	2684.3130	0.0000	1.0000
SYM	60	-1290.4227	2700.8453	0.0000	1.0000
K2P	56	-1294.5212	2701.0424	0.0000	1.0000
JC69G	56	-1295.8288	2703.6575	0.0000	1.0000
SYMI	61	-1292.3035	2706.6070	0.0000	1.0000
K2PI	57	-1296.4074	2706.8149	0.0000	1.0000
JC69IG	57	-1297.6836	2709.3671	0.0000	1.0000
JC69	55	-1313.0843	2736.1687	0.0000	1.0000
JC69I	56	-1314.9859	2741.9719	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1245.6700	2883.0174	0.9624	0.9624
HKYIG	61	-1245.6789	2889.5632	0.0365	0.9988
GTRG	64	-1239.3632	2896.5158	0.0011	1.0000
GTRIG	65	-1239.3723	2903.0618	0.0000	1.0000
HKY	59	-1262.4118	2909.9730	0.0000	1.0000
HKYI	60	-1262.4196	2916.5167	0.0000	1.0000
F81G	59	-1266.3519	2917.8532	0.0000	1.0000
GTR	63	-1255.7983	2922.8580	0.0000	1.0000
F81IG	60	-1266.3606	2924.3986	0.0000	1.0000
K2PG	57	-1276.8209	2925.7353	0.0000	1.0000
GTRI	64	-1255.8060	2929.4012	0.0000	1.0000
K2PIG	58	-1278.6081	2935.8378	0.0000	1.0000
F81	58	-1283.1488	2944.9191	0.0000	1.0000
SYMG	61	-1273.3893	2944.9840	0.0000	1.0000
F81I	59	-1283.1565	2951.4625	0.0000	1.0000
K2P	56	-1294.5212	2954.6081	0.0000	1.0000
SYMIG	62	-1275.1503	2955.0340	0.0000	1.0000
JC69G	56	-1295.8288	2957.2232	0.0000	1.0000
K2PI	57	-1296.4074	2964.9085	0.0000	1.0000
JC69IG	57	-1297.6836	2967.4607	0.0000	1.0000
SYM	60	-1290.4227	2972.5228	0.0000	1.0000
SYMI	61	-1292.3035	2982.8124	0.0000	1.0000
JC69	55	-1313.0843	2985.2064	0.0000	1.0000
JC69I	56	-1314.9859	2995.5375	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/9099_hg18.phylip = 1 - 684;
 Partition Dummy = 1:261-aggressive-phylip/9099_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/9099_hg18.phylip = 1 - 684;
 Partition Dummy = 1:261-aggressive-phylip/9099_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/9099_hg18.phylip = 1 - 684;
 Partition Dummy = 1:261-aggressive-phylip/9099_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

