## 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 <- drm(count~conc, data=S.capricornutum, fct=LN.3(), adjust="bc1")
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
s.cap.m3<-drm(count~conc, data=S.capricornutum, fct=LN.3(loge=FALSE), adjust="bc1")
ED(s.cap.m3, c(10, 20, 50), ci="delta")
Run the code above in your browser using DataLab