本文關鍵詞：R軟件nlme程序包在網狀Meta分析中的應用，由筆耕文化傳播整理發(fā)布。

中國循證醫(yī)學雜志 2014, 14(3): 355～360

論 著

方法學研究

R軟件nlme程序包在網狀Meta分析中的應用

張 超1 牛玉明2 曾憲濤1,2*

1. 湖北醫(yī)藥學院附屬太和醫(yī)院循證醫(yī)學中心（湖北十堰 442000）；2. 湖北醫(yī)藥學院附屬太和醫(yī)院口腔醫(yī)學中心（湖北十堰 442000）

摘要 nlme程序包是基于廣義最小二乘法和線性混合效應模型研發(fā)的、可通過R軟件實現廣義線性和非線性混合效應模型下的Meta分析。該程序包實現Meta分析時，需要對數據先行轉化為效應量的對數值才可進行。本文介紹了使用R軟件nlme程序包實現網狀Meta分析的過程，詳細呈現了如何轉化數據這一核心步驟。

關鍵詞 網狀Meta分析；廣義線性混合效應模型；R語言；nlme程序包

Implementation of Network Meta-Analysis with nlme Package in R Soft ware

ZHANG Chao1, NIU Yu-ming2, ZENG Xian-tao1,2*

1. Center for Evidence-Based Medicine and Clinical Research, Taihe Hospital, Hubei University of Medicine, Shiyan 442000, China;2. Department of Stomatology, Taihe Hospital, Hubei University of Medicine, Shiyan 442000, China

Abstract Th e nlme package is developed based on the generalized least squares (gls) and linear mixed-eff ects model (lme). It can perform meta-analysis based on linear and nonlinear mixed eff ects models in R language. When conducting meta-analysis using nlme package in R language, the fi rst step is to translate the data into its logarithm estimation. In this article, we introduce how to perform network meta-analysis using R language nlme package and show the core step of data translation in detail.

Key words Network meta-analysis; Generalized linear mixed model; R language; nlme package

隨著網狀Meta分析的發(fā)展及各種軟件的成功研發(fā)，運算模型的選擇也日益受到關注，廣義線性模型（generalized linear model，GLM）就是其中之一。GLM是線性模型的擴展，其特點是不強行改變數據的自然度量，數據可以具有非線性和非恒定方差結構；與線性模型相比，GLM模型中Y的分布可以是任何形式的指數分布（如高斯分布、泊松分布、二項式分布），聯結函數可以是任何單調可微函數（如對數函數logarithm或邏輯函數logit），這些優(yōu)點使得GLM模型可處理多種變量，日益受到關注（欲詳細了解GLM的讀者建議閱讀參考文獻 [1,2]）。

nlme是一款基于S語言在S-PLUS軟件中實現混合效應模型（mixed-eff ects models）分析的軟件。因此，其可以同時進行線性和非線性混合效應模型（linear and nonlinear mixed effects models，NLME）分析 [3,4]。當前，混合效應模型也可通過nlme程序

DOI: 10.7507/1672-2531.20140060基金項目：湖北省教育科學“十二五”規(guī)劃2012年度重點課題（編號：2012A050），湖北醫(yī)藥學院2011年度優(yōu)秀中青年科技創(chuàng)新團隊項目（編號：2011CZX01），湖北醫(yī)藥學院附屬太和醫(yī)院2013年循證醫(yī)學“苗圃基金”項目（編號：EBM2013004）。作者簡介：張超，男（1988年～），本科，，住院醫(yī)師，以循證醫(yī)學、系統(tǒng)評價與Meta分析為主要研究方向。Email: zhangchao0803@*通訊作者，Email: zengxiantao1128@

包在R語言中實現 [5,6]。nlme程序包亦可以實現網狀Meta分析 [6,7]，本文仍以《R軟件R2WinBUGS程序包在網狀Meta分析中的應用》一文中實例 [8,9]為例進行展示。

1 軟件及程序包的安裝/加載

R軟件的安裝已在前文講述，本文使用的是R-3.0.1 [8,10,11]。

此外，需要安裝nlme程序包，具體命令為：in-stall.packages(“nlme”)。在彈出的對話框中選擇某個鏡像（CRAN）安裝，安裝完成后再由library(“nlme”)命令完成加載。

2 數據的加載與預處理

2.1 數據加載

R軟件nlme程序包實現網狀Meta分析的基本思路是將因變量y設置為效應量，自變量x設置為干預措施，結合相關的干預措施之間的比較來計算得出線性關系 [7]。本文將y設為logOR，x設為13種藥物以及安慰劑，并對相關藥物之間的比較作一定標記。

因實例中含有3臂試驗，所以，首先將數據中

CJEBM

355

© 2014 中國循證醫(yī)學雜志編輯部

  本文關鍵詞：R軟件nlme程序包在網狀Meta分析中的應用，由筆耕文化傳播整理發(fā)布。