model1 <- multdrc(hormesis[,c(2,1)], fct=BC.5())
anova(model1)
plot(model1)
model2 <- multdrc(hormesis[,c(2,1)], fct=BC.4())
summary(model2)
ED(model2, c(50)) # compare the parameter estimate and
# its estimated standard error
# to the values in the paper by
# van Ewijk and Hoekstra
rm(model1, model2)
Run the code above in your browser using DataLab