
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:30 04/04/2011

Input data from file "591-tomodel-phylip/chr1_28849.phylip" (7 taxa, 629 characters)

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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
GTRI	20	-1972.1627	3985.7070	0.2617	0.2617
HKYI	16	-1976.6488	3986.1865	0.2059	0.4676
GTRG	20	-1972.4457	3986.2731	0.1972	0.6648
HKYG	16	-1976.8493	3986.5875	0.1685	0.8333
GTRIG	21	-1972.1352	3987.7925	0.0922	0.9256
HKYIG	17	-1976.6099	3988.2215	0.0744	1.0000
HKY	15	-1994.4809	4019.7448	0.0000	1.0000
GTR	19	-1990.3055	4019.8590	0.0000	1.0000
K2PI	13	-2004.8995	4036.3909	0.0000	1.0000
K2PG	13	-2005.0977	4036.7873	0.0000	1.0000
K2PIG	14	-2004.9539	4038.5918	0.0000	1.0000
SYMI	17	-2002.0367	4039.0751	0.0000	1.0000
SYMG	17	-2002.1709	4039.3435	0.0000	1.0000
SYMIG	18	-2002.0005	4041.1224	0.0000	1.0000
F81I	15	-2016.2567	4063.2965	0.0000	1.0000
F81G	15	-2016.5850	4063.9531	0.0000	1.0000
F81IG	16	-2016.2536	4065.3961	0.0000	1.0000
K2P	12	-2022.4677	4069.4419	0.0000	1.0000
SYM	16	-2019.6866	4072.2621	0.0000	1.0000
F81	14	-2032.4345	4093.5530	0.0000	1.0000
JC69I	12	-2039.5406	4103.5877	0.0000	1.0000
JC69G	12	-2039.8195	4104.1454	0.0000	1.0000
JC69IG	13	-2039.5359	4105.6637	0.0000	1.0000
JC69	11	-2055.8036	4134.0350	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1972.1627	3984.3254	0.2869	0.2869
GTRG	20	-1972.4457	3984.8915	0.2162	0.5031
HKYI	16	-1976.6488	3985.2976	0.1765	0.6796
HKYG	16	-1976.8493	3985.6986	0.1444	0.8240
GTRIG	21	-1972.1352	3986.2703	0.1085	0.9325
HKYIG	17	-1976.6099	3987.2199	0.0675	1.0000
GTR	19	-1990.3055	4018.6111	0.0000	1.0000
HKY	15	-1994.4809	4018.9617	0.0000	1.0000
K2PI	13	-2004.8995	4035.7990	0.0000	1.0000
K2PG	13	-2005.0977	4036.1955	0.0000	1.0000
K2PIG	14	-2004.9539	4037.9078	0.0000	1.0000
SYMI	17	-2002.0367	4038.0734	0.0000	1.0000
SYMG	17	-2002.1709	4038.3419	0.0000	1.0000
SYMIG	18	-2002.0005	4040.0011	0.0000	1.0000
F81I	15	-2016.2567	4062.5134	0.0000	1.0000
F81G	15	-2016.5850	4063.1700	0.0000	1.0000
F81IG	16	-2016.2536	4064.5072	0.0000	1.0000
K2P	12	-2022.4677	4068.9354	0.0000	1.0000
SYM	16	-2019.6866	4071.3732	0.0000	1.0000
F81	14	-2032.4345	4092.8689	0.0000	1.0000
JC69I	12	-2039.5406	4103.0812	0.0000	1.0000
JC69G	12	-2039.8195	4103.6389	0.0000	1.0000
JC69IG	13	-2039.5359	4105.0719	0.0000	1.0000
JC69	11	-2055.8036	4133.6071	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1976.6488	4056.4037	0.5376	0.5376
HKYG	16	-1976.8493	4056.8047	0.4399	0.9775
HKYIG	17	-1976.6099	4062.7701	0.0223	0.9998
GTRI	20	-1972.1627	4073.2080	0.0001	0.9999
GTRG	20	-1972.4457	4073.7741	0.0001	1.0000
GTRIG	21	-1972.1352	4079.5971	0.0000	1.0000
HKY	15	-1994.4809	4085.6237	0.0000	1.0000
K2PI	13	-2004.8995	4093.5727	0.0000	1.0000
K2PG	13	-2005.0977	4093.9692	0.0000	1.0000
K2PIG	14	-2004.9539	4100.1256	0.0000	1.0000
GTR	19	-1990.3055	4103.0496	0.0000	1.0000
SYMI	17	-2002.0367	4113.6237	0.0000	1.0000
SYMG	17	-2002.1709	4113.8921	0.0000	1.0000
SYMIG	18	-2002.0005	4119.9955	0.0000	1.0000
K2P	12	-2022.4677	4122.2650	0.0000	1.0000
F81I	15	-2016.2567	4129.1754	0.0000	1.0000
F81G	15	-2016.5850	4129.8320	0.0000	1.0000
F81IG	16	-2016.2536	4135.6133	0.0000	1.0000
SYM	16	-2019.6866	4142.4793	0.0000	1.0000
F81	14	-2032.4345	4155.0868	0.0000	1.0000
JC69I	12	-2039.5406	4156.4108	0.0000	1.0000
JC69G	12	-2039.8195	4156.9685	0.0000	1.0000
JC69IG	13	-2039.5359	4162.8456	0.0000	1.0000
JC69	11	-2055.8036	4182.4926	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_28849.phylip = 1 - 629;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28849.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (GTRI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_28849.phylip = 1 - 629;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28849.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_28849.phylip = 1 - 629;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28849.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

