R/clm__getSigmaMatrix__getSigmaExchangeable__rhoCalculation.R

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

rhoCalculation <- function(fit) {
  rho <- ((sum(apply(fit$residuals, 1, sum)^2) - sum(fit$residuals^2)) /
    (fit$info$N * fit$info$T * (fit$info$T - 1))) /
    (sum(fit$residuals^2) / (fit$info$N * fit$info$T))

  return(rho)
}