R/1.4.fitModelWithSigmaIdentity.R

In Mmcsd: Modeling Complex Longitudinal Data in a Quick and Easy Way

fitModelWithSigmaIdentity <- function(modelComponents) {
  explanatoryMatrix <- as.matrix(modelComponents$explanatory[, -1])
  fit <- mmcsd(modelComponents$response$value ~ explanatoryMatrix,
    waves = modelComponents$waves$value, ids = modelComponents$ids$value, weights = modelComponents$weights$value,
    stratum = modelComponents$stratum$value, cluster = modelComponents$cluster$value, sigma = "identity"
  )

  return(fit)
}