Difference between revisions of "Kausal Paths"

library(OpasnetUtils)

openv.setN(0)

objects.latest("Op_en2031", code_name="ERF2")
ERFchoice <- 1

ERF@formula <- function (...) 
{
  out <- OpasnetUtils::combine(ERF_env, ERF_omega3, ERF_mehg, 
                               ERF_diox, ERF_vit, ERF_micr, ERF_pfas)
#  out$Exposure_unit <- NULL
#  out$Exposure <- NULL
  if ("Iter" %in% colnames(out@output)) {
    out@output <- fillna(out@output, marginals = "Iter")
  }
  out <- unkeep(out * ERFchoice, cols = c("ERFchoiceSource", 
                                          "ERFchoiceResult"))
  marginals <- colnames(out@output)[out@marginal]
  coln <- colnames(out@output)
  out@output <- subgrouping(out@output, subg = "Subgroup")
  out@marginal <- colnames(out@output) %in% union(marginals, 
                                                  setdiff(colnames(out@output), coln))
  return(out)
}

ERF <- EvalOutput(ERF)

write.csv(ERF@output, "exposure_response_functions.csv", row.names = FALSE)