## Fitting a two-parameter logistic model
## to binomial responses (a logit model)
earthworms.m1 <- drm(number/total~dose, weights=total,
data = earthworms, fct = LL.2(), type = "binomial")
plot(earthworms.m1) # not fitting at the upper limit!
Run the code above in your browser using DataLab