基于Markov链的上海市枫林社区骨质疏松全人群防治工作的成本效果分析

苏瑾　寿涓　顾文钦　易春涛　徐莉苹　程莉莉　丁宏娟　周鹏　吴颖华　秦杰　薛斌　魏百川　王谦　彭燕　程毅敏　杨蓝　卫洋洋　王磊　祁瑨麟　邵迎　蔡立明

【摘要】 背景 隨着社区卫生服务综合改革的逐步深化，社区对全人群的单病种防治工作也在不断优化，但目前仍缺乏合适的方法来评价社区防治的成效。目的 通过分析2016—2022年度上海市徐汇区枫林街道社区卫生服务中心现有骨质疏松全人群防治措施的成本效果，探究社区单病种全人群防治的初步成效及可能存在的问题。方法 本研究以4 293例枫林社区2016—2022年度接受骨质疏松防治的社区居民为研究对象进行成本-效果分析。根据骨密度值将研究对象划分为三种状态：健康、低骨量、骨质疏松，以此构建Markov模型，分析不同状态转移的影响因素，并依据Markov模型预测值结合质量调整生命年（QALY）计算QALY增量，统计2016—2022年防治工作投入的所有成本，以成本/QALY增量比来评价防治效果，通过文献检索确定的健康效用值和期望寿命计算求得QALY。结果 枫林社区骨质疏松防治总成本为33 814 102.15元。第一次诊断时社区骨质疏松防治人群的总QALY为77 098.288 9，人均均值为17.959，标准差为9.34；第二次诊断时防治人群的总QALY为79 616.936 1，人均均值为18.546，标准差为9.342，两次诊断QALY的差值为2 518.647 2。增量成本效果比（ICER）为2 132.907 0元/QALY，ICER值小于1倍（66 965.10元）人均GDP，表明干预方案完全具有成本效果。结论 本研究结合Markov模型及QALY对枫林社区持续接受骨质疏松随访管理的人群的防治效果进行了初步评价，结果显示枫林社区的骨质疏松防治干预方案是值得投入的，为后续骨质疏松的精准防治决策提供了一定的依据。

【关键词】 骨质疏松；社区单病种；防治；成本效果分析；马尔科夫模型

【中图分类号】 R 681 【文献标识码】 A DOI：10.12114/j.issn.1007-9572.2023.0420

Cost Effectiveness Analysis of the Prevention and Treatment of Osteoporosis among the Entire Population in Fenglin Community，Shanghai Based on Markov Chain

SU Jin 1，SHOU Juan2*，GU Wenqin3，YI Chuntao4，XU Liping1，CHENG Lili1，DING Hongjuan1，ZHOU Peng1，

WU Yinghua1，QIN Jie1，XUE Bin1，WEI Baichuan1，WANG Qian1，PENG Yan1，CHENG Yimin1，YANG Lan1，

WEI Yangyang1，WANG Lei1，QI Jinlin1，SHAO Ying1，CAI Liming1

1.Shanghai Community Health Service Center，Fenglin Street，Xuhui District，Shanghai 200030，China

2.Zhongshan Hospital Affiliated to Fudan University，Shanghai 200032，China

3.Shanghai Xuhui District Health Management and Development Center，Shanghai 200237，China

4.Shanghai Xuhui District Health Commission Supervision Office，Shanghai 200030，China

*Corresponding author：SHOU Juan，Chief physician；E-mail：asyura55@126.com

【Abstract】 Background With the gradual deepening of the comprehensive reform of community health services，the prevention and control of single diseases for the entire population in the community is also constantly being optimized，However，there is still a lack of appropriate evaluation methods in the community to reflect its prevention and control effectiveness evaluation. Objective By analyzing and exploring the cost effect of the existing whole-population prevention and control measures of osteoporosis in Fenglin community，from 2016 to 2022，this paper answers the initial results and possible problems of the whole-population prevention and control of single disease in the community. Methods This study conducted a cost-effectiveness analysis on 4 293 community residents who received osteoporosis prevention and treatment in Fenglin community from 2016 to 2022. Divide the research object into three states based on bone density values：healthy，low bone mass，and osteoporosis，and construct a Markov model to analyze the influencing factors of different state transitions. And based on the predicted values of the Markov model combined with quality adjusted life years （QALY），the increment of quality adjusted life years （QALY） is calculated. All costs invested in prevention and control work from 2016 to 2022 are counted，and the effectiveness of prevention and control is evaluated by the cost/quality adjusted life year increment ratio. Quality adjusted life years （QALY） are obtained by calculating the health utility value and expected life determined through literature search. Results The total cost of osteoporosis prevention and treatment in Fenglin community was 33 814 102.15 yuan. The total quality-adjusted life years of 4 293 community osteoporosis prevention and treatment population were 77 098.288 9 at the first diagnosis，the average per capita was 17.959，and the standard deviation was 9.34. At the second diagnosis，the total quality-adjusted life years were 79 616.936 1，the mean per capita was 18.546，the standard deviation was 9.342，and the difference between the two diagnoses was 2 518.647 2.

