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

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

Minimum AIC  model: GTRG
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	-1171.7732	2486.2793	0.2952	0.2952
HKYI	60	-1171.8716	2486.4761	0.2675	0.5627
GTRG	64	-1166.6343	2487.4321	0.1659	0.7286
GTRI	64	-1166.9038	2487.9711	0.1267	0.8552
HKYIG	61	-1171.7729	2489.1097	0.0717	0.9269
GTRIG	65	-1166.6342	2490.3346	0.0389	0.9658
K2PG	57	-1178.9963	2492.3373	0.0143	0.9801
K2PI	57	-1179.0484	2492.4415	0.0136	0.9936
K2PIG	58	-1178.9975	2495.1185	0.0036	0.9972
SYMG	61	-1175.7623	2497.0885	0.0013	0.9985
SYMI	61	-1175.8759	2497.3157	0.0012	0.9997
SYMIG	62	-1175.7702	2499.9529	0.0003	1.0000
HKY	59	-1188.3707	2516.6610	0.0000	1.0000
GTR	63	-1183.5359	2518.3508	0.0000	1.0000
F81G	59	-1190.7273	2521.3742	0.0000	1.0000
F81I	59	-1191.0355	2521.9905	0.0000	1.0000
K2P	56	-1195.2181	2522.0190	0.0000	1.0000
F81IG	60	-1190.7234	2524.1797	0.0000	1.0000
JC69G	56	-1196.5215	2524.6258	0.0000	1.0000
JC69I	56	-1196.8055	2525.1938	0.0000	1.0000
SYM	60	-1191.5912	2525.9154	0.0000	1.0000
JC69IG	57	-1196.5203	2527.3852	0.0000	1.0000
F81	58	-1206.9355	2550.9944	0.0000	1.0000
JC69	55	-1212.6075	2554.0530	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1166.6343	2461.2686	0.3477	0.3477
GTRI	64	-1166.9038	2461.8075	0.2656	0.6134
GTRIG	65	-1166.6342	2463.2684	0.1279	0.7413
HKYG	60	-1171.7732	2463.5463	0.1113	0.8526
HKYI	60	-1171.8716	2463.7432	0.1009	0.9535
HKYIG	61	-1171.7729	2465.5459	0.0410	0.9945
K2PG	57	-1178.9963	2471.9927	0.0016	0.9961
K2PI	57	-1179.0484	2472.0969	0.0015	0.9977
SYMG	61	-1175.7623	2473.5246	0.0008	0.9984
SYMI	61	-1175.8759	2473.7519	0.0007	0.9991
K2PIG	58	-1178.9975	2473.9951	0.0006	0.9997
SYMIG	62	-1175.7702	2475.5404	0.0003	1.0000
GTR	63	-1183.5359	2493.0718	0.0000	1.0000
HKY	59	-1188.3707	2494.7415	0.0000	1.0000
F81G	59	-1190.7273	2499.4547	0.0000	1.0000
F81I	59	-1191.0355	2500.0710	0.0000	1.0000
F81IG	60	-1190.7234	2501.4468	0.0000	1.0000
K2P	56	-1195.2181	2502.4362	0.0000	1.0000
SYM	60	-1191.5912	2503.1825	0.0000	1.0000
JC69G	56	-1196.5215	2505.0430	0.0000	1.0000
JC69I	56	-1196.8055	2505.6110	0.0000	1.0000
JC69IG	57	-1196.5203	2507.0406	0.0000	1.0000
F81	58	-1206.9355	2529.8710	0.0000	1.0000
JC69	55	-1212.6075	2535.2150	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1178.9963	2697.0307	0.4239	0.4239
K2PI	57	-1179.0484	2697.1349	0.4024	0.8263
HKYG	60	-1171.7732	2700.4284	0.0775	0.9038
HKYI	60	-1171.8716	2700.6253	0.0703	0.9741
K2PIG	58	-1178.9975	2702.9811	0.0216	0.9957
HKYIG	61	-1171.7729	2706.3760	0.0040	0.9997
GTRG	64	-1166.6343	2713.9428	0.0001	0.9998
SYMG	61	-1175.7623	2714.3548	0.0001	0.9999
GTRI	64	-1166.9038	2714.4818	0.0001	0.9999
SYMI	61	-1175.8759	2714.5820	0.0001	1.0000
GTRIG	65	-1166.6342	2719.8907	0.0000	1.0000
SYMIG	62	-1175.7702	2720.3186	0.0000	1.0000
K2P	56	-1195.2181	2723.5261	0.0000	1.0000
JC69G	56	-1196.5215	2726.1329	0.0000	1.0000
JC69I	56	-1196.8055	2726.7010	0.0000	1.0000
HKY	59	-1188.3707	2727.6755	0.0000	1.0000
JC69IG	57	-1196.5203	2732.0786	0.0000	1.0000
F81G	59	-1190.7273	2732.3887	0.0000	1.0000
F81I	59	-1191.0355	2733.0050	0.0000	1.0000
F81IG	60	-1190.7234	2738.3289	0.0000	1.0000
SYM	60	-1191.5912	2740.0646	0.0000	1.0000
GTR	63	-1183.5359	2741.7980	0.0000	1.0000
JC69	55	-1212.6075	2752.3570	0.0000	1.0000
F81	58	-1206.9355	2758.8570	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/6665_hg18.phylip = 1 - 383;
 Partition Dummy = 1:261-aggressive-phylip/6665_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/6665_hg18.phylip = 1 - 383;
 Partition Dummy = 1:261-aggressive-phylip/6665_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/6665_hg18.phylip = 1 - 383;
 Partition Dummy = 1:261-aggressive-phylip/6665_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

