# Plot everything for the Puromycin data
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_profile(Pur.prof, plot.ellipse = TRUE)
# Show the corresponding plot from MASS/stats:
if (getRversion() < "4.4.0") {
loadNamespace("MASS")
} else
loadNamespace("stats")
graphics::pairs(Pur.prof)
Run the code above in your browser using DataLab