R/2.5.fittingULCS.R

In Mmcsd: Modeling Complex Longitudinal Data in a Quick and Easy Way

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

    sum(diag((fit$sigmaHat - sigmaTheta) %*% t(fit$sigmaHat - sigmaTheta))) / 2
  }
}