Title
stjmgraph -- Joint longitudinal and survival graph
Syntax
stjmgraph long_depvar [if] [in] [, options]
options Description ------------------------------------------------------------------------- panel(varname) panel identification variable censgraphopts(string) options to pass to the twoway plot for censored observations eventgraphopts(string) options to pass to the twoway plot for observations with the event combineopts(string) options to pass to the final graph combine draw displays all graphs used to create the final combined graph -------------------------------------------------------------------------
Description
stjmgraph creates a longitudinal trajectory plot, whereby the timescale is adjusted by taking away each patient's event/censoring time. This form of graph can be useful to display joint longitudinal and survival data, giving an indication of any association between the two processes. A separate plot is created for patients who were censored and for patients who experienced the event of interest. They are then combined using graph combine.
The dataset must be stset correctly into enter and exit times, using the enter option; see [ST] stset. stjmgraph uses _t0 to denote measurement times. For example, below we have 3 patients with 2, 5 and 3 measurements each, respectively.
--------------------------------- id _t0 _t _d long_resp --------------------------------- 1 0 0.2 0 0.93 1 0.2 0.7 0 1.32 2 0 0.5 0 1.15 2 0.5 1.2 0 1.67 2 1.2 1.6 0 1.92 2 1.6 1.9 0 2.65 2 1.9 2.6 1 3.15 3 0 2 0 0.25 3 2 2.3 0 0.21 3 2.3 2.4 1 0.31 ---------------------------------
See stjm for more details. Options
panel(varname) defines the panel identification variable.
censgraphopts(string) pass options to the twoway graph of censored observations. See twoway_options.
eventgraphopts(string) pass options to the twoway graph of observations who experienced the event of interest. See twoway_options.
combineopts(string) pass options to the final graph combine. See graph combine.
draw displays the intermediate twoway plots used to create the final graph.
Example
Load simulated example dataset: . use http://fmwww.bc.edu/repec/bocode/s/stjm_example
stset the data: . stset stop, enter(start) f(event=1) id(id)
Create joint plot: . stjmgraph long_response, panel(id)
Author
Michael J. Crowther, University of Leicester, UK. E-mail: mjc76@le.ac.uk.
Please report any errors you may find.
References
Crowther, M. J. & Abrams, K. R. & Lambert, P. C. Flexible parametric joint modelling of longitudinal and survival data. 2011. Submitted.
Wulfsohn, M. S. & Tsiatis, A. A. A joint model for survival and longitudinal data measured with error. Biometrics, 1997, 53, 330-339.
Also see
Online: stjm, stjm postestimation