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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	60	-1131.2713	2405.4894	0.4599	0.4599
HKYG	60	-1131.3651	2405.6769	0.4187	0.8786
HKYIG	61	-1131.3646	2408.5153	0.1013	0.9799
GTRI	64	-1129.3645	2413.1416	0.0100	0.9900
GTRG	64	-1129.5734	2413.5596	0.0081	0.9981
GTRIG	65	-1129.5728	2416.4705	0.0019	1.0000
F81IG	60	-1143.6223	2430.1913	0.0000	1.0000
F81I	59	-1146.1874	2432.4999	0.0000	1.0000
F81G	59	-1146.3647	2432.8545	0.0000	1.0000
HKY	59	-1149.7089	2439.5427	0.0000	1.0000
GTR	63	-1147.3533	2446.2256	0.0000	1.0000
K2PI	57	-1162.4216	2459.3773	0.0000	1.0000
K2PG	57	-1162.5786	2459.6913	0.0000	1.0000
SYMI	61	-1157.6282	2461.0427	0.0000	1.0000
SYMG	61	-1157.7986	2461.3834	0.0000	1.0000
K2PIG	58	-1162.5782	2462.4774	0.0000	1.0000
SYMIG	62	-1157.7980	2464.2396	0.0000	1.0000
F81	58	-1164.9822	2467.2852	0.0000	1.0000
JC69G	56	-1171.6267	2475.0181	0.0000	1.0000
JC69IG	57	-1171.6211	2477.7764	0.0000	1.0000
JC69I	56	-1174.4331	2480.6309	0.0000	1.0000
K2P	56	-1182.1331	2496.0309	0.0000	1.0000
SYM	60	-1177.7490	2498.4447	0.0000	1.0000
JC69	55	-1194.1282	2517.2687	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	60	-1131.2713	2382.5427	0.3991	0.3991
HKYG	60	-1131.3651	2382.7302	0.3634	0.7624
HKYIG	61	-1131.3646	2384.7292	0.1337	0.8962
GTRI	64	-1129.3645	2386.7289	0.0492	0.9454
GTRG	64	-1129.5734	2387.1469	0.0399	0.9853
GTRIG	65	-1129.5728	2389.1457	0.0147	1.0000
F81IG	60	-1143.6223	2407.2446	0.0000	1.0000
F81I	59	-1146.1874	2410.3749	0.0000	1.0000
F81G	59	-1146.3647	2410.7295	0.0000	1.0000
HKY	59	-1149.7089	2417.4177	0.0000	1.0000
GTR	63	-1147.3533	2420.7067	0.0000	1.0000
SYMI	61	-1157.6282	2437.2565	0.0000	1.0000
SYMG	61	-1157.7986	2437.5973	0.0000	1.0000
K2PI	57	-1162.4216	2438.8431	0.0000	1.0000
K2PG	57	-1162.5786	2439.1572	0.0000	1.0000
SYMIG	62	-1157.7980	2439.5960	0.0000	1.0000
K2PIG	58	-1162.5782	2441.1565	0.0000	1.0000
F81	58	-1164.9822	2445.9643	0.0000	1.0000
JC69G	56	-1171.6267	2455.2534	0.0000	1.0000
JC69IG	57	-1171.6211	2457.2422	0.0000	1.0000
JC69I	56	-1174.4331	2460.8662	0.0000	1.0000
SYM	60	-1177.7490	2475.4980	0.0000	1.0000
K2P	56	-1182.1331	2476.2662	0.0000	1.0000
JC69	55	-1194.1282	2498.2563	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	60	-1131.2713	2618.9529	0.5109	0.5109
HKYG	60	-1131.3651	2619.1405	0.4652	0.9761
HKYIG	61	-1131.3646	2625.0796	0.0239	0.9999
GTRI	64	-1129.3645	2638.8999	0.0000	1.0000
GTRG	64	-1129.5734	2639.3178	0.0000	1.0000
F81I	59	-1146.1874	2642.8450	0.0000	1.0000
F81G	59	-1146.3647	2643.1996	0.0000	1.0000
F81IG	60	-1143.6223	2643.6549	0.0000	1.0000
GTRIG	65	-1129.5728	2645.2568	0.0000	1.0000
HKY	59	-1149.7089	2649.8878	0.0000	1.0000
K2PI	57	-1162.4216	2663.4329	0.0000	1.0000
K2PG	57	-1162.5786	2663.7469	0.0000	1.0000
GTR	63	-1147.3533	2668.9374	0.0000	1.0000
K2PIG	58	-1162.5782	2669.6864	0.0000	1.0000
F81	58	-1164.9822	2674.4943	0.0000	1.0000
JC69G	56	-1171.6267	2675.9030	0.0000	1.0000
SYMI	61	-1157.6282	2677.6069	0.0000	1.0000
SYMG	61	-1157.7986	2677.9477	0.0000	1.0000
JC69I	56	-1174.4331	2681.5158	0.0000	1.0000
JC69IG	57	-1171.6211	2681.8320	0.0000	1.0000
SYMIG	62	-1157.7980	2683.8867	0.0000	1.0000
K2P	56	-1182.1331	2696.9158	0.0000	1.0000
SYM	60	-1177.7490	2711.9083	0.0000	1.0000
JC69	55	-1194.1282	2714.9657	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYI)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7466_hg18.phylip = 1 - 380;
 Partition Dummy = 1:261-aggressive-phylip/7466_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (HKYI)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7466_hg18.phylip = 1 - 380;
 Partition Dummy = 1:261-aggressive-phylip/7466_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7466_hg18.phylip = 1 - 380;
 Partition Dummy = 1:261-aggressive-phylip/7466_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

