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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRG
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
GTRG	64	-2475.4021	5095.3779	0.5237	0.5237
GTRIG	65	-2474.2237	5095.5731	0.4750	0.9987
HKYG	60	-2486.6442	5107.7549	0.0011	0.9998
HKYI	60	-2489.0691	5112.6046	0.0001	0.9999
HKYIG	61	-2487.8986	5112.7755	0.0001	1.0000
GTRI	64	-2496.2870	5137.1478	0.0000	1.0000
SYMIG	62	-2500.4292	5140.3584	0.0000	1.0000
SYMG	61	-2502.5722	5142.1226	0.0000	1.0000
K2PIG	58	-2511.2488	5151.9700	0.0000	1.0000
K2PG	57	-2513.6193	5154.2288	0.0000	1.0000
SYMI	61	-2521.3487	5179.6756	0.0000	1.0000
K2PI	57	-2530.8544	5188.6991	0.0000	1.0000
GTR	63	-2568.5718	5279.1755	0.0000	1.0000
HKY	59	-2579.7994	5291.5633	0.0000	1.0000
F81G	59	-2584.7807	5301.5260	0.0000	1.0000
F81IG	60	-2588.2507	5310.9677	0.0000	1.0000
SYM	60	-2599.7897	5334.0459	0.0000	1.0000
F81I	59	-2604.4746	5340.9136	0.0000	1.0000
JC69IG	57	-2607.6122	5342.2145	0.0000	1.0000
JC69G	56	-2609.7514	5344.0204	0.0000	1.0000
K2P	56	-2610.1596	5344.8368	0.0000	1.0000
JC69I	56	-2625.2709	5375.0595	0.0000	1.0000
F81	58	-2681.9123	5493.2971	0.0000	1.0000
JC69	55	-2706.9362	5535.9272	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-2474.2237	5078.4474	0.5443	0.5443
GTRG	64	-2475.4021	5078.8042	0.4553	0.9996
HKYG	60	-2486.6442	5093.2885	0.0003	0.9999
HKYIG	61	-2487.8986	5097.7972	0.0000	1.0000
HKYI	60	-2489.0691	5098.1382	0.0000	1.0000
GTRI	64	-2496.2870	5120.5741	0.0000	1.0000
SYMIG	62	-2500.4292	5124.8584	0.0000	1.0000
SYMG	61	-2502.5722	5127.1443	0.0000	1.0000
K2PIG	58	-2511.2488	5138.4975	0.0000	1.0000
K2PG	57	-2513.6193	5141.2386	0.0000	1.0000
SYMI	61	-2521.3487	5164.6974	0.0000	1.0000
K2PI	57	-2530.8544	5175.7089	0.0000	1.0000
GTR	63	-2568.5718	5263.1437	0.0000	1.0000
HKY	59	-2579.7994	5277.5988	0.0000	1.0000
F81G	59	-2584.7807	5287.5615	0.0000	1.0000
F81IG	60	-2588.2507	5296.5013	0.0000	1.0000
SYM	60	-2599.7897	5319.5795	0.0000	1.0000
F81I	59	-2604.4746	5326.9491	0.0000	1.0000
JC69IG	57	-2607.6122	5329.2243	0.0000	1.0000
JC69G	56	-2609.7514	5331.5028	0.0000	1.0000
K2P	56	-2610.1596	5332.3192	0.0000	1.0000
JC69I	56	-2625.2709	5362.5419	0.0000	1.0000
F81	58	-2681.9123	5479.8247	0.0000	1.0000
JC69	55	-2706.9362	5523.8724	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2486.6442	5353.7101	0.7299	0.7299
GTRG	64	-2475.4021	5356.5872	0.1732	0.9030
HKYI	60	-2489.0691	5358.5598	0.0646	0.9676
GTRIG	65	-2474.2237	5360.5707	0.0236	0.9913
HKYIG	61	-2487.8986	5362.5592	0.0087	1.0000
K2PG	57	-2513.6193	5388.6391	0.0000	1.0000
K2PIG	58	-2511.2488	5390.2384	0.0000	1.0000
SYMG	61	-2502.5722	5391.9063	0.0000	1.0000
SYMIG	62	-2500.4292	5393.9607	0.0000	1.0000
GTRI	64	-2496.2870	5398.3571	0.0000	1.0000
K2PI	57	-2530.8544	5423.1094	0.0000	1.0000
SYMI	61	-2521.3487	5429.4593	0.0000	1.0000
HKY	59	-2579.7994	5533.6800	0.0000	1.0000
GTR	63	-2568.5718	5536.5863	0.0000	1.0000
F81G	59	-2584.7807	5543.6427	0.0000	1.0000
F81IG	60	-2588.2507	5556.9229	0.0000	1.0000
JC69G	56	-2609.7514	5574.5629	0.0000	1.0000
K2P	56	-2610.1596	5575.3793	0.0000	1.0000
JC69IG	57	-2607.6122	5576.6248	0.0000	1.0000
SYM	60	-2599.7897	5580.0010	0.0000	1.0000
F81I	59	-2604.4746	5583.0303	0.0000	1.0000
JC69I	56	-2625.2709	5605.6020	0.0000	1.0000
F81	58	-2681.9123	5731.5655	0.0000	1.0000
JC69	55	-2706.9362	5762.5922	0.0000	1.0000



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

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

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


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/10845_hg18.phylip = 1 - 567;
 Partition Dummy = 1:261-aggressive-phylip/10845_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 (HKYG)]

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

