R/cov_clm__fitThetaModel__fitThetaCovariance__fitThetaCovarPML__evalDeriv2SigmaThetaExprList.R

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

evalDeriv2SigmaThetaExprList <- function(derivativesExprList, sigmaThetaExprList, paramsList) {
  names(paramsList) <- sigmaThetaExprList$paramsNames

  derivativesList <- lapply(derivativesExprList, function(expr) eval(expr, paramsList))
  deriv2SigmaThetaCombined <- sapply(derivativesList, function(expr) attr(expr, "hessian"))
  deriv2SigmaThetaList <- split(deriv2SigmaThetaCombined, seq(nrow(deriv2SigmaThetaCombined)))

  deriv2SigmaThetaMatrixList <- lapply(deriv2SigmaThetaList, function(expr) {
    matrix(expr, ncol = sqrt(length(derivativesList)), nrow = sqrt(length(derivativesList)))
  })

  return(deriv2SigmaThetaMatrixList)
}