在nlme中拟合混合效应非线性模型时的两个错误消息

Question

这是我的代码：

> pred4<-function(a1,b1,a3,b3,Pa,Temp,D) log(Pa)+a1+b1*(1000/(Temp+273.15) - 0.81)-log(1+Pa*exp(a3+b3*(1000/(Temp+273.15) - 0.81)))+D
> nlme.fit <- nlme(ln.R. ~ pred4(a1,b1,a3,b3,Pa,Temp,Day), fixed = a1+b1+a3+b3 ~1,random = D~1|Day, data =Dat4param, start=start0)

如果我删除随机日效果D，

nlme非常适合非线性模型。但是，上面的代码会导致以下错误和重复警告：

Error in chol.default((value + t(value))/2) : 
  the leading minor of order 1 is not positive definite
In addition: Warning messages:
1: In Ops.factor(log(Pa) + a1 + b1 * (1000/(Temp + 273.15) - 0.81) -  :
  + not meaningful for factors
2: In Ops.factor(log(Pa) + a1 + b1 * (1000/(Temp + 273.15) - 0.81) -  :
  + not meaningful for factors

等

有什么建议吗？