JMFit1 <- mvJointModelBayes(MixedModelFit1, CoxFit1, timeVar = "year")
summary(JMFit1)

## Output
## Call:
## mvJointModelBayes(mvglmerObject = MixedModelFit1, coxphObject = CoxFit1, 
##     timeVar = "year")
## 
## Data Descriptives:
## Number of Groups: 312		Number of events: 169 (54.2%)
## Number of Observations:
##   log(serBilir): 1945
## 
## Random-effects covariance matrix:
##       StdDev    Corr
## (I)1  0.9839  (Int)1
## yer1  0.1749  0.3902
## 
## Survival Outcome:
##                     PostMean  StDev  StErr    2.5%  97.5%     P
## drugD-penicil        -0.1966 0.1501 0.0047 -0.3742 0.2224 0.646
## age                   0.0374 0.0069 0.0002  0.0295 0.0564 0.000
## log(serBilir)_value   1.3523 0.0890 0.0029  1.1501 1.4991 0.000
## 
## Longitudinal Outcome: log(serBilir) (family = gaussian, link = identity)
##             PostMean  StDev  StErr   2.5%  97.5% P
## (Intercept)   0.6427 0.0589 0.0019 0.3851 0.6109 0
## year          0.1808 0.0133 0.0004 0.1519 0.2043 0
## sigma         0.3391 0.0068 0.0002 0.3370 0.3639 0
##
## [truncated]