Semiparametric estimation of the accelerated mean model with panel count data under informative examination times

Determination of correlations in multivariate count data with informative observation times

We consider there are various types of recurrent events and the total number of occurrences are collected at the random observation times. It has concerned that the observation process may not be independent to the multivariate event processes, hence the total counts and observation times may be correlated and the dependence may exist among different types of the event processes as well. Many methods have developed nonparametric models to accommodate such unknown structures; however, it is difficult to assess and directly quantify their correlation relationships. A multivariate frailty model is proposed to this study, in which the event and observation processes are linked by frailty variables whose joint distribution can be implicitly specified through the multivariate normal distribution with some unknown covariance matrix. The Bayesian inference method is conducted to obtain the estimates of the regression coefficients and correlation parameters. We use a form of trigonometric functions to represent the covariance matrix, so that it meets the positive-definiteness condition efficiently during the estimation schemes. The simulation studies demonstrate the utility of the proposed models. We apply the model to a skin cancer prevention study, and aim to determine the covariate and association effects. We found treatment is significant in determining the duration of examination times; prior-counts, age and gender are significant variables on the occurrence rates of tumor counts. Using the covariance matrix to access the underlying dependent structure, the mutual correlations among them are all positive, and the basal cell counts are more related to the examination times.

Conditional modeling of panel count data with partly interval-censored failure event

In longitudinal follow-up studies, panel count data arise from discrete observations on recurrent events. We investigate a more general situation where a partly interval-censored failure event is informative to recurrent events. The existing methods for the informative failure event are based on the latent variable model, which provides indirect interpretation for the effect of failure event. To solve this problem, we propose a failure-time-dependent proportional mean model with panel count data through an unspecified link function. For estimation of model parameters, we consider a conditional expectation of least squares function to overcome the challenges from partly interval-censoring, and develop a two-stage estimation procedure by treating the distribution function of the failure time as a functional nuisance parameter and using the B-spline functions to approximate unknown baseline mean and link functions. Furthermore, we derive the overall convergence rate of the proposed estimators and establish the asymptotic normality of finite-dimensional estimator and functionals of infinite-dimensional estimator. The proposed estimation procedure is evaluated by extensive simulation studies, in which the finite-sample performances coincide with the theoretical results. We further illustrate our method with a longitudinal healthy longevity study and draw some insightful conclusions.

Regression Modeling for Recurrent Events Possibly with an Informative Terminal Event Using R Package reReg

Recurrent event analyses have found a wide range of applications in biomedicine, public health, and engineering, among others, where study subjects may experience a sequence of event of interest during follow-up. The R package reReg offers a comprehensive collection of practical and easy-to-use tools for regression analysis of recurrent events, possibly with the presence of an informative terminal event. The regression framework is a general scale-change model which encompasses the popular Cox-type model, the accelerated rate model, and the accelerated mean model as special cases. Informative censoring is accommodated through a subject-specific frailty without any need for parametric specification. Different regression models are allowed for the recurrent event process and the terminal event. Also included are visualization and simulation tools.

On computation of semiparametric maximum likelihood estimators with shape constraints

Large sample theory of semiparametric models based on maximum likelihood estimation (MLE) with shape constraint on the nonparametric component is well studied. Relatively less attention has been paid to the computational aspect of semiparametric MLE. The computation of semiparametric MLE based on existing approaches such as the expectation-maximization (EM) algorithm can be computationally prohibitive when the missing rate is high. In this paper, we propose a computational framework for semiparametric MLE based on an inexact block coordinate ascent (BCA) algorithm. We show theoretically that the proposed algorithm converges. This computational framework can be applied to a wide range of data with different structures, such as panel count data, interval-censored data, and degradation data, among others. Simulation studies demonstrate favorable performance compared with existing algorithms in terms of accuracy and speed. Two data sets are used to illustrate the proposed computational method. We further implement the proposed computational method in R package BCA1SG, available at CRAN.

Semiparametric Regression Analysis of Panel Count Data: A Practical Review

Panel count data arise in many applications when the event history of a recurrent event process is only examined at a sequence of discrete time points. In spite of the recent methodological developments, the availability of their software implementations has been rather limited. Focusing on a practical setting where the effects of some time-independent covariates on the recurrent events are of primary interest, we review semiparametric regression modelling approaches for panel count data that have been implemented in R package spef. The methods are grouped into two categories depending on whether the examination times are associated with the recurrent event process after conditioning on covariates. The reviewed methods are illustrated with a subset of the data from a skin cancer clinical trial.