源自风暴统计网:一键统计分析与绘图的AI网站

  引言

今天这篇文章太有意思了,提供了一种从状态转移角度研究CHARLS纵向数据的新思路,再加上有高分指标肌少症傍身,一区Top文章不是手拿把掐!(附方法解读)
肌肉减少症是指骨骼肌质量、肌肉力量和身体功能全面下降的一种综合性病症,普遍存在于老年群体中。尽管在多种慢性病如阿尔茨海默病、糖尿病中已观察到病前状态存在双向转换,但肌少症的状态转换规律尚不明确。
近日,重庆医科大学附属医院用CHARLS数据库,在期刊Journal of Cachexia Sarcopenia And Muscle》(医学一区top,IF=9.4)发表了一篇题为:“Bidirectional transitions of sarcopenia states in older adults: The longitudinal evidence from CHARLS”的研究论文,旨在探讨我国老年人肌少症状态(包括非肌少症、可能肌少症、肌少症和死亡)的双向转变模式估计状态间的转移强度与概率,并识别影响转换的关键因素。

研究结果表明,我国老年人肌少症状态具有显著的双向转换特征,其中,年龄、BMI、身体功能损害程度是重要的影响因素。

研究基于中国健康与退休纵向研究(CHARLS)数据库2011~2015年的数据,经过纳排,最终纳入 4395 名 60 岁及以上、至少完成两次肌少症状态评估的个体。
研究方法主要采用连续时间多状态MSM模型,估计状态间的转移强度矩阵、1–5 年观察期内的转移概率、平均停留时间及各状态在死亡前的总停留时间,并利用风险比分析年龄、性别、BMI、身体功能等多种协变量对状态转换的影响。

主要研究结果

研究发现,对于可能肌少症患者,其逆转至非肌少状态的转移强度(0.252,95%CI 0.231–0.275)是恶化为肌少症的2.8倍(0.090,95%CI 0.080–0.100)

同时,MSM 模型结果显示,在一年观察期内,可能肌少症患者逆转为非肌少状态、进展为肌少症及死亡的概率估计值分别为0.181、0.066和0.035;而肌少症患者逆转为非肌少状态、改善为可能肌少状态及死亡的概率估计值则分别为0.016、0.125和0.075。
图2  1年、2年、3年和5年观察区间内的肌少症状态转换图
研究进一步指出,年龄、BMI 和身体机能等因素会显著影响肌肉减少症的双向转变的方向与强度。
图3 不同亚组中可能发生肌少症状态转变的估计可能性
此外,研究还预估了不同状态的肌少症的平均停留时间,并预测了这些状态在死亡前的总停留时间。
 图4 肌少症短暂状态的估计停留时间和预测总停留时间

综上所述,老年人肌少症状态具有显著的双向转换特征,尤其可能肌少症存在较高的自然恢复潜力。因此,未来可通过进行风险因素(如年龄、BMI、身体功能)筛查和加强干预可能促进恢复性转变

新方法–Markov模型

这篇文章除了应用了高分指标–肌少症外,另一个亮点就是出现了charls数据库中的研究新面孔——马尔可夫模型(MSM)。模型目前在其他领域使用较多,如机器学习研究领域,在医学领域中主要应用于疾病发展相关研究。

而多状态马尔可夫模型(MSM)是马尔可夫模型的扩展,它允许存在多个状态,并且能够更细致地描述复杂的健康状态变化过程。与简单的马尔可夫模型相比,MSM 可以更准确地反映现实中疾病的多样性和动态性。

为了方便理解,我们接下来结合本文研究,对该模型行拆解分析:

又是肌少症?CHARLS新思路踢开一区Top大门(IF 9.4)

1.马尔可夫性核心理念

即个体未来状态的转移风险仅取决于当前状态,而与历史状态路径无关。

eg:对于一个目前处于“可能肌少症”的个体,其下一步会好转、恶化还是死亡,只与他当前是“可能肌少症”这一状态有关,而与他之前是否已经在这个状态待了多久、或者是刚从“非肌少症”变过来无关。

2.状态空间

所有可能状态的集合。这些状态必须是互斥穷尽(即任何一个体在任一时刻必须且只能处于其中一个状态)。

eg:本文的状态空间为{非肌少症, 可能肌少症, 肌少症, 死亡}

3.转移强度

这是连续时间MSM的核心参数,也称为瞬时风险率。它表示在极短的时间区间内,从一个状态转移到另一个状态的速率。它描述了状态变化的“倾向”或“势头”。

eg:文章中的 q12=0.303表示,对于一个处于“非肌少症”的个体,他瞬时转移到“可能肌少症”的强度是0.303/人年。

4.转移概率

这是更具临床解释性的结果。转移强度是“速率”,而转移概率是在一段时间内发生转换的“可能性”。

eg:文中给出的“可能肌少症”者在1年内恢复为“非肌少症”的概率是0.181,这就是一个从转移强度推导出的、非常直观的临床指标。