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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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
GTRG	64	-1174.7207	2506.7373	0.9674	0.9674
GTRIG	65	-1176.6548	2513.6276	0.0309	0.9983
GTRI	64	-1181.2219	2519.7396	0.0015	0.9997
HKYG	60	-1188.9430	2523.3026	0.0002	0.9999
HKYIG	61	-1189.2069	2526.7692	0.0000	1.0000
HKYI	60	-1192.3970	2530.2107	0.0000	1.0000
GTR	63	-1190.7954	2535.8855	0.0000	1.0000
HKY	59	-1205.8030	2554.1042	0.0000	1.0000
K2PG	57	-1227.7261	2592.1739	0.0000	1.0000
F81G	59	-1225.5511	2593.6004	0.0000	1.0000
K2PIG	58	-1227.1335	2593.8670	0.0000	1.0000
F81IG	60	-1225.9692	2597.3550	0.0000	1.0000
SYMG	61	-1224.6102	2597.5757	0.0000	1.0000
K2PI	57	-1230.6557	2598.0330	0.0000	1.0000
SYMIG	62	-1224.1489	2599.6125	0.0000	1.0000
F81I	59	-1229.3742	2601.2466	0.0000	1.0000
SYMI	61	-1228.0301	2604.4156	0.0000	1.0000
F81	58	-1241.0310	2621.6620	0.0000	1.0000
K2P	56	-1246.4898	2626.8427	0.0000	1.0000
SYM	60	-1242.8010	2631.0186	0.0000	1.0000
JC69G	56	-1258.9692	2651.8015	0.0000	1.0000
JC69IG	57	-1258.3660	2653.4536	0.0000	1.0000
JC69I	56	-1261.9528	2657.7686	0.0000	1.0000
JC69	55	-1277.4032	2685.8303	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1174.7207	2477.4415	0.9481	0.9481
GTRIG	65	-1176.6548	2483.3096	0.0504	0.9985
GTRI	64	-1181.2219	2490.4438	0.0014	1.0000
HKYG	60	-1188.9430	2497.8860	0.0000	1.0000
HKYIG	61	-1189.2069	2500.4138	0.0000	1.0000
HKYI	60	-1192.3970	2504.7940	0.0000	1.0000
GTR	63	-1190.7954	2507.5908	0.0000	1.0000
HKY	59	-1205.8030	2529.6059	0.0000	1.0000
F81G	59	-1225.5511	2569.1021	0.0000	1.0000
K2PG	57	-1227.7261	2569.4523	0.0000	1.0000
K2PIG	58	-1227.1335	2570.2670	0.0000	1.0000
SYMG	61	-1224.6102	2571.2203	0.0000	1.0000
F81IG	60	-1225.9692	2571.9384	0.0000	1.0000
SYMIG	62	-1224.1489	2572.2978	0.0000	1.0000
K2PI	57	-1230.6557	2575.3113	0.0000	1.0000
F81I	59	-1229.3742	2576.7484	0.0000	1.0000
SYMI	61	-1228.0301	2578.0602	0.0000	1.0000
F81	58	-1241.0310	2598.0620	0.0000	1.0000
K2P	56	-1246.4898	2604.9797	0.0000	1.0000
SYM	60	-1242.8010	2605.6019	0.0000	1.0000
JC69G	56	-1258.9692	2629.9384	0.0000	1.0000
JC69IG	57	-1258.3660	2630.7320	0.0000	1.0000
JC69I	56	-1261.9528	2635.9056	0.0000	1.0000
JC69	55	-1277.4032	2664.8064	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-1174.7207	2724.1661	0.9122	0.9122
HKYG	60	-1188.9430	2729.1903	0.0740	0.9862
GTRIG	65	-1176.6548	2733.8892	0.0071	0.9932
HKYIG	61	-1189.2069	2735.5732	0.0030	0.9963
HKYI	60	-1192.3970	2736.0983	0.0023	0.9986
GTRI	64	-1181.2219	2737.1684	0.0014	1.0000
GTR	63	-1190.7954	2750.4603	0.0000	1.0000
HKY	59	-1205.8030	2757.0551	0.0000	1.0000
K2PG	57	-1227.7261	2789.1914	0.0000	1.0000
K2PIG	58	-1227.1335	2793.8611	0.0000	1.0000
K2PI	57	-1230.6557	2795.0504	0.0000	1.0000
F81G	59	-1225.5511	2796.5514	0.0000	1.0000
F81IG	60	-1225.9692	2803.2427	0.0000	1.0000
F81I	59	-1229.3742	2804.1976	0.0000	1.0000
SYMG	61	-1224.6102	2806.3797	0.0000	1.0000
SYMIG	62	-1224.1489	2811.3123	0.0000	1.0000
SYMI	61	-1228.0301	2813.2195	0.0000	1.0000
K2P	56	-1246.4898	2820.8637	0.0000	1.0000
F81	58	-1241.0310	2821.6561	0.0000	1.0000
SYM	60	-1242.8010	2836.9063	0.0000	1.0000
JC69G	56	-1258.9692	2845.8225	0.0000	1.0000
JC69IG	57	-1258.3660	2850.4711	0.0000	1.0000
JC69I	56	-1261.9528	2851.7896	0.0000	1.0000
JC69	55	-1277.4032	2876.8354	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/12508_hg18.phylip = 1 - 349;
 Partition Dummy = 1:261-aggressive-phylip/12508_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 (GTRG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/12508_hg18.phylip = 1 - 349;
 Partition Dummy = 1:261-aggressive-phylip/12508_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 BIC model (GTRG)]

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


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

