## Fitting a hormesis model
fit.ml3c <- multdrc(DW~Dose, data=leaflength, fct=ml3c(), boxcox=TRUE)
## Plottinf fitted curve and original data
plot(fit.ml3c, conLevel=0.001, type="points", legend=FALSE,
ylab="Produced leaf length(cm)", xlab="Metsulfuron-methyl(mg/l)",
main="Hormesis: leaf length of barley")
rm(fit.ml3c)
Run the code above in your browser using DataLab