
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:09 04/04/2011

Input data from file "683-tomodel-phylip/chr11_2599.phylip" (7 taxa, 830 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	16	-2027.9006	4088.4704	0.4700	0.4700
HKYIG	17	-2027.4207	4089.5952	0.2679	0.7379
HKYI	16	-2028.7986	4090.2663	0.1915	0.9294
GTRG	20	-2026.2702	4093.5787	0.0365	0.9659
GTRIG	21	-2025.7936	4094.7308	0.0205	0.9865
GTRI	20	-2027.2647	4095.5676	0.0135	1.0000
K2PG	13	-2043.9390	4114.3240	0.0000	1.0000
K2PIG	14	-2043.4305	4115.3764	0.0000	1.0000
K2PI	13	-2044.8167	4116.0796	0.0000	1.0000
SYMG	17	-2041.9818	4118.7174	0.0000	1.0000
SYMIG	18	-2041.4269	4119.6971	0.0000	1.0000
SYMI	17	-2042.9989	4120.7515	0.0000	1.0000
HKY	15	-2045.4862	4121.5621	0.0000	1.0000
F81G	15	-2047.2919	4125.1734	0.0000	1.0000
GTR	19	-2043.7553	4126.4489	0.0000	1.0000
F81IG	16	-2046.9067	4126.4826	0.0000	1.0000
F81I	15	-2047.9604	4126.5105	0.0000	1.0000
K2P	12	-2061.4693	4147.3205	0.0000	1.0000
JC69G	12	-2062.8386	4150.0591	0.0000	1.0000
JC69I	12	-2063.5083	4151.3986	0.0000	1.0000
JC69IG	13	-2062.5462	4151.5385	0.0000	1.0000
SYM	16	-2059.4566	4151.5822	0.0000	1.0000
F81	14	-2063.8131	4156.1416	0.0000	1.0000
JC69	11	-2079.2890	4180.9008	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2027.9006	4087.8013	0.4576	0.4576
HKYIG	17	-2027.4207	4088.8415	0.2720	0.7296
HKYI	16	-2028.7986	4089.5972	0.1864	0.9160
GTRG	20	-2026.2702	4092.5404	0.0428	0.9588
GTRIG	21	-2025.7936	4093.5872	0.0254	0.9842
GTRI	20	-2027.2647	4094.5293	0.0158	1.0000
K2PG	13	-2043.9390	4113.8779	0.0000	1.0000
K2PIG	14	-2043.4305	4114.8610	0.0000	1.0000
K2PI	13	-2044.8167	4115.6335	0.0000	1.0000
SYMG	17	-2041.9818	4117.9637	0.0000	1.0000
SYMIG	18	-2041.4269	4118.8537	0.0000	1.0000
SYMI	17	-2042.9989	4119.9978	0.0000	1.0000
HKY	15	-2045.4862	4120.9724	0.0000	1.0000
F81G	15	-2047.2919	4124.5838	0.0000	1.0000
GTR	19	-2043.7553	4125.5107	0.0000	1.0000
F81IG	16	-2046.9067	4125.8135	0.0000	1.0000
F81I	15	-2047.9604	4125.9208	0.0000	1.0000
K2P	12	-2061.4693	4146.9386	0.0000	1.0000
JC69G	12	-2062.8386	4149.6772	0.0000	1.0000
SYM	16	-2059.4566	4150.9131	0.0000	1.0000
JC69I	12	-2063.5083	4151.0167	0.0000	1.0000
JC69IG	13	-2062.5462	4151.0924	0.0000	1.0000
F81	14	-2063.8131	4155.6263	0.0000	1.0000
JC69	11	-2079.2890	4180.5781	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2027.9006	4163.3441	0.6815	0.6815
HKYI	16	-2028.7986	4165.1400	0.2776	0.9592
HKYIG	17	-2027.4207	4169.1057	0.0382	0.9974
K2PG	13	-2043.9390	4175.2565	0.0018	0.9992
K2PI	13	-2044.8167	4177.0120	0.0007	0.9999
K2PIG	14	-2043.4305	4180.9610	0.0001	1.0000
GTRG	20	-2026.2702	4186.9689	0.0000	1.0000
GTRI	20	-2027.2647	4188.9578	0.0000	1.0000
HKY	15	-2045.4862	4191.7938	0.0000	1.0000
GTRIG	21	-2025.7936	4192.7372	0.0000	1.0000
F81G	15	-2047.2919	4195.4051	0.0000	1.0000
F81I	15	-2047.9604	4196.7422	0.0000	1.0000
SYMG	17	-2041.9818	4198.2279	0.0000	1.0000
SYMI	17	-2042.9989	4200.2620	0.0000	1.0000
F81IG	16	-2046.9067	4201.3563	0.0000	1.0000
K2P	12	-2061.4693	4203.5957	0.0000	1.0000
SYMIG	18	-2041.4269	4203.8394	0.0000	1.0000
JC69G	12	-2062.8386	4206.3343	0.0000	1.0000
JC69I	12	-2063.5083	4207.6738	0.0000	1.0000
JC69IG	13	-2062.5462	4212.4709	0.0000	1.0000
GTR	19	-2043.7553	4215.2177	0.0000	1.0000
F81	14	-2063.8131	4221.7262	0.0000	1.0000
SYM	16	-2059.4566	4226.4559	0.0000	1.0000
JC69	11	-2079.2890	4232.5137	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 683-tomodel-phylip/chr11_2599.phylip = 1 - 830;
 Partition Dummy = 1:683-tomodel-phylip/chr11_2599.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 683-tomodel-phylip/chr11_2599.phylip = 1 - 830;
 Partition Dummy = 1:683-tomodel-phylip/chr11_2599.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 683-tomodel-phylip/chr11_2599.phylip = 1 - 830;
 Partition Dummy = 1:683-tomodel-phylip/chr11_2599.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

