Как исправить ошибку нелинейной регрессии: функция nls - PullRequest
1 голос
/ 27 мая 2020

Я использую следующие команды в 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==”"

Как я могу исправить эту ошибку?

1 Ответ

3 голосов
/ 27 мая 2020

Вот полное решение с нанесенной линией.

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")

enter image description here

...