Local partial likelihood estimation in proportional hazards regression

Abstract

Fan, Gijbels and King [Ann. Statist. 25 (1997) 1661–1690] considered the estimation of the risk function ψ(x) in the proportional hazards model. Their proposed estimator is based on integrating the estimated derivative function obtained through a local version of the partial likelihood. They proved the large sample properties of the derivative function, but the large sample properties of the estimator for the risk function itself were not established. In this paper, we consider direct estimation of the relative risk function ψ(x₂)−ψ(x₁) for any location normalization point x₁. The main novelty in our approach is that we select observations in shrinking neighborhoods of both x₁ and x₂ when constructing a local version of the partial likelihood, whereas Fan, Gijbels and King [Ann. Statist. 25 (1997) 1661–1690] only concentrated on a single neighborhood, resulting in the cancellation of the risk function in the local likelihood function. The asymptotic properties of our estimator are rigorously established and the variance of the estimator is easily estimated. The idea behind our approach is extended to estimate the differences between groups. A simulation study is carried out.