Abstract
In this paper, we review available methods for determination of the functional form of the relation between a covariate and the log hazard ratio for a Cox model. We pay special attention to the detection of influential observations to the extent that they influence the estimated functional form of the relation between a covariate and the log hazard ratio. Our paper is motivated by a data set from a cohort study of lung cancer and silica exposure, where the nonlinear shape of the estimated log hazard ratio for silica exposure plotted against cumulative exposure and hereafter referred to as the exposure–response curve was greatly affected by whether or not two individuals with the highest exposures were included in the analysis. Formal influence diagnostics did not identify these two individuals but did identify the three highest exposed cases. Removal of these three cases resulted in a biologically plausible exposure–response curve.
Acknowledgments
We are grateful to Prof. Harvey Checkoway for use of the silica data.
Disclosure statement
No potential conflict of interest was reported by the authors.
Funding
This research was supported by National Cancer Institute grant CA081345.