thanks very much for your model with the fixed-effects! Do you have any idea how to combine the splines with random effects instead of the fixed ones? You would be a great help.

Thanks in advance,

Joline

Y_hat = a0*to_vector(X) + to_vector(a*B);

and in bs( … , intercept = TRUE)

]]>http://jofrhwld.github.io/blog/2014/07/29/nonlinear_stan.html?

]]>There’s also a huge amount of work on spline linear/regression models, which helps. Fitting similar models using GPs would probably need A) custom programs, and B) lots of computer time.

]]>