The incremental cost-effectiveness ratio （ICER） is 2 132.907 0 yuan /QALY，and the incremental cost-effectiveness ratio （ICER） is less than 1 times （66 965.10 yuan） per capita GDP，indicating that the intervention scheme is fully cost-effective. Conclusion This study combined Markov model and quality-adjusted life years （QALY） to make a preliminary evaluation of the prevention and treatment effect of osteoporosis in Fenglin community. The results showed that the intervention program of osteoporosis prevention and treatment in Fenglin community was worth investing，and provided a certain basis for the accurate prevention and treatment decision of osteoporosis in the future.

【Key words】 Osteoporosis；Community single disease；Prevention and cure；Cost-effectiveness analysis；Markov model

作为国民健康的“守卫者”，基层社区卫生服务中心在慢性病管理和防治中发挥着至关重要的作用。骨质疏松（osteoporosis，OP）是严重危害老年人健康的常见慢性病，据估计截至2020年底我国罹患OP的人群增至 2.866 亿［1］，其中骨质疏松性骨折的患者占 10.0％以上［2］，且65岁以上老人的OP患病率高达32.0%。虽然基层社区正在逐步开展OP防治工作，但尚缺乏与社区人群OP罹患状态及未来发展趋势相关的研究，因此需要有效的工具来进一步探索社区防治现状和预测未来趋势。

多状态Markov模型［3］既能用于描述当下状态，也能用来预测未来状态，可作为处理慢性病多状态资料的工具，也可测量某种干预所产生的状态变化总量，并能够将状态变化进一步折算成效用或效益，计算投入产出比。基于此，本研究以上海市徐汇区枫林街道社区卫生服务中心2016—2022年度在自然状态下推进的社区OP全人群防治工作为例，通过构建社区OP防治人群的Markov模型，计算防治过程中投入的成本以及所带来的质量调整生命年（quality adjusted life years，QALY）［4］的产出，进行OP全人群防治的成本-效用分析，评价社区卫生某单病种防治服务的实际成效，为当前社区卫生服务单病种防治工作的开展提供参考和借鉴，为针对OP的精准防治提供改进方向。

1 对象与方法

1.1 研究对象

研究者选取2016年5月—2022年8月间接受上海市徐汇区枫林街道社区卫生服务中心OP防治的社区居民为研究对象。纳入标准：（1）符合WHO原发性OP的诊断标准［5］；（2）55～80岁；（3）所纳入病例与其主要照护者认知能力均正常；（4）自愿参与、接受治疗方案，并签署知情同意书。排除标准：（1）有严重心血管或呼吸系统疾病史，合并心、脑、肾等脏器急性疾病；（2）近期参与过其他类似研究；（3）有精神疾病或因故无法完成研究；（4）有严重躯体疾病如恶性肿瘤。根据以上纳排标准，最终共计筛选出4 293例研究对象。本研究已通过上海临床研究伦理委员会伦理学审查（伦审号SECCR/2021-35-01）。

1.2 研究方法

1.2.1 数据收集：社区全科医生及专病医生以社区为单位开展宣传，招募研究对象，告知注意事项并指导其签署知情同意书，进行专病管理数据整理汇总（在此期间进行第一次诊断），随后对其进行持续随访（在此期间进行第二次诊断）。研究数据来源主要分为两类，一类为医院信息系统（Hospital Information System，HIS）可详尽记录医疗服务及与之相对应的数据，据此可统计研究对象使用过的所有药物名称与数量、检查名称与数量、门诊转诊名称和数量等数据；一类为未在HIS中记录但确系向研究对象提供过的各类服务及与之相对应的数据（其数据出自上海市徐汇区枫林街道社区卫生服务中心OP服务记录册），据此可统计服务主题、服务对象、每次执行人、服务起止时间等数据［6］。

