gistfile1.txt

total$timef <-as.factor(total$time)
total$BMI<-total$Weight_at_consultation/((total$Length_at_consultation/100)^2)

## Order data and remove duplicates
head(total%>%group_by(Patientnr,time)%>%arrange(Patientnr,time));dim(total)
total<-total[!duplicated(total),];dim(total)

## The time and Survtime need to be on the same scale for joint modelling  - transform to weeks
total$SurvTimeWeek<-total$SurvTime/7
pancreas$SurvTimeWeek<-pancreas$SurvTime/7