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