1.2.2 模型构建：根据骨密度测量值划分骨密度状态［5］，将骨密度状态分为三种：健康（S1）、低骨量（S2）、OP（S3）［7-10］。通过比较2次诊断结果计算状态转移概率，构建各状态间的转移频数和概率矩阵，并结合初始分布（根据首次骨密度的例数合计4 293例进行计算）对未来30年枫林社区不同骨密度人群的概率分布进行预测及多状态轉移的影响因素分析（根据连续2次骨密度测量的例数合计362例数据进行计算）。

1.2.3 成本分析：成本参数采用医疗服务成本［11］来统计，“成本”指2016年5月—2022年8月间OP防治环节真正投入使用的直接成本。根据枫林社区实际防治情况，此次需要支出成本的项目有：（1）门诊总费用（挂号、处方等）；（2）门诊诊疗费（检测等）；（3）住院费；（4）化验费；（5）药品费；（6）枫骨俱乐部人工成本。其中门诊总费用、门诊诊疗费、住院费和化验费通过HIS获取；药品费总成本依据2016—2019年总药品费的均值（612 030.3元）×时间（2016年5月—2022年8月）获取；枫骨俱乐部人工成本为枫林社区进行OP防治健康宣教投入人力所花费成本，计算方式为：每年宣教次数×人数×平均工资（3 000元/人）。

1.2.4 效果分析：本研究效用参数采用QALY来统计，通过计算干预前后总的QALY产出变化，以评价社区卫生服务大面积、广范围、长周期的干预是否产生了可观的成效及尚存哪些问题。QALY的计算公式如下［4，7，12-14］：QALY=期望寿命×健康效用值。

本研究中的期望寿命引用了参考文献［14］中的期望寿命值数据，分别为20～（63.36）、25～（58.36）、30～（53.36）、40～（43.59）、45～（38.81）、50～（34.01）、55～（29.41）、60～（24.84）、65～（20.3）、70～（16.26）、75～（12.33）、80～（9.43）、85～（6.67），85岁及以上（占比7.2%）年龄段期望寿命均设置为6.67。

本研究的健康效用值通过参考文献［15-18］确定，选取文献中年龄段人数最多的效用值进行下一步计算，最终四种状态［包括健康、低骨量、OP（骨折）、死亡］的健康效用值分别为：1［15］、0.89［16］、0.84［17］、0［18］。骨密度状态评价：经过连续的双能X线吸收检测法测量骨密度值的362例研究对象使用实际测量值；而无连续测量值的3 931例研究对象，其骨密度状态是根据上述362例研究对象建立的Markov模型推算而来。

1.2.5 成本效果评价：对于成本效果计算结果，使用WHO的增量成本效果比（ICER）进行评判（2010 年

版）［19］，其表示在干预下每多获得单位效用所增加的成本。在本研究中，ICER=增量成本/QALY，反映延长患者一定质量的生存时间所需的额外支出。WHO推荐：当 ICER<1倍人均GDP时，新的治疗方式增加的成本是非常值得的；当ICER<3倍人均GDP时，增加的成本是可以接受的；当ICER>3倍人均GDP，成本投入过大，则防治工作是需要考虑的。

1.2.6 单变量敏感度分析：为探究不同参数对增量成本效用比值的影响，应用单变量敏感度分析的方法比较5个参数，用龙卷风图展示结果，当某个参数在其范围内波动时，保持其他参数不变。

1.3 统计学方法

调查对象的基线特征及诊断结果用频数及构成比进行描述，此部分的基础统计在SPSS 25.0、Excel 2016软件中实现。多状态Markov模型的构建、转移强度和转移概率、平均逗留时间均在R 4.1.2软件中的“msm” 软件包中完成，影响因素分析时对分类变量设置哑变量，各变量对状态间转移强度的影响采用HR（95%CI）表示。双侧检验，检验水准α=0.05。

2 结果

2.1 基本情况

