无随机效应结构规范的非线性混合模型

Question

我想用nlme软件包来拟合一个非线性模型，它与固定的结构规范是一致的。

model <- nlme(y ~ Asym/(1+exp((xmid-x)/scal)),
                      data = data,
                      fixed = list(Asym + xmid + scal ~ treatment))
                      #random =  Asym ~ 1|subject)

但是，我得到了以下错误：

Error in parse(text = paste("~", paste(nVal, collapse = "/"))) : 
  <text>:2:0: unexpected end of input
1: ~ 
   ^

有办法解决这个问题吗？任何建议都是受欢迎的。

Answer 1

我相信您希望gnls()函数(也来自nlme包)带有params=参数，而不是fixed=。试试这个：

model <- gnls(y ~ Asym/(1+exp((xmid-x)/scal)),
                  data = data,
                  params = list(Asym + xmid + scal ~ treatment),
                  start= ...)

FWIW，如果你真的在拟合一个逻辑(这不仅仅是你想要做的一个简单的例子)，用SSlogis()自启动函数代替你的显式公式可能更快/更健壮.