Difference between revisions of "ERFs of vitamins"

#This is code Op_en6866/ on page [[ERFs of vitamins]]

library(OpasnetUtils)
objects.latest("Op_en6866", code_name="initiate")
ERF_vit <- EvalOutput(ERF_vit)
threshold_vit <- EvalOutput(threshold_vit)

cat("ERF_vit\n")
oprint(summary(ERF_vit))
cat("threshold_vit\n")
oprint(summary(threshold_vit))

# This is code Op_en6866/initiate on page [[ERFs of vitamins]]

library(OpasnetUtils)

d <- opbase.data("Op_en6866")[-1] #Remove Obs
colnames(d) <- gsub(" ", "_", colnames(d))
d$Exposure_agent <- gsub(" ", "_", d$Exposure_agent)
d$Result <- ifelse(d$Result == "", "0", as.character(d$Result))

ERF_vit <- Ovariable("ERF_vit", data = d[d$Observation == "ERF", colnames(d) != "Observation"])

threshold_vit <- Ovariable("threshold_vit", data = d[d$Observation == "Threshold", colnames(d) != "Observation"])

objects.store(ERF_vit, threshold_vit)
cat("Ovariables ERF_vit, threshold_vit stored.\n")