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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
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
HKYI	16	-1957.2864	3947.5090	0.4110	0.4110
HKYG	16	-1957.4141	3947.7646	0.3617	0.7727
HKYIG	17	-1957.1972	3949.4495	0.1558	0.9284
GTRI	20	-1955.5192	3952.4942	0.0340	0.9624
GTRG	20	-1955.8087	3953.0732	0.0254	0.9879
GTRIG	21	-1955.4760	3954.5562	0.0121	1.0000
K2PI	13	-1976.6520	3979.9274	0.0000	1.0000
K2PG	13	-1976.8767	3980.3767	0.0000	1.0000
SYMI	17	-1973.0122	3981.0796	0.0000	1.0000
K2PIG	14	-1976.5930	3981.9065	0.0000	1.0000
SYMG	17	-1973.5598	3982.1748	0.0000	1.0000
SYMIG	18	-1973.0044	3983.1901	0.0000	1.0000
HKY	15	-1979.0392	3988.9032	0.0000	1.0000
GTR	19	-1977.6736	3994.6621	0.0000	1.0000
K2P	12	-1997.3225	4019.1783	0.0000	1.0000
SYM	16	-1994.3955	4021.7274	0.0000	1.0000
F81I	15	-2001.3408	4033.5063	0.0000	1.0000
F81G	15	-2001.4737	4033.7721	0.0000	1.0000
F81IG	16	-2001.2678	4035.4720	0.0000	1.0000
JC69I	12	-2016.2832	4057.0997	0.0000	1.0000
JC69G	12	-2016.4432	4057.4197	0.0000	1.0000
JC69IG	13	-2016.2080	4059.0393	0.0000	1.0000
F81	14	-2021.0524	4070.8253	0.0000	1.0000
JC69	11	-2035.7604	4093.9713	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1957.2864	3946.5727	0.3983	0.3983
HKYG	16	-1957.4141	3946.8282	0.3505	0.7487
HKYIG	17	-1957.1972	3948.3943	0.1602	0.9089
GTRI	20	-1955.5192	3951.0384	0.0427	0.9516
GTRG	20	-1955.8087	3951.6174	0.0320	0.9836
GTRIG	21	-1955.4760	3952.9520	0.0164	1.0000
K2PI	13	-1976.6520	3979.3041	0.0000	1.0000
K2PG	13	-1976.8767	3979.7534	0.0000	1.0000
SYMI	17	-1973.0122	3980.0244	0.0000	1.0000
SYMG	17	-1973.5598	3981.1196	0.0000	1.0000
K2PIG	14	-1976.5930	3981.1860	0.0000	1.0000
SYMIG	18	-1973.0044	3982.0088	0.0000	1.0000
HKY	15	-1979.0392	3988.0785	0.0000	1.0000
GTR	19	-1977.6736	3993.3472	0.0000	1.0000
K2P	12	-1997.3225	4018.6450	0.0000	1.0000
SYM	16	-1994.3955	4020.7911	0.0000	1.0000
F81I	15	-2001.3408	4032.6816	0.0000	1.0000
F81G	15	-2001.4737	4032.9473	0.0000	1.0000
F81IG	16	-2001.2678	4034.5356	0.0000	1.0000
JC69I	12	-2016.2832	4056.5664	0.0000	1.0000
JC69G	12	-2016.4432	4056.8864	0.0000	1.0000
JC69IG	13	-2016.2080	4058.4160	0.0000	1.0000
F81	14	-2021.0524	4070.1048	0.0000	1.0000
JC69	11	-2035.7604	4093.5207	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1957.2864	4016.8702	0.5195	0.5195
HKYG	16	-1957.4141	4017.1257	0.4572	0.9767
HKYIG	17	-1957.1972	4023.0854	0.0232	0.9999
K2PI	13	-1976.6520	4036.4208	0.0000	1.0000
K2PG	13	-1976.8767	4036.8701	0.0000	1.0000
GTRI	20	-1955.5192	4038.9102	0.0000	1.0000
GTRG	20	-1955.8087	4039.4892	0.0000	1.0000
K2PIG	14	-1976.5930	4042.6963	0.0000	1.0000
GTRIG	21	-1955.4760	4045.2174	0.0000	1.0000
HKY	15	-1979.0392	4053.9823	0.0000	1.0000
SYMI	17	-1973.0122	4054.7154	0.0000	1.0000
SYMG	17	-1973.5598	4055.8106	0.0000	1.0000
SYMIG	18	-1973.0044	4061.0934	0.0000	1.0000
K2P	12	-1997.3225	4071.3681	0.0000	1.0000
GTR	19	-1977.6736	4076.8255	0.0000	1.0000
SYM	16	-1994.3955	4091.0885	0.0000	1.0000
F81I	15	-2001.3408	4098.5854	0.0000	1.0000
F81G	15	-2001.4737	4098.8512	0.0000	1.0000
F81IG	16	-2001.2678	4104.8331	0.0000	1.0000
JC69I	12	-2016.2832	4109.2895	0.0000	1.0000
JC69G	12	-2016.4432	4109.6094	0.0000	1.0000
JC69IG	13	-2016.2080	4115.5327	0.0000	1.0000
F81	14	-2021.0524	4131.6151	0.0000	1.0000
JC69	11	-2035.7604	4141.8502	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr11_4723.phylip = 1 - 598;
 Partition Dummy = 1:591-tomodel-phylip/chr11_4723.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr11_4723.phylip = 1 - 598;
 Partition Dummy = 1:591-tomodel-phylip/chr11_4723.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr11_4723.phylip = 1 - 598;
 Partition Dummy = 1:591-tomodel-phylip/chr11_4723.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

