R/1.4.4.main.R

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

getSigmaExchangeable <- function(modelComponents) {
  fit <- fitModelWithSigmaIdentity(modelComponents)
  rho <- rhoCalculation(fit)

  sigma <- matrix(rho, ncol = fit$info$T, nrow = fit$info$T)
  diag(sigma) <- 1

  attr(sigma, "type") <- "exchangeable"
  attr(sigma, "rho") <- rho

  return(sigma)
}