R/2.6.2.getGFIgofULSC.R

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

getGFIgofULSC <- function(fit, fitTheta) {
  nominator <- fit$sigmaHat - fitTheta$sigmaTheta

  GFI <- 1 - sum(diag(nominator %*% nominator)) / sum(diag(fit$sigmaHat %*% fit$sigmaHat))

  return(GFI)
}