conc <- c(0, rev(unlist(lapply(1:18, function(x) 10*(2/3)^x))),10)
response <- c(1/(1+exp(-0.7*(4-conc[-20])))+rnorm(19)/50, 0)
TestRun <- growth.drFitModel(conc, response, drID = 'test')
print(summary(TestRun))
Run the code above in your browser using DataLab