
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:59 02/17/2011

Input data from file "183-phylip/chr7_5228.phylip" (29 taxa, 654 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	-1966.7784	4075.6824	0.4774	0.4774
GTRIG	65	-1965.8890	4076.3697	0.3386	0.8159
GTRI	64	-1968.1190	4078.3637	0.1249	0.9409
HKYG	60	-1974.3131	4080.9702	0.0339	0.9748
HKYIG	61	-1973.7647	4082.3063	0.0174	0.9922
HKYI	60	-1975.7836	4083.9111	0.0078	1.0000
K2PG	57	-1994.5344	4114.1627	0.0000	1.0000
SYMG	61	-1990.2240	4115.2250	0.0000	1.0000
K2PIG	58	-1993.9031	4115.3088	0.0000	1.0000
SYMIG	62	-1989.0846	4115.3876	0.0000	1.0000
SYMI	61	-1991.0796	4116.9363	0.0000	1.0000
K2PI	57	-1996.0367	4117.1673	0.0000	1.0000
GTR	63	-1990.2916	4120.2509	0.0000	1.0000
HKY	59	-1997.8386	4125.5963	0.0000	1.0000
F81G	59	-2005.7497	4141.4185	0.0000	1.0000
F81IG	60	-2005.2720	4142.8881	0.0000	1.0000
F81I	59	-2007.6614	4145.2420	0.0000	1.0000
K2P	56	-2017.6863	4158.0661	0.0000	1.0000
SYM	60	-2012.9464	4158.2368	0.0000	1.0000
JC69G	56	-2022.7098	4168.1131	0.0000	1.0000
JC69IG	57	-2022.2048	4169.5036	0.0000	1.0000
JC69I	56	-2024.2293	4171.1520	0.0000	1.0000
F81	58	-2027.9350	4183.3726	0.0000	1.0000
JC69	55	-2045.2031	4210.7073	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1966.7784	4061.5568	0.4522	0.4522
GTRIG	65	-1965.8890	4061.7779	0.4049	0.8571
GTRI	64	-1968.1190	4064.2381	0.1183	0.9754
HKYG	60	-1974.3131	4068.6262	0.0132	0.9886
HKYIG	61	-1973.7647	4069.5293	0.0084	0.9970
HKYI	60	-1975.7836	4071.5671	0.0030	1.0000
SYMIG	62	-1989.0846	4102.1693	0.0000	1.0000
SYMG	61	-1990.2240	4102.4480	0.0000	1.0000
K2PG	57	-1994.5344	4103.0688	0.0000	1.0000
K2PIG	58	-1993.9031	4103.8063	0.0000	1.0000
SYMI	61	-1991.0796	4104.1593	0.0000	1.0000
K2PI	57	-1996.0367	4106.0733	0.0000	1.0000
GTR	63	-1990.2916	4106.5831	0.0000	1.0000
HKY	59	-1997.8386	4113.6771	0.0000	1.0000
F81G	59	-2005.7497	4129.4993	0.0000	1.0000
F81IG	60	-2005.2720	4130.5441	0.0000	1.0000
F81I	59	-2007.6614	4133.3228	0.0000	1.0000
SYM	60	-2012.9464	4145.8928	0.0000	1.0000
K2P	56	-2017.6863	4147.3726	0.0000	1.0000
JC69G	56	-2022.7098	4157.4196	0.0000	1.0000
JC69IG	57	-2022.2048	4158.4096	0.0000	1.0000
JC69I	56	-2024.2293	4160.4585	0.0000	1.0000
F81	58	-2027.9350	4171.8701	0.0000	1.0000
JC69	55	-2045.2031	4200.4063	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1974.3131	4337.6126	0.7672	0.7672
HKYI	60	-1975.7836	4340.5536	0.1763	0.9435
HKYIG	61	-1973.7647	4342.9988	0.0519	0.9954
GTRG	64	-1966.7784	4348.4757	0.0034	0.9988
GTRI	64	-1968.1190	4351.1570	0.0009	0.9997
GTRIG	65	-1965.8890	4353.1799	0.0003	1.0000
K2PG	57	-1994.5344	4358.6059	0.0000	1.0000
K2PI	57	-1996.0367	4361.6104	0.0000	1.0000
K2PIG	58	-1993.9031	4363.8265	0.0000	1.0000
SYMG	61	-1990.2240	4375.9176	0.0000	1.0000
SYMI	61	-1991.0796	4377.6288	0.0000	1.0000
HKY	59	-1997.8386	4378.1805	0.0000	1.0000
SYMIG	62	-1989.0846	4380.1219	0.0000	1.0000
GTR	63	-1990.2916	4389.0189	0.0000	1.0000
F81G	59	-2005.7497	4394.0027	0.0000	1.0000
F81I	59	-2007.6614	4397.8262	0.0000	1.0000
K2P	56	-2017.6863	4398.4267	0.0000	1.0000
F81IG	60	-2005.2720	4399.5305	0.0000	1.0000
JC69G	56	-2022.7098	4408.4736	0.0000	1.0000
JC69I	56	-2024.2293	4411.5125	0.0000	1.0000
JC69IG	57	-2022.2048	4413.9467	0.0000	1.0000
SYM	60	-2012.9464	4414.8792	0.0000	1.0000
F81	58	-2027.9350	4431.8903	0.0000	1.0000
JC69	55	-2045.2031	4446.9772	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 183-phylip/chr7_5228.phylip = 1 - 654;
 Partition Dummy = 1:183-phylip/chr7_5228.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 183-phylip/chr7_5228.phylip = 1 - 654;
 Partition Dummy = 1:183-phylip/chr7_5228.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 183-phylip/chr7_5228.phylip = 1 - 654;
 Partition Dummy = 1:183-phylip/chr7_5228.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

