Skip to contents

Run one EM fit from a single initial responsibility matrix.

Source: R/em_driver.R
mixqr_em.Rd

Run one EM fit from a single initial responsibility matrix.

Usage

mixqr_em(
  engine,
  X,
  y,
  tau,
  m,
  p_init,
  error_density = "unequal",
  weights = NULL,
  control = mixqr_control()
)

Arguments

engine

A mixqr_engine (list of closures, frozen contract).

X, y

Design matrix and response.

tau

Quantile level.

m

Number of components.

p_init

Initial n x m responsibility matrix.

error_density

"unequal" or "equal" (recorded in state).

weights

Observation weights.

control

A mixqr_control() list.

Value

A list with state, posterior, iter, converged, objective, objective_path.