By Brian Cervantes Alvarez, Willa Van Liew, Hans Lehndorff
The study was developed from these concerns:
The global temperature increase might cause:
Increased transmission
Decreased transmission
Pesticides that are used to deal with those parasites are not being tested in different temperatures settings.
Does temperature affect the toxicity of pesticides?
Atrazine Herbicide
Azoxystrobin Fungicide
Carbaryl Insecticide
Methodology
Prepared a sample of 1,076 cercariae tadpoles.
Use fixed temperature of [TEMP]
Tested [NAME] pesticide
Apply various concentrations to tadpole groupings.
Repeat previous steps with next pesticide [NAME].
Record the Status and TTKD of each tadpole in group.
Repeat steps 1 & 2 with next [TEMP]
Key data
Id: unique identifier for each parasite
Temp: temperature of the water
Detected_conc: concentration of pesticide
TTKD : Time to ‘knockdown’
Status: Was tadpole knocked down?
Pesticide: Pesticide in test
at3 = ds %>%
filter(pesticide == "atrazine")
at3sf = survfit(Surv(at3$TTKD, at3$status)~at3$temp)
atrazine = ggsurvplot(fit = at3sf,
data = at3,
legend = "bottom",
palette = "Dark2",
risk.table = F,
legend.labs = c("18°C","24°C", "30°C"),
conf.int=F,
surv.median.line="hv") +
labs(title = "Survival Rate with Atrazine",
x = "Time in Hours")
atrazine
survdiff(Surv(at3$TTKD, at3$status)~at3$temp)
at1 = ds %>%
filter(pesticide == "azoxystrobin")
at1sf = survfit(Surv(at1$TTKD, at1$status)~at1$temp)
azoxystrobin = ggsurvplot(fit = at1sf,
data = at1,
palette = "Dark2",
legend = "bottom",
legend.labs = c("18°C","24°C", "30°C"),
risk.table = F,
conf.int=F,
surv.median.line="hv") +
labs(title = "Survival Rate with Pesticide Azoxystrobin",
x = "Time in Hours")
azoxystrobin
survdiff(Surv(at1$TTKD, at1$status)~at1$temp)
at2 = ds %>%
filter(pesticide == "carbaryl")
at2sf = survfit(Surv(at2$TTKD, at2$status)~at2$temp)
carbaryl = ggsurvplot(fit = at2sf,
data = at2,
palette = "Dark2",
legend = "bottom",
legend.labs = c("18°C","24°C", "30°C"),
risk.table = F,
conf.int=F,
surv.median.line="hv") +
labs(title = "Survival Rate with Carbaryl",
x = "Time in Hours")
carbaryl
survdiff(Surv(at2$TTKD, at2$status)~at2$temp)
Main findings were:
Original Study: https://purr.purdue.edu/publications/3933/1