如何在nlmer中为四参数逻辑模型添加固定效应

时间:2012-06-15 18:56:52

标签: r lme4

我正在尝试使用nlmerSSfpl来使用四参数逻辑函数来拟合某些数据。我可以使用以下方法很好地适应整体数据:

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),然后比较参数,但这看起来不是正确的方法。

1 个答案:

答案 0 :(得分: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))