R/cov_clm__getFittingFunction__fittingGLSC2.R

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

fittingGLSC2 <- function(sigmaThetaExprList, fit) {
  function(params) {
    sigmaTheta <- evalSigmaThetaExprList(sigmaThetaExprList, as.list(params))

    sum(diag((diag(1, nrow = fit$info$T, ncol = fit$info$T) - (sigmaTheta %*% solve(fit$sigmaHat))) %*%
      t(diag(1, nrow = fit$info$T, ncol = fit$info$T) - (sigmaTheta %*% solve(fit$sigmaHat))))) / 2
  }
}