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

Input data from file "261-aggressive-phylip/12659_hg18.phylip" (29 taxa, 998 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	-6581.7809	13300.4793	0.5119	0.5119
GTRIG	65	-6581.5083	13302.2225	0.2141	0.7259
HKYG	60	-6587.2986	13302.4093	0.1950	0.9209
HKYIG	61	-6587.0668	13304.2149	0.0791	1.0000
K2PG	57	-6612.6820	13346.3980	0.0000	1.0000
K2PIG	58	-6612.4777	13348.2439	0.0000	1.0000
SYMG	61	-6610.9409	13351.9629	0.0000	1.0000
SYMIG	62	-6610.8085	13353.9720	0.0000	1.0000
GTRI	64	-6646.9043	13430.7261	0.0000	1.0000
HKYI	60	-6651.6241	13431.0604	0.0000	1.0000
K2PI	57	-6679.1797	13479.3934	0.0000	1.0000
SYMI	61	-6678.2652	13486.6116	0.0000	1.0000
F81G	59	-6777.2019	13679.9518	0.0000	1.0000
F81IG	60	-6777.2693	13682.3507	0.0000	1.0000
JC69G	56	-6792.2229	13703.2300	0.0000	1.0000
JC69IG	57	-6792.2562	13705.5464	0.0000	1.0000
F81I	59	-6838.5716	13802.6912	0.0000	1.0000
JC69I	56	-6853.9682	13826.7207	0.0000	1.0000
HKY	59	-6855.5362	13836.6204	0.0000	1.0000
GTR	63	-6852.2970	13839.2279	0.0000	1.0000
K2P	56	-6879.7651	13878.3146	0.0000	1.0000
SYM	60	-6876.5099	13880.8320	0.0000	1.0000
F81	58	-7033.2922	14189.8730	0.0000	1.0000
JC69	55	-7047.0080	14210.5553	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-6581.7809	13291.5618	0.5543	0.5543
GTRIG	65	-6581.5083	13293.0165	0.2678	0.8221
HKYG	60	-6587.2986	13294.5972	0.1215	0.9436
HKYIG	61	-6587.0668	13296.1337	0.0564	1.0000
K2PG	57	-6612.6820	13339.3640	0.0000	1.0000
K2PIG	58	-6612.4777	13340.9553	0.0000	1.0000
SYMG	61	-6610.9409	13343.8817	0.0000	1.0000
SYMIG	62	-6610.8085	13345.6169	0.0000	1.0000
GTRI	64	-6646.9043	13421.8086	0.0000	1.0000
HKYI	60	-6651.6241	13423.2482	0.0000	1.0000
K2PI	57	-6679.1797	13472.3594	0.0000	1.0000
SYMI	61	-6678.2652	13478.5304	0.0000	1.0000
F81G	59	-6777.2019	13672.4038	0.0000	1.0000
F81IG	60	-6777.2693	13674.5385	0.0000	1.0000
JC69G	56	-6792.2229	13696.4458	0.0000	1.0000
JC69IG	57	-6792.2562	13698.5124	0.0000	1.0000
F81I	59	-6838.5716	13795.1432	0.0000	1.0000
JC69I	56	-6853.9682	13819.9364	0.0000	1.0000
HKY	59	-6855.5362	13829.0724	0.0000	1.0000
GTR	63	-6852.2970	13830.5941	0.0000	1.0000
K2P	56	-6879.7651	13871.5303	0.0000	1.0000
SYM	60	-6876.5099	13873.0198	0.0000	1.0000
F81	58	-7033.2922	14182.5844	0.0000	1.0000
JC69	55	-7047.0080	14204.0160	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-6587.2986	13588.9424	0.9614	0.9614
HKYIG	61	-6587.0668	13595.3846	0.0384	0.9997
GTRG	64	-6581.7809	13605.5300	0.0002	1.0000
GTRIG	65	-6581.5083	13611.8905	0.0000	1.0000
K2PG	57	-6612.6820	13618.9919	0.0000	1.0000
K2PIG	58	-6612.4777	13625.4890	0.0000	1.0000
SYMG	61	-6610.9409	13643.1327	0.0000	1.0000
SYMIG	62	-6610.8085	13649.7736	0.0000	1.0000
HKYI	60	-6651.6241	13717.5934	0.0000	1.0000
GTRI	64	-6646.9043	13735.7768	0.0000	1.0000
K2PI	57	-6679.1797	13751.9873	0.0000	1.0000
SYMI	61	-6678.2652	13777.7814	0.0000	1.0000
F81G	59	-6777.2019	13961.8432	0.0000	1.0000
F81IG	60	-6777.2693	13968.8837	0.0000	1.0000
JC69G	56	-6792.2229	13971.1679	0.0000	1.0000
JC69IG	57	-6792.2562	13978.1403	0.0000	1.0000
F81I	59	-6838.5716	14084.5826	0.0000	1.0000
JC69I	56	-6853.9682	14094.6586	0.0000	1.0000
HKY	59	-6855.5362	14118.5119	0.0000	1.0000
GTR	63	-6852.2970	14139.6565	0.0000	1.0000
K2P	56	-6879.7651	14146.2525	0.0000	1.0000
SYM	60	-6876.5099	14167.3650	0.0000	1.0000
F81	58	-7033.2922	14467.1181	0.0000	1.0000
JC69	55	-7047.0080	14473.8325	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/12659_hg18.phylip = 1 - 998;
 Partition Dummy = 1:261-aggressive-phylip/12659_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/12659_hg18.phylip = 1 - 998;
 Partition Dummy = 1:261-aggressive-phylip/12659_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/12659_hg18.phylip = 1 - 998;
 Partition Dummy = 1:261-aggressive-phylip/12659_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

