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

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

Minimum AIC  model: GTRG
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	-1246.8857	2642.8882	0.4586	0.4586
GTRI	64	-1247.0035	2643.1237	0.4076	0.8662
GTRIG	65	-1246.8913	2645.6146	0.1173	0.9835
HKYG	60	-1256.4041	2651.2002	0.0072	0.9907
HKYI	60	-1256.4666	2651.3252	0.0068	0.9975
HKYIG	61	-1256.4035	2653.8598	0.0019	0.9994
F81G	59	-1260.9998	2657.7440	0.0003	0.9996
F81I	59	-1261.0603	2657.8650	0.0003	0.9999
F81IG	60	-1261.0023	2660.3965	0.0001	1.0000
SYMG	61	-1261.2533	2663.5596	0.0000	1.0000
GTR	63	-1258.6157	2663.6466	0.0000	1.0000
SYMI	61	-1261.4257	2663.9043	0.0000	1.0000
SYMIG	62	-1261.2291	2666.1855	0.0000	1.0000
HKY	59	-1268.8648	2673.4739	0.0000	1.0000
K2PG	57	-1272.6862	2675.8612	0.0000	1.0000
K2PI	57	-1272.8039	2676.0967	0.0000	1.0000
K2PIG	58	-1272.6657	2678.4413	0.0000	1.0000
F81	58	-1273.4559	2680.0218	0.0000	1.0000
JC69G	56	-1276.5676	2681.0158	0.0000	1.0000
JC69I	56	-1276.6726	2681.2257	0.0000	1.0000
JC69IG	57	-1276.5532	2683.5952	0.0000	1.0000
SYM	60	-1273.6687	2685.7294	0.0000	1.0000
K2P	56	-1285.2361	2698.3528	0.0000	1.0000
JC69	55	-1289.1024	2703.4901	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1246.8857	2621.7715	0.4416	0.4416
GTRI	64	-1247.0035	2622.0069	0.3926	0.8342
GTRIG	65	-1246.8913	2623.7825	0.1616	0.9958
HKYG	60	-1256.4041	2632.8082	0.0018	0.9975
HKYI	60	-1256.4666	2632.9332	0.0017	0.9992
HKYIG	61	-1256.4035	2634.8069	0.0007	0.9999
F81G	59	-1260.9998	2639.9997	0.0000	0.9999
F81I	59	-1261.0603	2640.1206	0.0000	1.0000
F81IG	60	-1261.0023	2642.0045	0.0000	1.0000
GTR	63	-1258.6157	2643.2314	0.0000	1.0000
SYMG	61	-1261.2533	2644.5067	0.0000	1.0000
SYMI	61	-1261.4257	2644.8514	0.0000	1.0000
SYMIG	62	-1261.2291	2646.4582	0.0000	1.0000
HKY	59	-1268.8648	2655.7295	0.0000	1.0000
K2PG	57	-1272.6862	2659.3724	0.0000	1.0000
K2PI	57	-1272.8039	2659.6079	0.0000	1.0000
K2PIG	58	-1272.6657	2661.3313	0.0000	1.0000
F81	58	-1273.4559	2662.9118	0.0000	1.0000
JC69G	56	-1276.5676	2665.1352	0.0000	1.0000
JC69I	56	-1276.6726	2665.3451	0.0000	1.0000
JC69IG	57	-1276.5532	2667.1064	0.0000	1.0000
SYM	60	-1273.6687	2667.3375	0.0000	1.0000
K2P	56	-1285.2361	2682.4722	0.0000	1.0000
JC69	55	-1289.1024	2688.2047	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1256.4041	2880.5512	0.3931	0.3931
HKYI	60	-1256.4666	2880.6762	0.3692	0.7623
F81G	59	-1260.9998	2883.6136	0.0850	0.8473
F81I	59	-1261.0603	2883.7346	0.0800	0.9273
GTRG	64	-1246.8857	2886.0307	0.0254	0.9527
GTRI	64	-1247.0035	2886.2662	0.0226	0.9753
HKYIG	61	-1256.4035	2886.6790	0.0184	0.9936
F81IG	60	-1261.0023	2889.7475	0.0040	0.9976
GTRIG	65	-1246.8913	2892.1708	0.0012	0.9988
K2PG	57	-1272.6862	2894.7283	0.0003	0.9991
K2PI	57	-1272.8039	2894.9638	0.0003	0.9994
JC69G	56	-1276.5676	2896.3620	0.0001	0.9995
SYMG	61	-1261.2533	2896.3787	0.0001	0.9997
JC69I	56	-1276.6726	2896.5719	0.0001	0.9998
SYMI	61	-1261.4257	2896.7235	0.0001	0.9999
HKY	59	-1268.8648	2899.3435	0.0000	1.0000
K2PIG	58	-1272.6657	2900.8162	0.0000	1.0000
F81	58	-1273.4559	2902.3967	0.0000	1.0000
SYMIG	62	-1261.2291	2902.4593	0.0000	1.0000
JC69IG	57	-1276.5532	2902.4622	0.0000	1.0000
GTR	63	-1258.6157	2903.3616	0.0000	1.0000
K2P	56	-1285.2361	2913.6990	0.0000	1.0000
SYM	60	-1273.6687	2915.0805	0.0000	1.0000
JC69	55	-1289.1024	2915.3025	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/1831_hg18.phylip = 1 - 459;
 Partition Dummy = 1:261-aggressive-phylip/1831_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 AICc model (GTRG)]

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

