Regression analysis of multivariate current status data with dependent censoring

Multivariate current-status failure time data consist of several possibly related event times of interest, in which the status of each event is determined at a single examination time. If the examination time is intrinsically related to the event times, the examination is referred to as dependent censoring and needs to be taken into account. Such data often occur in clinical studies and animal carcinogenicity experiments. To accommodate for possible dependent censoring, this paper proposes a joint frailty model for event times and dependent censoring time. We develop a likelihood approach using Gaussian quadrature techniques for obtaining maximum likelihood estimates. We conduct extensive simulation studies for investigating finite-sample properties of the proposed method. We illustrate the proposed method with an analysis of patients with ankylosing spondylitis, where the examination time may be dependent on the event times of interest.