rpois.hsmm {mhsmm} | R Documentation |
This generates values from a Poisson distributed emission state j
given parameters in model
.
rpois.hsmm(j, model)
j |
An integer representing the state |
model |
A |
This is essentially a wrapper for rpois
. Users may build functions with the same
arguments and return values so they can use their own custom emission distributions.
A single value from the emission distribution.
Jared O'Connell jaredoconnell@gmail.com
J<-3 initial <- rep(1/J,J) P <- matrix(c(.8,.5,.1,0.05,.2,.5,.15,.3,.4),nrow=J) b <- list(lambda=c(1,3,6)) model <- hmmspec(init=initial, trans=P, parms.emission=b,dens.emission=dpois.hsmm) model train <- simulate(model, nsim=300, seed=1234, rand.emis=rpois.hsmm) plot(train,xlim=c(0,100)) h1 = hmmfit(train,model,mstep=mstep.pois)