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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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	-1274.0260	2710.0112	0.9594	0.9594
GTRIG	65	-1276.2054	2717.5747	0.0219	0.9812
GTRI	64	-1278.3871	2718.7334	0.0122	0.9935
HKYG	60	-1285.5704	2720.5384	0.0050	0.9985
HKYIG	61	-1285.5700	2723.6399	0.0011	0.9995
HKYI	60	-1287.9043	2725.2062	0.0005	1.0000
GTR	63	-1287.6054	2733.9914	0.0000	1.0000
HKY	59	-1296.1937	2738.7074	0.0000	1.0000
K2PG	57	-1306.9446	2754.1273	0.0000	1.0000
K2PIG	58	-1306.9440	2757.1550	0.0000	1.0000
K2PI	57	-1309.5909	2759.4198	0.0000	1.0000
SYMG	61	-1303.7654	2760.0309	0.0000	1.0000
SYMIG	62	-1303.7649	2763.1573	0.0000	1.0000
SYMI	61	-1306.3181	2765.1362	0.0000	1.0000
K2P	56	-1318.0939	2773.4210	0.0000	1.0000
SYM	60	-1314.7188	2778.8351	0.0000	1.0000
F81G	59	-1326.9623	2800.2445	0.0000	1.0000
F81IG	60	-1326.9630	2803.3236	0.0000	1.0000
F81I	59	-1329.3290	2804.9780	0.0000	1.0000
F81	58	-1337.5771	2818.4212	0.0000	1.0000
JC69G	56	-1342.5584	2822.3500	0.0000	1.0000
JC69IG	57	-1342.5613	2825.3607	0.0000	1.0000
JC69I	56	-1345.1542	2827.5417	0.0000	1.0000
JC69	55	-1352.5445	2839.3409	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1274.0260	2676.0520	0.9477	0.9477
GTRIG	65	-1276.2054	2682.4107	0.0394	0.9872
GTRI	64	-1278.3871	2684.7742	0.0121	0.9993
HKYG	60	-1285.5704	2691.1408	0.0005	0.9998
HKYIG	61	-1285.5700	2693.1399	0.0002	0.9999
HKYI	60	-1287.9043	2695.8086	0.0000	1.0000
GTR	63	-1287.6054	2701.2109	0.0000	1.0000
HKY	59	-1296.1937	2710.3874	0.0000	1.0000
K2PG	57	-1306.9446	2727.8892	0.0000	1.0000
SYMG	61	-1303.7654	2729.5309	0.0000	1.0000
K2PIG	58	-1306.9440	2729.8880	0.0000	1.0000
SYMIG	62	-1303.7649	2731.5298	0.0000	1.0000
K2PI	57	-1309.5909	2733.1817	0.0000	1.0000
SYMI	61	-1306.3181	2734.6362	0.0000	1.0000
K2P	56	-1318.0939	2748.1878	0.0000	1.0000
SYM	60	-1314.7188	2749.4375	0.0000	1.0000
F81G	59	-1326.9623	2771.9245	0.0000	1.0000
F81IG	60	-1326.9630	2773.9260	0.0000	1.0000
F81I	59	-1329.3290	2776.6580	0.0000	1.0000
F81	58	-1337.5771	2791.1542	0.0000	1.0000
JC69G	56	-1342.5584	2797.1168	0.0000	1.0000
JC69IG	57	-1342.5613	2799.1226	0.0000	1.0000
JC69I	56	-1345.1542	2802.3085	0.0000	1.0000
JC69	55	-1352.5445	2815.0889	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-1274.0260	2915.1927	0.4775	0.4775
HKYG	60	-1285.5704	2915.3351	0.4447	0.9223
HKYI	60	-1287.9043	2920.0029	0.0431	0.9654
HKYIG	61	-1285.5700	2921.0708	0.0253	0.9906
GTRI	64	-1278.3871	2923.9149	0.0061	0.9967
GTRIG	65	-1276.2054	2925.2879	0.0031	0.9998
HKY	59	-1296.1937	2930.8451	0.0002	1.0000
GTR	63	-1287.6054	2936.6149	0.0000	1.0000
K2PG	57	-1306.9446	2940.8738	0.0000	1.0000
K2PI	57	-1309.5909	2946.1663	0.0000	1.0000
K2PIG	58	-1306.9440	2946.6092	0.0000	1.0000
K2P	56	-1318.0939	2957.4358	0.0000	1.0000
SYMG	61	-1303.7654	2957.4618	0.0000	1.0000
SYMI	61	-1306.3181	2962.5671	0.0000	1.0000
SYMIG	62	-1303.7649	2963.1973	0.0000	1.0000
SYM	60	-1314.7188	2973.6319	0.0000	1.0000
F81G	59	-1326.9623	2992.3823	0.0000	1.0000
F81I	59	-1329.3290	2997.1157	0.0000	1.0000
F81IG	60	-1326.9630	2998.1203	0.0000	1.0000
JC69G	56	-1342.5584	3006.3648	0.0000	1.0000
F81	58	-1337.5771	3007.8754	0.0000	1.0000
JC69I	56	-1345.1542	3011.5565	0.0000	1.0000
JC69IG	57	-1342.5613	3012.1072	0.0000	1.0000
JC69	55	-1352.5445	3020.6004	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/9003_hg18.phylip = 1 - 310;
 Partition Dummy = 1:261-aggressive-phylip/9003_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/9003_hg18.phylip = 1 - 310;
 Partition Dummy = 1:261-aggressive-phylip/9003_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 (GTRG)]

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


-------------------------------------------------------------
End of Output

