药物经济学模型之Markov模型(四)——如何将事件参数纳入到模型中
姚佑楠

姚佑楠

广东省人民医院

擅长:心力衰竭管理与药物经济学研究
已关注
关注
2021-01-18 来源:医咖会

上一期(药物经济学模型之Markov模型(三):手把手教如何进行研究设计)我们完成了Markov模型骨架的完整搭建,其运行的内在逻辑是:不同NYHA心功能分级的65岁射血分数减低型心衰患者使用达格列净(10mg,qd)+标准化治疗或标准化治疗15年(预期寿命是80岁);每三个月作为一个周期,患者在每个周期内可无事件发生,也可以发生心衰住院或者死亡(心血管或非心血管死亡);发生心衰住院的患者在该周期内仍有再住院的风险,可能会出现无事件或者心衰再住院;在每个循环的最后,除了死亡的患者,其他患者会发生NYHA心功能分级的转换,然后进入下一个周期(如图)。

以上一期在TreeAge软件里构建的模型为基础,接下来我们把事件相关的参数纳入到模型中。首先,我们需要确定的是各类事件在三个月内的发生率,包括心衰住院率、心衰再住院率、心血管死亡率和非心血管死亡率(无事件的概率用1减去心衰住院率、心血管死亡率和非心血管死亡率三者之和就可以了),由于一共有两组(达格列净组和标准化治疗组),所以总计需要获得8个事件率的相关参数。

Markov模型的事件参数(尤其是主要终点事件及其组成成分)应来源于临床随机对照试验DAPA-HF研究的主要终点是心衰相关事件和心血管死亡,因此模型中关于心衰住院率、心衰再住院率、心血管死亡率我们就从DAPA-HF研究中提取。

我们假设在模型模拟的过程当中,心血管死亡率和心衰相关的住院率是固定不变的。根据DAPA-HF研究中随访18.2个月,达格列净组心血管死亡率为227/2373,标准化治疗组心血管死亡率为273/2371 [1]。我们通过公式:

三个月的心血管死亡率

= 1 -(1-心血管死亡率)3/18.2

计算得出达格列净组和标准化治疗组每个周期的心血管死亡率分别是1.644%1.996%

同理,达格列净组心衰住院率为231/2373,标准化治疗组心衰住院率为318/2371 [1],相应每个周期的心衰住院率为1.674%和2.346%。

在我们的模型中,所有的再住院都发生在心衰住院的患者中,根据DAPA-HF研究中的数据(文章中的Table 2)[1],我们假设两组发生心衰住院的患者其发生心衰再住院的风险一样,可以计算其风险为(567+742-227-273-231-318)/(231+318)=47.359%。

对于非心血管死亡,我们采用国家疾病控制中心的数据,模型模拟到不同年龄段时予以不同的非心血管死亡率[2],且假设组间非心血管死亡率相同。此处之所以不应用DAPA-HF研究的数据主要是考虑到:

(1)非心血管死亡不是DAPA-HF研究的主要终点,应用中国的数据可以更接近中国人口死亡的真实情况;

(2)DAPA-HF研究中非心血管死亡率通过换算,达格列净组约为每3个月0.343%,标准化治疗组约为0.393%,绝对值相对其他几个事件发生率而言小得多,且也在国家疾控中心数据范围之内;

(3)在后面的敏感性分析时,上述相应的数据也在分析范围之内。

上述参数总结在下表中:

接下来我们在TreeAge软件中展示一下如何将这些参数输入到模型中。如下图,点击最左侧的决策节,再点击Variable Definitions(红色方框)后,选择绿色加号即可弹出一个变量定义的对话框。

Name可以输入变量名称(CVd),Description可输入变量定义(cardiovascular mortality of dapagliflozin group),点击OK后弹出另一个对话框。

在“CVd=”中输入相应的数值(0.01644),点击OK。

有问题
找小咖
公众号
统计咨询
扫一扫添加小咖个人微信,立即咨询统计分析服务!
数据库搭建
想提升科研效率和数据质量?扫码添加小E,立即咨询数据采集与管理相关产品和服务!查看详细>>
意见反馈