Я пытался реализовать конкретную модель в 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.
Для кода шва я смог использовать следующую функцию для моделирования данных.
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? Если да, то как можно моделировать его в этих пакетах?