## Fitting 3-parameter log-normal model
s.cap.m1 <- drm(count ~ conc, data = S.capricornutum, fct = LN.3())
## Residual plot
plot(fitted(s.cap.m1), residuals(s.cap.m1))
## Fitting model with transform-both-sides approach
s.cap.m2 <- boxcox(s.cap.m1, method = "anova")
summary(s.cap.m2)
## Residual plot after transformation (looks better)
plot(fitted(s.cap.m2), residuals(s.cap.m2))
## Calculating ED values on log scale (not so interesting)
ED(s.cap.m2, c(10, 20, 50), ci="delta")
## Calculating ED values on the original
ED(s.cap.m2, c(10, 20, 50), ci="fls")
## Fitting model with ED50 as parameter
## (for comparison)
s.cap.m3 <- drm(count ~ conc, data = S.capricornutum, fct = LN.3(loge=FALSE))
ED(s.cap.m3, c(10, 20, 50), ci = "delta")
Run the code above in your browser using DataLab