隨機(jī)化檢驗(yàn)在藥物經(jīng)濟(jì)學(xué)評價中的應(yīng)用

盧 雙 薛 培 王 煊 朱文濤

隨機(jī)化檢驗(yàn)在藥物經(jīng)濟(jì)學(xué)評價中的應(yīng)用

盧 雙1薛 培2王 煊2朱文濤2

目前，國內(nèi)關(guān)于藥物經(jīng)濟(jì)學(xué)評價大多采用點(diǎn)估計(jì)的方法簡單對比不同用藥方案的評價指標(biāo)，但是基于統(tǒng)計(jì)推斷的理論，這種簡單處理并不能反映研究人群總體的比較結(jié)果。本文將隨機(jī)化檢驗(yàn)應(yīng)用于藥物經(jīng)濟(jì)學(xué)評價中，能在一定程度上解決藥物經(jīng)濟(jì)學(xué)療效評價的問題。

隨機(jī)化檢驗(yàn)；統(tǒng)計(jì)推斷；藥物經(jīng)濟(jì)學(xué)評價

藥物經(jīng)濟(jì)學(xué)評價是通過識別、測量和比較不同藥物、治療方案及衛(wèi)生服務(wù)項(xiàng)目的成本和社會經(jīng)濟(jì)效果，運(yùn)用統(tǒng)計(jì)方法，并結(jié)合經(jīng)濟(jì)模型對藥物經(jīng)濟(jì)學(xué)評價的結(jié)果進(jìn)行客觀、科學(xué)的分析，為藥品費(fèi)用控制的關(guān)鍵環(huán)節(jié)，如藥品報(bào)銷目錄及價格的制訂、臨床醫(yī)療指南的制訂和修訂提供依據(jù)。其評價方法有成本-效果分析、成本-效用分析等。針對成本-效用分析，傳統(tǒng)計(jì)算方法如下：A=C1/E1－C2/E2

其中A代表評價結(jié)果，C1、C2分別為兩組樣本的平均成本，E1、E2分別為兩組樣本的平均療效。相應(yīng)A代表兩組樣本最終的成本-效用指標(biāo)評價結(jié)果。一般國內(nèi)文獻(xiàn)均根據(jù)此樣本計(jì)算結(jié)果給出藥物經(jīng)濟(jì)學(xué)評價的結(jié)論。但是，根據(jù)統(tǒng)計(jì)推斷的原理，如果僅關(guān)心進(jìn)入研究的樣本人群評價結(jié)果，分析到此是可以接受的。想要給出樣本人群背后所代表的的總體人群的結(jié)論，這個結(jié)果是不科學(xué)的，因此需要根據(jù)樣本人群去推斷總體人群的結(jié)果。

常用的統(tǒng)計(jì)推斷包括參數(shù)檢驗(yàn)和非參數(shù)檢驗(yàn)。前者要求明確數(shù)據(jù)的分布、構(gòu)建概率分布、完成統(tǒng)計(jì)推斷；后者對數(shù)據(jù)的分布無任何要求。明確成本-效用指標(biāo)的統(tǒng)計(jì)分布比較麻煩，因此筆者檢驗(yàn)采用非參數(shù)檢驗(yàn)的方法完成統(tǒng)計(jì)推斷的工作。隨機(jī)化檢驗(yàn)是比較理想的選擇。

隨機(jī)化檢驗(yàn)，也稱為排列檢驗(yàn)[1]，是最早用以檢驗(yàn)兩個總體均值差的方法之一。雖然這種方法本質(zhì)上對數(shù)據(jù)的分布沒有任何要求，但因涉及大量的計(jì)算，所以直到最近才真正可行。它指的是從所有結(jié)果中隨機(jī)產(chǎn)生部分結(jié)果。而不是采用所有的可能結(jié)果進(jìn)行計(jì)算。具體操作如圖 1，最后根據(jù)結(jié)果是否落在拒絕域決定是否拒絕原假設(shè)。下面根據(jù)某臨床研究數(shù)據(jù)實(shí)現(xiàn)以上過程。

1 資料來源及數(shù)據(jù)特征

本研究數(shù)據(jù)來源于北京中醫(yī)藥大學(xué)藥物經(jīng)濟(jì)學(xué)研究所某關(guān)于心血管疾病的課題，從中隨機(jī)抽取102份樣本，其中使用A、B兩種治療方案的樣本分別為70份、32份，見表1。

圖1 隨機(jī)化抽樣流程

表1 人口學(xué)及疾病信息