研究对象平均年龄为（68.4±11.8）岁，以60～70岁居多（占53.5％），体质量（62.08±14.52）kg，身高（157.8±16.7）cm。BMI在正常范围的共1 986例（46.0%）；其次为超重，共1 585例（37.0%）；性别以女性为主，占总研究对象的66.3%，绝经年龄为（50.1±2.6）岁，具体人口统计学特征如表1所示。

2.2 多状态Markov模型

2.2.1 Markov模型转移概率矩阵：骨密度各状态的转移频数矩阵如表2所示，各状态的转移概率矩阵如表3所示。在一个周期（即以1年为循环周期）内，健康转换为低骨量的概率为8.5%，低骨量转化为健康概率为16.4%，低骨量转化为OP概率为17.9%，OP转化为低骨量的概率为24.9%；在一个周期内，健康状态和OP状态不会发生相互转换。社区OP防治人群不同骨密度状态的转换概率见图1。

2.2.2 未来30年不同骨密度状态概率分布预测：依据2.2.1中的不同状态的转移概率，结合初始分布（根据首次骨密度的例数合计4 293例进行计算）探究未来30年可能出现的各状态概率分布情况。表4体现出从首次诊断之后的30年中骨密度在各状态所占比例变化情况，其中健康状态概率为0.527，低骨量状态概率为0.274，OP状态概率为0.199。随着时间的推移，未來30年间低骨量和OP状态人群所占比例将有所减少，健康状态所占比例有所增加，但概率最终趋于稳定，骨密度的变化趋势见图2。

2.2.3 社区OP防治人群不同骨密度状态转移的影响因素分析：剔除BMI及疾病史空值5例，以BMI、年龄、性别、是否腰背疼痛或全身骨痛、有无脆性骨折史、是否明显变矮和是否驼背为自变量，以不同的状态转移（S1→S2，S2→S1，S2→S3，S3→S2）为因变量，以女、无腰背疼痛或全身骨痛、无脆性骨折史、无明显变矮、无驼背为各分类自变量的对照，拟合msm模型估计不同状态转换的HR从而进行影响因素分析。表5仅表示相邻骨密度状态转移，结果显示，BMI、年龄、性别和是否明显变矮是OP防治人群骨密度状态转移的影响因素（P<0.05）。与BMI较低的患者相比，BMI趋于正常的患者S2→S1的转移强度较高；与年龄较小的患者相比，年龄越大的患者S3→S2的转移强度较低；与女性相比，男性S2→S1的转移程度较高；与未明显变矮的患者相比，明显变矮的患者S1→S2的转移程度较高（表5）。

2.3 成本投入与QALY产出分析

2.3.1 投入的成本总额与分布：枫林社区OP防治总成本为33 814 102.15元，其中门诊总费用成本为

14 464 180.49元，门诊诊疗费用成本为6 681 264.20元，住院总成本为8 485 958.69元，化验费总成本为

129 575.00元，药品费总成本为3 828 123.77元，人工总成本为225 000.0元，其中门诊总费用成本占比最高，为42%，共管理研究对象4 293例，其中OP患者1 584例，低骨量人群2 038例，健康人群671例，详见表6。

2.3.2 产出的QALY：根据表4计算4 293例研究对象的QALY。第一次诊断出4 293例社区OP防治人群的总QALY为77 098.288 9，人均均值为17.959，标准差为9.34；根据第二次诊断收集到的362例研究对象数据建立Markov模型，预测其干预后的诊断结果，进行QALY的计算，计算得出总QALY为79 616.936 1，人均均值为18.546，标准差为9.342，两次诊断QALY的差值为2 518.6472。

2.3.3 投入产出比：成本部分以总成本的年平均值，按照2016—2022年6个年份计算，以确保对成本的充分估计，则社区OP防治人群每年在干预方案上的耗费为5 372 040.308元，一年跟踪期带来的社区人群OP防治结果是QALY增加2 518.647 2。

进一步测算ICER：ICER=ΔC/ΔE=5 372 040.308/2 518.647 2=2 132.907 0元/QALY。2016—2021年期间我国人均GDP平均值为66 965.10元（2016年53 783元；2017年59 592.25元；2018年65 533.74元；2019年70 077.69元；2020年71 828.15元；2021年80 975.79元），ICER值小于1倍（66 965.10元）人均GDP ，干预方案完全具有成本效果。

