R/cov_clm__getFittingFunction__fittingGLSC3.R

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

fittingGLSC3 <- function(sigmaThetaExpr, fit) {
  function(params) {
    sigmaTheta <- evalSigmaThetaExprList(sigmaThetaExpr, as.list(params))

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