生物等效性試驗中藥代動力學報表自動輸出的SAS實現

孟娜 張糈 胡建華

[摘要]目的為滿足藥物一致性評價中藥代動力學（PK）參數生物等效性分析計算以及結果報表規范化自動輸出。方法將 SAS 的數據處理功能與輸出交付系統相結合，編寫自動輸出 PK 參數統計分析報表的 SAS 宏程序%AbeTable。結果以加拿大衛生部發布有關文件中的 PK 數據為實例，程序計算結果準確，輸出的統計報表格式規范，符合臨床試驗統計分析報告的基本要求。結論宏%AbeTable具有較好的通用性，結果可靠，生成報表效率高，可以滿足生物等效性分析的統計需要。

[關鍵詞]生物等效性試驗;WinNonlin;藥代動力學參數; SAS 制表

[中圖分類號] R285 [文獻標識碼] A? [文章編號]2095-0616（2021）22-0216-04

Realization of automatic output of pharmacokinetic reports for bioequivalence test by SAS

MENG? Na?? ZHANG? Xu? ?HU? Jianhua

Department of Pharmaceutical Information and Biostatistics， School of Science of China Pharmaceutical University， Jiangsu， Nanjing 210000， China

[Abstract] Objective To satisfy the analysis and calculation of bioequivalence and pharmacokinetics （PK） parameters for drug consistency evaluation， and to realize the standardized automatic output of result reports. Methods The data processing function of Statistical Analysis System （SAS） was combined with the output delivery system to program the SAS macro program， %Abe Table， for automatic output of statistical analysis reports of PK parameters. Results Taking the PK data in the relevant documents released by Health Canada as an example， the output statistical reports was standardized in format if the results calculated by the program were accurate， which met the basic requirements of statistical analysis reports of clinical trials. Conclusion The macro program %Abe Table has good generality with reliable results， and high efficiency in generating reports， which can meet the statistical requirements of bioequivalence analysis.

[Key words] Bioequivalence test; WinNonlin; Pharmacokinetic parameters; SAS tabulation

截止至2017年11月底，我國生產的仿制藥占全部藥品數量的95%以上[1]。為保證仿制藥的療效與安全，監管部門對藥物的一致性評價工作愈發重視[2-5]。制藥企業使用WinNonlin進行生物等效性評價。WinNonlin輸出結果的格式固定，冗余內容多，報告撰寫人員需人工完成數據整理工作。此過程耗時長、錯誤率高，違背了數據結果的可重復性原則[6-7]。針對上述問題，編寫 SAS 宏程序%AbeTable，自動輸出規范化統計報表，提高統計工作效率與質量。

1生物等效性研究的藥代動力學（PK）分析

1.1生物等效性計算

本研究以2×2交叉設計展開研究。試驗的終點指標一般為峰濃度（Cmax ），血藥濃度0到 t 下面積（area under the curve from zero to t time point， AUC0-t ），血藥濃度0到∞下面積（area under the curve from zero to infinity， AUC0-∞）[8-9]。2×2試驗的生物等效性評價是在方差分析的基礎上進行置信區間計算與雙單側 t 檢驗。

1.1.1方差分析方差分析將給藥順序、藥物、周期作為固定效應，受試者（順序）作為隨機效應[10-11]建立模型。通過方差分析可以得出統計量與 P 值，但不能得出等效性結論，雙單側 t 檢驗與置信區間法實質性地解決了該問題。

1.1.2置信區間利用方差分析的結果，計算置信區間，其上下限計算公式分別為：

式中，S是方差分析誤差項均方的平方根; t1-α，v是檢驗水平為α，自由度為 v = n -2時 t 分布的界值;L T、LR 分別是兩制劑 P K 參數取對數后的算數均值。若計算出100（1-2α）%的置信區間介于80%～125%范圍，可判定為生物等效。

1.1.3雙單側 t 檢驗利用方差分解后獲得的均方誤差（mean square error，MSE）進行雙單側 t 檢驗計算。雙單側 t 檢驗的檢驗假設為：

H0：InμT - InμR≤ Inr1或 InμT - InμR ≥ Inr2 H1：Inr1< InμT - InμR < Inr2

其中， In μT? In μR 分別為兩制劑 P K 參數取對數后的算術均值，r 1和 r 2取值一般為0.8和1.25。

兩次單側檢驗的統計量分別為：

式中，InXT -InXR ～ N （In μT - In μR ，σInXT -InXR ），其標準誤為SInXT -InXR =。

若兩次檢驗的 P 值均<0.05，兩次檢驗均拒絕H0，可推斷兩種制劑具有生物等效性。

1.2生物等效性評價結果輸出

在等效性試驗中，統計分析報告包含方差分析、100（1-2α）%置信區間與雙單側 t 檢驗的 PK 分析結果。方差分析的目的是展示出 PK 參數的給藥順序、受試者（順序）、周期與制劑的變異信息。100（1-2α）%置信區間分析結果中需要提供各參數的幾何均值及比值、變異系數、置信區間和把握度。最后，列出雙單側 t 檢驗的統計量與 P 值。

2 SAS宏程序%AbeTable

SAS 宏程序%AbeTable包含數據預處理、統計分析和報表格式化輸出三個子宏。數據預處理通過調用子宏%PREPROC 確保輸入數據集的規范性。統計分析是進行統計量的計算及統計結果規范化，由宏%PK_CAC 實現。宏%PK_CAC 根據 PROC MIXED 過程，獲得方差分析差值估計，進而計算得到置信區間，幾何均值比，t 值，P 值，個體內變異及幾何均值比等結果。格式化輸出宏%REPORT 結合了 SAS 的輸出交付系統與 PROC REPORT 過程，以及 PROC TEMPLATE 過程定義的標準化統計三線表模板 style.pk，實現了報表的格式化輸出。宏%AbeTable的數據處理過程見圖1。

