[Bioperl-l] PAML::Codeml outputs unstable value, why?
Roy Chaudhuri
roy at colibase.bham.ac.uk
Tue May 29 14:05:12 UTC 2007
Hi Xianjun,
I'm not sure if it is the cause of your problem, but your sequences seem
to be quite short. This paper:
http://mbe.oxfordjournals.org/cgi/content/full/21/12/2290
suggests that the codeml method of calculating Ka and Ks may be
unreliable for sequences shorter than 300 codons.
Roy.
--
Dr. Roy Chaudhuri
Department of Veterinary Medicine
University of Cambridge, U.K.
More information about the Bioperl-l
mailing list