DNase1 <- subset(DNase, Run == 1)
DNase1$density <- sapply(DNase1$density, function(x) rnorm(1, x, 0.1 * x))
mod <- onls(density ~ Asym/(1 + exp((xmid - log(conc))/scal)),
data = DNase1, start = list(Asym = 3, xmid = 0, scal = 1))
logLik_o(mod)
## compare AIC of vertical versus orthogonal residuals.
AIC(mod)
AIC(logLik_o(mod))
Run the code above in your browser using DataLab