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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-1375.4439	2884.5444	0.4247	0.4247
HKYI	60	-1375.6652	2884.9870	0.3404	0.7651
HKYIG	61	-1375.4147	2886.9678	0.1264	0.8915
GTRG	64	-1372.7336	2889.1062	0.0434	0.9349
GTRI	64	-1372.9702	2889.5795	0.0343	0.9692
HKY	59	-1379.9020	2890.9884	0.0169	0.9861
GTRIG	65	-1372.6852	2891.5286	0.0129	0.9990
GTR	63	-1377.8029	2896.7352	0.0010	1.0000
K2PG	57	-1392.1249	2910.5169	0.0000	1.0000
K2PI	57	-1392.2991	2910.8654	0.0000	1.0000
K2PIG	58	-1392.1107	2912.9425	0.0000	1.0000
K2P	56	-1396.1907	2916.2037	0.0000	1.0000
SYMG	61	-1390.7685	2917.6753	0.0000	1.0000
SYMI	61	-1390.9290	2917.9964	0.0000	1.0000
SYMIG	62	-1390.7529	2920.1351	0.0000	1.0000
SYM	60	-1395.0293	2923.7152	0.0000	1.0000
F81G	59	-1398.0148	2927.2139	0.0000	1.0000
F81I	59	-1398.2396	2927.6636	0.0000	1.0000
F81IG	60	-1397.9975	2929.6517	0.0000	1.0000
F81	58	-1401.8897	2932.5006	0.0000	1.0000
JC69G	56	-1412.7914	2949.4050	0.0000	1.0000
JC69I	56	-1412.9898	2949.8018	0.0000	1.0000
JC69IG	57	-1412.7793	2951.8258	0.0000	1.0000
JC69	55	-1416.5413	2954.4689	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1375.4439	2870.8877	0.3552	0.3552
HKYI	60	-1375.6652	2871.3303	0.2847	0.6399
HKYIG	61	-1375.4147	2872.8295	0.1345	0.7744
GTRG	64	-1372.7336	2873.4671	0.0978	0.8722
GTRI	64	-1372.9702	2873.9404	0.0772	0.9494
GTRIG	65	-1372.6852	2875.3704	0.0378	0.9871
HKY	59	-1379.9020	2877.8040	0.0112	0.9983
GTR	63	-1377.8029	2881.6057	0.0017	1.0000
K2PG	57	-1392.1249	2898.2498	0.0000	1.0000
K2PI	57	-1392.2991	2898.5982	0.0000	1.0000
K2PIG	58	-1392.1107	2900.2213	0.0000	1.0000
SYMG	61	-1390.7685	2903.5370	0.0000	1.0000
SYMI	61	-1390.9290	2903.8581	0.0000	1.0000
K2P	56	-1396.1907	2904.3815	0.0000	1.0000
SYMIG	62	-1390.7529	2905.5058	0.0000	1.0000
SYM	60	-1395.0293	2910.0585	0.0000	1.0000
F81G	59	-1398.0148	2914.0295	0.0000	1.0000
F81I	59	-1398.2396	2914.4792	0.0000	1.0000
F81IG	60	-1397.9975	2915.9950	0.0000	1.0000
F81	58	-1401.8897	2919.7794	0.0000	1.0000
JC69G	56	-1412.7914	2937.5828	0.0000	1.0000
JC69I	56	-1412.9898	2937.9796	0.0000	1.0000
JC69IG	57	-1412.7793	2939.5586	0.0000	1.0000
JC69	55	-1416.5413	2943.0826	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1375.4439	3134.4027	0.4698	0.4698
HKYI	60	-1375.6652	3134.8453	0.3765	0.8463
HKY	59	-1379.9020	3136.9271	0.1330	0.9793
HKYIG	61	-1375.4147	3140.7364	0.0198	0.9991
K2PG	57	-1392.1249	3148.5890	0.0004	0.9995
K2PI	57	-1392.2991	3148.9375	0.0003	0.9998
K2P	56	-1396.1907	3150.3289	0.0002	0.9999
GTRG	64	-1372.7336	3154.5498	0.0000	1.0000
K2PIG	58	-1392.1107	3154.9525	0.0000	1.0000
GTRI	64	-1372.9702	3155.0231	0.0000	1.0000
GTR	63	-1377.8029	3158.2965	0.0000	1.0000
GTRIG	65	-1372.6852	3160.8450	0.0000	1.0000
SYMG	61	-1390.7685	3171.4439	0.0000	1.0000
SYMI	61	-1390.9290	3171.7650	0.0000	1.0000
F81G	59	-1398.0148	3173.1526	0.0000	1.0000
SYM	60	-1395.0293	3173.5735	0.0000	1.0000
F81I	59	-1398.2396	3173.6023	0.0000	1.0000
F81	58	-1401.8897	3174.5106	0.0000	1.0000
SYMIG	62	-1390.7529	3177.8047	0.0000	1.0000
F81IG	60	-1397.9975	3179.5100	0.0000	1.0000
JC69G	56	-1412.7914	3183.5301	0.0000	1.0000
JC69I	56	-1412.9898	3183.9270	0.0000	1.0000
JC69	55	-1416.5413	3184.6380	0.0000	1.0000
JC69IG	57	-1412.7793	3189.8979	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/4293_hg18.phylip = 1 - 597;
 Partition Dummy = 1:261-aggressive-phylip/4293_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/4293_hg18.phylip = 1 - 597;
 Partition Dummy = 1:261-aggressive-phylip/4293_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 (HKYG)]

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

