matrixdist (version 1.1.9)

MoE,bivdph-method: MoE method for bivdph Class

Description

MoE method for bivdph Class

Usage

# S4 method for bivdph
MoE(
  x,
  formula,
  y,
  data,
  alpha_vecs = NULL,
  weight = numeric(0),
  stepsEM = 1000,
  every = 10,
  rand_init = TRUE
)

Value

An object of class sph.

Arguments

x

An object of class bivdph.

formula

A regression formula.

y

A matrix of observations.

data

A data frame of covariates.

alpha_vecs

Matrix of initial probabilities.

weight

Vector of weights.

stepsEM

Number of EM steps to be performed.

every

Number of iterations between likelihood display updates.

rand_init

Random initiation in the R-step.

Examples

x <- bivdph(dimensions = c(3, 3))
n <- 100
responses <- cbind(rpois(n, 3) + 1, rbinom(n, 5, 0.5))
covariates <- data.frame(age = sample(18:65, n, replace = TRUE) / 100, income = runif(n, 0, 0.99))
f <- responses ~ age + income
MoE(x = x, formula = f, y = responses, data = covariates, stepsEM = 20)