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

Input data from file "261-aggressive-phylip/7445_hg18.phylip" (29 taxa, 310 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	-1086.2916	2321.9808	0.9685	0.9685
HKYIG	61	-1088.8400	2330.1801	0.0161	0.9846
HKYI	60	-1090.5629	2330.5233	0.0135	0.9981
GTRG	64	-1087.2399	2336.4391	0.0007	0.9988
GTRIG	65	-1085.7291	2336.6222	0.0006	0.9994
GTRI	64	-1088.0255	2338.0102	0.0003	0.9998
K2PG	57	-1100.0508	2340.3397	0.0001	0.9999
K2PIG	58	-1098.8831	2341.0331	0.0001	0.9999
K2PI	57	-1100.4696	2341.1774	0.0001	1.0000
SYMG	61	-1097.2227	2346.9455	0.0000	1.0000
SYMIG	62	-1095.7399	2347.1074	0.0000	1.0000
SYMI	61	-1097.9292	2348.3584	0.0000	1.0000
F81G	59	-1109.9003	2366.1207	0.0000	1.0000
F81I	59	-1110.2764	2366.8728	0.0000	1.0000
F81IG	60	-1108.9755	2367.3486	0.0000	1.0000
JC69G	56	-1120.3885	2378.0102	0.0000	1.0000
JC69I	56	-1120.7739	2378.7810	0.0000	1.0000
JC69IG	57	-1119.4587	2379.1554	0.0000	1.0000
HKY	59	-1129.7418	2405.8037	0.0000	1.0000
GTR	63	-1129.1691	2417.1188	0.0000	1.0000
K2P	56	-1142.0195	2421.2723	0.0000	1.0000
SYM	60	-1140.1479	2429.6933	0.0000	1.0000
F81	58	-1155.0101	2453.2871	0.0000	1.0000
JC69	55	-1165.2421	2464.7361	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1086.2916	2292.5832	0.9390	0.9390
HKYIG	61	-1088.8400	2299.6801	0.0270	0.9660
HKYI	60	-1090.5629	2301.1257	0.0131	0.9791
GTRIG	65	-1085.7291	2301.4583	0.0111	0.9902
GTRG	64	-1087.2399	2302.4799	0.0067	0.9969
GTRI	64	-1088.0255	2304.0510	0.0030	0.9999
K2PIG	58	-1098.8831	2313.7661	0.0000	0.9999
K2PG	57	-1100.0508	2314.1016	0.0000	1.0000
K2PI	57	-1100.4696	2314.9393	0.0000	1.0000
SYMIG	62	-1095.7399	2315.4798	0.0000	1.0000
SYMG	61	-1097.2227	2316.4455	0.0000	1.0000
SYMI	61	-1097.9292	2317.8584	0.0000	1.0000
F81G	59	-1109.9003	2337.8007	0.0000	1.0000
F81IG	60	-1108.9755	2337.9510	0.0000	1.0000
F81I	59	-1110.2764	2338.5528	0.0000	1.0000
JC69G	56	-1120.3885	2352.7770	0.0000	1.0000
JC69IG	57	-1119.4587	2352.9173	0.0000	1.0000
JC69I	56	-1120.7739	2353.5478	0.0000	1.0000
HKY	59	-1129.7418	2377.4837	0.0000	1.0000
GTR	63	-1129.1691	2384.3383	0.0000	1.0000
K2P	56	-1142.0195	2396.0391	0.0000	1.0000
SYM	60	-1140.1479	2400.2957	0.0000	1.0000
F81	58	-1155.0101	2426.0202	0.0000	1.0000
JC69	55	-1165.2421	2440.4841	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1086.2916	2516.7775	0.9718	0.9718
HKYI	60	-1090.5629	2525.3201	0.0136	0.9853
K2PG	57	-1100.0508	2527.0862	0.0056	0.9910
HKYIG	61	-1088.8400	2527.6110	0.0043	0.9953
K2PI	57	-1100.4696	2527.9239	0.0037	0.9990
K2PIG	58	-1098.8831	2530.4873	0.0010	1.0000
GTRG	64	-1087.2399	2541.6205	0.0000	1.0000
GTRI	64	-1088.0255	2543.1917	0.0000	1.0000
GTRIG	65	-1085.7291	2544.3355	0.0000	1.0000
SYMG	61	-1097.2227	2544.3764	0.0000	1.0000
SYMI	61	-1097.9292	2545.7894	0.0000	1.0000
SYMIG	62	-1095.7399	2547.1473	0.0000	1.0000
F81G	59	-1109.9003	2558.2585	0.0000	1.0000
F81I	59	-1110.2764	2559.0105	0.0000	1.0000
JC69G	56	-1120.3885	2562.0250	0.0000	1.0000
F81IG	60	-1108.9755	2562.1453	0.0000	1.0000
JC69I	56	-1120.7739	2562.7959	0.0000	1.0000
JC69IG	57	-1119.4587	2565.9019	0.0000	1.0000
HKY	59	-1129.7418	2597.9414	0.0000	1.0000
K2P	56	-1142.0195	2605.2871	0.0000	1.0000
GTR	63	-1129.1691	2619.7423	0.0000	1.0000
SYM	60	-1140.1479	2624.4900	0.0000	1.0000
F81	58	-1155.0101	2642.7414	0.0000	1.0000
JC69	55	-1165.2421	2645.9956	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/7445_hg18.phylip = 1 - 310;
 Partition Dummy = 1:261-aggressive-phylip/7445_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/7445_hg18.phylip = 1 - 310;
 Partition Dummy = 1:261-aggressive-phylip/7445_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/7445_hg18.phylip = 1 - 310;
 Partition Dummy = 1:261-aggressive-phylip/7445_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