3 SAS宏程序%AbeTable的應用

以加拿大衛生部發布的指導文件[12]中提供的 PK 數據為實例，展示%AbeTable在生物等效性試驗統計分析中的應用。臨床數據交換標準協會要求統計分析需在標準數據制表模型與分析數據模型上完成[13-15]。宏%AbeTable的適用對象為分析數據模型中的 A DPP 數據集，該數據集由 W i n N o nli n 計算并導出的 P K 參數衍生而來。

調用%AbeTable后，SAS 將進行 PK 參數的方差分析、置信區間的計算與雙單側檢驗，并將結果以三線表的形式輸出至 RTF 文檔，見表1～3。宏%AbeTable生成的表格基本符合統計分析報告的內容要求與格式規范，各統計量計算結果與WinNonlin軟件完全一致。

4討論

臨床試驗統計分析的核心體現于統計分析報告中，在保證質量的前提下提升工作效率是臨床試驗工作者長期追求的目標。吳正宇等用 R 語言計算藥動學參數和生物等效性分析[16]，張露等分別用 SAS 及stata實現2×2生物等效性檢驗，王建峰研究了Ⅰ期臨床中分析報告規范化結果輸出[17]。本研究根據醫藥公司現存的統計分析工作情況，結合該領域常用統計軟件的獨特優勢，編寫了一套自動生成 PK 參數統計分析報表的 SAS 宏程序，提升了工作效率。但仍存在一定不足：首先，本文僅研究了生物等效性 PK 參數的 SAS 宏生成，還未對統計分析報告中安全性結果報表生成進行探索;其次，本文還未實現用 SAS 對半重復或重復交叉設計等效性分析進行相應報告輸出。在今后的進一步研究中，將對生物等效性研究中安全性結果報表輸出及部分重復或完全重復的試驗類型展開討論。

[參考文獻]

[1]曾俊芬，宋金春.仿制藥質量和療效一致性評價思考[J].醫藥導報，2019，38（7）：884-887.

[2]張露，鄭倩雯，李亞文，等.生物等效性檢驗在 SAS 及stata中的實現及對比[J].中國衛生統計，2019，36（2）：306-308.

[3]王可瑩，劉平.我國仿制藥一致性評價現狀及展望[J].福建醫科大學學報（社會科學版），2018，19（3）：16-20.

[4]劉冬，哈莉莉，李芳，等.我國化學仿制藥一致性評價進展與展望[J].中國臨床藥理學雜志，2020，36（16）：2381-2385.

[5]胡宇，宗欣，于淼，等.我國仿制藥一致性評價政策環境現狀分析[J].中國藥物評價，2020，37（5）：321-326.

[6]STODDEN V，LEISCH F，PENG R D.Implementing Reproducible Research[M].Boca Raton：CRC Press，2014：1-29.

[7]EMMERT-STREIB F，DEHMER M，YLI-HARJA O.Ensuring Quality Standards and Reproducible Research for Data Analysis Services in Oncology： A Cooperative Service Model[J].Frontiers in Cell and Developmental Biology，2019，7（349）：1-5.

[8]隋鑫，王華偉，李曉斌，等.鹽酸特比萘芬片在中國健康受試者體內的生物等效性研究[J].中國臨床藥理學雜志，2021，37（2）：99-103.

[9]國家藥品監督管理局.總局關于發布藥物臨床試驗數據管理與統計分析的計劃和報告指導原則的通告（2016年第113號）[EB/OL].（2016）.[2021-04-22]. https：//www.nmpa.gov.cn/directory/web/nmpa/xxgk/ggtg/ qtggtg/20160729184001935.html

[10]王暉，崔暢，楊雙，等.琥珀酸索利那新片在中國健康受試者的生物等效性研究[J].中國臨床藥理學與治療學，2021，26（3）：299-304.

[11]魏敏吉，單愛蓮，趙藝璇.不依賴專業軟件的藥代動力學參數的計算方法 [J]. 中國臨床藥理學雜志，2016，32（23）：2199-2202.

[12]CANADA G O. Guidance Document： Conduct and Analysis of Comparative Bioavailability Studies[EB/OL].（2018）.[2021-04-18].https：//www.canada.ca/en/health-canada/services/drugs-health-products/drug-products/applicationssubmissions/guidance-documents/bioavailabilitybioequivalence/conduct-analysis-comparative.html.

[13] 沈健 .ADaM應用于臨床試驗的研究與探討 [D]. 南京：東南大學，2017：1-118.

[14] 王崇旭 . 臨床試驗數據管理電子化和 CDISC 數據自動化的實現 [D]. 南京：東南大學，2018：1-97.

[15] 朱天一，蔡晶晶，黃家俊，等 .ADaM在仿制藥生物等效性試驗安全性分析中的應用 [J]. 中國臨床藥理學與治療學，2020，25（6）：658-63.

[16] 吳正宇，方文韜，欒迪，等 .R 語言計算藥動學參數及生物等效性分析 [J]. 中國現代應用藥學，2020，37（3）：308-313.

[17] 王建峰 . Ⅰ期臨床藥物試驗統計分析報告規范化圖表輸出 [D]. 南京：東南大學，2018：1-8.

（收稿日期：2021-06-21）