R/cov_clm__getGradientFunction__fittingGradientPML__evalDerivSigmaThetaExprList.R

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

evalDerivSigmaThetaExprList <- function(derivativesExprList, sigmaThetaExprList, paramsList) {
  if (is.null(names(paramsList))) {
    names(paramsList) <- sigmaThetaExprList$paramsNames
  }

  derivativesList <- lapply(derivativesExprList, function(expr) eval(expr, paramsList))
  derivSigmaThetaCombined <- sapply(derivativesList, function(expr) attr(expr, "gradient"))
  derivSigmaThetaList <- stats::setNames(
    split(derivSigmaThetaCombined, seq(nrow(derivSigmaThetaCombined))),
    sigmaThetaExprList$paramsNames
  )
  derivSigmaThetaMatrixList <- lapply(derivSigmaThetaList, function(expr) {
    matrix(expr, ncol = sqrt(length(derivativesList)), nrow = sqrt(length(derivativesList)))
  })

  return(derivSigmaThetaMatrixList)
}