我在R中使用以下命令。
attach(Puromycin)
Puromycin
plot(Puromycin$conc,Puromycin$rate)
mm=function(conc,vmax,k) vmax*conc/(k+conc)
mm
mm1=nls(rate~mm(conc,vmax,k),data=Puromycin,
start=c(vmax=50,k=0.05),subset=state=="treated”)
mm1
不幸的是,我收到了下面的错误。
Error: unexpected input in "mm1=nls(rate~mm(conc,vmax,k),data=Puromycin,start=c(vmax=50,k=0.05),subset=state==”"
我怎样才能解决这个错误?
这是一个完整的解决方案,并绘制了拟合线。
data(Puromycin)
mm <- function(conc, vmax, k) vmax*conc/(k + conc)
mm1 <- nls(rate ~ mm(conc,vmax,k), data = Puromycin,
start = list(vmax = 50, k = 0.05),
subset = state == "treated")
summary(mm1)
with(Puromycin, plot(conc, rate))
x <- with(Puromycin, seq(min(conc), max(conc), length.out = 100))
y <- predict(mm1, newdata = data.frame(conc = x), se = FALSE)
lines(x, y, col = "blue")