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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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	-1900.8295	3952.2065	0.6594	0.6594
GTRIG	65	-1900.1219	3953.5590	0.3353	0.9947
GTRI	64	-1905.6597	3961.8669	0.0053	1.0000
SYMG	61	-1915.9597	3974.2527	0.0000	1.0000
SYMIG	62	-1915.2405	3975.5377	0.0000	1.0000
HKYG	60	-1920.1151	3979.8549	0.0000	1.0000
HKYIG	61	-1919.3462	3981.0258	0.0000	1.0000
SYMI	61	-1921.3314	3984.9961	0.0000	1.0000
K2PG	57	-1927.7206	3987.0264	0.0000	1.0000
K2PIG	58	-1927.0481	3988.3468	0.0000	1.0000
HKYI	60	-1925.0992	3989.8231	0.0000	1.0000
K2PI	57	-1933.1695	3997.9240	0.0000	1.0000
GTR	63	-1946.6046	4041.0038	0.0000	1.0000
SYM	60	-1958.7763	4057.1773	0.0000	1.0000
HKY	59	-1965.4079	4067.7462	0.0000	1.0000
K2P	56	-1972.0277	4072.9892	0.0000	1.0000
F81G	59	-1971.2696	4079.4698	0.0000	1.0000
JC69G	56	-1975.3416	4079.6168	0.0000	1.0000
F81IG	60	-1970.5865	4080.7976	0.0000	1.0000
JC69IG	57	-1974.6877	4080.9604	0.0000	1.0000
F81I	59	-1976.3454	4089.6213	0.0000	1.0000
JC69I	56	-1980.6573	4090.2484	0.0000	1.0000
F81	58	-2014.5279	4163.3065	0.0000	1.0000
JC69	55	-2018.6773	4163.6509	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1900.8295	3929.6591	0.5700	0.5700
GTRIG	65	-1900.1219	3930.2438	0.4255	0.9954
GTRI	64	-1905.6597	3939.3195	0.0046	1.0000
SYMG	61	-1915.9597	3953.9194	0.0000	1.0000
SYMIG	62	-1915.2405	3954.4811	0.0000	1.0000
HKYG	60	-1920.1151	3960.2302	0.0000	1.0000
HKYIG	61	-1919.3462	3960.6925	0.0000	1.0000
SYMI	61	-1921.3314	3964.6628	0.0000	1.0000
K2PG	57	-1927.7206	3969.4413	0.0000	1.0000
K2PIG	58	-1927.0481	3970.0961	0.0000	1.0000
HKYI	60	-1925.0992	3970.1984	0.0000	1.0000
K2PI	57	-1933.1695	3980.3389	0.0000	1.0000
GTR	63	-1946.6046	4019.2092	0.0000	1.0000
SYM	60	-1958.7763	4037.5526	0.0000	1.0000
HKY	59	-1965.4079	4048.8157	0.0000	1.0000
K2P	56	-1972.0277	4056.0555	0.0000	1.0000
F81G	59	-1971.2696	4060.5393	0.0000	1.0000
F81IG	60	-1970.5865	4061.1730	0.0000	1.0000
JC69G	56	-1975.3416	4062.6831	0.0000	1.0000
JC69IG	57	-1974.6877	4063.3753	0.0000	1.0000
F81I	59	-1976.3454	4070.6908	0.0000	1.0000
JC69I	56	-1980.6573	4073.3147	0.0000	1.0000
F81	58	-2014.5279	4145.0558	0.0000	1.0000
JC69	55	-2018.6773	4147.3546	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-1900.8295	4190.3339	0.8986	0.8986
GTRIG	65	-1900.1219	4194.9917	0.0875	0.9861
GTRI	64	-1905.6597	4199.9943	0.0072	0.9933
K2PG	57	-1927.7206	4201.6048	0.0032	0.9965
SYMG	61	-1915.9597	4202.3751	0.0022	0.9987
HKYG	60	-1920.1151	4204.6129	0.0007	0.9994
K2PIG	58	-1927.0481	4206.3327	0.0003	0.9997
SYMIG	62	-1915.2405	4207.0098	0.0002	0.9999
HKYIG	61	-1919.3462	4209.1482	0.0001	1.0000
K2PI	57	-1933.1695	4212.5025	0.0000	1.0000
SYMI	61	-1921.3314	4213.1185	0.0000	1.0000
HKYI	60	-1925.0992	4214.5811	0.0000	1.0000
GTR	63	-1946.6046	4275.8110	0.0000	1.0000
SYM	60	-1958.7763	4281.9353	0.0000	1.0000
K2P	56	-1972.0277	4284.1460	0.0000	1.0000
HKY	59	-1965.4079	4289.1253	0.0000	1.0000
JC69G	56	-1975.3416	4290.7736	0.0000	1.0000
JC69IG	57	-1974.6877	4295.5388	0.0000	1.0000
F81G	59	-1971.2696	4300.8489	0.0000	1.0000
JC69I	56	-1980.6573	4301.4052	0.0000	1.0000
F81IG	60	-1970.5865	4305.5557	0.0000	1.0000
F81I	59	-1976.3454	4311.0004	0.0000	1.0000
JC69	55	-2018.6773	4371.3721	0.0000	1.0000
F81	58	-2014.5279	4381.2924	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/12234_hg18.phylip = 1 - 434;
 Partition Dummy = 1:261-aggressive-phylip/12234_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/12234_hg18.phylip = 1 - 434;
 Partition Dummy = 1:261-aggressive-phylip/12234_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 (GTRG)]

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


-------------------------------------------------------------
End of Output

