R/cov_clm__getFittingFunction__fittingPML____main.R

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

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

    sum(diag(fit$sigmaHat %*% solve(sigmaTheta))) -
      unclass(determinant(fit$sigmaHat %*% solve(sigmaTheta), logarithm = TRUE)$modulus) - fit$info$T
  }
}