2 數(shù)據(jù)分析及結(jié)果

以上程序中g(shù)roup為分組變量，f1、f2分別代表效用及成本，col1、col2分別代表兩種治療方案的成本-效用比，f3代表兩種治療方案的成本-效用比差值。原始數(shù)據(jù)效用、成本信息見表2、3，原始數(shù)據(jù)成本-效用比分析結(jié)果見表4，隨機(jī)抽樣500次計(jì)算結(jié)果見表5。

根據(jù)表4和表5的計(jì)算結(jié)果，原始研究數(shù)據(jù)的成本-效用比組間差值在成本-效用比差值95% CI范圍內(nèi)。設(shè)定可接受一類錯誤為雙側(cè)0.05，據(jù)此分析結(jié)果，尚不能認(rèn)為兩種不同用藥方案成本-效用比存在差異。

表2 兩種治療方案視覺模擬量表(VAS)評分結(jié)果(分)

表3 兩種治療方案成本信息(元)

表4 根據(jù)試驗(yàn)原始數(shù)據(jù)成本-效用比計(jì)算結(jié)果

表5 隨機(jī)抽樣成本-效用比計(jì)算結(jié)果

3 討論

本研究假設(shè)原始數(shù)據(jù)來源于隨機(jī)試驗(yàn)，另外計(jì)算過程中未考慮到極值問題，其計(jì)算過程并不完美。但為藥物經(jīng)濟(jì)學(xué)統(tǒng)計(jì)推斷提供了一個思考的方向。對于非隨機(jī)化的觀察性研究數(shù)據(jù)，一些與評價指標(biāo)相關(guān)的因素可能在組間分布不均衡，其因果推斷不像隨機(jī)化試驗(yàn)可信度高，此時建議在完成傾向評分匹配[2-3]的基礎(chǔ)上，再進(jìn)行隨機(jī)化檢驗(yàn)不失為一種可以研究的方向，比如本研究的數(shù)據(jù)中患病時間存在組間差異，可以先根據(jù)傾向值匹配的原理進(jìn)行組間匹配，之后在完成隨機(jī)化檢驗(yàn)。有學(xué)者提出[4]采用自主法估算區(qū)間，與本文中提出的方法具有一定的相似性，讀者可自行選擇合適的方法。另有學(xué)者[5]研究提出其他區(qū)間估計(jì)方法，如BOX法、橢圓法、Taylor級數(shù)擴(kuò)展方法。但在計(jì)算機(jī)廣泛應(yīng)用的今天，采用重復(fù)抽樣的方法計(jì)算區(qū)間相對比較簡單。

[1]William Navidi.統(tǒng)計(jì)學(xué)—科學(xué)與工程應(yīng)用[M].北京∶清華大學(xué)出版社,2007∶449-449.

[2]郭申陽.傾向值分析∶統(tǒng)計(jì)方法與應(yīng)用[M].重慶∶重慶大學(xué)出版社,2012∶84-109.

[3]李智文,任愛國.傾向性評分在 SAS軟件中的實(shí)現(xiàn)[J].中國生育健康雜志,2010,21(5)∶320-322.

[4]謝穎,宗欣,孫麗華.非參數(shù) Bootstrap法在成本-效果置信區(qū)間估計(jì)中的應(yīng)用[J].中國藥房,2010,21(22)∶2027-2029.

[5]楊莉.藥物經(jīng)濟(jì)學(xué)評價需要求、方法和指南研究[D].復(fù)旦大學(xué)博士論文,上海,2010.

The Application of Randomization Test in Pharmacoeconomics Evaluation

Lu Shuang Xue Pei Wang Xuan Zhu Wentao

At present,the method of point estimation is used generally to assess drug economics between different regimen.But based on the theory of statistical inference,this simple process can't reflect the overall study population. By randomization test,we can solve the problems.

Randomization test;Statistical Inference;Pharmacoeconomics evaluation

R956

A

1673-5846(2015)01-0012-03

1北京愛恩康臨床醫(yī)學(xué)研究（北京）有限公司統(tǒng)計(jì)部，北京 100013

2北京中醫(yī)藥大學(xué)，北京 100013

朱文濤，北京中醫(yī)藥大學(xué)管理院教授、碩士生導(dǎo)師，北京中醫(yī)藥大學(xué)中藥藥物經(jīng)濟(jì)學(xué)評價研究所所長。研究方向：藥物經(jīng)濟(jì)學(xué)、中藥質(zhì)量風(fēng)險管理。E-mail：wentao67@126.com