stateprobs: Calculate conditional local state probabilities in HMMs

In LaMa: Fast Numerical Maximum Likelihood Estimation for Latent Markov Models

Calculate conditional local state probabilities in HMMs

Description

Computes

\Pr(\text{State}_t = j \mid X_1, ..., X_T)

for a given HMM.

Usage

stateprobs(
  delta,
  Gamma,
  allprobs,
  trackID = NULL,
  mod = NULL,
  forecast = FALSE
)

Arguments

delta	initial distribution; either a vector of length nStates, or a matrix of dimension c(nTracks, nStates) if trackID is provided
Gamma	transition probability matrix; either a matrix of dimension c(nStates, nStates), an array of dimension c(nStates, nStates, nTracks) if trackID is provided, or an array of dimension c(nStates, nStates, nObs) for time-varying transition probabilities, in which case stateprobs_g is called internally
allprobs	matrix of state-dependent probabilities or density values of dimension c(nObs, nStates)
trackID	optional vector of length nObs containing nTracks unique IDs that separate tracks
mod	optional model object containing delta, Gamma, allprobs, and optionally trackID. When using RTMB::MakeADFun or qreml with forward in the likelihood, these are reported automatically after model fitting and the object returned by RTMB::report() or qreml can be passed directly.
forecast	logical, indicating if forecast probabilities \Pr(\text{State}_t = j \mid X_1, ..., X_{t-1}) should be calculated instead.

Value

matrix of conditional state probabilities of dimension c(nObs, nStates)

See Also

Other decoding functions: stateprobs_g(), stateprobs_p(), viterbi(), viterbi_g(), viterbi_p()

Examples

Gamma = tpm(c(-1,-2))
delta = stationary(Gamma)
allprobs = matrix(runif(10), nrow = 10, ncol = 2)

probs = stateprobs(delta, Gamma, allprobs)

LaMa documentation built on July 25, 2026, 1:08 a.m.