bivrecReg {BivRec}  R Documentation 
Semiparametric Regression Analysis of Bivariate Alternating Recurrent Event Gap Time Data
Description
This function allows the user to evaluate covariate effects on two alternating recurrent events gap times (referred as Type I and Type II gap times) under the assumption that the two gap times follow accelerated failure time (AFT) models. See details for the estimation methods provided.
Usage
bivrecReg(formula, data, method)
Arguments
formula 
A formula with a 
data 
A data frame that includes the vectors needed for the 
method 
A string indicating which method to use to estimate effects of the covariates. See details. 
Details
Two different estimation methods are available:

method = "Lee.et.al"
(default) is a Ustatisticsbased smooth estimating function approach. See Lee, Huang, Xu, Luo (2018) for further details. 
method = "Chang"
is a rankbased estimating function approach. See Chang (2004) for further details. Note that following the Chang method, the variances of the estimated regression coefficients are approximated using the resampling techniques developed by Parzen, Wei, Ying (1994). This approximation requires extensive computing time for a relatively small sample size. In addition, using the Chang method does not guarantee convergence for the estimation of the coefficients and user may get the message, "Error: Max iterations reached. Did not converge.".
Related methods: coef.bivrecReg
, confint.bivrecReg
, plot
, print
, summary.bivrecReg
, vcov.bivrecReg
.
Value
A bivrecReg object that contains:

call

lee_fit
orchang_fit

formula

data
References
Chang SH. (2004). Estimating marginal effects in accelerated failure time models for serial sojourn times among repeated events. Lifetime Data Analysis, 10: 175190. doi: 10.1023/B:LIDA.0000030202.20842.c9
Lee CH, Huang CY, Xu G, Luo X. (2018). Semiparametric regression analysis for alternating recurrent event data. Statistics in Medicine, 37: 9961008. doi: 10.1002/sim.7563
Parzen MI, Wei LJ, Ying Z. (1994). A resampling method based on pivotal estimating functions. Biometrika, 81: 341350. http://www.people.fas.harvard.edu/~mparzen/published/parzen1.pdf
Examples
## Not run:
library(BivRec)
set.seed(28)
# Simulate bivariate alternating recurrent event data
bivrec_data < simBivRec(nsize=100, beta1=c(0.5,0.5), beta2=c(0,0.5),
tau_c=63, set=1.1)
# Apply Lee, Huang, Xu, Luo (2018) method using two covariates
lee_reg < bivrecReg(bivrecSurv(id, epi, xij, yij, d1, d2) ~ a1 + a2,
bivrec_data, "Lee.et.al")
summary(lee_reg)
## End(Not run)