I know how to fit a hidden markov model to a data sequence, using the matlab-implementation of the baum-welch algorithm.
But what should I do if I do not have one data sequence, but a bunch of them? How can I fit to a set of sequences?
I think what I want to do is to maximize the likelihood of the HMM to output any of my sequences.
Is there a best practice for that kind of problem?
R
, to my understanding, thedepmixS4
package may be able to handle this. See here for the paper: jstatsoft.org/v36/i07 $\endgroup$