Marce
1/20/2020
A way to asses public exposure to pesticides is to measure the concentration of the compound in urine.
One can convert the biomonitoring data into an estimated dose using reverse dosimetry.Reverse dosimetry can range from a deterministic one-compartment PK model to a complex probabilistic PBPK model.
In this case, the formula behind the App is a deterministic model that uses bmi, age and pregnancy to adjust the estimation the daily intake correspondent to the concentration found in urine.
```{r}
shinyUI(fluidPage(
titlePanel(“App to estimate the Daily Intake of Ethylenethiourea in Women”),
sidebarLayout(
sidebarPanel(
h3(“Introduce your data”),
numericInput(“Agenum”,“Age”,value=25,min=20, max=50, step=1),
numericInput(“ETUnum”,“ETU in urine (ppm)”,value=1,min=0, max=5, step=0.01),
sliderInput(“Heightslider”,“Height (m)”,1.40,1.90,1.40),
sliderInput(“Weightslider”,“Weight (kg)”,30,150,30),
selectInput(“Pregnant”,“Are you Pregnant?”, choices = c(“Yes”,“No”)),
submitButton(“Apply”)),
mainPanel(
h4(“Your Result”),
plotOutput(“plot”)
) ) ))`
```
```{r}
shinyServer(
function(input,output){
etu_er <- 0.76 # Average humans
etu_EPA <- 0.2 #EPA
EDI.calculation <- reactive({
mage <- input\(Agenum mweight <- input\)Weightslider
mheight <- input\(Heightslider mbmi <- mweight/(mheight^2) ETU <- input\)ETUnum
CER =ifelse(mage < 18 & mbmi< 30, 1.008mheight2.045exp(0.01552(mheight-90)),
ifelse(mage < 18 & mbmi> 30, 1.008mheight2.045exp(0.01552(mheight-90))* (1.429- 0.0198mbmi),
ifelse(mage > 18 & mbmi< 30, (0.9931.64(140-mage) ((mweight1.5)* (mheight0.5)))/1000,
ifelse(mage > 18 & mbmi> 30, ((0.9931.64(140-mage)* ((mweight1.5)* (mheight0.5)))/1000)* (1.429- 0.0198*mbmi), 0))))
CER <- ifelse(input$Pregnant==“Yes”, CER+(0.3*CER),CER)
CER.gkgd= (CER/mweight)/1000
EDI = (ETU*CER.gkgd)/ etu_er
print(EDI)})
output$plot <- renderPlot({
plot(EDI.calculation() ,xlab=“ETU (ug)”, ylab= “Estimated Daily Intake”, ylim=c(0, 1), xaxt=“n”, col=“#8FDA40”, pch=16,size=4)
abline(h=etu_EPA,col=“#40D5DA”,lty=“dashed”)
text(x=0.7,y=0.25, labels = “EPA Reference Dose”, col = “#40D5DA”)
text(x=1,y=EDI.calculation()+0.05, labels = paste(“Your EDI is here”, round(EDI.calculation(),digits = 4 ) ), col = “#8FDA40”)})})`
```
Visit my app to check your results ETU-EDI Calculator