Department of Health Statistics, Faculty of Preventive Medicine, Air Force Medical University, No.169 Changle West Road, Xi'an, 710032, Shaanxi, China.
BMC Med Res Methodol. 2021 Sep 21;21(1):192. doi: 10.1186/s12874-021-01378-8.
In follow-up studies, the occurrence of the intermediate event may influence the risk of the outcome of interest. Existing methods estimate the effect of the intermediate event by including a time-varying covariate in the outcome model. However, the insusceptible fraction to the intermediate event in the study population has not been considered in the literature, leading to effect estimation bias due to the inaccurate dataset.
In this paper, we propose a new effect estimation method, in which the susceptible subpopulation is identified firstly so that the estimation could be conducted in the right population. Then, the effect is estimated via the extended Cox regression and landmark methods in the identified susceptible subpopulation. For susceptibility identification, patients with observed intermediate event time are classified as susceptible. Based on the mixture cure model fitted the incidence and time of the intermediate event, the susceptibility of the patient with censored intermediate event time is predicted by the residual intermediate event time imputation. The effect estimation performance of the new method was investigated in various scenarios via Monte-Carlo simulations with the performance of existing methods serving as the comparison. The application of the proposed method to mycosis fungoides data has been reported as an example.
The simulation results show that the estimation bias of the proposed method is smaller than that of the existing methods, especially in the case of a large insusceptible fraction. The results hold for small sample sizes. Besides, the estimation bias of the new method decreases with the increase of the covariates, especially continuous covariates, in the mixture cure model. The heterogeneity of the effect of covariates on the outcome in the insusceptible and susceptible subpopulation, as well as the landmark time, does not affect the estimation performance of the new method.
Based on the pre-identification of the susceptible, the proposed new method could improve the effect estimation accuracy of the intermediate event on the outcome when there is an insusceptible fraction to the intermediate event in the study population.
在随访研究中,中间事件的发生可能会影响感兴趣结局的风险。现有的方法通过在结局模型中纳入时变协变量来估计中间事件的影响。然而,研究人群中对中间事件无抵抗力的部分在文献中尚未得到考虑,导致由于数据集不准确而导致估计偏差。
在本文中,我们提出了一种新的效应估计方法,首先识别易感亚组,以便在正确的人群中进行估计。然后,在识别出的易感亚组中,通过扩展的 Cox 回归和 landmark 方法来估计效应。对于易感性识别,将观察到中间事件时间的患者分类为易感。基于拟合中间事件发生率和时间的混合治愈模型,通过残差中间事件时间插补来预测censored 中间事件时间患者的易感性。通过 Monte-Carlo 模拟研究了新方法在各种情况下的效果,将现有方法的效果作为比较。作为一个例子,报告了该方法在蕈样真菌病数据中的应用。
模拟结果表明,与现有方法相比,该方法的估计偏差较小,尤其是在无抵抗力部分较大的情况下。该结果适用于小样本量。此外,新方法的估计偏差随着混合治愈模型中协变量的增加而减小,特别是连续协变量。易感和易感亚组中协变量对结局的影响以及 landmark 时间的异质性不影响新方法的估计性能。
基于易感的预先识别,当研究人群中对中间事件无抵抗力时,该新方法可以提高中间事件对结局的影响的估计准确性。