R/clm__fitModel__fitSigmaHat.R

In phmpacheco-ufjf/clm: A package to deal with complex longitudinal survey data

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)
}