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