綜合利用兩個數(shù)據(jù)庫進行藥品不良反應風險信號挖掘方法的探索

王芳，李永輝，郭瑞鋒，劉欣欣，張晶晶

目前，我國藥品的質(zhì)量評價所面臨的一個難題就是仿制藥品眾多，同一種藥品由多家企業(yè)生產(chǎn)，其生物效價、臨床療效和不良反應不盡相同[1]，如何選用和比較同種藥品的質(zhì)量、臨床療效以及安全風險是當前所必須面對、亟待解決的問題，也是藥學工作者推進合理用藥、減少藥源性損害發(fā)生的一項重要工作[2]。研究[3-9]表明藥品不良反應的發(fā)生和藥品質(zhì)量具有相關性，通過比較不同廠家生產(chǎn)的同通用名藥品不良反應發(fā)生率探究藥品安全風險與質(zhì)量的相關性具有一定的挑戰(zhàn)性和新穎性。藥品不良反應監(jiān)測數(shù)據(jù)是做好藥品監(jiān)管工作的技術基礎和科學依據(jù)，充分分析、有效利用這些數(shù)據(jù)，檢測和質(zhì)量相關的風險信號，是風險管理和科學監(jiān)管的重要環(huán)節(jié)和前提。

國家食品藥品監(jiān)督管理總局發(fā)布的《國家藥品不良反應監(jiān)測年度報告》[10-15]中連續(xù)六年(2011--2016年)抗感染藥不良反應病例報告數(shù)量居首位，排名首位的是頭孢菌素類抗生素，其嚴重不良反應以過敏反應為主。

本研究將國家藥品不良反應監(jiān)測系統(tǒng)的不良反應數(shù)據(jù)庫與河北省食品藥品監(jiān)督管理局藥品電子檢驗報告監(jiān)管部門核查系統(tǒng)中醫(yī)療機構藥品購進數(shù)據(jù)相結合進行分析，選取頭孢菌素類中的3個品種為樣本研究不同藥品生產(chǎn)企業(yè)間的同通用名藥品不良反應報告率的差異，差異有統(tǒng)計學意義的企業(yè)及其品種作為風險信號，為進一步的風險信號篩選、分析、確認、分析和風險管理找準方向。

1 資料與方法

1.1研究假設為研究不同企業(yè)生產(chǎn)的同通用名藥品之間安全性差異，可通過比較不同企業(yè)生產(chǎn)的該藥品不良反應發(fā)生率來實現(xiàn)。藥品的不良反應實際發(fā)生率可通過如下公式計算：

Ki=νi/ηi

其中Ki代表某企業(yè)生產(chǎn)的某類藥品的不良反應發(fā)生率，νi為該企業(yè)該藥品的不良反應發(fā)生次數(shù)，ηi為該企業(yè)的藥品使用人次。但由于實際工作中，νi與ηi均無法直接獲得，故可以考慮利用藥品不良反應上報數(shù)及藥品采購量來反映藥品不良反應發(fā)生率。

由于用藥人次ηi無法直接測量，可采用文獻中的辦法來估算用藥人數(shù)[16]。

用藥人數(shù)≈藥品采購量/(用藥日數(shù)×定義日劑量)

假設病人感染程度不同時對所用藥品的生產(chǎn)企業(yè)沒有傾向性，則可認為用藥日數(shù)×定義日劑量為一個常數(shù)，直接使用采購量來代替用藥人數(shù)，藥品不良反應報告率的絕對數(shù)值都會偏小，但ADR的發(fā)生強度的比較仍可進行。

不良反應自發(fā)報告系統(tǒng)統(tǒng)計數(shù)據(jù)與藥品的真實不良反應發(fā)生率之間存在差異，自發(fā)報告上報數(shù)ν′與真實不良反應發(fā)生數(shù)ν之間存在如下關系：

ν′=αν{|0≤α≤1}

