rpois.hsmm {mhsmm}

Random number generation from a Poisson distributed emission distribution

Description

This generates values from a Poisson distributed emission state j given parameters in model.

Usage

rpois.hsmm(j, model)

Arguments

j

An integer representing the state

model

A hmmspec or hsmmspec object

Details

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.

Value

A single value from the emission distribution.

Author(s)

Jared O'Connell [e-mail]

See Also

mstep.pois, dpois.hsmm

Examples

  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)