2.3.4 单因素敏感性分析：单变量敏感性分析结果见图3。研究指标ICER对于健康到低骨量的转移概率最为敏感，接下来依次为低骨量到健康的转移概率和低骨量到OP的转移概率，CPI取值的敏感性最弱。同时，不管参数如何变化，整体始终小于1倍（66 965.10元）人均GDP，即社区OP防治始终具有成本效用，由此可见确定性分析结果较为稳健。

3 讨论

3.1 基层社区防治评价工具现状

我国现尚缺乏合适的评价基层社区防治成果的有效工具，其发展受到了多种因素的制约：（1）干预模式的不同，与传统的干预模式相比，社区针对全人群防治，改变的是整个人群的状态，传统的对照试验难以测量出社区防治的效果；（2）结构不同，社区防治工作覆盖范围较为宽泛，包括社区动员、筛查、诊断、干预、随访等环节，需形成闭环管理体系，其过程中的人力、物力和财力结构与传统医疗不同，不少环节没有设置明确的收费标准［20］；（3）面对的社区人群初始疾病状态、过程干预方案、干预进程中停留时间和临床结局指标皆不同，因此传统的临床疗效分析思维也不再适合社区防治工作［6］。

Markov模型是描述一类随机动态系统的模型，具有无后效性和稳定性，其在临床疗效评价、筛查试验的评价、疾病预后因素的研究等领域广泛应用。20世纪80年代，临床医学开始将Markov模型用于模拟慢性疾病的发展过程，20世纪90年代后，逐渐应用到卫生经济学领域。冯鑫等［8］建立Markov经济学模型分析了阿伦磷酸钠预防绝经后妇女OP性骨折的成本-效用；SI等［11］利用Markov经济学模型评估了采用不同手段对处于不同骨折状态下的虚拟患者进行干预的效果，都描述了Markov模型在骨骼疾病状态分析方面的应用情况。

目前社区所用技术大部分都是被验证过疗效和安全的成熟技术，当社区卫生机构将已被验证过的成熟技术结合社区具体情况进行创新性组合应用时，面向社区人群的干预、管理的成效和评价指标应需要重新考量［21］，以此对社区防治效果进行精准评价。利用Markov模型则可以动态、准确地评价社区人群经社区卫生机构慢性病管理后的疾病状态，预测系统变量的现在状态及其变化趋势，预测其在未来某一特定时间可能出现的状态，从而为决策提供依据，引导社区卫生服务精细化发展和竞争式发展，避免粗放地盲目推行或长期依赖于大量的过程性考核指标来驱动管理。

3.2 社区OP防治人群各状态及其影响因素分析

各状态的转移概率分析结果提示，接受OP防治的骨密度状态为健康的人群健康状态可维持在较高概率，相比纠正OP患者（转换为低骨量），防止健康人群和低骨量人群地向更差的状态转化能阻断更远期危害的发生，对低骨量人群进行合理干预能有效控制OP，而健康人群也无法直接进入OP状态，其提示低骨量人群是社区OP防治工作的重点，但同时也应做好健康人群的管理。

本研究构建的Markov链模型及未来30年间各状态的构成比结果显示，在干预的状态下，社区人群中OP患者和低骨量患者越来越少，健康人群持续增长，意味着从长远期看相较于及时干预，不及时干预或不干预将会消耗更多医疗资源。因此未来社区全科医生及其团队应更早、更积极地投入到慢性病高危人群的防治上来。

本研究结果表明，BMI、年龄、性别和是否明显变矮是OP患者不同状态转移的影响因素，与BMI趋于正常范围的患者相比，BMI偏低的患者向健康状态转移的风险较低，向OP状态转移的风险较高；年长的患者从低骨量转移成OP状态的风险较高；女性从健康向低骨量状态转移的风险较男性高；明显变矮的患者从低骨量转健康较未明显变矮的患者轉移风险较低，这提示社区卫生人员在开展OP防治工作时，应对偏瘦、年长的女性患者有所侧重。

3.3 投入成本效用分析

从2.3.3的投入产出比来看，本研究所在机构正在运行的社区OP防治干预方案对人群产生的QALY增加明显，相当于投入5 372 040.308元，产出2 518.647 2个QALY，所产生的ICER 2 132.907 0元/QALY小于1倍人均GDP，其可被判定为枫林社区的OP防治干预方案完全具有成本效果，也即当前的OP防治干预方案是值得投入的。

