# Plot everything for the Puromycin data
library(nls)
data(Puromycin)
Purboth <- nls(formula = rate ~ ((Vm + delV * (state == "treated"))* conc)/(K + conc),
data = Puromycin, start = list(Vm = 160, delV = 40, K = 0.05))
Pur.prof <- profile(Purboth)
pairs(Pur.prof, plot.ellipse = TRUE)
Run the code above in your browser using DataLab