cdghmm {CDGHMM}

Hidden Markov Models for Multivariate Panel Data

Description

Estimates hidden Markov models from the CDGHMM family under various missingness schemes.

Usage

cdghmm(x,m,id,mu=NULL,sigma=NULL,gamma=NULL,delta=NULL,alpha=NULL,beta=NULL,
                        maxiter=10000,tol=1e-6,type="s",covtype="VVA")

Arguments

x

Data frame or matrix to perform variable selection on

m

Number to indicate the number of states to fit.

id

A vector of indicators to indicate observational unit.

mu

An m-length list of matrices to be used as an initial estimate for mu. If no initial mu is provided the algorithm will initialize via k-means.

sigma

An m-length list of matrices to be used as an initial estimate for sigma. If no initial sigma is provided the algorithm will initialize.

gamma

A mxm matrix to be used as an initial estimate for gamma. If no initial gamma is provided the algorithm will initialize.

delta

A vector to be used as an initial estimate for delta. If no initial delta is provided the algorithm will initialize.

alpha

A mxpxt array to be used as an initial estimate for alpha. If no initial alpha is provided the algorithm will initialize.

beta

A mxpxt array to be used as an initial estimate for beta. If no initial beta is provided the algorithm will initialize.

maxiter

A number to indicate the maximum number of iterations allowed, default is 10000.

tol

A number to indicate the tolerance value, default is 1e-6.

type

A character to indicate which type of missingness mechanism to use. The allowed values are: "mar" (missing at random), "s" (3.2.1 in cited paper), "sv" (3.2.2), "st" (3.2.3 model without beta), "svt" (3.2.4 model without beta), "st2" (3.2.3 model with beta), "svt2" (3.2.4 model with beta). The default is type="s".

covtype

A string to indicate which covariance estimate to use. The allowed values are: EEA, VVA, VEA, EVA, VVI, VEI, EVI, EEI. The default is covtype="VVA".

Value

mu

Estimated mean matrices.

sigma

Estimated covariance matrices.

gamma

Estimated gamma matrix.

delta

Estimated delta vector.

alpha

Estimated alpha missingness parameters.

beta

Estimated beta missingness parameters.

llk

Estimated log-likelihood.

AIC

The value of the Akaike information criterion.

BIC

The value of the Bayes information criterion.

ICL

The value of the integrated completed likelihood.

Avg_Silouette

The value of the average silhouette score.

probs

A matrix whose entries correspond to the probability of belonging to a state.

states

Estimated states via map(probs).

mod

The CDGHMM family member fit.

Author(s)

Mackenzie R. Neal, Alexa A. Sochaniwsky, Paul D. McNicholas

References

See citation("CDGHMM").

Examples

data("simulated_data")
id=simulated_data$V5
x <- simulated_data[,1:4]
EEI_mod=cdghmm(x,2,id=id,covtype="EEI",tol=1e-4)
table(simulated_data$V7,EEI_mod$states)