本研究局限性：在构建模型方面，本研究基于现有数据初步构建了本社区OP防治人群的Markov模型，但因初次骨密度值的4 293例研究对象中，仅362例有第二次骨密度测量数据，虽然在模拟自然状态下持续接受随访管理人群的防治效果时，有一定代表性和参考价值，但较少的样本量无疑影响了模拟全社区人群的精度。在计算投入成本方面，目前是将2016—2022年间所发生的所有成本投入总量进行测算，而未分解到每年，同时，由于4 293例研究对象是从2016年5月开始接受社区防治干预，为确保对成本的充分估计，研究根据5月后的数据计算出了平均年度投入，导致了数据的不准确。在进行效用分析方面，由于缺乏连续的数据，首先是仅计算到第二次诊断的状态，而不是与投入测算保持相同的时间段，因此在最终的处理上研究将投入和产出分别以年均来计量。按照目前的改善趋势，研究所用算法可能过于高估了产出部分，也即如果测算到2022年的人群健康效用值，按照2017年的趋势，枫林社区OP防治人群的QALY或许不会达到文中呈现的结果。其次，第二次诊断时年龄以第一次诊断的年龄情况进行计算，没有考虑人群在治疗期间的QALY可能发生的变化，即使是对第二次诊断值的测算，本研究所用到的第二次健康效用值也不全是基于健康状态测量值，有3 931例研究对象的健康状态使用了有第二次骨密度测量值的362例研究对象建立的Markov模型推算而来，一定程度上与真实世界存在着误差。未来更为合适的做法是做好患者的数据管理并基于完全相同的时间段来测算当年的投入和产出比，使得测算结果更加精确，更利于指导OP防治的后续决策。

4 小结

虽然通过构建Markov模型，引入QALY初步得到了社区开展OP防治的投入产出比结果。虽然本研究数据精度还有待进一步提升，数据管理还有待加强，但可以得出结论：对于社区OP防治，值得投入成本进行长期的推进和开展。

作者贡献：苏瑾提出研究目标，負责研究的构思与设计，研究的实施，论文起草与撰写；寿涓负责文章的质量与审校，对文章整体负责；顾文钦、易春涛负责监督管理骨质疏松防治工作的开展；徐莉苹、程莉莉、周鹏负责文献收集与整理；丁宏娟、卫洋洋负责数据统计分析，绘制图表；吴颖华、秦杰、薛斌、魏百川、王谦、彭燕、程毅敏、杨蓝负责研究的实施，数据收集与采集；王磊、祁瑨麟、邵迎、蔡立明负责数据清洗。

本文无利益冲突。

参考文献

《原发性骨质疏松症诊疗社区指导原则》编写组. 原发性骨质疏松症社区诊疗指导原则［J］. 中国全科医学，2019，22（10）：1125-1132. DOI：10.12114/j.issn.1007-9572.2019.00.116.

胡军，张华，牟青. 骨质疏松症的流行病学趋势与防治进展［J］. 临床荟萃，2011，26（8）：729-731.

刘迅，凌莉，王成，等. 多状态Markov模型在慢性肾脏病分级预后研究中的应用［J］. 第二军医大学学报，2009，30（7）：804-807. DOI：10.3724/SP.J.1008.2009.00804.

蒋新军，江华，李明子，等. 糖尿病患者教育成本效果分析［J］. 中国糖尿病杂志，2018，26（11）：964-968. DOI：10.3969/j.issn.1006-6187.2018.11.018.

中华医学会骨质疏松和骨矿盐疾病分会. 原发性骨质疏松症诊疗指南（2017）［J］. 中国骨质疏松杂志，2019，25（3）：281-309. DOI：10.3969/j.issn.1006-7108.2019.03.001.

周鹏，杨蓝，薛斌，等. 基于双向转诊模式构建社区骨质疏松症防治路径及其意义研究［J］. 中国全科医学，2020，23（5）：585-592. DOI：10.12114/j.issn.1007-9572.2020.00.043.

李戈，张春林，蔡承骅，等. 骨质疏松症药物及健康教育干预策略的成本-效益研究［J］. 中国骨质疏松杂志，2006，12（3）：310-312，217. DOI：10.3969/j.issn.1006-7108.2006.03.030.

