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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-1773.6797	3678.0924	0.6978	0.6978
HKYIG	61	-1773.6977	3680.5025	0.2091	0.9069
GTRG	64	-1771.1831	3682.6376	0.0719	0.9788
GTRIG	65	-1771.2012	3685.0760	0.0212	1.0000
K2PG	57	-1800.9200	3725.4925	0.0000	1.0000
K2PIG	58	-1800.9382	3727.8822	0.0000	1.0000
SYMG	61	-1798.7590	3730.6251	0.0000	1.0000
SYMIG	62	-1798.7771	3733.0424	0.0000	1.0000
HKY	59	-1822.4612	3773.2883	0.0000	1.0000
HKYI	60	-1822.4733	3775.6798	0.0000	1.0000
GTR	63	-1820.3421	3778.5605	0.0000	1.0000
GTRI	64	-1820.3546	3780.9806	0.0000	1.0000
F81G	59	-1826.5346	3781.4353	0.0000	1.0000
F81IG	60	-1826.5529	3783.8388	0.0000	1.0000
JC69G	56	-1847.5148	3816.3358	0.0000	1.0000
JC69IG	57	-1847.5331	3818.7187	0.0000	1.0000
K2P	56	-1849.8133	3820.9328	0.0000	1.0000
K2PI	57	-1849.8259	3823.3044	0.0000	1.0000
SYM	60	-1847.6376	3826.0084	0.0000	1.0000
SYMI	61	-1847.6502	3828.4075	0.0000	1.0000
F81	58	-1874.1886	3874.3830	0.0000	1.0000
F81I	59	-1874.2010	3876.7681	0.0000	1.0000
JC69	55	-1895.4673	3909.9011	0.0000	1.0000
JC69I	56	-1895.4797	3912.2656	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1773.6797	3667.3593	0.6010	0.6010
HKYIG	61	-1773.6977	3669.3953	0.2171	0.8181
GTRG	64	-1771.1831	3670.3662	0.1336	0.9517
GTRIG	65	-1771.2012	3672.4025	0.0483	1.0000
K2PG	57	-1800.9200	3715.8399	0.0000	1.0000
K2PIG	58	-1800.9382	3717.8763	0.0000	1.0000
SYMG	61	-1798.7590	3719.5179	0.0000	1.0000
SYMIG	62	-1798.7771	3721.5542	0.0000	1.0000
HKY	59	-1822.4612	3762.9223	0.0000	1.0000
HKYI	60	-1822.4733	3764.9466	0.0000	1.0000
GTR	63	-1820.3421	3766.6843	0.0000	1.0000
GTRI	64	-1820.3546	3768.7093	0.0000	1.0000
F81G	59	-1826.5346	3771.0693	0.0000	1.0000
F81IG	60	-1826.5529	3773.1057	0.0000	1.0000
JC69G	56	-1847.5148	3807.0296	0.0000	1.0000
JC69IG	57	-1847.5331	3809.0661	0.0000	1.0000
K2P	56	-1849.8133	3811.6267	0.0000	1.0000
K2PI	57	-1849.8259	3813.6518	0.0000	1.0000
SYM	60	-1847.6376	3815.2752	0.0000	1.0000
SYMI	61	-1847.6502	3817.3003	0.0000	1.0000
F81	58	-1874.1886	3864.3771	0.0000	1.0000
F81I	59	-1874.2010	3866.4021	0.0000	1.0000
JC69	55	-1895.4673	3900.9346	0.0000	1.0000
JC69I	56	-1895.4797	3902.9595	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1773.6797	3944.0011	0.9652	0.9652
HKYIG	61	-1773.6977	3950.6478	0.0348	1.0000
GTRG	64	-1771.1831	3965.4507	0.0000	1.0000
GTRIG	65	-1771.2012	3972.0977	0.0000	1.0000
K2PG	57	-1800.9200	3978.6496	0.0000	1.0000
K2PIG	58	-1800.9382	3985.2967	0.0000	1.0000
SYMG	61	-1798.7590	4000.7704	0.0000	1.0000
SYMIG	62	-1798.7771	4007.4173	0.0000	1.0000
HKY	59	-1822.4612	4034.9534	0.0000	1.0000
HKYI	60	-1822.4733	4041.5884	0.0000	1.0000
F81G	59	-1826.5346	4043.1003	0.0000	1.0000
F81IG	60	-1826.5529	4049.7475	0.0000	1.0000
GTR	63	-1820.3421	4057.1581	0.0000	1.0000
GTRI	64	-1820.3546	4063.7938	0.0000	1.0000
JC69G	56	-1847.5148	4065.2286	0.0000	1.0000
K2P	56	-1849.8133	4069.8257	0.0000	1.0000
JC69IG	57	-1847.5331	4071.8758	0.0000	1.0000
K2PI	57	-1849.8259	4076.4615	0.0000	1.0000
SYM	60	-1847.6376	4091.9170	0.0000	1.0000
SYMI	61	-1847.6502	4098.5528	0.0000	1.0000
F81	58	-1874.1886	4131.7975	0.0000	1.0000
F81I	59	-1874.2010	4138.4331	0.0000	1.0000
JC69	55	-1895.4673	4154.5229	0.0000	1.0000
JC69I	56	-1895.4797	4161.1585	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/9406_hg18.phylip = 1 - 743;
 Partition Dummy = 1:261-aggressive-phylip/9406_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 AICc model (HKYG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/9406_hg18.phylip = 1 - 743;
 Partition Dummy = 1:261-aggressive-phylip/9406_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 (HKYG)]

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

