hmmTMB

This R package implements flexible hidden Markov models, based on Template Model Builder (TMB): flexible state-dependent distributions, transition probability structures, random effects, and smoothing splines.

Preprint

The statistical background, as well as details about the implementation of the package, and several example analyses, are presented in the following preprint.

Michelot, T. (2022). hmmTMB: Hidden Markov models with flexible covariate effects in R. arXiv:2211.14139.

Package installation

The package is available on CRAN, and the stable version can therefore be installed using

install.packages("hmmTMB")

The development version of the package can be installed from Github using devtools,

devtools::install_github("TheoMichelot/hmmTMB")

Package documentation

To find help files for the methods implemented in the package, search for help using the name of the corresponding class, e.g.,

?MarkovChain
?Observation
?HMM

We describe functionalities of the package in several vignettes, which are listed below.

Main vignette

• 'Analysing time series data with hidden Markov models in hmmTMB': Overview of package workflow, using detailed example based on analysis of energy prices. This is a good starting point to learn how to use the package.

Other features

• 'Bayesian inference in hmmTMB': Description of workflow for Bayesian analysis in hmmTMB, including specifying priors, and extracting posterior samples.

• 'Advanced features of hmmTMB': Description of some other useful functionalities, including (semi-)supervised learning, parameter constraints, selection of initial parameter values, etc.

• 'General dependence structures in hmmTMB': Implementation details for hidden Markov models (HMMs) with non-standard dependence structures, including hidden semi-Markov models, higher-order HMMs, autoregressive HMMs, and coupled HMMs.

Case studies

• 'Flexible animal movement modelling using hmmTMB': Description of wild haggis movement analysis, illustrating how non-parametric covariate effects can be included. This includes two different types of movement models: (1) correlated random walks based on step lengths and turning angles, and (2) correlated random walks based on locations directly.

• 'Occupancy modelling using hmmTMB': Analysis of occupancy data set of crossbill from Kéry et al. (2013).

• 'The categorical distribution in hmmTMB': The categorical distribution is slightly different from others, in that the number of parameters depends on the data (specifically, on the number of categories in the data). This vignette shows an example analysis for a data set on patient-therapist interactions from the mHMMbayes package.

  • 'The multivariate normal distribution in hmmTMB': We illustrate the state-dependent multivariate normal distribution for modelling correlated variables, using time series of log-returns. The vignette describes the formulation in terms of standard deviations and correlation parameters, and compares this approach to using multiple univariate normal distributions.

Reference vignettes

• 'List of distributions in hmmTMB': List of observation distributions currently available in hmmTMB.

• 'Adding new distributions in hmmTMB': Technical vignette about adding new observation distributions to hmmTMB, with detailed instructions to add the required R and C++ code.

• 'Syntax of model formulas in hmmTMB': Flexible covariate dependence can be specified using the syntax of mgcv formulas, and this vignette describes a few special cases, including random effects (intercepts or slopes), smooth non-linear relationships, factor-smooth interactions, multidimensional smooths, varying-coefficient models, cyclical effects, etc.

• 'HMM package comparison': Comparison of syntax between several R packages (hmmTMB, momentuHMM, depmixS4, LMest) and a Python package (hmmlearn) for a simple 2-state Gaussian HMM.