# NOT RUN {
obj.dr4pl <- dr4pl(Response ~ Dose, data = sample_data_1) # Fit a 4PL model to data
confint(obj.dr4pl) # Print conventional 95% confidence intervals
confint(obj.dr4pl, level = 0.99) # Print 99%confidence intervals
theta <- coef(obj.dr4pl)
theta[4] <- 0 # Set the lower asymptote to be zero
confint(obj.dr4pl, parm = theta) # Use our dr4pl object but different parameter estimates
# }
Run the code above in your browser using DataLab