# \donttest{
##
library(survival)
library(nlme)
library(dplyr)
jmcs1<-jmbayesBig(dtlong=long2,
dtsurv = surv2 ,
longm=y~ x7+visit,
survm=Surv(time,status)~x1+visit,
rd= ~ visit|id,
timeVar='visit',
nchain=1,
samplesize=200,
id='id')
mod3<-jmcs1
ydt<-long2%>%filter(id%in%c(900))
names(ydt)
cdt<-surv2[,'id']%>%filter(id%in%c(900))
names(cdt)
newdata<-full_join(ydt,cdt,by='id')
P2<-predJMbayes(model<-mod3,ids<-c(900),newdata=newdata,process = 'event')
plot(P2$p1[[1]])
##
# }
Run the code above in your browser using DataLab