The episodes from one patient share treatment, biological susceptibility and prior history, so they are not independent observations. Treating every row as an unrelated subject can substantially understate uncertainty. Two approaches address this within-patient recurrent-event dependence:
A patient-level bootstrap is another way to retain dependence: resample complete histories, not individual episode rows. It is not necessary to assume that robust marginal/working and frailty-conditional effect estimates target exactly the same parameter.
Shared frailty model 2026-10-07
A shared frailty model assigns the same latent positive hazard multiplier to correlated event histories or group members. Gamma or lognormal frailty distributions are common examples. Modeling dependence through this latent effect yields conditional covariate effects; integrating out frailty can change the form of marginal hazards. Eligibility and event history still have to be encoded in the at-risk indicator.
Repeated episodes from one person share susceptibility, treatment and history, so row-wise independence is generally inappropriate. Patient-clustered sandwich covariance matrices estimate uncertainty using whole-patient score contributions under a suitable working model. Shared frailty models instead model persistent latent heterogeneity. These methods require independent sampling units at the patient level and do not necessarily estimate the same marginal versus conditional effect.