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

Input data from file "591-tomodel-phylip/chr3_23724.phylip" (7 taxa, 609 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	-1887.4262	3807.7713	0.6123	0.6123
HKYIG	17	-1887.2139	3809.4633	0.2628	0.8751
GTRG	20	-1885.6602	3812.7491	0.0508	0.9259
HKYI	16	-1890.0260	3812.9709	0.0455	0.9714
GTRIG	21	-1885.4010	3814.3761	0.0225	0.9939
GTRI	20	-1887.7909	3817.0104	0.0060	1.0000
K2PG	13	-1901.0442	3828.7002	0.0000	1.0000
K2PIG	14	-1900.8082	3830.3235	0.0000	1.0000
SYMG	17	-1898.7862	3832.6079	0.0000	1.0000
K2PI	13	-1903.5622	3833.7362	0.0000	1.0000
SYMIG	18	-1898.5441	3834.2476	0.0000	1.0000
SYMI	17	-1901.4098	3837.8551	0.0000	1.0000
F81G	15	-1910.3165	3851.4424	0.0000	1.0000
F81IG	16	-1910.8110	3854.5410	0.0000	1.0000
F81I	15	-1912.3818	3855.5731	0.0000	1.0000
HKY	15	-1912.6641	3856.1376	0.0000	1.0000
GTR	19	-1910.3245	3859.9394	0.0000	1.0000
JC69G	12	-1925.6743	3875.8721	0.0000	1.0000
K2P	12	-1925.7363	3875.9962	0.0000	1.0000
JC69IG	13	-1925.5575	3877.7267	0.0000	1.0000
JC69I	12	-1927.3468	3879.2170	0.0000	1.0000
SYM	16	-1923.2204	3879.3597	0.0000	1.0000
F81	14	-1933.3142	3895.3355	0.0000	1.0000
JC69	11	-1947.4782	3917.3987	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1887.4262	3806.8524	0.5882	0.5882
HKYIG	17	-1887.2139	3808.4277	0.2676	0.8558
GTRG	20	-1885.6602	3811.3205	0.0630	0.9188
HKYI	16	-1890.0260	3812.0519	0.0437	0.9625
GTRIG	21	-1885.4010	3812.8020	0.0300	0.9925
GTRI	20	-1887.7909	3815.5818	0.0075	1.0000
K2PG	13	-1901.0442	3828.0885	0.0000	1.0000
K2PIG	14	-1900.8082	3829.6164	0.0000	1.0000
SYMG	17	-1898.7862	3831.5724	0.0000	1.0000
SYMIG	18	-1898.5441	3833.0882	0.0000	1.0000
K2PI	13	-1903.5622	3833.1245	0.0000	1.0000
SYMI	17	-1901.4098	3836.8195	0.0000	1.0000
F81G	15	-1910.3165	3850.6330	0.0000	1.0000
F81IG	16	-1910.8110	3853.6220	0.0000	1.0000
F81I	15	-1912.3818	3854.7637	0.0000	1.0000
HKY	15	-1912.6641	3855.3282	0.0000	1.0000
GTR	19	-1910.3245	3858.6490	0.0000	1.0000
JC69G	12	-1925.6743	3875.3486	0.0000	1.0000
K2P	12	-1925.7363	3875.4727	0.0000	1.0000
JC69IG	13	-1925.5575	3877.1149	0.0000	1.0000
SYM	16	-1923.2204	3878.4408	0.0000	1.0000
JC69I	12	-1927.3468	3878.6936	0.0000	1.0000
F81	14	-1933.3142	3894.6284	0.0000	1.0000
JC69	11	-1947.4782	3916.9565	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1887.4262	3877.4415	0.8733	0.8733
HKYI	16	-1890.0260	3882.6410	0.0649	0.9381
HKYIG	17	-1887.2139	3883.4286	0.0438	0.9819
K2PG	13	-1901.0442	3885.4421	0.0160	0.9979
K2PI	13	-1903.5622	3890.4781	0.0013	0.9992
K2PIG	14	-1900.8082	3891.3819	0.0008	1.0000
GTRG	20	-1885.6602	3899.5568	0.0000	1.0000
GTRI	20	-1887.7909	3903.8181	0.0000	1.0000
GTRIG	21	-1885.4010	3905.4502	0.0000	1.0000
SYMG	17	-1898.7862	3906.5733	0.0000	1.0000
SYMI	17	-1901.4098	3911.8205	0.0000	1.0000
SYMIG	18	-1898.5441	3912.5010	0.0000	1.0000
F81G	15	-1910.3165	3916.8103	0.0000	1.0000
F81I	15	-1912.3818	3920.9409	0.0000	1.0000
HKY	15	-1912.6641	3921.5054	0.0000	1.0000
F81IG	16	-1910.8110	3924.2111	0.0000	1.0000
JC69G	12	-1925.6743	3928.2904	0.0000	1.0000
K2P	12	-1925.7363	3928.4145	0.0000	1.0000
JC69I	12	-1927.3468	3931.6354	0.0000	1.0000
JC69IG	13	-1925.5575	3934.4686	0.0000	1.0000
GTR	19	-1910.3245	3942.4736	0.0000	1.0000
SYM	16	-1923.2204	3949.0299	0.0000	1.0000
F81	14	-1933.3142	3956.3939	0.0000	1.0000
JC69	11	-1947.4782	3965.4865	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 591-tomodel-phylip/chr3_23724.phylip = 1 - 609;
 Partition Dummy = 1:591-tomodel-phylip/chr3_23724.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 591-tomodel-phylip/chr3_23724.phylip = 1 - 609;
 Partition Dummy = 1:591-tomodel-phylip/chr3_23724.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 591-tomodel-phylip/chr3_23724.phylip = 1 - 609;
 Partition Dummy = 1:591-tomodel-phylip/chr3_23724.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

