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

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

Minimum AIC  model: GTRG
Minimum AICc model: HKYG
Minimum BIC  model: K2PG


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
HKYG	60	-1891.7180	3915.4164	0.6405	0.6405
GTRG	64	-1887.7031	3917.1129	0.2742	0.9147
HKYIG	61	-1893.5280	3921.4559	0.0313	0.9460
HKYI	60	-1895.0463	3922.0729	0.0230	0.9690
GTRIG	65	-1889.3467	3922.8519	0.0156	0.9845
GTRI	64	-1891.0178	3923.7424	0.0100	0.9945
K2PG	57	-1900.4685	3925.7058	0.0037	0.9982
K2PI	57	-1901.4282	3927.6252	0.0014	0.9997
SYMG	61	-1898.5625	3931.5250	0.0002	0.9999
K2PIG	58	-1902.7215	3932.6077	0.0001	1.0000
SYMIG	62	-1900.7501	3938.3278	0.0000	1.0000
SYMI	61	-1902.2457	3938.8913	0.0000	1.0000
HKY	59	-1907.6485	3944.8657	0.0000	1.0000
GTR	63	-1903.4910	3946.2451	0.0000	1.0000
K2P	56	-1916.5739	3955.5284	0.0000	1.0000
SYM	60	-1914.3403	3960.6610	0.0000	1.0000
F81G	59	-1923.0458	3975.6602	0.0000	1.0000
F81IG	60	-1925.5611	3983.1026	0.0000	1.0000
JC69G	56	-1930.3998	3983.1801	0.0000	1.0000
F81I	59	-1927.0463	3983.6611	0.0000	1.0000
JC69IG	57	-1930.1135	3984.9957	0.0000	1.0000
JC69I	56	-1931.3684	3985.1172	0.0000	1.0000
F81	58	-1938.7620	4004.6887	0.0000	1.0000
JC69	55	-1946.2887	4012.5773	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1887.7031	3903.4062	0.4562	0.4562
HKYG	60	-1891.7180	3903.4360	0.4494	0.9056
GTRIG	65	-1889.3467	3908.6935	0.0324	0.9380
HKYIG	61	-1893.5280	3909.0559	0.0271	0.9651
GTRI	64	-1891.0178	3910.0356	0.0166	0.9817
HKYI	60	-1895.0463	3910.0925	0.0161	0.9978
K2PG	57	-1900.4685	3914.9371	0.0014	0.9992
K2PI	57	-1901.4282	3916.8565	0.0005	0.9998
SYMG	61	-1898.5625	3919.1250	0.0002	0.9999
K2PIG	58	-1902.7215	3921.4430	0.0001	1.0000
SYMIG	62	-1900.7501	3925.5002	0.0000	1.0000
SYMI	61	-1902.2457	3926.4913	0.0000	1.0000
GTR	63	-1903.4910	3932.9819	0.0000	1.0000
HKY	59	-1907.6485	3933.2971	0.0000	1.0000
K2P	56	-1916.5739	3945.1479	0.0000	1.0000
SYM	60	-1914.3403	3948.6807	0.0000	1.0000
F81G	59	-1923.0458	3964.0916	0.0000	1.0000
F81IG	60	-1925.5611	3971.1222	0.0000	1.0000
F81I	59	-1927.0463	3972.0925	0.0000	1.0000
JC69G	56	-1930.3998	3972.7996	0.0000	1.0000
JC69IG	57	-1930.1135	3974.2270	0.0000	1.0000
JC69I	56	-1931.3684	3974.7367	0.0000	1.0000
F81	58	-1938.7620	3993.5240	0.0000	1.0000
JC69	55	-1946.2887	4002.5773	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1900.4685	4172.0218	0.5666	0.5666
K2PI	57	-1901.4282	4173.9412	0.2170	0.7836
HKYG	60	-1891.7180	4174.0515	0.2054	0.9890
HKYI	60	-1895.0463	4180.7080	0.0074	0.9964
K2PIG	58	-1902.7215	4183.0379	0.0023	0.9987
HKYIG	61	-1893.5280	4184.1817	0.0013	1.0000
GTRG	64	-1887.7031	4192.0627	0.0000	1.0000
SYMG	61	-1898.5625	4194.2507	0.0000	1.0000
K2P	56	-1916.5739	4197.7223	0.0000	1.0000
GTRI	64	-1891.0178	4198.6922	0.0000	1.0000
HKY	59	-1907.6485	4199.4023	0.0000	1.0000
SYMI	61	-1902.2457	4201.6171	0.0000	1.0000
GTRIG	65	-1889.3467	4201.8603	0.0000	1.0000
SYMIG	62	-1900.7501	4205.1362	0.0000	1.0000
GTR	63	-1903.4910	4217.1282	0.0000	1.0000
SYM	60	-1914.3403	4219.2962	0.0000	1.0000
JC69G	56	-1930.3998	4225.3741	0.0000	1.0000
JC69I	56	-1931.3684	4227.3112	0.0000	1.0000
F81G	59	-1923.0458	4230.1968	0.0000	1.0000
JC69IG	57	-1930.1135	4231.3117	0.0000	1.0000
F81I	59	-1927.0463	4238.1977	0.0000	1.0000
F81IG	60	-1925.5611	4241.7377	0.0000	1.0000
JC69	55	-1946.2887	4250.6415	0.0000	1.0000
F81	58	-1938.7620	4255.1190	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/819_hg18.phylip = 1 - 672;
 Partition Dummy = 1:261-aggressive-phylip/819_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 (HKYG)]

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


[Mrbayes block for the best BIC model (K2PG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/819_hg18.phylip = 1 - 672;
 Partition Dummy = 1:261-aggressive-phylip/819_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

