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

Input data from file "261-aggressive-phylip/5142_hg18.phylip" (29 taxa, 655 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	-1906.6665	3955.4347	0.7940	0.7940
GTRIG	65	-1906.7988	3958.1646	0.2028	0.9967
HKYG	60	-1917.0740	3966.4712	0.0032	0.9999
SYMG	61	-1920.2048	3975.1650	0.0000	1.0000
SYMIG	62	-1920.2138	3977.6235	0.0000	1.0000
HKYIG	61	-1921.4685	3977.6925	0.0000	1.0000
K2PG	57	-1927.5055	3980.0863	0.0000	1.0000
K2PIG	58	-1931.6311	3990.7453	0.0000	1.0000
F81G	59	-1949.7933	4029.4857	0.0000	1.0000
GTR	63	-1944.9888	4029.6222	0.0000	1.0000
GTRI	64	-1945.0577	4032.2171	0.0000	1.0000
JC69G	56	-1958.3334	4039.3424	0.0000	1.0000
F81IG	60	-1955.6929	4043.7091	0.0000	1.0000
HKY	59	-1957.0974	4044.0939	0.0000	1.0000
SYM	60	-1959.4005	4051.1242	0.0000	1.0000
JC69IG	57	-1964.2229	4053.5211	0.0000	1.0000
SYMI	61	-1959.4033	4053.5621	0.0000	1.0000
K2P	56	-1965.7664	4054.2083	0.0000	1.0000
HKYI	60	-1961.7079	4055.7390	0.0000	1.0000
K2PI	57	-1971.9336	4068.9425	0.0000	1.0000
F81	58	-1989.9826	4107.4485	0.0000	1.0000
F81I	59	-1992.9158	4115.7307	0.0000	1.0000
JC69	55	-1999.2081	4118.7001	0.0000	1.0000
JC69I	56	-2002.1827	4127.0410	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1906.6665	3941.3330	0.7553	0.7553
GTRIG	65	-1906.7988	3943.5975	0.2434	0.9987
HKYG	60	-1917.0740	3954.1479	0.0012	1.0000
SYMG	61	-1920.2048	3962.4095	0.0000	1.0000
SYMIG	62	-1920.2138	3964.4276	0.0000	1.0000
HKYIG	61	-1921.4685	3964.9370	0.0000	1.0000
K2PG	57	-1927.5055	3969.0110	0.0000	1.0000
K2PIG	58	-1931.6311	3979.2621	0.0000	1.0000
GTR	63	-1944.9888	4015.9776	0.0000	1.0000
F81G	59	-1949.7933	4017.5865	0.0000	1.0000
GTRI	64	-1945.0577	4018.1154	0.0000	1.0000
JC69G	56	-1958.3334	4028.6668	0.0000	1.0000
F81IG	60	-1955.6929	4031.3858	0.0000	1.0000
HKY	59	-1957.0974	4032.1948	0.0000	1.0000
SYM	60	-1959.4005	4038.8010	0.0000	1.0000
SYMI	61	-1959.4033	4040.8066	0.0000	1.0000
JC69IG	57	-1964.2229	4042.4458	0.0000	1.0000
HKYI	60	-1961.7079	4043.4157	0.0000	1.0000
K2P	56	-1965.7664	4043.5328	0.0000	1.0000
K2PI	57	-1971.9336	4057.8671	0.0000	1.0000
F81	58	-1989.9826	4095.9653	0.0000	1.0000
F81I	59	-1992.9158	4103.8316	0.0000	1.0000
JC69	55	-1999.2081	4108.4163	0.0000	1.0000
JC69I	56	-2002.1827	4116.3654	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1917.0740	4223.2261	0.6342	0.6342
K2PG	57	-1927.5055	4224.6352	0.3135	0.9478
GTRG	64	-1906.6665	4228.3496	0.0489	0.9967
GTRIG	65	-1906.7988	4235.0988	0.0017	0.9984
SYMG	61	-1920.2048	4235.9723	0.0011	0.9995
HKYIG	61	-1921.4685	4238.4998	0.0003	0.9998
K2PIG	58	-1931.6311	4239.3710	0.0002	1.0000
SYMIG	62	-1920.2138	4242.4750	0.0000	1.0000
JC69G	56	-1958.3334	4279.8064	0.0000	1.0000
F81G	59	-1949.7933	4282.1800	0.0000	1.0000
K2P	56	-1965.7664	4294.6723	0.0000	1.0000
HKY	59	-1957.0974	4296.7882	0.0000	1.0000
JC69IG	57	-1964.2229	4298.0700	0.0000	1.0000
GTR	63	-1944.9888	4298.5096	0.0000	1.0000
F81IG	60	-1955.6929	4300.4639	0.0000	1.0000
GTRI	64	-1945.0577	4305.1321	0.0000	1.0000
SYM	60	-1959.4005	4307.8791	0.0000	1.0000
HKYI	60	-1961.7079	4312.4938	0.0000	1.0000
K2PI	57	-1971.9336	4313.4913	0.0000	1.0000
SYMI	61	-1959.4033	4314.3693	0.0000	1.0000
JC69	55	-1999.2081	4355.0712	0.0000	1.0000
F81	58	-1989.9826	4356.0741	0.0000	1.0000
JC69I	56	-2002.1827	4367.5050	0.0000	1.0000
F81I	59	-1992.9158	4368.4250	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/5142_hg18.phylip = 1 - 655;
 Partition Dummy = 1:261-aggressive-phylip/5142_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/5142_hg18.phylip = 1 - 655;
 Partition Dummy = 1:261-aggressive-phylip/5142_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/5142_hg18.phylip = 1 - 655;
 Partition Dummy = 1:261-aggressive-phylip/5142_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

