model1 <- multdrc(hormesis[,c(2,1)], fct=bcl3())
summary(model1)
ED(model1, c(50)) # compare the parameter estimate and
# its estimated standard error
# to the values in the above paper
rm(model1)
Run the code above in your browser using DataLab