simPH: An R package for illustrating estimates from cox proportional hazard models including for interactive and nonlinear effects

The R package simPH provides tools for effectively communicating results from Cox proportional hazard (PH) models, including models with interactive and nonlinear effects. The Cox (PH) model is a popular tool for examining event data. However, previously available computational tools have not made it easy to explore and communicate quantities of interest and associated uncertainty estimated from them. This is especially true when the effects are interactions or nonlinear transformations of continuous variables. These transformations are especially useful with Cox PH models because they can be employed to correctly specifying models that would otherwise violate the nonproportional hazards assumption. Package simPH makes it easy to simulate and then plot quantities of interest for a variety of effects estimated from Cox PH models including interactive effects, nonlinear effects, as well as standard linear effects. Package simPH employs visual weighting in order to effectively communicate estimation uncertainty. There are options to show either the standard central interval of the simulation's distribution or the shortest probability interval - which can be useful for asymmetrically distributed estimates. This paper uses hypothetical and empirical examples to illustrate package simPH 's syntax and capabilities.

[1]  Ying Liu,et al.  Simulation-efficient shortest probability intervals , 2013, Stat. Comput..

[2]  Thomas Brambor,et al.  Understanding Interaction Models: Improving Empirical Analyses , 2006, Political Analysis.

[3]  May,et al.  [Wiley Series in Probability and Statistics] Applied Survival Analysis (Regression Modeling of Time-to-Event Data) || Extensions of the Proportional Hazards Model , 2008 .

[4]  P. V. Rao,et al.  Applied Survival Analysis: Regression Modeling of Time to Event Data , 2000 .

[5]  R. Michael Alvarez,et al.  Event History Modeling: A Guide for Social Scientists , 2004 .

[6]  P. Grambsch,et al.  Martingale-based residuals for survival models , 1990 .

[7]  L. Keele Semiparametric Regression for the Social Sciences , 2008 .

[8]  Jason Wittenberg,et al.  Making the Most Of Statistical Analyses: Improving Interpretation and Presentation , 2000 .

[9]  A. Licht Change Comes with Time: Substantive Interpretation of Nonproportional Hazards in Event History Analysis , 2011, Political Analysis.

[10]  J. Golub,et al.  How Time Affects EU Decision-Making , 2007 .

[11]  Rob J Hyndman,et al.  Computing and Graphing Highest Density Regions , 1996 .

[12]  P. Grambsch,et al.  Proportional hazards tests and diagnostics based on weighted residuals , 1994 .

[13]  D.,et al.  Regression Models and Life-Tables , 2022 .

[14]  Hadley Wickham,et al.  An Implementation of the Grammar of Graphics , 2015 .

[15]  Paul H. C. Eilers,et al.  Flexible smoothing with B-splines and penalties , 1996 .

[16]  Gary King,et al.  Zelig: Everyone's Statistical Software , 2006 .

[17]  Chandler Stolp,et al.  The Visual Display of Quantitative Information , 1983 .

[18]  Cedric E. Ginestet ggplot2: Elegant Graphics for Data Analysis , 2011 .

[19]  Luke Keele,et al.  Proportionally Difficult: Testing for Nonproportional Hazards in Cox Models , 2010, Political Analysis.

[20]  J. Klenk,et al.  Use of penalized splines in extended Cox-type additive hazard regression to flexibly estimate the effect of time-varying serum uric acid on risk of cancer incidence: a prospective, population-based study in 78,850 men. , 2009, Annals of epidemiology.

[21]  Janet M. Box-Steffensmeier,et al.  Duration models and proportional hazards in political science , 2001 .

[22]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[23]  Solomon Hsiang,et al.  Visually-Weighted Regression , 2013 .

[24]  Daniel Carpenter,et al.  Groups, the Media, Agency Waiting Costs, and FDA Drug Approval , 2002 .

[25]  G. C. Tiao,et al.  Bayesian inference in statistical analysis , 1973 .

[26]  David W. Hosmer,et al.  Applied Survival Analysis: Regression Modeling of Time-to-Event Data , 2008 .

[27]  A. Dreher Modeling Survival Data Extending The Cox Model , 2016 .

[28]  Christopher Zorn,et al.  Nonproportional Hazards and Event History Analysis in International Relations , 2003 .

[29]  P. Grambsch,et al.  Modeling Survival Data: Extending the Cox Model , 2000 .