Marce Project Week 4 DDP

Marce

1/20/2020

Introduction

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.

The code behind the App… UI



```{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”)
) ) ))`

```

The code behind the App… Server



```{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.993
1.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”)})})`

```

Some Instructions

Slide with Plot

Visit my app to check your results ETU-EDI Calculator