5

我正在尝试使用nlmerwithSSfpl来拟合一些具有四参数逻辑函数的数据。我可以使用以下方法很好地拟合整体数据:

nm.fpl <- nlmer(meanFix ~ SSfpl(Time, A, B, xmid, scal) ~ (scal | Subject), 
          data = dataSubset, start = c(A = 0.2, B = 0.7, xmid = 600, scal = 100))

现在我想添加一个固定的效果Condition,它有 2 个主题内级别。我想评估这两个条件是否在 4 个参数(A、B、xmid、scal)中的任何一个方面有所不同,但我不知道如何在这个公式中指定它。我可以将模型分别拟合到两个子集(条件 A 和条件 B),然后比较参数,但这似乎不是正确的方法。

4

1 回答 1

0

您是否尝试过添加交互?

nm.fpl <- nlmer(meanFix ~ SSfpl(Time, A * Condition, B, xmid, scal) ~ (scal | Subject), 
          data = dataSubset, start = c(A = 0.2, B = 0.7, xmid = 600, scal = 100))
nm.fpl <- nlmer(meanFix ~ SSfpl(Time, A : Condition, B, xmid, scal) ~ (scal | Subject), 
          data = dataSubset, start = c(A = 0.2, B = 0.7, xmid = 600, scal = 100))
于 2012-06-15T22:33:48.177 回答