R/model-hetnormal-gam.R

In IQSS/ZeligGAM: Zelig for Generalized Additive Models

zhetnormalgam <- setRefClass("Zelig-hetnormal-gam",
                             contains = "Zelig")


zhetnormalgam$methods(
  initialize = function() {
    callSuper()
    .self$name <- "hetnormal.gam"
    # .self$fn <- zlavaan
    # .self$authors <- "Christine Choirat"
    # .self$description <- "Structural Equation Model"
    .self$year <- 2015
    .self$category <- "continuous"
    .self$wrapper <- "hetnormal.gam"
  }
)

zhetnormalgam$methods(
  zelig = function(formula, data, ..., weights = NULL, by = NULL) {
    .self$zelig.call <- match.call(expand.dots = TRUE)
    .self$model.call <- .self$zelig.call
    callSuper(formula = formula, data = data, ...,
              weights = NULL, by = by)
  }
)

zhetnormalgam$methods(
  set = function(z.out) {
    return()
  }
)

##----- QI's need to be defined

zhetnormalgam$methods(
  param = function(z.out) {
    return(mvrnorm(.self$num, coef(z.out), vcov(z.out)))
  }
)

zhetnormalgam$methods(
  qi = function(simparam, mm) {
    ev <- simparam %*% t(mm)
    pv <- ev
    return(list(ev = ev, pv = pv))
  }
)