在lme4或nlme中建模特定模型

Question

我一直试图在nlme中使用特定模型：

  

y_ij = a_i-b_i * e ^ {-g_i（A_ {ij} -12）} + d_i * l_ {ij}（+ e_ {ij}），   其中A_ {ij}和l_ {ij}来自数据，而a，b，d和g是混合效应。

我已经成功地在saemix包中实现了该模型，但是没有找到任何方法可以在lme4或nlme中提供我自己的模型。

对于seamix代码，我能够使用以下函数对数据进行建模。

Model<-function(psi,id,x){
  t <- x[,1]
  l <- x[,2]
  a <- psi[id,1]
  b <- psi[id,2]
  g <- psi[id,3]
  d <- psi[id,4]
  fpred <- a - b*exp(-g*(t-12))+d*l
  return(fpred)
}

是否可以在nlme或lme4中使用此模型？如果是这样，如何在这些软件包中对其建模？