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

Input data from file "591-tomodel-phylip/chrZ_6362.phylip" (7 taxa, 621 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	-1896.7818	3826.4643	0.4267	0.4267
HKYG	16	-1897.2643	3827.4293	0.2634	0.6901
GTRI	20	-1893.6671	3828.7341	0.1372	0.8272
HKYIG	17	-1897.3431	3829.7012	0.0846	0.9118
GTRIG	21	-1893.7073	3830.9571	0.0451	0.9569
GTRG	20	-1894.8256	3831.0512	0.0431	1.0000
HKY	15	-1915.9601	3862.7136	0.0000	1.0000
GTR	19	-1913.3476	3865.9598	0.0000	1.0000
F81I	15	-1925.2361	3881.2655	0.0000	1.0000
F81G	15	-1925.7509	3882.2952	0.0000	1.0000
F81IG	16	-1925.2697	3883.4401	0.0000	1.0000
K2PI	13	-1942.0813	3910.7622	0.0000	1.0000
K2PG	13	-1942.0898	3910.7793	0.0000	1.0000
K2PIG	14	-1942.1411	3912.9752	0.0000	1.0000
SYMI	17	-1940.4698	3915.9545	0.0000	1.0000
SYMG	17	-1940.7446	3916.5041	0.0000	1.0000
F81	14	-1943.9942	3916.6815	0.0000	1.0000
SYMIG	18	-1940.7957	3918.7277	0.0000	1.0000
K2P	12	-1962.5304	3949.5740	0.0000	1.0000
JC69I	12	-1964.6077	3953.7286	0.0000	1.0000
JC69G	12	-1964.7806	3954.0744	0.0000	1.0000
SYM	16	-1961.2162	3955.3330	0.0000	1.0000
JC69IG	13	-1964.8313	3956.2624	0.0000	1.0000
JC69	11	-1984.5256	3991.4846	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1896.7818	3825.5637	0.3976	0.3976
HKYG	16	-1897.2643	3826.5286	0.2454	0.6430
GTRI	20	-1893.6671	3827.3341	0.1641	0.8071
HKYIG	17	-1897.3431	3828.6863	0.0834	0.8905
GTRIG	21	-1893.7073	3829.4145	0.0580	0.9485
GTRG	20	-1894.8256	3829.6512	0.0515	1.0000
HKY	15	-1915.9601	3861.9203	0.0000	1.0000
GTR	19	-1913.3476	3864.6953	0.0000	1.0000
F81I	15	-1925.2361	3880.4722	0.0000	1.0000
F81G	15	-1925.7509	3881.5018	0.0000	1.0000
F81IG	16	-1925.2697	3882.5394	0.0000	1.0000
K2PI	13	-1942.0813	3910.1625	0.0000	1.0000
K2PG	13	-1942.0898	3910.1796	0.0000	1.0000
K2PIG	14	-1942.1411	3912.2821	0.0000	1.0000
SYMI	17	-1940.4698	3914.9396	0.0000	1.0000
SYMG	17	-1940.7446	3915.4892	0.0000	1.0000
F81	14	-1943.9942	3915.9884	0.0000	1.0000
SYMIG	18	-1940.7957	3917.5914	0.0000	1.0000
K2P	12	-1962.5304	3949.0609	0.0000	1.0000
JC69I	12	-1964.6077	3953.2155	0.0000	1.0000
JC69G	12	-1964.7806	3953.5613	0.0000	1.0000
SYM	16	-1961.2162	3954.4324	0.0000	1.0000
JC69IG	13	-1964.8313	3955.6627	0.0000	1.0000
JC69	11	-1984.5256	3991.0511	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1896.7818	3896.4650	0.6097	0.6097
HKYG	16	-1897.2643	3897.4299	0.3763	0.9860
HKYIG	17	-1897.3431	3904.0189	0.0140	1.0000
GTRI	20	-1893.6671	3915.9607	0.0000	1.0000
GTRG	20	-1894.8256	3918.2778	0.0000	1.0000
GTRIG	21	-1893.7073	3922.4725	0.0000	1.0000
HKY	15	-1915.9601	3928.3902	0.0000	1.0000
F81I	15	-1925.2361	3946.9421	0.0000	1.0000
F81G	15	-1925.7509	3947.9718	0.0000	1.0000
GTR	19	-1913.3476	3948.8906	0.0000	1.0000
F81IG	16	-1925.2697	3953.4407	0.0000	1.0000
K2PI	13	-1942.0813	3967.7698	0.0000	1.0000
K2PG	13	-1942.0898	3967.7869	0.0000	1.0000
K2PIG	14	-1942.1411	3974.3208	0.0000	1.0000
F81	14	-1943.9942	3978.0271	0.0000	1.0000
SYMI	17	-1940.4698	3990.2722	0.0000	1.0000
SYMG	17	-1940.7446	3990.8218	0.0000	1.0000
SYMIG	18	-1940.7957	3997.3554	0.0000	1.0000
K2P	12	-1962.5304	4002.2369	0.0000	1.0000
JC69I	12	-1964.6077	4006.3915	0.0000	1.0000
JC69G	12	-1964.7806	4006.7373	0.0000	1.0000
JC69IG	13	-1964.8313	4013.2700	0.0000	1.0000
SYM	16	-1961.2162	4025.3337	0.0000	1.0000
JC69	11	-1984.5256	4039.7957	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/chrZ_6362.phylip = 1 - 621;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_6362.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/chrZ_6362.phylip = 1 - 621;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_6362.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/chrZ_6362.phylip = 1 - 621;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_6362.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

