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]