R/2.8.generatePreOrder1SigmaThetaExprList.R

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

generatePreOrder1SigmaThetaExprList <- function(fit) {
  wavesCombination <- expand.grid(t = seq(fit$info$T), tlinha = seq(fit$info$T))

  sigmaThetaExpr <- mapply(function(t, tlinha) {
    if (t == tlinha) {
      parse(text = paste0("sigma2u + sigmav", t, "^2"))
    } else if (tlinha < t) {
      parse(text = paste0("sigma2u + sigmav", tlinha, " * sigmav", t, Reduce("paste0", paste0(" * rho", seq(t - 1, from = tlinha)))))
    } else {
      parse(text = paste0("sigma2u + sigmav", t, " * sigmav", tlinha, Reduce("paste0", paste0(" * rho", seq(tlinha - 1, from = t)))))
    }
  }, wavesCombination$t, wavesCombination$tlinha, SIMPLIFY = FALSE)

  sigmaThetaExprList <- list(
    expr = sigmaThetaExpr,
    paramsNames = c("sigma2u", paste0("sigmav", seq(fit$info$T)), paste0("rho", seq(fit$info$T - 1)))
  )

  return(sigmaThetaExprList)
}