mfit <- vglm(velocity ~ 1, micmen, data = enzyme, trace = TRUE,
crit = "coef", form2 = ~ conc - 1)
summary(mfit)
if (FALSE) {
plot(velocity ~ conc, enzyme, xlab = "concentration", las = 1,
col = "blue",
main = "Michaelis-Menten equation for the enzyme data",
ylim = c(0, max(velocity)), xlim = c(0, max(conc)))
points(fitted(mfit) ~ conc, enzyme, col = 2, pch = "+", cex = 2)
# This predicts the response at a finer grid:
newenzyme <- data.frame(conc = seq(0, max(with(enzyme, conc)),
len = 200))
mfit@extra$Xm2 <- newenzyme$conc # This is needed for prediction
lines(predict(mfit, newenzyme, "response") ~ conc, newenzyme,
col = "red") }
Run the code above in your browser using DataLab