冯鑫，高颖，潘雪梅，等. 阿仑膦酸钠预防绝经后妇女骨质疏松性骨折的成本-效用分析［J］. 中国药房，2017，28（17）：2313-2318. DOI：10.6039/j.issn.1001-0408.2017.17.03.

苏晓清，黄乐松，杨洲，等. 阿仑膦酸钠的成本-效果分析［J］. 中国老年学杂志，2013，33（15）：3625-3628. DOI：10.3969/j.issn.1005-9202.2013.15.035.

IGLESIAS C P，TORGERSON D J，BEARNE A，et al. The cost utility of bisphosphonate treatment in established osteoporosis［J］. QJM，2002，95（5）：305-311. DOI：10.1093/qjmed/95.5.305.

SI L，WINZENBERG T M，JIANG Q，et al. Screening for and treatment of osteoporosis：construction and validation of a state-transition microsimulation cost-effectiveness model［J］. Osteoporos Int，2015，26（5）：1477-1489. DOI：10.1007/s00198-014-2999-4.

杨旦红，张伟东. 上海城郊社区高血压患者质量调整生命年及其影响因素研究［J］. 中国全科医学，2014，17（10）：1158-1160. DOI：10.3969/j.issn.1007-9572.2014.10.019.

徐迅，李凡，朱云霞，等. 社区规范化管理对慢性阻塞性肺疾病患者质量调整生命年的影响研究［J］. 中国全科医学，2013，16（9）：798-801. DOI：10.3969/j.issn.1007-9572.2013.07.026.

卞宏毅，郑英杰，曹佳燕. 2013年上海市塘桥社区居民寿命表及去死因寿命表分析［J］. 健康教育与健康促进，2016，11（1）：52-57. DOI：10.16117/j.cnki.31-1974/r.201601017.

官海静，徐菲，刘国恩. 基于EQ-5D量表的质量调整生命年计算方法探讨［J］. 中国卫生经济，2015，34（10）：5-8. DOI：10.7664/CHE20151001.

王俊. 基于Markov模型的糖尿病患者行为干预经济学评价［D］. 合肥：安徽医科大学，2018.

周挺，官海静，刘国恩，等. 中国普通人群生命质量系统综述［J］. 中国卫生事业管理，2016，33（8）：621-623，630.

BRAZIER J E，GREEN C，KANIS J A，et al. A systematic review of health state utility values for osteoporosis-related conditions［J］. Osteoporos Int，2002，13（10）：768-776. DOI：10.1007/s001980200107.

《中国药物经济学评价指南》课题组，刘国恩，胡善联，等.中国药物经济学评价指南（2011版）［J］.中国药物经济学，2011（3）：6-9，11-48.

张丽芳，刘中元，林春梅，等. 我国社区卫生服务中心专科服务开展现状及发展策略［J］. 中国全科医学，2019，22（16）：1900-1903. DOI：10.12114/j.issn.1007-9572.2019.00.284.

沈福来，蔡雨阳，吴建平，等. 社区医疗服务产品属性分析与产品开发路径探索［J］. 中华全科医师杂志，2017，16（3）：174-178. DOI：10.3760/cma.j.issn.1671-7368.2017.03.002.

（收稿日期：2023-08-03；修回日期：2024-03-26）

（本文编辑：程圣）

基金项目：上海市徐汇区医学科技项目（SHXH202019）；上海市徐汇区高原学科建設（SHXHZDXK202325）；上海市徐汇区科委重大项目（SHXH202203）

引用本文：苏瑾，寿涓，顾文钦，等. 基于Markov链的上海市枫林社区骨质疏松全人群防治工作的成本效果分析［J］. 中国全科医学，2024，27（22）：2789-2796. DOI：10.12114/j.issn.1007-9572.2023.0420. ［www.chinagp.net］

SU J，SHOU J，GU W Q，et al. Cost effectiveness analysis of the prevention and treatment of osteoporosis among the entire population in Fenglin community，Shanghai based on markov chain［J］. Chinese General Practice，2024，27（22）：2789-2796.

? Editorial Office of Chinese General Practice. This is an open access article under the CC BY-NC-ND 4.0 license.

*通信作者：寿涓，主任医师；E-mail：asyura55@126.com