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

Input data from file "261-aggressive-phylip/7829_hg18.phylip" (29 taxa, 413 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	-1448.5332	3037.8619	0.5585	0.5585
HKYIG	61	-1448.1017	3039.7533	0.2169	0.7754
GTRG	64	-1444.6136	3041.1353	0.1087	0.8841
HKYI	60	-1450.7103	3042.2160	0.0633	0.9475
GTRIG	65	-1444.0608	3042.8478	0.0462	0.9936
GTRI	64	-1447.4489	3046.8058	0.0064	1.0000
K2PG	57	-1467.1677	3066.9607	0.0000	1.0000
K2PIG	58	-1466.6991	3068.7315	0.0000	1.0000
K2PI	57	-1469.3106	3071.2465	0.0000	1.0000
SYMG	61	-1464.1169	3071.7837	0.0000	1.0000
SYMIG	62	-1463.6437	3073.6074	0.0000	1.0000
SYMI	61	-1466.4345	3076.4190	0.0000	1.0000
HKY	59	-1473.5911	3085.2388	0.0000	1.0000
GTR	63	-1469.9772	3089.0604	0.0000	1.0000
F81G	59	-1482.9470	3103.9507	0.0000	1.0000
F81IG	60	-1482.5083	3105.8120	0.0000	1.0000
F81I	59	-1485.1387	3108.3341	0.0000	1.0000
K2P	56	-1493.2706	3116.4738	0.0000	1.0000
SYM	60	-1490.9986	3122.7927	0.0000	1.0000
JC69G	56	-1497.7516	3125.4358	0.0000	1.0000
JC69IG	57	-1497.2521	3127.1295	0.0000	1.0000
JC69I	56	-1500.1111	3130.1547	0.0000	1.0000
F81	58	-1507.7857	3150.9047	0.0000	1.0000
JC69	55	-1523.4514	3174.1577	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1448.5332	3017.0665	0.3080	0.3080
GTRG	64	-1444.6136	3017.2273	0.2842	0.5922
GTRIG	65	-1444.0608	3018.1216	0.1817	0.7739
HKYIG	61	-1448.1017	3018.2034	0.1744	0.9484
HKYI	60	-1450.7103	3021.4205	0.0349	0.9833
GTRI	64	-1447.4489	3022.8978	0.0167	1.0000
K2PG	57	-1467.1677	3048.3354	0.0000	1.0000
K2PIG	58	-1466.6991	3049.3982	0.0000	1.0000
SYMG	61	-1464.1169	3050.2339	0.0000	1.0000
SYMIG	62	-1463.6437	3051.2874	0.0000	1.0000
K2PI	57	-1469.3106	3052.6211	0.0000	1.0000
SYMI	61	-1466.4345	3054.8691	0.0000	1.0000
HKY	59	-1473.5911	3065.1821	0.0000	1.0000
GTR	63	-1469.9772	3065.9544	0.0000	1.0000
F81G	59	-1482.9470	3083.8941	0.0000	1.0000
F81IG	60	-1482.5083	3085.0165	0.0000	1.0000
F81I	59	-1485.1387	3088.2775	0.0000	1.0000
K2P	56	-1493.2706	3098.5412	0.0000	1.0000
SYM	60	-1490.9986	3101.9972	0.0000	1.0000
JC69G	56	-1497.7516	3107.5032	0.0000	1.0000
JC69IG	57	-1497.2521	3108.5041	0.0000	1.0000
JC69I	56	-1500.1111	3112.2221	0.0000	1.0000
F81	58	-1507.7857	3131.5714	0.0000	1.0000
JC69	55	-1523.4514	3156.9028	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1448.5332	3258.4733	0.8406	0.8406
HKYI	60	-1450.7103	3262.8274	0.0953	0.9360
HKYIG	61	-1448.1017	3263.6337	0.0637	0.9996
GTRG	64	-1444.6136	3274.7279	0.0002	0.9999
K2PG	57	-1467.1677	3277.6719	0.0001	1.0000
GTRIG	65	-1444.0608	3279.6457	0.0000	1.0000
GTRI	64	-1447.4489	3280.3984	0.0000	1.0000
K2PI	57	-1469.3106	3281.9577	0.0000	1.0000
K2PIG	58	-1466.6991	3282.7581	0.0000	1.0000
SYMG	61	-1464.1169	3295.6642	0.0000	1.0000
SYMI	61	-1466.4345	3300.2994	0.0000	1.0000
SYMIG	62	-1463.6437	3300.7411	0.0000	1.0000
HKY	59	-1473.5911	3302.5655	0.0000	1.0000
GTR	63	-1469.9772	3319.4316	0.0000	1.0000
F81G	59	-1482.9470	3321.2775	0.0000	1.0000
K2P	56	-1493.2706	3323.8543	0.0000	1.0000
F81I	59	-1485.1387	3325.6609	0.0000	1.0000
F81IG	60	-1482.5083	3326.4234	0.0000	1.0000
JC69G	56	-1497.7516	3332.8162	0.0000	1.0000
JC69I	56	-1500.1111	3337.5352	0.0000	1.0000
JC69IG	57	-1497.2521	3337.8406	0.0000	1.0000
SYM	60	-1490.9986	3343.4041	0.0000	1.0000
F81	58	-1507.7857	3364.9313	0.0000	1.0000
JC69	55	-1523.4514	3378.1924	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/7829_hg18.phylip = 1 - 413;
 Partition Dummy = 1:261-aggressive-phylip/7829_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/7829_hg18.phylip = 1 - 413;
 Partition Dummy = 1:261-aggressive-phylip/7829_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/7829_hg18.phylip = 1 - 413;
 Partition Dummy = 1:261-aggressive-phylip/7829_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

