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

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

Minimum AIC  model: GTRI
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	-1108.2038	2354.6165	0.3253	0.3253
HKYG	60	-1108.2283	2354.6655	0.3174	0.6427
GTRI	64	-1103.7761	2356.4568	0.1296	0.7724
GTRG	64	-1103.8967	2356.6980	0.1149	0.8873
HKYIG	61	-1108.2384	2357.3397	0.0834	0.9706
GTRIG	65	-1103.9093	2359.4306	0.0293	0.9999
HKY	59	-1118.6846	2372.9375	0.0000	1.0000
GTR	63	-1113.6164	2373.4432	0.0000	1.0000
SYMI	61	-1118.0342	2376.9312	0.0000	1.0000
F81I	59	-1120.7041	2376.9764	0.0000	1.0000
F81G	59	-1120.7154	2376.9991	0.0000	1.0000
SYMG	61	-1118.1585	2377.1798	0.0000	1.0000
F81IG	60	-1120.7657	2379.7403	0.0000	1.0000
SYMIG	62	-1118.1701	2379.8703	0.0000	1.0000
K2PI	57	-1125.2258	2380.7776	0.0000	1.0000
K2PG	57	-1125.2387	2380.8034	0.0000	1.0000
K2PIG	58	-1125.2832	2383.5069	0.0000	1.0000
SYM	60	-1127.6761	2393.5611	0.0000	1.0000
F81	58	-1131.2502	2395.4411	0.0000	1.0000
K2P	56	-1135.6854	2399.0950	0.0000	1.0000
JC69I	56	-1137.0299	2401.7839	0.0000	1.0000
JC69G	56	-1137.0576	2401.8394	0.0000	1.0000
JC69IG	57	-1137.0674	2404.4607	0.0000	1.0000
JC69	55	-1147.3774	2419.8900	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	64	-1103.7761	2335.5523	0.2682	0.2682
GTRG	64	-1103.8967	2335.7934	0.2377	0.5059
HKYI	60	-1108.2038	2336.4076	0.1749	0.6808
HKYG	60	-1108.2283	2336.4565	0.1706	0.8514
GTRIG	65	-1103.9093	2337.8185	0.0864	0.9378
HKYIG	61	-1108.2384	2338.4768	0.0621	0.9999
GTR	63	-1113.6164	2353.2327	0.0000	1.0000
HKY	59	-1118.6846	2355.3692	0.0000	1.0000
SYMI	61	-1118.0342	2358.0683	0.0000	1.0000
SYMG	61	-1118.1585	2358.3170	0.0000	1.0000
F81I	59	-1120.7041	2359.4082	0.0000	1.0000
F81G	59	-1120.7154	2359.4309	0.0000	1.0000
SYMIG	62	-1118.1701	2360.3403	0.0000	1.0000
F81IG	60	-1120.7657	2361.5313	0.0000	1.0000
K2PI	57	-1125.2258	2364.4517	0.0000	1.0000
K2PG	57	-1125.2387	2364.4775	0.0000	1.0000
K2PIG	58	-1125.2832	2366.5663	0.0000	1.0000
SYM	60	-1127.6761	2375.3522	0.0000	1.0000
F81	58	-1131.2502	2378.5005	0.0000	1.0000
K2P	56	-1135.6854	2383.3709	0.0000	1.0000
JC69I	56	-1137.0299	2386.0598	0.0000	1.0000
JC69G	56	-1137.0576	2386.1153	0.0000	1.0000
JC69IG	57	-1137.0674	2388.1348	0.0000	1.0000
JC69	55	-1147.3774	2404.7549	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	60	-1108.2038	2584.6712	0.4943	0.4943
HKYG	60	-1108.2283	2584.7202	0.4823	0.9767
HKYIG	61	-1108.2384	2590.8782	0.0222	0.9988
HKY	59	-1118.6846	2599.4951	0.0003	0.9991
K2PI	57	-1125.2258	2600.3021	0.0002	0.9993
K2PG	57	-1125.2387	2600.3279	0.0002	0.9995
GTRI	64	-1103.7761	2600.3668	0.0002	0.9997
GTRG	64	-1103.8967	2600.6080	0.0002	0.9999
F81I	59	-1120.7041	2603.5341	0.0000	0.9999
F81G	59	-1120.7154	2603.5568	0.0000	1.0000
K2PIG	58	-1125.2832	2606.5545	0.0000	1.0000
GTRIG	65	-1103.9093	2606.7708	0.0000	1.0000
F81IG	60	-1120.7657	2609.7950	0.0000	1.0000
SYMI	61	-1118.0342	2610.4697	0.0000	1.0000
SYMG	61	-1118.1585	2610.7183	0.0000	1.0000
GTR	63	-1113.6164	2613.9095	0.0000	1.0000
K2P	56	-1135.6854	2615.0836	0.0000	1.0000
SYMIG	62	-1118.1701	2616.8793	0.0000	1.0000
JC69I	56	-1137.0299	2617.7725	0.0000	1.0000
JC69G	56	-1137.0576	2617.8280	0.0000	1.0000
F81	58	-1131.2502	2618.4886	0.0000	1.0000
SYM	60	-1127.6761	2623.6158	0.0000	1.0000
JC69IG	57	-1137.0674	2623.9853	0.0000	1.0000
JC69	55	-1147.3774	2632.3299	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/10298_hg18.phylip = 1 - 463;
 Partition Dummy = 1:261-aggressive-phylip/10298_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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

