Mandel Micha, Gauthier Susan A, Guttmann Charles R G, Weiner Howard L, Betensky Rebecca A
Micha Mandel is a postdoctoral fellow, Department of Biostatistics, Harvard School of Public Health, Boston, MA 02115 (
J Am Stat Assoc. 2007 Dec;102(480):1254-1266. doi: 10.1198/016214507000000059.
The expanded disability status scale (EDSS) is an ordinal score that measures progression in multiple sclerosis (MS). Progression is defined as reaching EDSS of a certain level (absolute progression) or increasing of one point of EDSS (relative progression). Survival methods for time to progression are not adequate for such data since they do not exploit the EDSS level at the end of follow-up. Instead, we suggest a Markov transitional model applicable for repeated categorical or ordinal data. This approach enables derivation of covariate-specific survival curves, obtained after estimation of the regression coefficients and manipulations of the resulting transition matrix. Large sample theory and resampling methods are employed to derive pointwise confidence intervals, which perform well in simulation. Methods for generating survival curves for time to EDSS of a certain level, time to increase of EDSS of at least one point, and time to two consecutive visits with EDSS greater than three are described explicitly. The regression models described are easily implemented using standard software packages. Survival curves are obtained from the regression results using packages that support simple matrix calculation. We present and demonstrate our method on data collected at the Partners MS center in Boston, MA. We apply our approach to progression defined by time to two consecutive visits with EDSS greater than three, and calculate crude (without covariates) and covariate-specific curves.
扩展残疾状态量表(EDSS)是一种用于衡量多发性硬化症(MS)病情进展的序贯评分。病情进展定义为达到特定EDSS水平(绝对进展)或EDSS增加1分(相对进展)。针对此类数据,传统的进展时间生存分析方法并不适用,因为它们没有利用随访结束时的EDSS水平。相反,我们建议使用适用于重复分类或序贯数据的马尔可夫转换模型。这种方法能够在估计回归系数并对所得转换矩阵进行操作后,推导出特定协变量的生存曲线。我们采用大样本理论和重采样方法来推导逐点置信区间,该方法在模拟中表现良好。文中详细描述了针对达到特定EDSS水平的时间、EDSS至少增加1分的时间以及连续两次就诊时EDSS大于3的时间生成生存曲线的方法。所描述的回归模型使用标准软件包即可轻松实现。生存曲线可通过支持简单矩阵计算的软件包从回归结果中获得。我们在马萨诸塞州波士顿的 Partners MS 中心收集的数据上展示并演示了我们的方法。我们将我们的方法应用于由连续两次就诊时EDSS大于3所定义的病情进展,并计算了粗生存率(无协变量)和特定协变量的曲线。