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

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

Minimum AIC  model: GTRIG
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	-1170.9660	2481.7157	0.3979	0.3979
HKYIG	61	-1170.4502	2483.3991	0.1715	0.5694
GTRG	64	-1166.4900	2483.7123	0.1466	0.7160
HKYI	60	-1172.0214	2483.8265	0.1385	0.8545
GTRIG	65	-1165.3301	2484.1671	0.1168	0.9713
GTRI	64	-1168.1214	2486.9751	0.0287	1.0000
HKY	59	-1186.7021	2510.4877	0.0000	1.0000
GTR	63	-1183.5143	2515.0013	0.0000	1.0000
F81G	59	-1194.8927	2526.8689	0.0000	1.0000
F81IG	60	-1194.3054	2528.3946	0.0000	1.0000
F81I	59	-1195.8041	2528.6918	0.0000	1.0000
F81	58	-1210.4967	2555.3912	0.0000	1.0000
K2PG	57	-1220.1919	2572.1104	0.0000	1.0000
K2PIG	58	-1219.0540	2572.5059	0.0000	1.0000
K2PI	57	-1221.1209	2573.9684	0.0000	1.0000
SYMG	61	-1219.0535	2580.6056	0.0000	1.0000
SYMIG	62	-1217.8593	2580.9469	0.0000	1.0000
SYMI	61	-1220.0536	2582.6058	0.0000	1.0000
K2P	56	-1236.4745	2602.0186	0.0000	1.0000
JC69G	56	-1236.9692	2603.0080	0.0000	1.0000
JC69IG	57	-1236.0932	2603.9130	0.0000	1.0000
JC69I	56	-1237.8592	2604.7879	0.0000	1.0000
SYM	60	-1235.4331	2610.6500	0.0000	1.0000
JC69	55	-1253.0854	2632.5975	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-1165.3301	2460.6602	0.3269	0.3269
GTRG	64	-1166.4900	2460.9801	0.2786	0.6055
HKYG	60	-1170.9660	2461.9319	0.1731	0.7786
HKYIG	61	-1170.4502	2462.9004	0.1067	0.8852
HKYI	60	-1172.0214	2464.0427	0.0602	0.9455
GTRI	64	-1168.1214	2464.2429	0.0545	1.0000
HKY	59	-1186.7021	2491.4042	0.0000	1.0000
GTR	63	-1183.5143	2493.0285	0.0000	1.0000
F81G	59	-1194.8927	2507.7854	0.0000	1.0000
F81IG	60	-1194.3054	2508.6108	0.0000	1.0000
F81I	59	-1195.8041	2509.6083	0.0000	1.0000
F81	58	-1210.4967	2536.9933	0.0000	1.0000
K2PIG	58	-1219.0540	2554.1081	0.0000	1.0000
K2PG	57	-1220.1919	2554.3839	0.0000	1.0000
K2PI	57	-1221.1209	2556.2419	0.0000	1.0000
SYMIG	62	-1217.8593	2559.7187	0.0000	1.0000
SYMG	61	-1219.0535	2560.1070	0.0000	1.0000
SYMI	61	-1220.0536	2562.1071	0.0000	1.0000
K2P	56	-1236.4745	2584.9491	0.0000	1.0000
JC69G	56	-1236.9692	2585.9384	0.0000	1.0000
JC69IG	57	-1236.0932	2586.1864	0.0000	1.0000
JC69I	56	-1237.8592	2587.7183	0.0000	1.0000
SYM	60	-1235.4331	2590.8662	0.0000	1.0000
JC69	55	-1253.0854	2616.1708	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1170.9660	2705.8984	0.6996	0.6996
HKYI	60	-1172.0214	2708.0092	0.2435	0.9431
HKYIG	61	-1170.4502	2710.9330	0.0564	0.9996
GTRG	64	-1166.4900	2721.2110	0.0003	0.9999
GTRI	64	-1168.1214	2724.4738	0.0001	0.9999
GTRIG	65	-1165.3301	2724.9572	0.0001	1.0000
HKY	59	-1186.7021	2731.3045	0.0000	1.0000
F81G	59	-1194.8927	2747.6857	0.0000	1.0000
GTR	63	-1183.5143	2749.1933	0.0000	1.0000
F81I	59	-1195.8041	2749.5086	0.0000	1.0000
F81IG	60	-1194.3054	2752.5773	0.0000	1.0000
F81	58	-1210.4967	2772.8276	0.0000	1.0000
K2PG	57	-1220.1919	2786.1520	0.0000	1.0000
K2PI	57	-1221.1209	2788.0101	0.0000	1.0000
K2PIG	58	-1219.0540	2789.9423	0.0000	1.0000
SYMG	61	-1219.0535	2808.1396	0.0000	1.0000
SYMI	61	-1220.0536	2810.1397	0.0000	1.0000
SYMIG	62	-1217.8593	2811.8174	0.0000	1.0000
K2P	56	-1236.4745	2812.6511	0.0000	1.0000
JC69G	56	-1236.9692	2813.6405	0.0000	1.0000
JC69I	56	-1237.8592	2815.4204	0.0000	1.0000
JC69IG	57	-1236.0932	2817.9546	0.0000	1.0000
SYM	60	-1235.4331	2834.8327	0.0000	1.0000
JC69	55	-1253.0854	2839.8067	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/3430_hg18.phylip = 1 - 431;
 Partition Dummy = 1:261-aggressive-phylip/3430_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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

