--- title: "Parkinson's Data: Joint Model for Informative Dropout" output: rmarkdown::html_vignette vignette: > %\VignetteIndexEntry{Informative Dropout} %\VignetteEngine{knitr::rmarkdown} %\VignetteEncoding{UTF-8} --- ```{r, include = FALSE} knitr::opts_chunk$set(collapse = TRUE, comment = "#>") ``` Patients with faster UPDRS progression drop out earlier, so a naive longitudinal model is biased. `pmrm_joint()` links the two processes through a shared association structure. ```{r setup} library(progme) data(parkinsons_simulated) fit <- tryCatch( pmrm_joint(updrs ~ time * trt, Surv(time_to_dropout, dropout) ~ trt, data = parkinsons_simulated, association = "current_value", surv_model = "weibull"), error = function(e) e ) if (!inherits(fit, "error")) { s <- summary(fit) str(s$association_coefficients) } else { message("Demo fit skipped: ", conditionMessage(fit)) } ```