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

Input data from file "261-aggressive-phylip/8760_hg18.phylip" (29 taxa, 582 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	-2800.2364	5744.5656	0.7312	0.7312
GTRIG	65	-2799.9747	5746.5774	0.2674	0.9987
HKYG	60	-2811.9109	5757.8717	0.0009	0.9996
HKYIG	61	-2811.5427	5759.6315	0.0004	1.0000
GTRI	64	-2811.7246	5767.5421	0.0000	1.0000
HKYI	60	-2821.8756	5777.8010	0.0000	1.0000
K2PG	57	-2840.7943	5808.2070	0.0000	1.0000
SYMG	61	-2836.2051	5808.9563	0.0000	1.0000
K2PIG	58	-2840.2826	5809.6513	0.0000	1.0000
SYMIG	62	-2837.4612	5813.9745	0.0000	1.0000
K2PI	57	-2851.0610	5828.7403	0.0000	1.0000
SYMI	61	-2848.7492	5834.0445	0.0000	1.0000
F81G	59	-2868.6822	5868.9277	0.0000	1.0000
F81IG	60	-2868.5866	5871.2231	0.0000	1.0000
GTR	63	-2867.9886	5877.5449	0.0000	1.0000
F81I	59	-2879.0767	5889.7166	0.0000	1.0000
HKY	59	-2880.2190	5892.0011	0.0000	1.0000
JC69G	56	-2884.3327	5892.8254	0.0000	1.0000
JC69IG	57	-2889.7375	5906.0934	0.0000	1.0000
JC69I	56	-2900.6583	5925.4766	0.0000	1.0000
K2P	56	-2911.3892	5946.9385	0.0000	1.0000
SYM	60	-2906.9704	5947.9907	0.0000	1.0000
F81	58	-2929.1382	5987.3623	0.0000	1.0000
JC69	55	-2951.6959	6025.1029	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2800.2364	5728.4728	0.6763	0.6763
GTRIG	65	-2799.9747	5729.9495	0.3232	0.9995
HKYG	60	-2811.9109	5743.8218	0.0003	0.9998
HKYIG	61	-2811.5427	5745.0854	0.0002	1.0000
GTRI	64	-2811.7246	5751.4493	0.0000	1.0000
HKYI	60	-2821.8756	5763.7511	0.0000	1.0000
SYMG	61	-2836.2051	5794.4101	0.0000	1.0000
K2PG	57	-2840.7943	5795.5887	0.0000	1.0000
K2PIG	58	-2840.2826	5796.5653	0.0000	1.0000
SYMIG	62	-2837.4612	5798.9225	0.0000	1.0000
K2PI	57	-2851.0610	5816.1220	0.0000	1.0000
SYMI	61	-2848.7492	5819.4984	0.0000	1.0000
F81G	59	-2868.6822	5855.3645	0.0000	1.0000
F81IG	60	-2868.5866	5857.1732	0.0000	1.0000
GTR	63	-2867.9886	5861.9773	0.0000	1.0000
F81I	59	-2879.0767	5876.1534	0.0000	1.0000
HKY	59	-2880.2190	5878.4379	0.0000	1.0000
JC69G	56	-2884.3327	5880.6654	0.0000	1.0000
JC69IG	57	-2889.7375	5893.4751	0.0000	1.0000
JC69I	56	-2900.6583	5913.3166	0.0000	1.0000
SYM	60	-2906.9704	5933.9408	0.0000	1.0000
K2P	56	-2911.3892	5934.7785	0.0000	1.0000
F81	58	-2929.1382	5974.2763	0.0000	1.0000
JC69	55	-2951.6959	6013.3919	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2811.9109	6005.8100	0.7014	0.7014
GTRG	64	-2800.2364	6007.9269	0.2434	0.9448
HKYIG	61	-2811.5427	6011.4401	0.0420	0.9869
GTRIG	65	-2799.9747	6013.7701	0.0131	1.0000
HKYI	60	-2821.8756	6025.7394	0.0000	1.0000
GTRI	64	-2811.7246	6030.9034	0.0000	1.0000
K2PG	57	-2840.7943	6044.4775	0.0000	1.0000
K2PIG	58	-2840.2826	6049.8205	0.0000	1.0000
SYMG	61	-2836.2051	6060.7648	0.0000	1.0000
K2PI	57	-2851.0610	6065.0108	0.0000	1.0000
SYMIG	62	-2837.4612	6069.6436	0.0000	1.0000
SYMI	61	-2848.7492	6085.8531	0.0000	1.0000
F81G	59	-2868.6822	6112.9863	0.0000	1.0000
F81IG	60	-2868.5866	6119.1615	0.0000	1.0000
JC69G	56	-2884.3327	6125.1878	0.0000	1.0000
F81I	59	-2879.0767	6133.7751	0.0000	1.0000
HKY	59	-2880.2190	6136.0597	0.0000	1.0000
GTR	63	-2867.9886	6137.0649	0.0000	1.0000
JC69IG	57	-2889.7375	6142.3639	0.0000	1.0000
JC69I	56	-2900.6583	6157.8389	0.0000	1.0000
K2P	56	-2911.3892	6179.3008	0.0000	1.0000
SYM	60	-2906.9704	6195.9290	0.0000	1.0000
F81	58	-2929.1382	6227.5316	0.0000	1.0000
JC69	55	-2951.6959	6253.5478	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/8760_hg18.phylip = 1 - 582;
 Partition Dummy = 1:261-aggressive-phylip/8760_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/8760_hg18.phylip = 1 - 582;
 Partition Dummy = 1:261-aggressive-phylip/8760_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/8760_hg18.phylip = 1 - 582;
 Partition Dummy = 1:261-aggressive-phylip/8760_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

