cov.mix.wt: weighted covariance

In hhsmm: Hidden Hybrid Markov/Semi-Markov Model Fitting

weighted covariance

Description

The weighted means and variances using the observation matrix and the estimated weight vectors

Usage

cov.mix.wt(
  x,
  wt1 = rep(1/nrow(x), nrow(x)),
  wt2 = rep(1/nrow(x), nrow(x)),
  cor = FALSE,
  center = TRUE,
  method = c("unbiased", "ML")
)

Arguments

x	the observation matrix
wt1	the state probabilities matrix (number of observations times number of states)
wt2	the mixture components probabilities list (of length nstate) of matrices (number of observations times number of mixture components)
cor	logical. if TRUE the weighted correlation is also given
center	logical. if TRUE the weighted mean is also given
method	with two possible entries: "unbiased" the unbiased estimator is given "ML" the maximum likelihood estimator is given

Value

list containing the following items:

• center the weighted mean of x

• cov the weighted covariance of x

• n.obs the number of observations in x

• cor the weighted correlation of x, if the parameter cor is TRUE

• wt1 the state weighs wt1

• wt2 the mixture component weights wt2

• pmix the estimated mixture proportions

Author(s)

Morteza Amini, [e-mail], Afarin Bayat, [e-mail]

Examples

data(CMAPSS)
n = nrow(CMAPSS$train$x)
wt1 = runif(n)
wt2 = runif(n)
cov.mix.wt(CMAPSS$train$x, wt1, wt2)

hhsmm documentation built on Aug. 8, 2023, 9:06 a.m.