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

Input data from file "261-aggressive-phylip/10088_hg18.phylip" (29 taxa, 604 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	-1361.5902	2866.6164	0.7253	0.7253
GTRIG	65	-1361.6045	2869.1570	0.2036	0.9290
GTR	63	-1365.5907	2872.1148	0.0464	0.9754
GTRI	64	-1365.6035	2874.6429	0.0131	0.9885
HKYG	60	-1371.0441	2875.5688	0.0083	0.9967
HKYIG	61	-1371.0585	2878.0728	0.0024	0.9991
HKY	59	-1374.7345	2880.4838	0.0007	0.9998
HKYI	60	-1374.7476	2882.9758	0.0002	1.0000
F81G	59	-1394.4540	2919.9226	0.0000	1.0000
K2PG	57	-1397.5663	2921.2425	0.0000	1.0000
F81IG	60	-1394.4684	2922.4174	0.0000	1.0000
K2PIG	58	-1397.5809	2923.7196	0.0000	1.0000
SYMG	61	-1394.3439	2924.6435	0.0000	1.0000
F81	58	-1398.1945	2924.9469	0.0000	1.0000
SYMIG	62	-1394.3584	2927.1567	0.0000	1.0000
F81I	59	-1398.2074	2927.4295	0.0000	1.0000
K2P	56	-1401.9067	2927.4844	0.0000	1.0000
SYM	60	-1398.1951	2929.8708	0.0000	1.0000
K2PI	57	-1401.9198	2929.9495	0.0000	1.0000
SYMI	61	-1398.2080	2932.3718	0.0000	1.0000
JC69G	56	-1417.2969	2958.2648	0.0000	1.0000
JC69IG	57	-1417.3112	2960.7324	0.0000	1.0000
JC69	55	-1421.4784	2964.1977	0.0000	1.0000
JC69I	56	-1421.4912	2966.6532	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1361.5902	2851.1804	0.6960	0.6960
GTRIG	65	-1361.6045	2853.2091	0.2524	0.9484
GTR	63	-1365.5907	2857.1814	0.0346	0.9831
GTRI	64	-1365.6035	2859.2069	0.0126	0.9957
HKYG	60	-1371.0441	2862.0881	0.0030	0.9986
HKYIG	61	-1371.0585	2864.1170	0.0011	0.9997
HKY	59	-1374.7345	2867.4691	0.0002	0.9999
HKYI	60	-1374.7476	2869.4951	0.0001	1.0000
F81G	59	-1394.4540	2906.9079	0.0000	1.0000
F81IG	60	-1394.4684	2908.9367	0.0000	1.0000
K2PG	57	-1397.5663	2909.1326	0.0000	1.0000
SYMG	61	-1394.3439	2910.6878	0.0000	1.0000
K2PIG	58	-1397.5809	2911.1618	0.0000	1.0000
F81	58	-1398.1945	2912.3891	0.0000	1.0000
SYMIG	62	-1394.3584	2912.7168	0.0000	1.0000
F81I	59	-1398.2074	2914.4148	0.0000	1.0000
K2P	56	-1401.9067	2915.8135	0.0000	1.0000
SYM	60	-1398.1951	2916.3901	0.0000	1.0000
K2PI	57	-1401.9198	2917.8396	0.0000	1.0000
SYMI	61	-1398.2080	2918.4161	0.0000	1.0000
JC69G	56	-1417.2969	2946.5939	0.0000	1.0000
JC69IG	57	-1417.3112	2948.6225	0.0000	1.0000
JC69	55	-1421.4784	2952.9568	0.0000	1.0000
JC69I	56	-1421.4912	2954.9823	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1371.0441	3126.3026	0.5777	0.5777
HKY	59	-1374.7345	3127.2800	0.3544	0.9321
HKYIG	61	-1371.0585	3132.7351	0.0232	0.9553
GTRG	64	-1361.5902	3133.0091	0.0202	0.9755
HKYI	60	-1374.7476	3133.7096	0.0142	0.9897
GTR	63	-1365.5907	3134.6066	0.0091	0.9988
GTRIG	65	-1361.6045	3139.4414	0.0008	0.9996
GTRI	64	-1365.6035	3141.0357	0.0004	1.0000
K2PG	57	-1397.5663	3160.1363	0.0000	1.0000
K2P	56	-1401.9067	3162.4136	0.0000	1.0000
K2PIG	58	-1397.5809	3166.5691	0.0000	1.0000
F81G	59	-1394.4540	3166.7188	0.0000	1.0000
F81	58	-1398.1945	3167.7964	0.0000	1.0000
K2PI	57	-1401.9198	3168.8433	0.0000	1.0000
F81IG	60	-1394.4684	3173.1512	0.0000	1.0000
F81I	59	-1398.2074	3174.2257	0.0000	1.0000
SYMG	61	-1394.3439	3179.3058	0.0000	1.0000
SYM	60	-1398.1951	3180.6046	0.0000	1.0000
SYMIG	62	-1394.3584	3185.7384	0.0000	1.0000
SYMI	61	-1398.2080	3187.0341	0.0000	1.0000
JC69G	56	-1417.2969	3193.1941	0.0000	1.0000
JC69	55	-1421.4784	3195.1534	0.0000	1.0000
JC69IG	57	-1417.3112	3199.6262	0.0000	1.0000
JC69I	56	-1421.4912	3201.5825	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/10088_hg18.phylip = 1 - 604;
 Partition Dummy = 1:261-aggressive-phylip/10088_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/10088_hg18.phylip = 1 - 604;
 Partition Dummy = 1:261-aggressive-phylip/10088_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/10088_hg18.phylip = 1 - 604;
 Partition Dummy = 1:261-aggressive-phylip/10088_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

