head(O.mykiss)
## Fitting exponential model
O.mykiss.m1<-drm(weight~conc, data=O.mykiss, fct=EXD.2(), na.action=na.omit)
summary(O.mykiss.m1)
anova(O.mykiss.m1)
## Fitting same model with transform-both-sides approach
O.mykiss.m2<-drm(weight~conc, data=O.mykiss, fct=EXD.2(), na.action=na.omit, adjust="bc1")
summary(O.mykiss.m2)
# no need for a transformation
## Plotting the fit
plot(O.mykiss.m1, type = "points", xlim = c(0, 500), ylim = c(0,4),
xlab = "Concentration (mg/l)", ylab = "Weight (g)", broken = TRUE)
Run the code above in your browser using DataLab