為應對這種差異，本研究假設醫(yī)療機構對于不同企業(yè)生產(chǎn)的同類藥品的不良反應上報情況不存在傾向性，各企業(yè)產(chǎn)品的α值近似相等，使用ν′值計算的不良反應事件差異可以反映產(chǎn)品ν值的差異。

將不良反應上報數(shù)與用藥頻度作為替代指標衡量藥品安全性，則：

本研究使用藥品不良反應上報數(shù)與藥品采購量的比值作為藥品不良反應報告率，通過比較同種藥品不同生產(chǎn)企業(yè)間的藥品不良反應報告率，反映不同企業(yè)間藥品安全性差異情況。

1.2數(shù)據(jù)來源及研究方法藥品不良反應報告資料來自國家藥品不良反應監(jiān)測系統(tǒng)河北省數(shù)據(jù)庫，不良反應報告時間為2014年1月至2015年12月，根據(jù)藥品不良反應上報情況，選取抗感染類藥品中不良反應報告占比較高的頭孢菌素類中三個品種作為研究對象。納入研究范圍的通用名藥品包括注射用頭孢曲松鈉、注射用頭孢呋辛鈉和注射用頭孢哌酮鈉舒巴坦鈉;對不良反應數(shù)據(jù)進行標準化處理后，按通用名、生產(chǎn)企業(yè)分類做描述統(tǒng)計分析。研究進一步采集河北省食品藥品監(jiān)督管理局藥品電子檢驗報告監(jiān)管部門核查系統(tǒng)中醫(yī)療機構藥品采購數(shù)據(jù)。以研究藥品的ADR例數(shù)與對應的采購量比值來表示不良反應報告率。對同通用品不同企業(yè)的藥品不良反應報告率差異進行比較，以χ2檢驗判別結果，P<0.05為差異有統(tǒng)計學意義，并計算其相對危險度(Relative risk，RR)[17]。

由于無法獲得藥品實際使用量信息，本研究采用藥品采購數(shù)量替代藥品使用量信息計算不良反應報告率。藥品采購數(shù)量來源是醫(yī)療機構通過河北省食品藥品監(jiān)督管理局藥品電子檢驗報告監(jiān)管部門核查系統(tǒng)查驗購進藥品電子檢驗報告時所產(chǎn)生的購進數(shù)量。此系統(tǒng)是由河北省食品藥品監(jiān)督管理局開發(fā)并在全省范圍內(nèi)使用的。省內(nèi)各涉藥單位在采購藥品時必須索要或查驗所購進藥品的檢驗報告并如實填報采購數(shù)量。管理人員可通過登錄系統(tǒng)查詢省內(nèi)各類涉藥單位(包括藥品生產(chǎn)企業(yè)、經(jīng)營企業(yè)、藥品使用單位等)藥品供應或采購數(shù)量。

2 結果

2.1數(shù)據(jù)清洗本研究對不良反應數(shù)據(jù)中生產(chǎn)企業(yè)、劑型或給藥方式等信息缺失或明顯錯誤的條目進行人工清洗，并與藥品采購數(shù)據(jù)中的生產(chǎn)企業(yè)進行匹配，篩除不合格數(shù)據(jù)及未匹配數(shù)據(jù)。并以清洗匹配后的不良反應例數(shù)(n)與該藥品年度報告的不良反應例數(shù)(N)的比值計算匹配度。

3種藥品連續(xù)2年不良反應報告率及供應量匹配情況見表1。

2.2不良反應報告率比較為研究各企業(yè)生產(chǎn)的藥品與同品種其他企業(yè)產(chǎn)品的安全性，本研究按下式計算各企業(yè)與其他產(chǎn)品不良反應報告率的相對危險度(RR)，并以χ2檢驗結果判別其統(tǒng)計顯著性。

本研究共涉及三種同通用名藥品，企業(yè)采購量信息180條，其中41條(22.78%)企業(yè)不良反應報告率顯著高于同品種其他企業(yè);63條(35.00%)企業(yè)不良反應報告率顯著低于同品種其他企業(yè);34條(18.89%)與同品種其他企業(yè)相比差異無統(tǒng)計學意義;42條(23.33%)只有采購量數(shù)據(jù)，未見不良反應報告，見表2。

