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