如何在nlmer的四参数logistic模型中加入固定效应

我正在尝试使用

nlmer

和

SSfpl

来拟合一些具有四参数逻辑函数的数据。我可以使用以下方法对总体数据进行精确拟合：

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

---

nm.fpl您是否尝试添加交互

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

nm.fpl也会产生错误：nlmer中的错误（meanFix~SSfpl（时间，A，B，xmid*条件，scal）~：评估模型的梯度属性必须是数字矩阵。此外：警告消息：在运算因子中（xmid，条件）：*对Factors没有意义我遇到了同样的问题，到目前为止没有成功。我的尝试在上面的更新：见（中的修订材料）链接中进行了描述。我已经取得了一些成功，目前这个过程太复杂了，我无法在这里完整地复制它（但欢迎您将我的材料带到这里，并在这里归结为一个答案。）底线是，您必须手动构建自己的分组变量/评估（简单地指定@Thierry答案中的交互公式是行不通的），但这是可以做到的。您可能会发现我的答案很有用。