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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	21	-2159.0957	4361.3364	0.3370	0.3370
GTRG	20	-2160.3709	4361.7813	0.2698	0.6067
HKYIG	17	-2164.2029	4363.1605	0.1354	0.7421
GTRI	20	-2161.1940	4363.4275	0.1184	0.8605
HKYG	16	-2165.6034	4363.8767	0.0946	0.9552
HKYI	16	-2166.3502	4365.3704	0.0448	1.0000
K2PIG	14	-2181.9267	4392.3694	0.0000	1.0000
K2PG	13	-2183.2144	4392.8755	0.0000	1.0000
K2PI	13	-2183.9345	4394.3157	0.0000	1.0000
SYMG	17	-2180.5607	4395.8760	0.0000	1.0000
SYMIG	18	-2179.5937	4396.0318	0.0000	1.0000
SYMI	17	-2181.1895	4397.1337	0.0000	1.0000
F81IG	16	-2192.2987	4417.2673	0.0000	1.0000
F81G	15	-2193.6090	4417.8085	0.0000	1.0000
F81I	15	-2194.4233	4419.4370	0.0000	1.0000
JC69IG	13	-2207.3227	4441.0919	0.0000	1.0000
JC69G	12	-2208.5216	4441.4255	0.0000	1.0000
JC69I	12	-2209.3004	4442.9832	0.0000	1.0000
GTR	19	-2205.0769	4449.0933	0.0000	1.0000
HKY	15	-2209.7007	4449.9918	0.0000	1.0000
K2P	12	-2226.2226	4476.8276	0.0000	1.0000
SYM	16	-2223.4533	4479.5765	0.0000	1.0000
F81	14	-2236.2648	4501.0455	0.0000	1.0000
JC69	11	-2250.6227	4523.5686	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	21	-2159.0957	4360.1915	0.3637	0.3637
GTRG	20	-2160.3709	4360.7417	0.2762	0.6399
GTRI	20	-2161.1940	4362.3879	0.1213	0.7611
HKYIG	17	-2164.2029	4362.4059	0.1202	0.8813
HKYG	16	-2165.6034	4363.2068	0.0805	0.9618
HKYI	16	-2166.3502	4364.7004	0.0382	1.0000
K2PIG	14	-2181.9267	4391.8534	0.0000	1.0000
K2PG	13	-2183.2144	4392.4289	0.0000	1.0000
K2PI	13	-2183.9345	4393.8691	0.0000	1.0000
SYMG	17	-2180.5607	4395.1214	0.0000	1.0000
SYMIG	18	-2179.5937	4395.1873	0.0000	1.0000
SYMI	17	-2181.1895	4396.3791	0.0000	1.0000
F81IG	16	-2192.2987	4416.5973	0.0000	1.0000
F81G	15	-2193.6090	4417.2181	0.0000	1.0000
F81I	15	-2194.4233	4418.8466	0.0000	1.0000
JC69IG	13	-2207.3227	4440.6453	0.0000	1.0000
JC69G	12	-2208.5216	4441.0431	0.0000	1.0000
JC69I	12	-2209.3004	4442.6008	0.0000	1.0000
GTR	19	-2205.0769	4448.1538	0.0000	1.0000
HKY	15	-2209.7007	4449.4014	0.0000	1.0000
K2P	12	-2226.2226	4476.4452	0.0000	1.0000
SYM	16	-2223.4533	4478.9065	0.0000	1.0000
F81	14	-2236.2648	4500.5295	0.0000	1.0000
JC69	11	-2250.6227	4523.2454	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2165.6034	4438.7303	0.6188	0.6188
HKYI	16	-2166.3502	4440.2240	0.2932	0.9120
HKYIG	17	-2164.2029	4442.6496	0.0872	0.9992
K2PG	13	-2183.2144	4453.7917	0.0003	0.9995
GTRG	20	-2160.3709	4455.1461	0.0002	0.9997
K2PI	13	-2183.9345	4455.2319	0.0002	0.9999
GTRI	20	-2161.1940	4456.7923	0.0001	0.9999
K2PIG	14	-2181.9267	4457.9365	0.0000	1.0000
GTRIG	21	-2159.0957	4459.3161	0.0000	1.0000
SYMG	17	-2180.5607	4475.3651	0.0000	1.0000
SYMI	17	-2181.1895	4476.6228	0.0000	1.0000
SYMIG	18	-2179.5937	4480.1513	0.0000	1.0000
F81G	15	-2193.6090	4488.0214	0.0000	1.0000
F81I	15	-2194.4233	4489.6499	0.0000	1.0000
F81IG	16	-2192.2987	4492.1208	0.0000	1.0000
JC69G	12	-2208.5216	4497.6858	0.0000	1.0000
JC69I	12	-2209.3004	4499.2435	0.0000	1.0000
JC69IG	13	-2207.3227	4502.0082	0.0000	1.0000
HKY	15	-2209.7007	4520.2047	0.0000	1.0000
K2P	12	-2226.2226	4533.0879	0.0000	1.0000
GTR	19	-2205.0769	4537.8380	0.0000	1.0000
SYM	16	-2223.4533	4554.4300	0.0000	1.0000
F81	14	-2236.2648	4566.6126	0.0000	1.0000
JC69	11	-2250.6227	4575.1678	0.0000	1.0000



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

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

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


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

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


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

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

