Hi: I can't give you a direct answer, but I'd suggest trying
library(sos) findFn('hidden Markov model') which should give you a pretty good idea of where to start looking. If you don't have the package, install it first from CRAN. HTH, Dennis On Tue, Mar 22, 2011 at 4:07 PM, Emmanuel Levy <emmanuel.l...@gmail.com>wrote: > Dear All, > > I would like to generate random protein sequences using a HMM model. > Has anybody done that before, or would you have any idea which package > is likely to be best for that? > > The important facts are that the HMM will be fitted on ~3 million > sequential observations, with 20 different states (one for each amino > acid). I guess that 2-5 hidden states should be enough, and an order > of 3 would also be enough I would say. > > Thanks a lot for any suggestion! > > All the best, > > Emmanuel > > ______________________________________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.