R/1.2.1.fitSigmaHat.R

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

fitSigmaHat <- function(fit) {
  Sw <- fit$individuals %>%
    rlist::list.map((sqrt(wi) * (yi - mui)) %*% t(sqrt(wi) * (yi - mui))) %>%
    Reduce("+", .)
  Sw <- (1 / rlist::list.map(fit$individuals, wi) %>% Reduce("+", .)) * Sw

  return(Sw)
}