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

Input data from file "261-aggressive-phylip/3719_hg18.phylip" (29 taxa, 752 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	-1972.1718	4084.4541	0.2416	0.2416
GTRI	64	-1972.2981	4084.7069	0.2129	0.4545
HKYG	60	-1977.0646	4084.7225	0.2112	0.6657
HKYI	60	-1977.1335	4084.8603	0.1972	0.8629
GTRIG	65	-1972.1772	4086.8617	0.0725	0.9354
HKYIG	61	-1977.0782	4087.1188	0.0637	0.9992
K2PG	57	-1987.2982	4098.1238	0.0003	0.9994
K2PI	57	-1987.3028	4098.1329	0.0003	0.9997
SYMG	61	-1983.5537	4100.0698	0.0001	0.9998
SYMI	61	-1983.5651	4100.0925	0.0001	0.9999
K2PIG	58	-1987.3006	4100.4772	0.0001	1.0000
SYMIG	62	-1983.5565	4102.4512	0.0000	1.0000
HKY	59	-1988.2110	4104.6533	0.0000	1.0000
GTR	63	-1983.7258	4105.1724	0.0000	1.0000
K2P	56	-1998.0786	4117.3429	0.0000	1.0000
SYM	60	-1994.2632	4119.1197	0.0000	1.0000
F81I	59	-1997.8657	4123.9626	0.0000	1.0000
F81G	59	-1997.9755	4124.1823	0.0000	1.0000
F81IG	60	-1997.9859	4126.5652	0.0000	1.0000
JC69I	56	-2006.6557	4134.4969	0.0000	1.0000
JC69G	56	-2006.7233	4134.6321	0.0000	1.0000
JC69IG	57	-2006.7244	4136.9762	0.0000	1.0000
F81	58	-2008.5338	4142.9436	0.0000	1.0000
JC69	55	-2017.2790	4153.4086	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1972.1718	4072.3435	0.3136	0.3136
GTRI	64	-1972.2981	4072.5962	0.2764	0.5900
HKYG	60	-1977.0646	4074.1292	0.1284	0.7184
HKYI	60	-1977.1335	4074.2669	0.1199	0.8383
GTRIG	65	-1972.1772	4074.3544	0.1147	0.9530
HKYIG	61	-1977.0782	4076.1565	0.0466	0.9996
K2PG	57	-1987.2982	4088.5964	0.0001	0.9997
K2PI	57	-1987.3028	4088.6056	0.0001	0.9998
SYMG	61	-1983.5537	4089.1075	0.0001	0.9999
SYMI	61	-1983.5651	4089.1302	0.0001	0.9999
K2PIG	58	-1987.3006	4090.6013	0.0000	1.0000
SYMIG	62	-1983.5565	4091.1131	0.0000	1.0000
GTR	63	-1983.7258	4093.4515	0.0000	1.0000
HKY	59	-1988.2110	4094.4221	0.0000	1.0000
K2P	56	-1998.0786	4108.1572	0.0000	1.0000
SYM	60	-1994.2632	4108.5264	0.0000	1.0000
F81I	59	-1997.8657	4113.7314	0.0000	1.0000
F81G	59	-1997.9755	4113.9511	0.0000	1.0000
F81IG	60	-1997.9859	4115.9719	0.0000	1.0000
JC69I	56	-2006.6557	4125.3113	0.0000	1.0000
JC69G	56	-2006.7233	4125.4465	0.0000	1.0000
JC69IG	57	-2006.7244	4127.4489	0.0000	1.0000
F81	58	-2008.5338	4133.0677	0.0000	1.0000
JC69	55	-2017.2790	4144.5580	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1977.0646	4351.4934	0.2876	0.2876
HKYI	60	-1977.1335	4351.6311	0.2685	0.5561
K2PG	57	-1987.2982	4352.0924	0.2132	0.7693
K2PI	57	-1987.3028	4352.1015	0.2122	0.9815
HKYIG	61	-1977.0782	4358.1434	0.0103	0.9918
K2PIG	58	-1987.3006	4358.7200	0.0078	0.9996
K2P	56	-1998.0786	4367.0305	0.0001	0.9997
HKY	59	-1988.2110	4367.1635	0.0001	0.9998
GTRG	64	-1972.1718	4368.1986	0.0001	0.9999
GTRI	64	-1972.2981	4368.4514	0.0001	1.0000
SYMG	61	-1983.5537	4371.0944	0.0000	1.0000
SYMI	61	-1983.5651	4371.1171	0.0000	1.0000
GTRIG	65	-1972.1772	4374.8323	0.0000	1.0000
SYMIG	62	-1983.5565	4377.7227	0.0000	1.0000
JC69I	56	-2006.6557	4384.1846	0.0000	1.0000
JC69G	56	-2006.7233	4384.3197	0.0000	1.0000
GTR	63	-1983.7258	4384.6839	0.0000	1.0000
SYM	60	-1994.2632	4385.8906	0.0000	1.0000
F81I	59	-1997.8657	4386.4728	0.0000	1.0000
F81G	59	-1997.9755	4386.6925	0.0000	1.0000
JC69IG	57	-2006.7244	4390.9448	0.0000	1.0000
F81IG	60	-1997.9859	4393.3360	0.0000	1.0000
JC69	55	-2017.2790	4398.8085	0.0000	1.0000
F81	58	-2008.5338	4401.1864	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/3719_hg18.phylip = 1 - 752;
 Partition Dummy = 1:261-aggressive-phylip/3719_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/3719_hg18.phylip = 1 - 752;
 Partition Dummy = 1:261-aggressive-phylip/3719_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/3719_hg18.phylip = 1 - 752;
 Partition Dummy = 1:261-aggressive-phylip/3719_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

