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