3 討論

3.1不良反應報告率可提示藥品生產(chǎn)企業(yè)質(zhì)量安全風險已有研究表明同品種藥品的不良反應發(fā)生率與不同的藥品生產(chǎn)企業(yè)存在關聯(lián)性：阿奇霉素國產(chǎn)較進口合資的不良反應發(fā)生率高[3];華蘭生物的流感疫苗較葛蘭素史克的ADR發(fā)生率高[4];國產(chǎn)甲氨蝶呤較進口的致肝損害的發(fā)生率高[5]。不同廠家生產(chǎn)的頭孢拉定引起血尿的發(fā)生率存在差異[6]。

此外有研究表明不良反應發(fā)生率與藥品質(zhì)量相關：抗生素的藥品不良反應發(fā)生率偏高的生產(chǎn)企業(yè)其雜質(zhì)總量的數(shù)據(jù)偏高[7];頭孢拉定中不溶性微粒增多會導致血尿不良反應發(fā)生[8]。藥品流通環(huán)節(jié)采取針對性風險管理措施能夠顯著降低藥品的不良發(fā)生發(fā)生率[9]。

藥品不良反應的發(fā)生率能夠快速反映藥品的質(zhì)量情況，揭示不同企業(yè)間的質(zhì)量差異性，是風險在臨床用藥安全的一個窗口，可作為藥品安全性評價的直接證據(jù)。由于實際工作中不良反應發(fā)生率難以獲得，本研究通過不良反應報告率間接反映不良反應發(fā)生率，來快速提取風險信號，識別風險和確認風險，為后續(xù)的信號分析評價研究及利用找準方向，適應當前形勢、經(jīng)濟、實用、可操作。

3.2兩個數(shù)據(jù)庫(不良反應數(shù)據(jù)與藥品采購數(shù)據(jù))匹配度較高河北省每年通過國家藥品不良反應監(jiān)測收集的報告數(shù)4萬余份，但僅能反映發(fā)生頻數(shù)，因沒有使用人次數(shù)，不能反映不良反應發(fā)生率的情況。藥品電子檢驗報告監(jiān)管部門核查系統(tǒng)僅能提供藥品采購量的情況，不能反映不良反應發(fā)生頻數(shù)的情況。綜合利用該兩個數(shù)據(jù)庫可獲得不良反應報告率，可作為挖掘不良反應風險信號的切入點。

國家藥品不良反應監(jiān)測系統(tǒng)的不良反應上報數(shù)據(jù)與藥品電子檢驗報告監(jiān)管部門核查系統(tǒng)中藥品采購量信息匹配度在73.04%～95.24%之間，總體匹配度為86.39%。總體來看，ADR數(shù)據(jù)與藥品采購數(shù)據(jù)匹配情況較理想，具有一定的分析利用價值，可綜合利用該兩個系統(tǒng)進行數(shù)據(jù)挖掘。

3.3不同企業(yè)生產(chǎn)的同種藥品不良反應報告率之間存在差異由本研究結果來看，不同企業(yè)生產(chǎn)的同品種藥品不良反應報告率之間確實存在一定的差異性，22.78%的藥品不良反應報告率高于同類其他企業(yè)生產(chǎn)的藥品，差異有統(tǒng)計學意義，18.89%的藥品不良反應報告率低于同類其他企業(yè)生產(chǎn)的藥品，差異有統(tǒng)計學意義。同品種之間的不良反應報告率差異較大，相對危險度最高達431.10，風險信號指示強烈。不良反應報告率高于同類其他藥品，差異有統(tǒng)計學意義的的條目可視為產(chǎn)生一個不良反應風險信號，提示我們該企業(yè)生產(chǎn)的該類藥品可能存在安全性隱患，應加強監(jiān)測并結合其他相關方法找出風險原因并判定其風險等級。

