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

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

Minimum AIC  model: GTRG
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	-1054.5339	2256.2796	0.4663	0.4663
GTRG	64	-1049.2314	2257.8590	0.2117	0.6780
HKYI	60	-1055.7254	2258.6627	0.1416	0.8197
GTRIG	65	-1048.8378	2260.1756	0.0665	0.8861
GTRI	64	-1050.5036	2260.4034	0.0593	0.9455
HKYIG	61	-1055.3695	2260.9629	0.0448	0.9903
SYMG	61	-1057.5149	2265.2536	0.0052	0.9956
K2PG	57	-1063.9762	2266.2612	0.0032	0.9987
K2PI	57	-1065.2201	2268.7490	0.0009	0.9996
K2PIG	58	-1065.2276	2271.7098	0.0002	0.9998
SYMI	61	-1061.5598	2273.3435	0.0001	0.9999
SYMIG	62	-1060.4167	2274.0918	0.0001	1.0000
HKY	59	-1077.9856	2300.1934	0.0000	1.0000
F81G	59	-1077.9960	2300.2142	0.0000	1.0000
JC69G	56	-1082.8613	2301.1073	0.0000	1.0000
GTR	63	-1072.6300	2301.5758	0.0000	1.0000
F81I	59	-1079.0603	2302.3427	0.0000	1.0000
JC69I	56	-1083.9205	2303.2256	0.0000	1.0000
JC69IG	57	-1082.8613	2304.0315	0.0000	1.0000
SYM	60	-1081.3038	2309.8196	0.0000	1.0000
K2P	56	-1088.1665	2311.7177	0.0000	1.0000
F81IG	60	-1082.4437	2312.0992	0.0000	1.0000
F81	58	-1101.0020	2343.2586	0.0000	1.0000
JC69	55	-1106.3620	2345.2058	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1049.2314	2226.4628	0.4486	0.4486
GTRIG	65	-1048.8378	2227.6756	0.2446	0.6932
GTRI	64	-1050.5036	2229.0072	0.1257	0.8189
HKYG	60	-1054.5339	2229.0677	0.1220	0.9409
HKYI	60	-1055.7254	2231.4508	0.0370	0.9779
HKYIG	61	-1055.3695	2232.7391	0.0195	0.9974
SYMG	61	-1057.5149	2237.0298	0.0023	0.9996
K2PG	57	-1063.9762	2241.9524	0.0002	0.9998
K2PI	57	-1065.2201	2244.4402	0.0001	0.9999
SYMIG	62	-1060.4167	2244.8333	0.0000	0.9999
SYMI	61	-1061.5598	2245.1196	0.0000	1.0000
K2PIG	58	-1065.2276	2246.4551	0.0000	1.0000
GTR	63	-1072.6300	2271.2601	0.0000	1.0000
HKY	59	-1077.9856	2273.9712	0.0000	1.0000
F81G	59	-1077.9960	2273.9920	0.0000	1.0000
F81I	59	-1079.0603	2276.1205	0.0000	1.0000
JC69G	56	-1082.8613	2277.7227	0.0000	1.0000
JC69IG	57	-1082.8613	2279.7226	0.0000	1.0000
JC69I	56	-1083.9205	2279.8409	0.0000	1.0000
SYM	60	-1081.3038	2282.6077	0.0000	1.0000
F81IG	60	-1082.4437	2284.8873	0.0000	1.0000
K2P	56	-1088.1665	2288.3331	0.0000	1.0000
F81	58	-1101.0020	2318.0040	0.0000	1.0000
JC69	55	-1106.3620	2322.7241	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1054.5339	2457.0133	0.5121	0.5121
K2PG	57	-1063.9762	2458.5007	0.2434	0.7555
HKYI	60	-1055.7254	2459.3963	0.1555	0.9110
K2PI	57	-1065.2201	2460.9885	0.0702	0.9812
HKYIG	61	-1055.3695	2464.4837	0.0122	0.9934
K2PIG	58	-1065.2276	2466.8025	0.0038	0.9973
SYMG	61	-1057.5149	2468.7744	0.0014	0.9987
GTRG	64	-1049.2314	2469.6047	0.0009	0.9996
GTRI	64	-1050.5036	2472.1491	0.0003	0.9999
GTRIG	65	-1048.8378	2474.6166	0.0001	1.0000
SYMI	61	-1061.5598	2476.8643	0.0000	1.0000
SYMIG	62	-1060.4167	2480.3771	0.0000	1.0000
JC69G	56	-1082.8613	2490.4718	0.0000	1.0000
JC69I	56	-1083.9205	2492.5901	0.0000	1.0000
JC69IG	57	-1082.8613	2496.2709	0.0000	1.0000
HKY	59	-1077.9856	2498.1177	0.0000	1.0000
F81G	59	-1077.9960	2498.1384	0.0000	1.0000
F81I	59	-1079.0603	2500.2670	0.0000	1.0000
K2P	56	-1088.1665	2501.0823	0.0000	1.0000
SYM	60	-1081.3038	2510.5532	0.0000	1.0000
GTR	63	-1072.6300	2510.6029	0.0000	1.0000
F81IG	60	-1082.4437	2512.8329	0.0000	1.0000
JC69	55	-1106.3620	2531.6742	0.0000	1.0000
F81	58	-1101.0020	2538.3514	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/11540_hg18.phylip = 1 - 330;
 Partition Dummy = 1:261-aggressive-phylip/11540_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 (HKYG)]

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

