# 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