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

Input data from file "261-aggressive-phylip/4023_hg18.phylip" (29 taxa, 335 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	-1141.3128	2429.3408	0.2230	0.2230
K2PG	57	-1145.7994	2429.4688	0.2092	0.4321
HKYI	60	-1141.9065	2430.5283	0.1231	0.5553
K2PI	57	-1146.3891	2430.6483	0.1160	0.6712
SYMG	61	-1140.5758	2430.8585	0.1044	0.7756
HKYIG	61	-1141.2007	2432.1083	0.0559	0.8315
K2PIG	58	-1145.7183	2432.2338	0.0525	0.8840
SYMI	61	-1141.3487	2432.4043	0.0482	0.9322
SYMIG	62	-1140.4384	2433.5973	0.0265	0.9587
GTRG	64	-1137.6429	2434.1007	0.0206	0.9794
GTRI	64	-1137.9778	2434.7704	0.0148	0.9941
GTRIG	65	-1137.3623	2436.6205	0.0059	1.0000
JC69G	56	-1167.6613	2470.2865	0.0000	1.0000
JC69I	56	-1168.2739	2471.5118	0.0000	1.0000
F81G	59	-1163.9444	2471.6344	0.0000	1.0000
F81I	59	-1164.5426	2472.8306	0.0000	1.0000
JC69IG	57	-1167.5908	2473.0516	0.0000	1.0000
F81IG	60	-1163.8714	2474.4582	0.0000	1.0000
K2P	56	-1174.0634	2483.0909	0.0000	1.0000
HKY	59	-1169.8111	2483.3678	0.0000	1.0000
SYM	60	-1170.6228	2487.9610	0.0000	1.0000
GTR	63	-1167.4758	2490.7081	0.0000	1.0000
JC69	55	-1195.3799	2522.8387	0.0000	1.0000
F81	58	-1191.6953	2524.1878	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1141.3128	2402.6255	0.1839	0.1839
SYMG	61	-1140.5758	2403.1516	0.1414	0.3253
GTRG	64	-1137.6429	2403.2859	0.1322	0.4575
HKYI	60	-1141.9065	2403.8130	0.1016	0.5591
GTRI	64	-1137.9778	2403.9556	0.0946	0.6537
HKYIG	61	-1141.2007	2404.4013	0.0757	0.7294
SYMI	61	-1141.3487	2404.6974	0.0653	0.7947
GTRIG	65	-1137.3623	2404.7246	0.0644	0.8591
SYMIG	62	-1140.4384	2404.8767	0.0597	0.9187
K2PG	57	-1145.7994	2405.5987	0.0416	0.9603
K2PI	57	-1146.3891	2406.7782	0.0231	0.9834
K2PIG	58	-1145.7183	2407.4367	0.0166	1.0000
F81G	59	-1163.9444	2445.8889	0.0000	1.0000
F81I	59	-1164.5426	2447.0851	0.0000	1.0000
JC69G	56	-1167.6613	2447.3225	0.0000	1.0000
F81IG	60	-1163.8714	2447.7429	0.0000	1.0000
JC69I	56	-1168.2739	2448.5478	0.0000	1.0000
JC69IG	57	-1167.5908	2449.1816	0.0000	1.0000
HKY	59	-1169.8111	2457.6223	0.0000	1.0000
K2P	56	-1174.0634	2460.1269	0.0000	1.0000
GTR	63	-1167.4758	2460.9517	0.0000	1.0000
SYM	60	-1170.6228	2461.2457	0.0000	1.0000
F81	58	-1191.6953	2499.3907	0.0000	1.0000
JC69	55	-1195.3799	2500.7598	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1145.7994	2623.0042	0.6099	0.6099
K2PI	57	-1146.3891	2624.1837	0.3382	0.9481
K2PIG	58	-1145.7183	2628.6562	0.0361	0.9842
HKYG	60	-1141.3128	2631.4734	0.0088	0.9930
HKYI	60	-1141.9065	2632.6608	0.0049	0.9979
SYMG	61	-1140.5758	2635.8135	0.0010	0.9989
HKYIG	61	-1141.2007	2637.0633	0.0005	0.9995
SYMI	61	-1141.3487	2637.3593	0.0005	0.9999
SYMIG	62	-1140.4384	2641.3528	0.0001	1.0000
GTRG	64	-1137.6429	2647.3902	0.0000	1.0000
GTRI	64	-1137.9778	2648.0600	0.0000	1.0000
GTRIG	65	-1137.3623	2652.6431	0.0000	1.0000
JC69G	56	-1167.6613	2660.9138	0.0000	1.0000
JC69I	56	-1168.2739	2662.1391	0.0000	1.0000
JC69IG	57	-1167.5908	2666.5870	0.0000	1.0000
F81G	59	-1163.9444	2670.9226	0.0000	1.0000
F81I	59	-1164.5426	2672.1188	0.0000	1.0000
K2P	56	-1174.0634	2673.7182	0.0000	1.0000
F81IG	60	-1163.8714	2676.5907	0.0000	1.0000
HKY	59	-1169.8111	2682.6560	0.0000	1.0000
SYM	60	-1170.6228	2690.0935	0.0000	1.0000
GTR	63	-1167.4758	2701.2419	0.0000	1.0000
JC69	55	-1195.3799	2710.5370	0.0000	1.0000
F81	58	-1191.6953	2720.6102	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/4023_hg18.phylip = 1 - 335;
 Partition Dummy = 1:261-aggressive-phylip/4023_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/4023_hg18.phylip = 1 - 335;
 Partition Dummy = 1:261-aggressive-phylip/4023_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/4023_hg18.phylip = 1 - 335;
 Partition Dummy = 1:261-aggressive-phylip/4023_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

