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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: K2PG


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	-1121.4864	2396.3974	0.6878	0.6878
HKYIG	61	-1121.1023	2398.9019	0.1966	0.8844
HKYI	60	-1124.2590	2401.9428	0.0430	0.9274
K2PG	57	-1129.3190	2402.4217	0.0338	0.9612
K2PIG	58	-1128.4227	2403.8136	0.0169	0.9781
GTRIG	65	-1117.4563	2405.0061	0.0093	0.9873
GTRG	64	-1119.2754	2405.2485	0.0082	0.9956
K2PI	57	-1131.9550	2407.6939	0.0024	0.9980
GTRI	64	-1121.6343	2409.9663	0.0008	0.9988
SYMG	61	-1126.7392	2410.1756	0.0007	0.9995
SYMIG	62	-1125.4403	2410.8806	0.0005	1.0000
SYMI	61	-1130.0207	2416.7387	0.0000	1.0000
F81G	59	-1172.3171	2494.8159	0.0000	1.0000
HKY	59	-1173.2590	2496.6998	0.0000	1.0000
F81IG	60	-1172.2523	2497.9292	0.0000	1.0000
JC69G	56	-1178.6955	2498.0189	0.0000	1.0000
F81I	59	-1174.9251	2500.0319	0.0000	1.0000
JC69IG	57	-1178.6058	2500.9954	0.0000	1.0000
JC69I	56	-1181.1644	2502.9567	0.0000	1.0000
GTR	63	-1169.8896	2503.1126	0.0000	1.0000
K2P	56	-1182.1205	2504.8688	0.0000	1.0000
SYM	60	-1177.7041	2508.8328	0.0000	1.0000
F81	58	-1221.6571	2590.2826	0.0000	1.0000
JC69	55	-1227.7961	2593.0922	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1121.4864	2362.9727	0.4577	0.4577
HKYIG	61	-1121.1023	2364.2047	0.2472	0.7049
GTRIG	65	-1117.4563	2364.9126	0.1735	0.8784
GTRG	64	-1119.2754	2366.5508	0.0765	0.9549
HKYI	60	-1124.2590	2368.5181	0.0286	0.9835
GTRI	64	-1121.6343	2371.2686	0.0072	0.9907
K2PG	57	-1129.3190	2372.6379	0.0036	0.9944
K2PIG	58	-1128.4227	2372.8453	0.0033	0.9976
SYMIG	62	-1125.4403	2374.8806	0.0012	0.9988
SYMG	61	-1126.7392	2375.4784	0.0009	0.9997
K2PI	57	-1131.9550	2377.9101	0.0003	1.0000
SYMI	61	-1130.0207	2382.0415	0.0000	1.0000
F81G	59	-1172.3171	2462.6341	0.0000	1.0000
F81IG	60	-1172.2523	2464.5046	0.0000	1.0000
HKY	59	-1173.2590	2464.5180	0.0000	1.0000
GTR	63	-1169.8896	2465.7793	0.0000	1.0000
F81I	59	-1174.9251	2467.8501	0.0000	1.0000
JC69G	56	-1178.6955	2469.3911	0.0000	1.0000
JC69IG	57	-1178.6058	2471.2117	0.0000	1.0000
JC69I	56	-1181.1644	2474.3289	0.0000	1.0000
SYM	60	-1177.7041	2475.4081	0.0000	1.0000
K2P	56	-1182.1205	2476.2410	0.0000	1.0000
F81	58	-1221.6571	2559.3143	0.0000	1.0000
JC69	55	-1227.7961	2565.5922	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1129.3190	2579.8209	0.5442	0.5442
HKYG	60	-1121.4864	2581.0601	0.2929	0.8371
K2PIG	58	-1128.4227	2583.6631	0.0797	0.9168
K2PI	57	-1131.9550	2585.0931	0.0390	0.9558
HKYIG	61	-1121.1023	2585.9269	0.0257	0.9815
HKYI	60	-1124.2590	2586.6055	0.0183	0.9998
SYMG	61	-1126.7392	2597.2006	0.0001	0.9999
GTRG	64	-1119.2754	2599.1774	0.0000	1.0000
SYMIG	62	-1125.4403	2600.2376	0.0000	1.0000
GTRIG	65	-1117.4563	2601.1740	0.0000	1.0000
SYMI	61	-1130.0207	2603.7637	0.0000	1.0000
GTRI	64	-1121.6343	2603.8952	0.0000	1.0000
JC69G	56	-1178.6955	2672.9393	0.0000	1.0000
F81G	59	-1172.3171	2677.0867	0.0000	1.0000
JC69I	56	-1181.1644	2677.8771	0.0000	1.0000
JC69IG	57	-1178.6058	2678.3947	0.0000	1.0000
HKY	59	-1173.2590	2678.9706	0.0000	1.0000
K2P	56	-1182.1205	2679.7892	0.0000	1.0000
F81I	59	-1174.9251	2682.3027	0.0000	1.0000
F81IG	60	-1172.2523	2682.5920	0.0000	1.0000
SYM	60	-1177.7041	2693.4955	0.0000	1.0000
GTR	63	-1169.8896	2694.7710	0.0000	1.0000
JC69	55	-1227.7961	2765.5056	0.0000	1.0000
F81	58	-1221.6571	2770.1321	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/1446_hg18.phylip = 1 - 280;
 Partition Dummy = 1:261-aggressive-phylip/1446_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/1446_hg18.phylip = 1 - 280;
 Partition Dummy = 1:261-aggressive-phylip/1446_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 (K2PG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/1446_hg18.phylip = 1 - 280;
 Partition Dummy = 1:261-aggressive-phylip/1446_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

