6.1 Analysis of Mortality and Immobility:
• Calculate percentage mortality and immobility for each concentration and time point, separately for males and females.
• Use the drc package in R (with the drm function) to determine LCx and ECx values for each sex.
Example R code: 1. model1 <- drm(Immobile/Exposed ~ Conc, Sex, weights = Exposed, data = Anthracene,
type = "binomial", fct = LL.2())
ED(model1, c(50), interval = "delta") # for EC50 and confidence intervals
2. model2 <- drm(Dead/Exposed ~ Conc, Sex, weights = Exposed, data = Anthracene,
type = "binomial", fct = LL.2())
ED(model1, c(50), interval = "delta") # for determining the LC50 and 95% confidence intervals
6.2 Dose-Response Curves:
• Generate dose-response curves for each PAH, separately for males and females.
Example R code: plot(model1, broken = TRUE, ylim = c(0, 1), xlim = c(0,5000))
6.3 Statistical Analysis:
• Use generalised linear models (GLM) to investigate interactions between more variables.
Example R code: model3 <- glm(Dead/Exposed ~ Concentration + PAH + Sex + Concentration*PAH + Concentration*Sex, family = binomial,
data = PAH_Data, weights = Exposed)