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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	65	-2863.4776	5866.2813	0.9028	0.9028
GTRG	64	-2866.8866	5870.8068	0.0940	0.9968
GTRI	64	-2870.5239	5878.0814	0.0025	0.9993
HKYIG	61	-2875.3490	5880.8842	0.0006	0.9999
HKYG	60	-2878.2204	5884.3544	0.0001	1.0000
HKYI	60	-2881.2956	5890.5047	0.0000	1.0000
SYMIG	62	-2886.3712	5905.2060	0.0000	1.0000
SYMG	61	-2889.1876	5908.5614	0.0000	1.0000
K2PIG	58	-2894.9944	5913.3718	0.0000	1.0000
SYMI	61	-2892.2971	5914.7804	0.0000	1.0000
K2PG	57	-2897.6179	5916.3608	0.0000	1.0000
K2PI	57	-2900.3254	5921.7758	0.0000	1.0000
F81IG	60	-2907.4581	5942.8296	0.0000	1.0000
F81G	59	-2909.9626	5945.5710	0.0000	1.0000
F81I	59	-2913.0161	5951.6780	0.0000	1.0000
JC69IG	57	-2924.0003	5969.1257	0.0000	1.0000
JC69G	56	-2926.4045	5971.6809	0.0000	1.0000
JC69I	56	-2929.3770	5977.6258	0.0000	1.0000
GTR	63	-2921.5020	5977.7502	0.0000	1.0000
HKY	59	-2930.4079	5986.4616	0.0000	1.0000
SYM	60	-2941.4906	6010.8946	0.0000	1.0000
K2P	56	-2948.6389	6016.1497	0.0000	1.0000
F81	58	-2959.6738	6042.7306	0.0000	1.0000
JC69	55	-2975.7457	6068.1151	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-2863.4776	5856.9553	0.9151	0.9151
GTRG	64	-2866.8866	5861.7731	0.0823	0.9974
GTRI	64	-2870.5239	5869.0477	0.0022	0.9996
HKYIG	61	-2875.3490	5872.6981	0.0003	0.9999
HKYG	60	-2878.2204	5876.4409	0.0001	1.0000
HKYI	60	-2881.2956	5882.5912	0.0000	1.0000
SYMIG	62	-2886.3712	5896.7423	0.0000	1.0000
SYMG	61	-2889.1876	5900.3752	0.0000	1.0000
K2PIG	58	-2894.9944	5905.9889	0.0000	1.0000
SYMI	61	-2892.2971	5906.5943	0.0000	1.0000
K2PG	57	-2897.6179	5909.2358	0.0000	1.0000
K2PI	57	-2900.3254	5914.6508	0.0000	1.0000
F81IG	60	-2907.4581	5934.9161	0.0000	1.0000
F81G	59	-2909.9626	5937.9252	0.0000	1.0000
F81I	59	-2913.0161	5944.0322	0.0000	1.0000
JC69IG	57	-2924.0003	5962.0007	0.0000	1.0000
JC69G	56	-2926.4045	5964.8090	0.0000	1.0000
GTR	63	-2921.5020	5969.0040	0.0000	1.0000
JC69I	56	-2929.3770	5970.7539	0.0000	1.0000
HKY	59	-2930.4079	5978.8158	0.0000	1.0000
SYM	60	-2941.4906	6002.9811	0.0000	1.0000
K2P	56	-2948.6389	6009.2778	0.0000	1.0000
F81	58	-2959.6738	6035.3477	0.0000	1.0000
JC69	55	-2975.7457	6061.4915	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2878.2204	6170.0603	0.5618	0.5618
HKYIG	61	-2875.3490	6171.2111	0.3160	0.8778
GTRG	64	-2866.8866	6174.9671	0.0483	0.9262
GTRIG	65	-2863.4776	6175.0429	0.0465	0.9727
HKYI	60	-2881.2956	6176.2106	0.0259	0.9986
GTRI	64	-2870.5239	6182.2417	0.0013	0.9999
K2PG	57	-2897.6179	6188.1742	0.0001	1.0000
K2PIG	58	-2894.9944	6189.8209	0.0000	1.0000
K2PI	57	-2900.3254	6193.5892	0.0000	1.0000
SYMG	61	-2889.1876	6198.8883	0.0000	1.0000
SYMIG	62	-2886.3712	6200.1490	0.0000	1.0000
SYMI	61	-2892.2971	6205.1073	0.0000	1.0000
F81G	59	-2909.9626	6226.6509	0.0000	1.0000
F81IG	60	-2907.4581	6228.5355	0.0000	1.0000
F81I	59	-2913.0161	6232.7579	0.0000	1.0000
JC69G	56	-2926.4045	6238.8537	0.0000	1.0000
JC69IG	57	-2924.0003	6240.9391	0.0000	1.0000
JC69I	56	-2929.3770	6244.7987	0.0000	1.0000
HKY	59	-2930.4079	6267.5416	0.0000	1.0000
GTR	63	-2921.5020	6277.3044	0.0000	1.0000
K2P	56	-2948.6389	6283.3226	0.0000	1.0000
SYM	60	-2941.4906	6296.6005	0.0000	1.0000
F81	58	-2959.6738	6319.1797	0.0000	1.0000
JC69	55	-2975.7457	6330.6426	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/3695_hg18.phylip = 1 - 986;
 Partition Dummy = 1:261-aggressive-phylip/3695_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 (GTRIG)]

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

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

