基于隊列模型生存率數(shù)據(jù)的跨年度生存患者干預(yù)方案醫(yī)?；鹬С龅挠嬎惴椒?/h1>
2021-09-22 06:55:13黃元楷陳磊李瑩席曉宇
中國藥房 2021年17期
黃元楷 陳磊 李瑩 席曉宇
中圖分類號 R956 文獻標(biāo)志碼 A 文章編號 1001-0408（2021）17-2097-06
DOI 10.6039/j.issn.1001-0408.2021.17.09
摘 要 目的：構(gòu)建一種能夠較為準(zhǔn)確地反映跨年度生存患者干預(yù)方案醫(yī)?；鹬С龅挠嬎惴椒?，為醫(yī)保預(yù)算影響分析（BIA）研究提供參考。方法：基于隊列模型的生存率數(shù)據(jù)，以各年份中每個周期確診的患者作為一個隊列，計算各隊列患者在研究年份內(nèi)各狀態(tài)下的人均生存周期數(shù)，即平均生存時間;在此基礎(chǔ)上，再根據(jù)各隊列的人數(shù)和各狀態(tài)下的周期人均費用計算研究年份內(nèi)所有隊列的患者產(chǎn)生的總費用。以某癌癥的干預(yù)方案作為示例，按所構(gòu)建的計算方法進行演算，并將計算結(jié)果與目前幾種常見算法結(jié)果進行比較;同時，對所構(gòu)建的算法在特殊情形下的拓展情況提出應(yīng)用建議。結(jié)果與結(jié)論：所構(gòu)建的算法相比目前幾種常見算法，其計算過程更符合真實世界中藥物干預(yù)方案相關(guān)醫(yī)保基金支出的過程，且其能夠靈活適應(yīng)多種特殊情形下的計算需求。該算法可較為準(zhǔn)確地計算特定年份下某一干預(yù)方案的醫(yī)?；鹬С?，在一定程度上解決了原有算法未充分考慮跨年度生存患者或計算過程簡易、粗略而導(dǎo)致醫(yī)?；鹬С鲱A(yù)測不準(zhǔn)確的問題，可為我國醫(yī)保BIA相關(guān)研究提供更準(zhǔn)確的方法選擇。
關(guān)鍵詞 預(yù)算影響分析;跨年度生存;醫(yī)?；鹬С?計算方法;生存率
Calculation Method of Medical Insurance Fund Expenditure of Intervention Scheme for Cross-year Survival Patients Based on Survival Data of Cohort Model
HUANG Yuankai，CHEN Lei，LI Ying，XI Xiaoyu（The Research Center of National Drug Policy&Ecosystem， China Pharmaceutical University， Nanjing 211198， China）
ABSTRACT OBJECTIVE： To construct a calculation method which can accurately reflect the medical insurance fund expenditure of intervention scheme for cross-year survival patients， and to provide reference for the research of medical insurance budget impact analysis （BIA）. METHODS： Based on survival data of cohort model， taking the patients diagnosed in each cycle in each year as a cohort， the number of per capita survival cycle of cohort patients in each state in the study year was calculated， i.e. the average survival time; on this basis， the total cost of patients in all cohorts in the study year was calculated according to the number of people in each cohort and the per capita cost each cycle in each state. Taking the intervention scheme of a cancer as an example， the calculation was carried out by the established algorithm， and the calculation results were compared with the results of several common algorithms; at the same time， the application suggestions were put forward for the expansion of the constructed algorithm in special cases. RESULTS & CONCLUSIONS： Compared with the several common algorithms， the calculation process of the constructed algorithm is more in line with the process of medical insurance fund expenditure related to drug intervention scheme in the real world， and it can flexibly adapt to the calculation needs in a variety of special situations. This algorithm can more accurately calculate the medical insurance fund expenditure of a intervention scheme in a specific year， and to a certain extent solve the problem of inaccurate prediction of medical insurance fund expenditure due to insufficient consideration of cross-year survival patients or simple and rough calculation process. It can provide a more accurate method choice for the research of medical insurance BIA in China.
KEYWORDS Budget impact analysis; Cross-year survival; Medical insurance fund expenditure; Calculation method; Survival rate
藥品的醫(yī)保預(yù)算影響分析（budget impact analysis，BIA）可評價醫(yī)保目錄外藥品進入醫(yī)保目錄對醫(yī)保基金支出的影響，是藥物經(jīng)濟學(xué)評價的重要補充。《2020年國家醫(yī)保藥品目錄調(diào)整工作方案》要求在藥品醫(yī)保準(zhǔn)入工作中，由地方醫(yī)保部門和相關(guān)學(xué)術(shù)團體推薦醫(yī)保管理、藥物經(jīng)濟學(xué)等方面的專家組成評審專家，分為基金測算組和藥物經(jīng)濟學(xué)組，分別從醫(yī)?；鹩绊懞退幤方?jīng)濟學(xué)評價兩方面針對談判藥品和競價藥品提出評估意見[1]，可見藥品醫(yī)保BIA日益受到我國醫(yī)保部門的重視。
在藥品醫(yī)保BIA中，藥品干預(yù)方案目標(biāo)人群的確定是計算醫(yī)保基金支出的重要環(huán)節(jié)，該目標(biāo)人群指的是符合待評價干預(yù)方案的申請醫(yī)保準(zhǔn)入適用情形的所有患者[2]。長期疾病的患者可能在1年之后病情仍未進展而繼續(xù)使用原干預(yù)方案，或生存期跨自然年，即可能存在跨年度生存患者，因此該類疾病的藥品干預(yù)方案的目標(biāo)人群不僅包括研究年份內(nèi)的新確診患者，也包括往年確診且存活至研究年份的患者，而且不同時間段確診的患者在研究年份內(nèi)的平均生存時間很可能存在差異[3-4]。但是，目前關(guān)于跨年度生存患者在某一特定年份醫(yī)?；鹬С龅挠嬎惴椒ㄉ形唇y(tǒng)一[5-7]，缺乏公認(rèn)的計算方法，且有部分研究甚至未考慮跨年度生存患者所產(chǎn)生的醫(yī)?；鹬С銮闆r[8-10]。
《中國藥物經(jīng)濟學(xué)評價指南（2020中英雙語版）》（以下簡稱“《指南》”）為藥品醫(yī)保BIA的框架構(gòu)建提供了參考，但對于包括上述問題在內(nèi)的BIA具體算法未提出詳細指導(dǎo)意見[2]。因此，從理論角度而言，確定一種能較為準(zhǔn)確反映跨年度生存患者在研究年份內(nèi)醫(yī)保基金支出的計算方法對《指南》的完善具有一定意義。就實踐角度而言，考慮跨年度生存患者產(chǎn)生的醫(yī)?；鹬С瞿茌^為真實地反映干預(yù)措施的預(yù)期醫(yī)?；鹬С銮闆r，有利于提高醫(yī)?；鹬С鲂?，促進我國醫(yī)療衛(wèi)生資源的合理配置。
如果目標(biāo)人群存在跨年度生存現(xiàn)象，無法準(zhǔn)確計算年份內(nèi)的平均生存時間，則可以考慮引入基于隊列模型或基于個體模型的模擬方法來計算不同市場環(huán)境下目標(biāo)人群大小的值[2]。因此，本研究基于隊列模型，構(gòu)建了一種根據(jù)隊列的生存數(shù)據(jù)計算使用某一干預(yù)方案的跨年度生存人群在特定年份下平均生存時間的方法，再將該平均生存時間和單位時間的干預(yù)方案醫(yī)?；鹬С鱿喑耍梢杂嬎愠鲈撎囟攴輧?nèi)的平均醫(yī)?；鹬С觯缓笤賹⑵涑艘愿深A(yù)方案的使用人數(shù)就可以得到該特定年份內(nèi)干預(yù)方案的醫(yī)保基金總支出;同時，本研究還以某癌癥的干預(yù)方案為例，對本算法的計算過程進行基本演示，可得出本算法適用于存在跨年度生存的長期疾病患者的干預(yù)方案的醫(yī)保基金支出計算。
1 平均生存時間計算原理
本算法與馬爾可夫（Markov）模型、分區(qū)生存模型等藥物經(jīng)濟學(xué)常用模型類似，可以考慮多種生存狀態(tài)，且以循環(huán)周期為單位對狀態(tài)的持續(xù)時間進行計算[11-12]。算法的基本理念是以各年份中每個周期確診的患者作為一個隊列，計算各隊列患者在研究年份內(nèi)各狀態(tài)下的人均生存周期數(shù)，即平均生存時間;在此基礎(chǔ)上，再根據(jù)各隊列的人數(shù)和各狀態(tài)下的周期人均費用計算研究年份內(nèi)所有隊列的患者產(chǎn)生的總費用?？梢?，患者在研究年份內(nèi)各狀態(tài)下的平均生存時間是本算法的關(guān)鍵內(nèi)容。
本算法需要下列資料作為支撐：（1）從研究的干預(yù)方案上市年份至研究年份期間，研究疾病的發(fā)病人數(shù)資料——由于大多數(shù)流行病學(xué)資料僅能提供研究疾病的年發(fā)病率，故本研究將以各年份發(fā)病率為例進行闡述;（2）各年份發(fā)病人數(shù)中使用該干預(yù)方案的比例;（3）上述患者人群的長期隨訪生存率數(shù)據(jù)。
本算法有如下假設(shè)：（1）研究適應(yīng)證下，同一年份內(nèi)不同周期確診的患者隊列人數(shù)相等;（2）研究適應(yīng)證下，確診后使用相同干預(yù)方案的患者在各狀態(tài)下各周期的生存率相同。
本算法的思路如圖1所示。圖中僅考慮生存和死亡2個狀態(tài)，在時間軸上羅列某干預(yù)方案下各周期確診的患者隊列的生存率，并截取某特定年份內(nèi)各周期（圖1中以季度為周期）下各隊列的生存率，將其視為該隊列在該周期下的平均生存周期數(shù)，以此計算各隊列在該年份下的平均生存時間（例如某隊列患者在2021年第1季度的生存率是70%，則視為該隊列在這個季度的平均生存時間是0.7個季度，后續(xù)3個季度以此類推;各季度下的平均生存時間相加即為該隊列在2021年的平均生存時間），并按照各隊列的總?cè)藬?shù)對各隊列在該年份下的平均生存時間做加權(quán)求和計算，求得使用該干預(yù)方案的全部人群在該年份下的總生存期數(shù)（實際運算時，一般會先計算出各年份確診患者隊列在研究年份內(nèi)的平均生存時間，然后用各年份發(fā)病人數(shù)中使用該干預(yù)方案的人數(shù)進行加權(quán)計算，而不是按照隊列計算，因為按年份計算便于理解和編制BIA模型）。多狀態(tài)模型在各狀態(tài)下的平均生存時間計算思路亦同理。在此基礎(chǔ)上，即可計算出使用該干預(yù)方案的人群在研究年份內(nèi)的總生存周期數(shù)，進而計算出該干預(yù)方案在該年份的醫(yī)?；鹂傊С?。
2 平均生存時間和醫(yī)保基金支出的計算方法
筆者現(xiàn)以“患者生存”這一單個狀態(tài)為視角，基于患者的生存數(shù)據(jù)，闡述上述原理的詳細計算方法如下。
2.1 定義基本信息和參數(shù)
本算法的基本設(shè)定是：各年份內(nèi)每個周期的確診患者人數(shù)相等，且每個周期產(chǎn)生1個患者隊列，并單獨計算其生存情況。本算法同時使用2套時間系統(tǒng)，分別為患者的確診時間和研究關(guān)注的時間段，兩者均以周期為單位進行計算。以圖2為例，各年份被平均劃分為4個周期，各列分別為各周期確診的患者隊列的生存情況，各行分別為不同周期確診的患者隊列在實際時間中同一周期的生存情況。定義基線年（即研究的年份）為y0，基線年前一年為y0+1，以此類推從基線年開始向前n年（包括基線年）為yn，研究時間內(nèi)生存隊列最早產(chǎn)生的年份為ymax;pt為患者在干預(yù)方案下第t個周期的生存率，c為每年的周期數(shù)。圖2中，y2為研究時間內(nèi)生存隊列最早產(chǎn)生的年份（ymax），c的取值為4。
2.2 計算某周期確診患者在研究年份內(nèi)的平均生存周期數(shù)
定義年內(nèi)生存周期數(shù)（survival cycles of cohort，csct）為在某一周期確診的患者隊列從確診第t周期開始1年內(nèi)各周期的生存率的數(shù)值之和，其實際含義是接受該干預(yù)方案治療的患者從第t周期開始1年內(nèi)的平均生存周期數(shù)，計算公式為：
csct=pt+pt+1+pt+2+……+pt+c-1…（公式1）
在研究年份的csct中，t的取值為患者隊列確診周期開始到研究年份第1周期為止所經(jīng)過的周期數(shù)。以圖2為例，y2年第2周期確診的患者在y0年第1周期時共生存了8個周期，所以其csct為csc8，其取值為第8、9、10和11周期的生存率的數(shù)值之和。
需要注意的是，為方便計算，csct中的周期數(shù)t可以取負值或0，且該2個取值下pt及其之后數(shù)個周期的生存率可能為0。以圖2右下角為例，csc-1為p1和p2之和，即y0年第3周期確診的患者在y0年的平均生存周期數(shù)在計算中依然視為其從y0年第1周期至第4周期的生存周期數(shù)之和，但因前2個周期中該隊列患者尚未確診，因此這2個周期的生存率視為0。
2.3 計算某年份確診患者在研究年份內(nèi)的平均生存周期數(shù)
定義單年產(chǎn)生隊列年內(nèi)平均生存周期數(shù)（average survival cycles of cohorts generated in one year，ascn）為yn年各周期確診的患者在y0年的csct的均值，其實際含義是在假定yn年各周期確診患者人數(shù)相等的情況下，yn年確診的患者在y0年的平均生存周期數(shù)，計算公式為：
ascn=[∑c×（n-1）+1
t=c×（n-2）+2][csct][c] …（公式2）
公式2中的c可理解為yn年確診的隊列數(shù)量，則公式2的計算思路為yn年確診的各隊列在y0年的生存周期數(shù)之和除以yn年產(chǎn)生的隊列數(shù)量，從而求得各隊列的平均生存周期數(shù)。
以圖2為例，y1年第1周期（周期5）確診的患者在y0年各周期的生存率是患者確診后第5～8周期的生存率，其總和為csc5;y1年第2、3、4周期確診的患者以此類推，分別為csc4、csc3、csc2;這4個隊列的生存周期數(shù)之和即為圖中asc2對應(yīng)的4個隊列的cscn之和，將其除以y1年的確診患者隊列數(shù)4即為asc2。
2.4 計算干預(yù)方案全部使用人群在研究年份內(nèi)的總生存周期數(shù)和年醫(yī)?；鹂傊С?/p>
求得y0至ymax各年份確診患者在y0年的平均生存周期數(shù)asc0、asc0+1……ascmax后，根據(jù)上述各年份使用該干預(yù)方案的人數(shù)對其做加權(quán)求和運算，即可得出使用該干預(yù)方案的人群在研究年份內(nèi)的總生存周期數(shù);再將其與該干預(yù)方案的單周期人均醫(yī)保費用相乘，即可得出該干預(yù)方案的年醫(yī)?；鹂傊С?。
3 本研究算法的實例演算及與目前常見算法的比較
本研究以某癌癥的干預(yù)方案為例，演算上述算法。該實例中，設(shè)定模型周期為1個月，則每年有12個周期。以2021年為研究年份，且該方案藥物于2019年初上市，則生存隊列從2019年1月開始，每月產(chǎn)生隊列，直至2021年12月。假設(shè)使用該干預(yù)方案治療的癌癥患者隊列的各周期無進展生存（progression-free survival，PFS）概率和總生存（overall survival，OS）概率符合Weibull分布，則用上述PFS期和OS期相關(guān)數(shù)據(jù)擬合PFS和OS曲線，并直接以各周期的時間中點作為周期數(shù)，以進行半周期校正，可得各周期PFS和OS狀態(tài)下的生存率表達式分別為：Pt=e-0.06×（t-0.5）1.44、Pt=e-0.02×（t-0.5）1.41-e-0.06×（t-0.5）1.44;而各周期疾病進展（progressive disease，PD）概率為相應(yīng)周期的OS概率減去PFS概率。
由于本實例從2019年1月確診患者隊列開始模擬，而研究年份為2021年，則2019年患者隊列的平均生存周期數(shù)csct相對應(yīng)的t值為2019年1月至2021年1月間的周期數(shù)（即25），其他隊列以此類推，從而可用公式1計算出各隊列在2021年的PFS期和PD期的csct;再據(jù)此計算2019、2020、2021年確診患者在2021年的PFS期和PD期的ascn，結(jié)果如表1所示（表中2020年和2021年的部分計算細節(jié)略去）。
按照上述計算結(jié)果，可根據(jù)該干預(yù)方案的PFS期和PD期的單周期人均醫(yī)?；鹬С鲆约?019年至2021年各年份新確診患者中使用該干預(yù)方案的人數(shù)，計算出2021年該干預(yù)方案產(chǎn)生的醫(yī)?；鹂傊С觥?/p>
本研究還將3種目前該領(lǐng)域的常見算法與本算法進行比較，結(jié)果如表2所示。表2中各算法分別為：（1）本研究構(gòu)建的算法;（2）僅考慮2021年確診患者在2021年產(chǎn)生的費用，按平均終生PFS期和終生PD期計算;（3）僅考慮2021年確診患者產(chǎn)生的費用，按1月確診后12個月（包括1月）的平均PFS期和PD期計算;（4）僅考慮2021年確診患者產(chǎn)生的費用，按7月確診后6個月（包括7月）的平均PFS期和PD期計算。各算法均使用上文計算實例中使用的生存數(shù)據(jù)來計算患者的PFS期和PD期，并假定2019年至2021年研究疾病的發(fā)病人數(shù)分別為10萬人、11萬人和12萬人，以及該干預(yù)方案的PFS期費用為5 000元/月、PD期費用為2 000元/月。
由表2可見，如果醫(yī)保BIA中的研究疾病和干預(yù)方案下存在跨年度生存的患者，目前常見的3種算法對干預(yù)方案在2021年產(chǎn)生的醫(yī)保基金總支出的計算結(jié)果與本研究構(gòu)建的算法結(jié)果比較差異較大，且大多低于本研究算法結(jié)果。從理論上講，相比目前常見算法，本研究算法能夠?qū)IA目標(biāo)人群在研究年份的生存情況進行較為貼近現(xiàn)實的模擬測算，因而能夠為醫(yī)保BIA提供更為準(zhǔn)確的人均生存或用藥時長參數(shù)。
4 討論
本研究基于患者生存率構(gòu)建了一種針對跨年度生存患者的醫(yī)?；鹬С龅挠嬎惴椒ǎ啾饶壳霸擃I(lǐng)域中常見的算法，本算法可較為準(zhǔn)確地計算特定年份下某一干預(yù)方案的醫(yī)?；鹬С觯谝欢ǔ潭壬辖鉀Q了常見算法未充分考慮跨年度生存患者或計算過程簡易、粗略而導(dǎo)致醫(yī)?；鹬С鲱A(yù)測不準(zhǔn)確的問題，可為我國醫(yī)保BIA相關(guān)研究以及《指南》的完善提供參考，為研究者和決策者提供建議，促進我國醫(yī)療衛(wèi)生資源的合理配置。
4.1 本算法關(guān)鍵問題的解釋
首先，本算法在構(gòu)建過程中引入了某周期確診患者在研究年份內(nèi)的平均生存周期數(shù)csct和某年份確診患者在研究年份內(nèi)的平均生存周期數(shù)ascn這2個概念以便于計算，其實際意義理解起來較為晦澀，但其基本含義均是指某一時間段內(nèi)確診的患者在研究年份或研究年份內(nèi)某一周期的平均生存時間。需要注意的是，這與傳統(tǒng)意義上的患者平均生存時間（患者從確診或開始干預(yù)直至死亡的時間跨度）的概念有所差異。
其次，本研究引入了負生存周期數(shù)的概念（公式1），該設(shè)計并無實際含義，僅為方便統(tǒng)一和簡化公式而設(shè)計，因而無需從現(xiàn)實意義角度進行理解。
第三，本算法可基于多狀態(tài)模型計算，其計算思路與上文中所述單個狀態(tài)的計算思路完全一致。例如在“3”項下實例中，本研究以存在PFS和PD 2個健康狀態(tài)的癌癥為例計算了各狀態(tài)下的患者平均生存時間。
4.2 本算法的使用條件和關(guān)鍵設(shè)置
因為本算法是基于患者生存率構(gòu)建的，因此使用本算法的基本條件是可以獲得符合研究需求的生存數(shù)據(jù);同時，為計算研究年份內(nèi)的醫(yī)?；鹬С觯枰@取干預(yù)方案在各年份新確診患者中的使用人數(shù)，該參數(shù)可通過專家咨詢、臨床調(diào)研等方式獲得。
在將本算法應(yīng)用到醫(yī)保BIA的過程中，有2個關(guān)鍵參數(shù)需要進行設(shè)置：1）研究范圍內(nèi)生存隊列的最早產(chǎn)生年份ymax，可以通過設(shè)置模型模擬結(jié)束的生存率閾值來確定，且要注意該值必須為整數(shù);2）每年的周期數(shù)，可以根據(jù)疾病特征來設(shè)計，但為方便計算，須盡量使每年周期數(shù)為整數(shù)。
4.3 特殊情形下本算法的拓展
由于實際市場和臨床情況比較復(fù)雜，上述算法需根據(jù)實際用藥情況靈活改變，以下列舉了本算法在一些特殊情形下的拓展情況。
4.3.1 新藥上市或適應(yīng)證獲批時間晚于理論上的ymax 對于一些新藥或獲批新適應(yīng)證的藥品，在進行醫(yī)保BIA過程中，若按照生存率閾值確定ymax，則可能出現(xiàn)新藥上市年份或新適應(yīng)證獲批年份晚于ymax的情況，此時計算所得的醫(yī)保基金支出有悖于實際情況。
因此，對于新藥或新獲批適應(yīng)證的藥品，ymax的設(shè)計應(yīng)考慮藥品上市或適應(yīng)證獲批的年份——若新藥上市或新適應(yīng)證獲批的年份早于理論上的ymax，則按照生存率閾值確定ymax;若新藥上市或新適應(yīng)證獲批的年份晚于理論上的ymax，則以新藥上市或新適應(yīng)證獲批年份作為ymax。在實際中，新藥上市時間或適應(yīng)證獲批時間可能在年末，導(dǎo)致該年的銷售量較低，此時可根據(jù)藥品銷售情況調(diào)整算式中的起始年份為藥品上市或適應(yīng)證獲批年份的后一年，或僅將獲批年份下產(chǎn)生的隊列作為該年度的全部患者隊列來計算ascn。
4.3.2 治療過程中僅接受干預(yù)方案的一次性治療 在部分研究中，可能會出現(xiàn)患者僅接受該研究關(guān)注的干預(yù)方案的一次性治療，后續(xù)治療周期中接受其他干預(yù)方案治療的情況[13]。該情況下，由于對醫(yī)保基金產(chǎn)生影響的變動僅發(fā)生在研究關(guān)注的年份內(nèi)，故ymax可直接設(shè)置為研究關(guān)注的多個年份中最早的一年。例如癌癥放療方案的BIA中，一般假設(shè)干預(yù)方案本身的費用在患者確診周期已全部發(fā)生，但后續(xù)的患者生存情況會因選用不同的放療技術(shù)而存在差異。在這樣的情況下，實際上干預(yù)方案調(diào)入或調(diào)出報銷范圍對該方案選擇的影響僅出現(xiàn)在報銷范圍調(diào)整生效后的年份，因此僅需計算報銷范圍調(diào)整后年份的人均生存或用藥時間。
4.3.3 實際用藥時長與理論時長有差異 這類情形主要分為2種：①患者中止治療。本算法中所使用的生存數(shù)據(jù)來源于臨床試驗和已公開發(fā)表的文獻，研究對象在用藥時應(yīng)處于嚴(yán)格的臨床試驗條件下或其用藥依從性高，其生存數(shù)據(jù)反映了患者用藥后理想的生存情況。但在實際治療過程中，用藥依從性較差、藥物可及性較低等原因可導(dǎo)致患者突然停止治療[14]，此時若仍使用文獻報道的生存數(shù)據(jù)可導(dǎo)致生存率虛高。在該情況下，應(yīng)考慮將ymax設(shè)置為絕大部分患者的最長用藥時長，例如發(fā)現(xiàn)研究對象中95%的患者都會在2.5年內(nèi)因依從性差而停止用藥，但理論上的ymax為5年，則ymax實際應(yīng)設(shè)置為3年。此處建議不以患者的平均用藥時長來估算ymax，否則會導(dǎo)致部分患者在超過該時長后的用藥被忽略，繼而低估醫(yī)?；鹬С觥＂诨颊咄Ｖ褂盟幒?，該研究所關(guān)注的干預(yù)方案的藥效仍持續(xù)了一段時間。生存率反映的是干預(yù)方案的臨床治療效果，但某些長期疾病的患者在停止用藥后，之前服用的藥物仍可維持病情的穩(wěn)定，此時基于生存率計算而得到的平均生存時間并不等于實際用藥時間，導(dǎo)致醫(yī)?；鹬С鲱A(yù)測不準(zhǔn)確。在該情況下，可設(shè)置用藥時長閾值，或根據(jù)實際用藥時長設(shè)置生存率閾值，使計算而得的患者隊列的平均用藥時長接近實際用藥時長。
4.3.4 年份內(nèi)各周期確診人數(shù)存在明顯差異 在常見的流行病學(xué)資料中，研究者往往僅能獲得疾病的年發(fā)病率數(shù)據(jù)，而沒有各周期的發(fā)病率數(shù)據(jù)。本研究為貼近實際，且為方便解釋計算原理，故假設(shè)一年中各周期的發(fā)病率相等。但在實際研究過程中，研究者也許能夠獲得或計算出更為詳盡的發(fā)病率數(shù)據(jù)，這種情況下，可以為模型中每個隊列單獨賦予確診人數(shù)，并單獨計算每個隊列的生存時間，以更充分地利用高質(zhì)量的發(fā)病率數(shù)據(jù)。
4.3.5 新確診患者和跨年度生存患者的費用結(jié)構(gòu)存在差異 部分研究中，新確診患者與跨年度生存患者的年均總費用構(gòu)成可能存在差異。例如放療、手術(shù)等干預(yù)方案成本在BIA中一般僅會被納入新確診患者的費用構(gòu)成中。對此，可以分開計算本年度新確診患者隊列和跨年度生存患者隊列的年均總費用，為本年度新確診患者單獨增加成本項目。
4.4 局限性
本研究構(gòu)建的算法存在一定局限性：第一，本研究構(gòu)建的算法需要使用該干預(yù)方案的患者的生存數(shù)據(jù)和各年份使用該干預(yù)方案的患者人數(shù)，如果這2項數(shù)據(jù)無法獲取，則無法使用本算法;第二，由于實際中難以追溯年份內(nèi)不同周期確診的患者人數(shù)，而本算法假設(shè)不同周期確診的患者人數(shù)相等，因此計算得到的平均生存率可能與實際存在一定偏差;第三，由于醫(yī)保BIA關(guān)注的實際市場情況和臨床情況復(fù)雜，而本算法尚未經(jīng)過大量研究的應(yīng)用和實證，因此本研究對一部分特殊情形下算法的拓展提出了建議，但在研究一些具體BIA問題時，本算法的適用性仍可能較低，有待繼續(xù)改進或開發(fā)更為合理的計算方法。
5 結(jié)語
長期疾病藥物干預(yù)方案的醫(yī)保BIA中，目前常見的計算方法常未考慮跨年度生存的患者群體，或計算方法較為粗略，對醫(yī)保基金支出估算的準(zhǔn)確性產(chǎn)生了影響。本研究以完善醫(yī)保BIA建模思路為目的，基于患者生存率構(gòu)建了一種存在跨年度生存患者情形下的醫(yī)?；鹬С鲇嬎惴椒?。本方法的基本思路是計算研究年份及其以往各年份內(nèi)使用某干預(yù)方案（如藥物治療方案）的確診患者在研究年份內(nèi)的平均生存周期數(shù)，再根據(jù)各年份的確診患者中使用該干預(yù)方案的人數(shù)和人均周期成本計算該干預(yù)方案在研究年份的總醫(yī)?；鹬С?。本算法能夠在一定程度上解決目前常見計算方法中對跨年度生存情形下的醫(yī)保基金支出預(yù)測不準(zhǔn)確的問題。但由于實際市場情況和臨床情況復(fù)雜，本算法需根據(jù)實際用藥情況靈活改變。
參考文獻
[ 1 ] 國家醫(yī)療保障局. 2020年國家醫(yī)保藥品目錄調(diào)整工作方案[EB/OL].（2020-08-17）[2021-08-06].http：//www.nhsa.gov.cn/art/2020/8/17/art_37_3407.html.
[ 2 ] 劉國恩.中國藥物經(jīng)濟學(xué)評價指南：2020中英雙語版[M]. 北京：中國市場出版社，2020：91-92.
[ 3 ] WU Y L，CHU D T，HAN B，et al. Phase Ⅲ，randomized，open-label，first-line study in Asia of gefitinib versus carboplatin/paclitaxel in clinically selected patients with advanced non-small-cell lung cancer：evaluation of patients recruited from mainland China[J]. Asia Pac J Clin Oncol，2012，8（3）：232-243.
[ 4 ] MOVAHEDI M，HEPWORTH E，Mirza R，et al. Discontinuation of biologic therapy due to lack/loss of response and adverse events is similar between TNFi and non-TNFi class：results from a real-world rheumatoid arthritis cohort[J]. Semin Arthritis Rheum，2020，50（5）：915-922.
[ 5 ] 謝麗，夏茹意，莊貴華.伴代償期肝硬化乙肝相關(guān)肝癌患者根治性切除術(shù)后抗病毒治療的預(yù)算影響分析[J].中國衛(wèi)生經(jīng)濟，2019，38（9）：64-68.
[ 6 ] 殷瀟，方欣，胡江藺，等.舒尼替尼治療轉(zhuǎn)移性腎細胞癌的預(yù)算影響分析[J].中國衛(wèi)生經(jīng)濟，2016，35（12）：67-69.
[ 7 ] 張欲曉，殷瀟，田夢媛，等.赫賽汀治療HER-2陽性乳腺癌的醫(yī)保預(yù)算影響分析[J].中國衛(wèi)生經(jīng)濟，2016，35（12）：63-66.
[ 8 ] 陳斌斌，范長生.注射用紫杉醇（白蛋白結(jié)合型）治療晚期乳腺癌的預(yù)算影響分析[J].中國醫(yī)療保險，2016（12）：60-63.
[ 9 ] 范長生，何敏媚，吳久鴻.復(fù)方黃黛片治療急性早幼粒細胞白血病預(yù)算影響分析[J].中國新藥雜志，2017，26（10）：1219-1224.
[10] 張強，崔穎，馬銳，等.規(guī)范使用唑來膦酸治療國內(nèi)晚期肺癌或乳腺癌合并骨轉(zhuǎn)移預(yù)算影響分析[J].中國醫(yī)療保險，2018（12）：59-67.
[11] 張潔，黃泰康. 藥物經(jīng)濟學(xué)評價中馬爾可夫模型的優(yōu)化：馬爾可夫決策模型[J].中國藥房，2008，19（17）：1284- 1286.
[12] 邵榮杰，唐文熙，馬愛霞.分區(qū)生存模型在藥物經(jīng)濟學(xué)評價中的應(yīng)用[J].中國衛(wèi)生經(jīng)濟，2019，38（9）：60-63.
[13] 汪紅艷，孔令玲，王凡，等.老年食管癌放射治療療效及預(yù)后影響因素分析[J].安徽醫(yī)科大學(xué)學(xué)報，2016，51（8）：1188-1192.
[14] 洪梅花，徐東曉，童素梅，等.急性心肌梗死支架術(shù)后患者不同時期服藥依從性及其影響因素的調(diào)查與分析[J].中國護理管理，2014，14（12）：1261-1265.
（收稿日期：2021-01-27 修回日期：2021-08-06）
（編輯：胡曉霖）
猜你喜歡
浮力計算方法匯集中學(xué)生數(shù)理化·八年級物理人教版(2022年5期)2022-06-05 06:57:38 “五年生存率”不等于只能活五年祝您健康·文摘版(2020年7期)2020-07-13 09:35:19 人工智能助力卵巢癌生存率預(yù)測中國生殖健康(2019年2期)2019-08-23 08:11:52 日本首次公布本國居民癌癥三年生存率世界知識(2018年19期)2018-11-21 03:01:12 “五年生存率”≠只能活五年老年教育(老年大學(xué))(2018年12期)2018-01-28 14:15:50 隨機振動試驗包絡(luò)計算方法現(xiàn)代工業(yè)經(jīng)濟和信息化(2016年2期)2016-05-17 05:34:16 HER2 表達強度對三陰性乳腺癌無病生存率的影響癌癥進展(2016年12期)2016-03-20 13:16:14 基于能量的球壓痕硬度的計算方法材料科學(xué)與工程學(xué)報(2016年5期)2016-02-27 07:11:22 不同應(yīng)變率比值計算方法在甲狀腺惡性腫瘤診斷中的應(yīng)用中國醫(yī)學(xué)影像學(xué)雜志(2015年9期)2015-12-15 11:03:26 中藥治療對三陰性乳腺癌患者生存率的影響西部中醫(yī)藥(2015年9期)2015-02-02 01:19:24
黃元楷 陳磊 李瑩 席曉宇
中圖分類號 R956 文獻標(biāo)志碼 A 文章編號 1001-0408（2021）17-2097-06
DOI 10.6039/j.issn.1001-0408.2021.17.09
摘 要 目的：構(gòu)建一種能夠較為準(zhǔn)確地反映跨年度生存患者干預(yù)方案醫(yī)?；鹬С龅挠嬎惴椒?，為醫(yī)保預(yù)算影響分析（BIA）研究提供參考。方法：基于隊列模型的生存率數(shù)據(jù)，以各年份中每個周期確診的患者作為一個隊列，計算各隊列患者在研究年份內(nèi)各狀態(tài)下的人均生存周期數(shù)，即平均生存時間;在此基礎(chǔ)上，再根據(jù)各隊列的人數(shù)和各狀態(tài)下的周期人均費用計算研究年份內(nèi)所有隊列的患者產(chǎn)生的總費用。以某癌癥的干預(yù)方案作為示例，按所構(gòu)建的計算方法進行演算，并將計算結(jié)果與目前幾種常見算法結(jié)果進行比較;同時，對所構(gòu)建的算法在特殊情形下的拓展情況提出應(yīng)用建議。結(jié)果與結(jié)論：所構(gòu)建的算法相比目前幾種常見算法，其計算過程更符合真實世界中藥物干預(yù)方案相關(guān)醫(yī)保基金支出的過程，且其能夠靈活適應(yīng)多種特殊情形下的計算需求。該算法可較為準(zhǔn)確地計算特定年份下某一干預(yù)方案的醫(yī)?；鹬С?，在一定程度上解決了原有算法未充分考慮跨年度生存患者或計算過程簡易、粗略而導(dǎo)致醫(yī)?；鹬С鲱A(yù)測不準(zhǔn)確的問題，可為我國醫(yī)保BIA相關(guān)研究提供更準(zhǔn)確的方法選擇。
關(guān)鍵詞 預(yù)算影響分析;跨年度生存;醫(yī)?；鹬С?計算方法;生存率
Calculation Method of Medical Insurance Fund Expenditure of Intervention Scheme for Cross-year Survival Patients Based on Survival Data of Cohort Model
HUANG Yuankai，CHEN Lei，LI Ying，XI Xiaoyu（The Research Center of National Drug Policy&Ecosystem， China Pharmaceutical University， Nanjing 211198， China）
ABSTRACT OBJECTIVE： To construct a calculation method which can accurately reflect the medical insurance fund expenditure of intervention scheme for cross-year survival patients， and to provide reference for the research of medical insurance budget impact analysis （BIA）. METHODS： Based on survival data of cohort model， taking the patients diagnosed in each cycle in each year as a cohort， the number of per capita survival cycle of cohort patients in each state in the study year was calculated， i.e. the average survival time; on this basis， the total cost of patients in all cohorts in the study year was calculated according to the number of people in each cohort and the per capita cost each cycle in each state. Taking the intervention scheme of a cancer as an example， the calculation was carried out by the established algorithm， and the calculation results were compared with the results of several common algorithms; at the same time， the application suggestions were put forward for the expansion of the constructed algorithm in special cases. RESULTS & CONCLUSIONS： Compared with the several common algorithms， the calculation process of the constructed algorithm is more in line with the process of medical insurance fund expenditure related to drug intervention scheme in the real world， and it can flexibly adapt to the calculation needs in a variety of special situations. This algorithm can more accurately calculate the medical insurance fund expenditure of a intervention scheme in a specific year， and to a certain extent solve the problem of inaccurate prediction of medical insurance fund expenditure due to insufficient consideration of cross-year survival patients or simple and rough calculation process. It can provide a more accurate method choice for the research of medical insurance BIA in China.
KEYWORDS Budget impact analysis; Cross-year survival; Medical insurance fund expenditure; Calculation method; Survival rate
藥品的醫(yī)保預(yù)算影響分析（budget impact analysis，BIA）可評價醫(yī)保目錄外藥品進入醫(yī)保目錄對醫(yī)保基金支出的影響，是藥物經(jīng)濟學(xué)評價的重要補充。《2020年國家醫(yī)保藥品目錄調(diào)整工作方案》要求在藥品醫(yī)保準(zhǔn)入工作中，由地方醫(yī)保部門和相關(guān)學(xué)術(shù)團體推薦醫(yī)保管理、藥物經(jīng)濟學(xué)等方面的專家組成評審專家，分為基金測算組和藥物經(jīng)濟學(xué)組，分別從醫(yī)?；鹩绊懞退幤方?jīng)濟學(xué)評價兩方面針對談判藥品和競價藥品提出評估意見[1]，可見藥品醫(yī)保BIA日益受到我國醫(yī)保部門的重視。
在藥品醫(yī)保BIA中，藥品干預(yù)方案目標(biāo)人群的確定是計算醫(yī)保基金支出的重要環(huán)節(jié)，該目標(biāo)人群指的是符合待評價干預(yù)方案的申請醫(yī)保準(zhǔn)入適用情形的所有患者[2]。長期疾病的患者可能在1年之后病情仍未進展而繼續(xù)使用原干預(yù)方案，或生存期跨自然年，即可能存在跨年度生存患者，因此該類疾病的藥品干預(yù)方案的目標(biāo)人群不僅包括研究年份內(nèi)的新確診患者，也包括往年確診且存活至研究年份的患者，而且不同時間段確診的患者在研究年份內(nèi)的平均生存時間很可能存在差異[3-4]。但是，目前關(guān)于跨年度生存患者在某一特定年份醫(yī)?；鹬С龅挠嬎惴椒ㄉ形唇y(tǒng)一[5-7]，缺乏公認(rèn)的計算方法，且有部分研究甚至未考慮跨年度生存患者所產(chǎn)生的醫(yī)?；鹬С銮闆r[8-10]。
《中國藥物經(jīng)濟學(xué)評價指南（2020中英雙語版）》（以下簡稱“《指南》”）為藥品醫(yī)保BIA的框架構(gòu)建提供了參考，但對于包括上述問題在內(nèi)的BIA具體算法未提出詳細指導(dǎo)意見[2]。因此，從理論角度而言，確定一種能較為準(zhǔn)確反映跨年度生存患者在研究年份內(nèi)醫(yī)保基金支出的計算方法對《指南》的完善具有一定意義。就實踐角度而言，考慮跨年度生存患者產(chǎn)生的醫(yī)?；鹬С瞿茌^為真實地反映干預(yù)措施的預(yù)期醫(yī)?；鹬С銮闆r，有利于提高醫(yī)?；鹬С鲂?，促進我國醫(yī)療衛(wèi)生資源的合理配置。
如果目標(biāo)人群存在跨年度生存現(xiàn)象，無法準(zhǔn)確計算年份內(nèi)的平均生存時間，則可以考慮引入基于隊列模型或基于個體模型的模擬方法來計算不同市場環(huán)境下目標(biāo)人群大小的值[2]。因此，本研究基于隊列模型，構(gòu)建了一種根據(jù)隊列的生存數(shù)據(jù)計算使用某一干預(yù)方案的跨年度生存人群在特定年份下平均生存時間的方法，再將該平均生存時間和單位時間的干預(yù)方案醫(yī)?；鹬С鱿喑耍梢杂嬎愠鲈撎囟攴輧?nèi)的平均醫(yī)?；鹬С觯缓笤賹⑵涑艘愿深A(yù)方案的使用人數(shù)就可以得到該特定年份內(nèi)干預(yù)方案的醫(yī)保基金總支出;同時，本研究還以某癌癥的干預(yù)方案為例，對本算法的計算過程進行基本演示，可得出本算法適用于存在跨年度生存的長期疾病患者的干預(yù)方案的醫(yī)保基金支出計算。
1 平均生存時間計算原理
本算法與馬爾可夫（Markov）模型、分區(qū)生存模型等藥物經(jīng)濟學(xué)常用模型類似，可以考慮多種生存狀態(tài)，且以循環(huán)周期為單位對狀態(tài)的持續(xù)時間進行計算[11-12]。算法的基本理念是以各年份中每個周期確診的患者作為一個隊列，計算各隊列患者在研究年份內(nèi)各狀態(tài)下的人均生存周期數(shù)，即平均生存時間;在此基礎(chǔ)上，再根據(jù)各隊列的人數(shù)和各狀態(tài)下的周期人均費用計算研究年份內(nèi)所有隊列的患者產(chǎn)生的總費用?？梢?，患者在研究年份內(nèi)各狀態(tài)下的平均生存時間是本算法的關(guān)鍵內(nèi)容。
本算法需要下列資料作為支撐：（1）從研究的干預(yù)方案上市年份至研究年份期間，研究疾病的發(fā)病人數(shù)資料——由于大多數(shù)流行病學(xué)資料僅能提供研究疾病的年發(fā)病率，故本研究將以各年份發(fā)病率為例進行闡述;（2）各年份發(fā)病人數(shù)中使用該干預(yù)方案的比例;（3）上述患者人群的長期隨訪生存率數(shù)據(jù)。
本算法有如下假設(shè)：（1）研究適應(yīng)證下，同一年份內(nèi)不同周期確診的患者隊列人數(shù)相等;（2）研究適應(yīng)證下，確診后使用相同干預(yù)方案的患者在各狀態(tài)下各周期的生存率相同。
本算法的思路如圖1所示。圖中僅考慮生存和死亡2個狀態(tài)，在時間軸上羅列某干預(yù)方案下各周期確診的患者隊列的生存率，并截取某特定年份內(nèi)各周期（圖1中以季度為周期）下各隊列的生存率，將其視為該隊列在該周期下的平均生存周期數(shù)，以此計算各隊列在該年份下的平均生存時間（例如某隊列患者在2021年第1季度的生存率是70%，則視為該隊列在這個季度的平均生存時間是0.7個季度，后續(xù)3個季度以此類推;各季度下的平均生存時間相加即為該隊列在2021年的平均生存時間），并按照各隊列的總?cè)藬?shù)對各隊列在該年份下的平均生存時間做加權(quán)求和計算，求得使用該干預(yù)方案的全部人群在該年份下的總生存期數(shù)（實際運算時，一般會先計算出各年份確診患者隊列在研究年份內(nèi)的平均生存時間，然后用各年份發(fā)病人數(shù)中使用該干預(yù)方案的人數(shù)進行加權(quán)計算，而不是按照隊列計算，因為按年份計算便于理解和編制BIA模型）。多狀態(tài)模型在各狀態(tài)下的平均生存時間計算思路亦同理。在此基礎(chǔ)上，即可計算出使用該干預(yù)方案的人群在研究年份內(nèi)的總生存周期數(shù)，進而計算出該干預(yù)方案在該年份的醫(yī)?；鹂傊С?。
2 平均生存時間和醫(yī)保基金支出的計算方法
筆者現(xiàn)以“患者生存”這一單個狀態(tài)為視角，基于患者的生存數(shù)據(jù)，闡述上述原理的詳細計算方法如下。
2.1 定義基本信息和參數(shù)
本算法的基本設(shè)定是：各年份內(nèi)每個周期的確診患者人數(shù)相等，且每個周期產(chǎn)生1個患者隊列，并單獨計算其生存情況。本算法同時使用2套時間系統(tǒng)，分別為患者的確診時間和研究關(guān)注的時間段，兩者均以周期為單位進行計算。以圖2為例，各年份被平均劃分為4個周期，各列分別為各周期確診的患者隊列的生存情況，各行分別為不同周期確診的患者隊列在實際時間中同一周期的生存情況。定義基線年（即研究的年份）為y0，基線年前一年為y0+1，以此類推從基線年開始向前n年（包括基線年）為yn，研究時間內(nèi)生存隊列最早產(chǎn)生的年份為ymax;pt為患者在干預(yù)方案下第t個周期的生存率，c為每年的周期數(shù)。圖2中，y2為研究時間內(nèi)生存隊列最早產(chǎn)生的年份（ymax），c的取值為4。
2.2 計算某周期確診患者在研究年份內(nèi)的平均生存周期數(shù)
定義年內(nèi)生存周期數(shù)（survival cycles of cohort，csct）為在某一周期確診的患者隊列從確診第t周期開始1年內(nèi)各周期的生存率的數(shù)值之和，其實際含義是接受該干預(yù)方案治療的患者從第t周期開始1年內(nèi)的平均生存周期數(shù)，計算公式為：
csct=pt+pt+1+pt+2+……+pt+c-1…（公式1）
在研究年份的csct中，t的取值為患者隊列確診周期開始到研究年份第1周期為止所經(jīng)過的周期數(shù)。以圖2為例，y2年第2周期確診的患者在y0年第1周期時共生存了8個周期，所以其csct為csc8，其取值為第8、9、10和11周期的生存率的數(shù)值之和。
需要注意的是，為方便計算，csct中的周期數(shù)t可以取負值或0，且該2個取值下pt及其之后數(shù)個周期的生存率可能為0。以圖2右下角為例，csc-1為p1和p2之和，即y0年第3周期確診的患者在y0年的平均生存周期數(shù)在計算中依然視為其從y0年第1周期至第4周期的生存周期數(shù)之和，但因前2個周期中該隊列患者尚未確診，因此這2個周期的生存率視為0。
2.3 計算某年份確診患者在研究年份內(nèi)的平均生存周期數(shù)
定義單年產(chǎn)生隊列年內(nèi)平均生存周期數(shù)（average survival cycles of cohorts generated in one year，ascn）為yn年各周期確診的患者在y0年的csct的均值，其實際含義是在假定yn年各周期確診患者人數(shù)相等的情況下，yn年確診的患者在y0年的平均生存周期數(shù)，計算公式為：
ascn=[∑c×（n-1）+1
t=c×（n-2）+2][csct][c] …（公式2）
公式2中的c可理解為yn年確診的隊列數(shù)量，則公式2的計算思路為yn年確診的各隊列在y0年的生存周期數(shù)之和除以yn年產(chǎn)生的隊列數(shù)量，從而求得各隊列的平均生存周期數(shù)。
以圖2為例，y1年第1周期（周期5）確診的患者在y0年各周期的生存率是患者確診后第5～8周期的生存率，其總和為csc5;y1年第2、3、4周期確診的患者以此類推，分別為csc4、csc3、csc2;這4個隊列的生存周期數(shù)之和即為圖中asc2對應(yīng)的4個隊列的cscn之和，將其除以y1年的確診患者隊列數(shù)4即為asc2。
2.4 計算干預(yù)方案全部使用人群在研究年份內(nèi)的總生存周期數(shù)和年醫(yī)?；鹂傊С?/p>
求得y0至ymax各年份確診患者在y0年的平均生存周期數(shù)asc0、asc0+1……ascmax后，根據(jù)上述各年份使用該干預(yù)方案的人數(shù)對其做加權(quán)求和運算，即可得出使用該干預(yù)方案的人群在研究年份內(nèi)的總生存周期數(shù);再將其與該干預(yù)方案的單周期人均醫(yī)保費用相乘，即可得出該干預(yù)方案的年醫(yī)?；鹂傊С?。
3 本研究算法的實例演算及與目前常見算法的比較
本研究以某癌癥的干預(yù)方案為例，演算上述算法。該實例中，設(shè)定模型周期為1個月，則每年有12個周期。以2021年為研究年份，且該方案藥物于2019年初上市，則生存隊列從2019年1月開始，每月產(chǎn)生隊列，直至2021年12月。假設(shè)使用該干預(yù)方案治療的癌癥患者隊列的各周期無進展生存（progression-free survival，PFS）概率和總生存（overall survival，OS）概率符合Weibull分布，則用上述PFS期和OS期相關(guān)數(shù)據(jù)擬合PFS和OS曲線，并直接以各周期的時間中點作為周期數(shù)，以進行半周期校正，可得各周期PFS和OS狀態(tài)下的生存率表達式分別為：Pt=e-0.06×（t-0.5）1.44、Pt=e-0.02×（t-0.5）1.41-e-0.06×（t-0.5）1.44;而各周期疾病進展（progressive disease，PD）概率為相應(yīng)周期的OS概率減去PFS概率。
由于本實例從2019年1月確診患者隊列開始模擬，而研究年份為2021年，則2019年患者隊列的平均生存周期數(shù)csct相對應(yīng)的t值為2019年1月至2021年1月間的周期數(shù)（即25），其他隊列以此類推，從而可用公式1計算出各隊列在2021年的PFS期和PD期的csct;再據(jù)此計算2019、2020、2021年確診患者在2021年的PFS期和PD期的ascn，結(jié)果如表1所示（表中2020年和2021年的部分計算細節(jié)略去）。
按照上述計算結(jié)果，可根據(jù)該干預(yù)方案的PFS期和PD期的單周期人均醫(yī)?；鹬С鲆约?019年至2021年各年份新確診患者中使用該干預(yù)方案的人數(shù)，計算出2021年該干預(yù)方案產(chǎn)生的醫(yī)?；鹂傊С觥?/p>
本研究還將3種目前該領(lǐng)域的常見算法與本算法進行比較，結(jié)果如表2所示。表2中各算法分別為：（1）本研究構(gòu)建的算法;（2）僅考慮2021年確診患者在2021年產(chǎn)生的費用，按平均終生PFS期和終生PD期計算;（3）僅考慮2021年確診患者產(chǎn)生的費用，按1月確診后12個月（包括1月）的平均PFS期和PD期計算;（4）僅考慮2021年確診患者產(chǎn)生的費用，按7月確診后6個月（包括7月）的平均PFS期和PD期計算。各算法均使用上文計算實例中使用的生存數(shù)據(jù)來計算患者的PFS期和PD期，并假定2019年至2021年研究疾病的發(fā)病人數(shù)分別為10萬人、11萬人和12萬人，以及該干預(yù)方案的PFS期費用為5 000元/月、PD期費用為2 000元/月。
由表2可見，如果醫(yī)保BIA中的研究疾病和干預(yù)方案下存在跨年度生存的患者，目前常見的3種算法對干預(yù)方案在2021年產(chǎn)生的醫(yī)保基金總支出的計算結(jié)果與本研究構(gòu)建的算法結(jié)果比較差異較大，且大多低于本研究算法結(jié)果。從理論上講，相比目前常見算法，本研究算法能夠?qū)IA目標(biāo)人群在研究年份的生存情況進行較為貼近現(xiàn)實的模擬測算，因而能夠為醫(yī)保BIA提供更為準(zhǔn)確的人均生存或用藥時長參數(shù)。
4 討論
本研究基于患者生存率構(gòu)建了一種針對跨年度生存患者的醫(yī)?；鹬С龅挠嬎惴椒ǎ啾饶壳霸擃I(lǐng)域中常見的算法，本算法可較為準(zhǔn)確地計算特定年份下某一干預(yù)方案的醫(yī)?；鹬С觯谝欢ǔ潭壬辖鉀Q了常見算法未充分考慮跨年度生存患者或計算過程簡易、粗略而導(dǎo)致醫(yī)?；鹬С鲱A(yù)測不準(zhǔn)確的問題，可為我國醫(yī)保BIA相關(guān)研究以及《指南》的完善提供參考，為研究者和決策者提供建議，促進我國醫(yī)療衛(wèi)生資源的合理配置。
4.1 本算法關(guān)鍵問題的解釋
首先，本算法在構(gòu)建過程中引入了某周期確診患者在研究年份內(nèi)的平均生存周期數(shù)csct和某年份確診患者在研究年份內(nèi)的平均生存周期數(shù)ascn這2個概念以便于計算，其實際意義理解起來較為晦澀，但其基本含義均是指某一時間段內(nèi)確診的患者在研究年份或研究年份內(nèi)某一周期的平均生存時間。需要注意的是，這與傳統(tǒng)意義上的患者平均生存時間（患者從確診或開始干預(yù)直至死亡的時間跨度）的概念有所差異。
其次，本研究引入了負生存周期數(shù)的概念（公式1），該設(shè)計并無實際含義，僅為方便統(tǒng)一和簡化公式而設(shè)計，因而無需從現(xiàn)實意義角度進行理解。
第三，本算法可基于多狀態(tài)模型計算，其計算思路與上文中所述單個狀態(tài)的計算思路完全一致。例如在“3”項下實例中，本研究以存在PFS和PD 2個健康狀態(tài)的癌癥為例計算了各狀態(tài)下的患者平均生存時間。
4.2 本算法的使用條件和關(guān)鍵設(shè)置
因為本算法是基于患者生存率構(gòu)建的，因此使用本算法的基本條件是可以獲得符合研究需求的生存數(shù)據(jù);同時，為計算研究年份內(nèi)的醫(yī)?；鹬С觯枰@取干預(yù)方案在各年份新確診患者中的使用人數(shù)，該參數(shù)可通過專家咨詢、臨床調(diào)研等方式獲得。
在將本算法應(yīng)用到醫(yī)保BIA的過程中，有2個關(guān)鍵參數(shù)需要進行設(shè)置：1）研究范圍內(nèi)生存隊列的最早產(chǎn)生年份ymax，可以通過設(shè)置模型模擬結(jié)束的生存率閾值來確定，且要注意該值必須為整數(shù);2）每年的周期數(shù)，可以根據(jù)疾病特征來設(shè)計，但為方便計算，須盡量使每年周期數(shù)為整數(shù)。
4.3 特殊情形下本算法的拓展
由于實際市場和臨床情況比較復(fù)雜，上述算法需根據(jù)實際用藥情況靈活改變，以下列舉了本算法在一些特殊情形下的拓展情況。
4.3.1 新藥上市或適應(yīng)證獲批時間晚于理論上的ymax 對于一些新藥或獲批新適應(yīng)證的藥品，在進行醫(yī)保BIA過程中，若按照生存率閾值確定ymax，則可能出現(xiàn)新藥上市年份或新適應(yīng)證獲批年份晚于ymax的情況，此時計算所得的醫(yī)保基金支出有悖于實際情況。
因此，對于新藥或新獲批適應(yīng)證的藥品，ymax的設(shè)計應(yīng)考慮藥品上市或適應(yīng)證獲批的年份——若新藥上市或新適應(yīng)證獲批的年份早于理論上的ymax，則按照生存率閾值確定ymax;若新藥上市或新適應(yīng)證獲批的年份晚于理論上的ymax，則以新藥上市或新適應(yīng)證獲批年份作為ymax。在實際中，新藥上市時間或適應(yīng)證獲批時間可能在年末，導(dǎo)致該年的銷售量較低，此時可根據(jù)藥品銷售情況調(diào)整算式中的起始年份為藥品上市或適應(yīng)證獲批年份的后一年，或僅將獲批年份下產(chǎn)生的隊列作為該年度的全部患者隊列來計算ascn。
4.3.2 治療過程中僅接受干預(yù)方案的一次性治療 在部分研究中，可能會出現(xiàn)患者僅接受該研究關(guān)注的干預(yù)方案的一次性治療，后續(xù)治療周期中接受其他干預(yù)方案治療的情況[13]。該情況下，由于對醫(yī)保基金產(chǎn)生影響的變動僅發(fā)生在研究關(guān)注的年份內(nèi)，故ymax可直接設(shè)置為研究關(guān)注的多個年份中最早的一年。例如癌癥放療方案的BIA中，一般假設(shè)干預(yù)方案本身的費用在患者確診周期已全部發(fā)生，但后續(xù)的患者生存情況會因選用不同的放療技術(shù)而存在差異。在這樣的情況下，實際上干預(yù)方案調(diào)入或調(diào)出報銷范圍對該方案選擇的影響僅出現(xiàn)在報銷范圍調(diào)整生效后的年份，因此僅需計算報銷范圍調(diào)整后年份的人均生存或用藥時間。
4.3.3 實際用藥時長與理論時長有差異 這類情形主要分為2種：①患者中止治療。本算法中所使用的生存數(shù)據(jù)來源于臨床試驗和已公開發(fā)表的文獻，研究對象在用藥時應(yīng)處于嚴(yán)格的臨床試驗條件下或其用藥依從性高，其生存數(shù)據(jù)反映了患者用藥后理想的生存情況。但在實際治療過程中，用藥依從性較差、藥物可及性較低等原因可導(dǎo)致患者突然停止治療[14]，此時若仍使用文獻報道的生存數(shù)據(jù)可導(dǎo)致生存率虛高。在該情況下，應(yīng)考慮將ymax設(shè)置為絕大部分患者的最長用藥時長，例如發(fā)現(xiàn)研究對象中95%的患者都會在2.5年內(nèi)因依從性差而停止用藥，但理論上的ymax為5年，則ymax實際應(yīng)設(shè)置為3年。此處建議不以患者的平均用藥時長來估算ymax，否則會導(dǎo)致部分患者在超過該時長后的用藥被忽略，繼而低估醫(yī)?；鹬С觥＂诨颊咄Ｖ褂盟幒?，該研究所關(guān)注的干預(yù)方案的藥效仍持續(xù)了一段時間。生存率反映的是干預(yù)方案的臨床治療效果，但某些長期疾病的患者在停止用藥后，之前服用的藥物仍可維持病情的穩(wěn)定，此時基于生存率計算而得到的平均生存時間并不等于實際用藥時間，導(dǎo)致醫(yī)?；鹬С鲱A(yù)測不準(zhǔn)確。在該情況下，可設(shè)置用藥時長閾值，或根據(jù)實際用藥時長設(shè)置生存率閾值，使計算而得的患者隊列的平均用藥時長接近實際用藥時長。
4.3.4 年份內(nèi)各周期確診人數(shù)存在明顯差異 在常見的流行病學(xué)資料中，研究者往往僅能獲得疾病的年發(fā)病率數(shù)據(jù)，而沒有各周期的發(fā)病率數(shù)據(jù)。本研究為貼近實際，且為方便解釋計算原理，故假設(shè)一年中各周期的發(fā)病率相等。但在實際研究過程中，研究者也許能夠獲得或計算出更為詳盡的發(fā)病率數(shù)據(jù)，這種情況下，可以為模型中每個隊列單獨賦予確診人數(shù)，并單獨計算每個隊列的生存時間，以更充分地利用高質(zhì)量的發(fā)病率數(shù)據(jù)。
4.3.5 新確診患者和跨年度生存患者的費用結(jié)構(gòu)存在差異 部分研究中，新確診患者與跨年度生存患者的年均總費用構(gòu)成可能存在差異。例如放療、手術(shù)等干預(yù)方案成本在BIA中一般僅會被納入新確診患者的費用構(gòu)成中。對此，可以分開計算本年度新確診患者隊列和跨年度生存患者隊列的年均總費用，為本年度新確診患者單獨增加成本項目。
4.4 局限性
本研究構(gòu)建的算法存在一定局限性：第一，本研究構(gòu)建的算法需要使用該干預(yù)方案的患者的生存數(shù)據(jù)和各年份使用該干預(yù)方案的患者人數(shù)，如果這2項數(shù)據(jù)無法獲取，則無法使用本算法;第二，由于實際中難以追溯年份內(nèi)不同周期確診的患者人數(shù)，而本算法假設(shè)不同周期確診的患者人數(shù)相等，因此計算得到的平均生存率可能與實際存在一定偏差;第三，由于醫(yī)保BIA關(guān)注的實際市場情況和臨床情況復(fù)雜，而本算法尚未經(jīng)過大量研究的應(yīng)用和實證，因此本研究對一部分特殊情形下算法的拓展提出了建議，但在研究一些具體BIA問題時，本算法的適用性仍可能較低，有待繼續(xù)改進或開發(fā)更為合理的計算方法。
5 結(jié)語
長期疾病藥物干預(yù)方案的醫(yī)保BIA中，目前常見的計算方法常未考慮跨年度生存的患者群體，或計算方法較為粗略，對醫(yī)保基金支出估算的準(zhǔn)確性產(chǎn)生了影響。本研究以完善醫(yī)保BIA建模思路為目的，基于患者生存率構(gòu)建了一種存在跨年度生存患者情形下的醫(yī)?；鹬С鲇嬎惴椒?。本方法的基本思路是計算研究年份及其以往各年份內(nèi)使用某干預(yù)方案（如藥物治療方案）的確診患者在研究年份內(nèi)的平均生存周期數(shù)，再根據(jù)各年份的確診患者中使用該干預(yù)方案的人數(shù)和人均周期成本計算該干預(yù)方案在研究年份的總醫(yī)?；鹬С?。本算法能夠在一定程度上解決目前常見計算方法中對跨年度生存情形下的醫(yī)保基金支出預(yù)測不準(zhǔn)確的問題。但由于實際市場情況和臨床情況復(fù)雜，本算法需根據(jù)實際用藥情況靈活改變。
參考文獻
[ 1 ] 國家醫(yī)療保障局. 2020年國家醫(yī)保藥品目錄調(diào)整工作方案[EB/OL].（2020-08-17）[2021-08-06].http：//www.nhsa.gov.cn/art/2020/8/17/art_37_3407.html.
[ 2 ] 劉國恩.中國藥物經(jīng)濟學(xué)評價指南：2020中英雙語版[M]. 北京：中國市場出版社，2020：91-92.
[ 3 ] WU Y L，CHU D T，HAN B，et al. Phase Ⅲ，randomized，open-label，first-line study in Asia of gefitinib versus carboplatin/paclitaxel in clinically selected patients with advanced non-small-cell lung cancer：evaluation of patients recruited from mainland China[J]. Asia Pac J Clin Oncol，2012，8（3）：232-243.
[ 4 ] MOVAHEDI M，HEPWORTH E，Mirza R，et al. Discontinuation of biologic therapy due to lack/loss of response and adverse events is similar between TNFi and non-TNFi class：results from a real-world rheumatoid arthritis cohort[J]. Semin Arthritis Rheum，2020，50（5）：915-922.
[ 5 ] 謝麗，夏茹意，莊貴華.伴代償期肝硬化乙肝相關(guān)肝癌患者根治性切除術(shù)后抗病毒治療的預(yù)算影響分析[J].中國衛(wèi)生經(jīng)濟，2019，38（9）：64-68.
[ 6 ] 殷瀟，方欣，胡江藺，等.舒尼替尼治療轉(zhuǎn)移性腎細胞癌的預(yù)算影響分析[J].中國衛(wèi)生經(jīng)濟，2016，35（12）：67-69.
[ 7 ] 張欲曉，殷瀟，田夢媛，等.赫賽汀治療HER-2陽性乳腺癌的醫(yī)保預(yù)算影響分析[J].中國衛(wèi)生經(jīng)濟，2016，35（12）：63-66.
[ 8 ] 陳斌斌，范長生.注射用紫杉醇（白蛋白結(jié)合型）治療晚期乳腺癌的預(yù)算影響分析[J].中國醫(yī)療保險，2016（12）：60-63.
[ 9 ] 范長生，何敏媚，吳久鴻.復(fù)方黃黛片治療急性早幼粒細胞白血病預(yù)算影響分析[J].中國新藥雜志，2017，26（10）：1219-1224.
[10] 張強，崔穎，馬銳，等.規(guī)范使用唑來膦酸治療國內(nèi)晚期肺癌或乳腺癌合并骨轉(zhuǎn)移預(yù)算影響分析[J].中國醫(yī)療保險，2018（12）：59-67.
[11] 張潔，黃泰康. 藥物經(jīng)濟學(xué)評價中馬爾可夫模型的優(yōu)化：馬爾可夫決策模型[J].中國藥房，2008，19（17）：1284- 1286.
[12] 邵榮杰，唐文熙，馬愛霞.分區(qū)生存模型在藥物經(jīng)濟學(xué)評價中的應(yīng)用[J].中國衛(wèi)生經(jīng)濟，2019，38（9）：60-63.
[13] 汪紅艷，孔令玲，王凡，等.老年食管癌放射治療療效及預(yù)后影響因素分析[J].安徽醫(yī)科大學(xué)學(xué)報，2016，51（8）：1188-1192.
[14] 洪梅花，徐東曉，童素梅，等.急性心肌梗死支架術(shù)后患者不同時期服藥依從性及其影響因素的調(diào)查與分析[J].中國護理管理，2014，14（12）：1261-1265.
（收稿日期：2021-01-27 修回日期：2021-08-06）
（編輯：胡曉霖）
猜你喜歡
中學(xué)生數(shù)理化·八年級物理人教版(2022年5期)2022-06-05 06:57:38
祝您健康·文摘版(2020年7期)2020-07-13 09:35:19
中國生殖健康(2019年2期)2019-08-23 08:11:52
世界知識(2018年19期)2018-11-21 03:01:12
老年教育(老年大學(xué))(2018年12期)2018-01-28 14:15:50
現(xiàn)代工業(yè)經(jīng)濟和信息化(2016年2期)2016-05-17 05:34:16
癌癥進展(2016年12期)2016-03-20 13:16:14
材料科學(xué)與工程學(xué)報(2016年5期)2016-02-27 07:11:22
中國醫(yī)學(xué)影像學(xué)雜志(2015年9期)2015-12-15 11:03:26
西部中醫(yī)藥(2015年9期)2015-02-02 01:19:24
