[QUESTION] `posterior_traj()` for `stan_mvmer` about rstanarm

[sambrilleman]

Hi Francesco

Yeah sorry I imagine it wasn't implemented as a bit of an oversight more than anything. Because I developed stan_jm as part of a Phd in joint modelling, stan_jm was treated as somewhat of a first class citizen, whilst stan_mvmer was more of a building block for stan_jm and so probably didn't get implemented as comprehensively. ☹️

I'm a bit too long out of these details to remember exactly what limitations there might or might not be, but I think you are right in the sense that you should be able to take the fitted model, and then some newly observed data for a new individual, and draw some random effects for the new individual, and predict conditional on those. In terms of implementation, the only suggestion I'd have would be to fork the repo and try to implement it on your own. Sorry! And I'm not sure how easy or realistic that is... e.g. I think you'd have to look at this code and see if you need to modify to draw new random effects conditional only on the longitudinal data and not anything to do with the survival submodel. It might actually be that it doesn't require conditioning on the survival endpoint itself, in which case it might not be too hard to modify. But sorry my head has been too long out of those details to really provide much useful context without diving back into the code...

[frbrz]

Hi Sam,

Thanks a lot for the very detailed reply. And no need to apologize, I think all packages I have seen (JMBayes, JoineREML) only focus on the joint model case, so I must be the one making weird requests ;)

Thanks also for the lead on where to look in the codes. Unfortunately, I do not have much time atm... but if I ever have the time to try what you suggest, I will update this post.