help.search("Data")
.Data
q()
.
help.search("enzyme")
?Puromycin
?SSmicmen
Puromycin
data(Puromycin)
Puromycin
ls()
help("puromycin")
help("Puromycin")
     plot(rate ~ conc, data = Puromycin, las = 1,
          xlab = "Substrate concentration (ppm)",
          ylab = "Reaction velocity (counts/min/min)",
          pch = as.integer(Puromycin$state),
          col = as.integer(Puromycin$state),
          main = "Puromycin data and fitted Michaelis-Menten curves")
     fm4 <- nls(rate ~ SSmicmen(conc, Vm, K), data = Puromycin,
                subset = state == "treated")
     summary(fm4)
     ## add fitted lines to the plot
     conc <- seq(0, 1.2, len = 101)
     lines(conc, predict(fm1, list(conc = conc)), lty = 1, col = 1)
     lines(conc, predict(fm2, list(conc = conc)), lty = 2, col = 2)
     legend(0.8, 120, levels(Puromycin$state),
            col = 1:2, lty = 1:2, pch = 1:2)
            col = 1:2, lty = 1:2, pch = 1:2)
graphics.off()
     fm4 <- nls(rate ~ SSmicmen(conc, Vm, K), data = Puromycin,
                subset = state == "treated")
     summary(fm4)
     ## add fitted lines to the plot
     conc <- seq(0, 1.2, len = 101)
     lines(conc, predict(fm1, list(conc = conc)), lty = 1, col = 1)
     lines(conc, predict(fm2, list(conc = conc)), lty = 2, col = 2)
     legend(0.8, 120, levels(Puromycin$state),
            col = 1:2, lty = 1:2, pch = 1:2)
graphics.off()
     plot(rate ~ conc, data = Puromycin, las = 1,
          xlab = "Substrate concentration (ppm)",
          ylab = "Reaction velocity (counts/min/min)",
          pch = as.integer(Puromycin$state),
          col = as.integer(Puromycin$state),
          main = "Puromycin data and fitted Michaelis-Menten curves")
     ## simplest form of fitting the Michaelis-Menten model to these data
     fm1 <- nls(rate ~ Vm * conc/(K + conc), data = Puromycin,
                subset = state == "treated",
                start = c(Vm = 200, K = 0.05), trace = TRUE)
     fm2 <- nls(rate ~ Vm * conc/(K + conc), data = Puromycin,
                subset = state == "untreated",
                start = c(Vm = 160, K = 0.05), trace = TRUE)
     ## add fitted lines to the plot
     conc <- seq(0, 1.2, len = 101)
     lines(conc, predict(fm1, list(conc = conc)), lty = 1, col = 1)
     lines(conc, predict(fm2, list(conc = conc)), lty = 2, col = 2)
     legend(0.8, 120, levels(Puromycin$state),
            col = 1:2, lty = 1:2, pch = 1:2)
q()