表1 目標藥品不良反應上報及用量匹配情況

注：采購量=規(guī)格×支數(shù)

表2 不同生產(chǎn)企業(yè)不良反應報告率比較情況/家

3.4綜合利用兩個數(shù)據(jù)庫進行藥品不良反應風險信號挖掘方法的特點

3.4.1快速、準確、經(jīng)濟甄別風險信號 獲取準確的不良反應發(fā)生率研究成本高，需投入大量的人力、財力、物力，比較眾多不同廠家的ADR發(fā)生率需要的樣本量大、工作量巨大、周期長，難以實現(xiàn)。該方法快速、準確、經(jīng)濟的提取不良反應風險信號，為進一步深入探究找準方向和目標。

使用不良反應報告率替代ADR發(fā)生率，有產(chǎn)生偏倚的可能性，為減少假陽性信號的干擾，提高方法準確性，可在差異有統(tǒng)計學意義的前提下為相對危險度設定高閾值，閾值過高會遺漏信號，閾值可根據(jù)實際情況進行調(diào)整。

3.4.2與既有信號檢測方法的不同之處 當前藥品不良反應比較經(jīng)典的信號檢測方法有頻數(shù)法(PRR、ROR、MHRA)和貝葉斯法(BCPNN、MGPS)等，數(shù)據(jù)來源為單個數(shù)據(jù)庫(如國家藥品不良反應監(jiān)測系統(tǒng)數(shù)據(jù)庫)，針對某藥品和某特定不良反應提取風險信號。本研究信號檢測方法綜合利用兩個數(shù)據(jù)庫，在確定為同品種一藥品的情況下，重在區(qū)別企業(yè)間差異性，針對藥品不良反應的發(fā)生情況與企業(yè)之間的關聯(lián)性提取風險信號。

3.4.3信號的應用 提取到的風險信號，可提供給臨床加強監(jiān)測，對風險進行進一步分析評估;可反饋給相關藥品生產(chǎn)企業(yè)，企業(yè)對風險信號進行分析評價，進一步探究其內(nèi)在原因加強藥品生產(chǎn)質(zhì)量控制;可作為藥品安全性風險信號供藥品監(jiān)管部門或相關研究機構做進一步研究與參考。

該方法可進一步拓展應用，如比較同一生產(chǎn)企業(yè)同一品種不同生產(chǎn)批次的ADR報告率差異，本研究進行過程中，曾發(fā)現(xiàn)某企業(yè)A批次的ADR報告率高于其他批次(χ2=6.85，P=0.008 9)，差異有統(tǒng)計學意義，在對企業(yè)進行藥品生產(chǎn)質(zhì)量管理規(guī)范(GMP)認證檢查時重點查閱了該兩個批次的生產(chǎn)批記錄，A批次生產(chǎn)期間懸浮粒子量有1 h超標，信號得到了初步確認。

3.5展望將藥品不良反應數(shù)據(jù)庫與藥品電子檢驗報告監(jiān)管部門核查系統(tǒng)數(shù)據(jù)進行聯(lián)合分析，能夠快速、準確甄別出同通用名下不同企業(yè)生產(chǎn)的藥品之間的安全性可能存在的差異，但受限于數(shù)據(jù)來源的局限性，該研究目前尚不能夠作為確定性證據(jù)反映藥品的安全等級，但可作為風險信號提示，縮小研究范圍，鎖定目標企業(yè)，據(jù)此進行藥物流行病學研究或相關臨床試驗研究以對相關藥品的安全性差異進行進一步探討和確認。待藥品不良反應監(jiān)測工作達到較成熟的階段，不良反應報告數(shù)接近甚至近似等于不良反應發(fā)生次數(shù)和藥品電子檢驗報告電子平臺更加全面完善時，不良反應報告率會更接近真實的不良反應發(fā)生率，風險信號提取的準確性將大大提升，方法進一步確立。