孔 亮 朱文瀟 劉 蓓 陳龍龍 曹向陽* 崔宏勛 王向陽

1 湖南中醫(yī)藥大學(xué) （湖南 長沙 410208）

2 河南省洛陽正骨醫(yī)院（河南省骨科醫(yī)院） （河南 洛陽 471000）

3 河南中醫(yī)藥大學(xué) （河南 鄭州 450000）

內(nèi)容提要： R軟件是一種功能強(qiáng)大且免費(fèi)的統(tǒng)計(jì)分析及制圖軟件，其為統(tǒng)計(jì)計(jì)算提供了靈活而強(qiáng)大的基礎(chǔ)，并且免費(fèi)開放給更多的使用者?；赗軟件對統(tǒng)計(jì)數(shù)據(jù)進(jìn)行Meta分析的程序包種類多樣，不僅可以適應(yīng)各種統(tǒng)計(jì)計(jì)算需求，更能滿足個(gè)性化的制圖需求，以高標(biāo)準(zhǔn)地完成Meta分析。

R軟件是一種功能強(qiáng)大且免費(fèi)的統(tǒng)計(jì)分析及制圖軟件，它類似于編程語言可自由地實(shí)現(xiàn)各種數(shù)據(jù)處理需求，而又區(qū)別于編程語言繁復(fù)的實(shí)現(xiàn)方式，其通過加載程序包的方式達(dá)到種類良多的擴(kuò)展功能，在自由及簡便之間尋找微妙的平衡[1]。隨著循證醫(yī)學(xué)的不斷發(fā)展與進(jìn)步，Meta分析也逐漸被公眾所接納認(rèn)可，其通過定量合成分析同類研究證據(jù)指標(biāo)的方法來客觀評價(jià)某一問題的解決方法，成為循證決策的最佳依據(jù)。通過R軟件對統(tǒng)計(jì)數(shù)據(jù)進(jìn)行Meta分析的程序包有Metafor、Meta等，不同程序包的操作方式不盡相同，本次以闡述Meta程序包的應(yīng)用方法為主。軟件環(huán)境以R（v.3.4.1）、R studio（1.1.463）版本為例進(jìn)行講解。

1.R軟件、運(yùn)行環(huán)境、Meta程序包的下載、安裝及加載

R軟件及其運(yùn)行環(huán)境（R studio）建議一并安裝使用，R studio的安裝可提供更為簡潔美觀的窗口展示、更高效的操作錄入方式等功能[2]。

R軟件的擴(kuò)展功能豐富，但為保證軟件本身的輕便、快速、更新等優(yōu)勢，使用者需在首次使用時(shí)下載相關(guān)程序包進(jìn)行安裝，并在每次使用時(shí)進(jìn)行加載操作。得益于遍布全球各地的數(shù)十個(gè)CRAN鏡像服務(wù)站，用戶可通過鍵入少許指令，從而快捷下載安裝各類擴(kuò)展程序包。

安裝完成后，雙擊R studio便可來到R的交互式窗口，通過引導(dǎo)界面創(chuàng)建新的項(xiàng)目。R的窗口整體可分為4大部分，最上方為選項(xiàng)菜單欄，可完成程序本身的各項(xiàng)基本操作和設(shè)置；左側(cè)為Console，為軟件運(yùn)行的主要操控臺(tái)；右上方包括Environment、History，可檢閱軟件運(yùn)行過程中定義的各類數(shù)據(jù)及歷史操作記錄；右下方為File、Plots、Packages、Help，可實(shí)現(xiàn)文件管理、圖像查看、加載及卸載已安裝程序包、查詢函數(shù)幫助的功能。

R軟件的操作以鍵入命令行的形式完成，于命令提示符“＞”后書寫操作指令install.packages(“meta”)即可完成Meta程序包的安裝，安裝完成后需鍵入library(“meta”)以完成Meta程序包的加載。完成上述操作后，即完成使用R語言Meta程序包行二變量薈萃分析的基本軟件設(shè)置。

2.數(shù)據(jù)的加載與處理

待使用的數(shù)據(jù)應(yīng)先進(jìn)行版式及格式調(diào)整，轉(zhuǎn)換為R軟件容易調(diào)取的格式。本文以“明膠海綿在經(jīng)皮椎體成形術(shù)中預(yù)防骨水泥滲漏的Meta分析”為例講解利用R軟件進(jìn)行Meta分析的基本操作方法。

在科研工作中，最為廣泛使用的數(shù)據(jù)統(tǒng)計(jì)及整理軟件是Excel，但其xlsx格式并不能直接為R所用。因此，需預(yù)先將xlsx格式的數(shù)據(jù)文件另存為文本文件（制表符分隔）（*.txt）格式。

在R-Console中鍵入x=read.delim("h:/leakage.txt")，由此在R中建立了一個(gè)名為“x”的9行6列的矩陣。

接下來鍵入metaleakage＜-metabin(Events.E,Total.E,Events.C,Total.C,data = x,method = "MH",sm="OR",comb.fixed = T,comb.random = T,study=paste(x$Author(s), x$Year)，以計(jì)算二變量資料的合成結(jié)果（metabin函數(shù)）并賦值給metaleakage，其中Events.E, Total.E, Events.C, Total.C分別為實(shí)驗(yàn)組的發(fā)生次數(shù)/樣本總量，對照組的發(fā)生次數(shù)/樣本總量，data為metabin函數(shù)使用數(shù)據(jù)的來源，x為已經(jīng)導(dǎo)入的矩陣數(shù)據(jù)，sm為合并統(tǒng)計(jì)量的形式{"RR", "OR", "RD", or "ASD" }，comb. fixed及comb.random為效應(yīng)模型的選擇{"固定"，"隨機(jī)"}，studlab為指向矩陣中標(biāo)簽信息的向量，當(dāng)指向多個(gè)標(biāo)簽信息時(shí)應(yīng)借用paste函數(shù)。若不同標(biāo)簽信息需要間隔時(shí)，可在paste中內(nèi)嵌sep函數(shù)，以逗號為例：sep=“,”。

3.森林圖的繪制與編輯

forest.meta(metaleakage,family=("sans"),fontsize=9.5,lab.e="實(shí)驗(yàn)組",lab.c="對照組",lwd=2,col.diamond.fixed="lightslategray",col.diamond.lines.fixed="lightslategray",col.diamond.random="maroon",col.diamond.lines.random="maroon",col.square="skyblue",col.study="lightslategray",lty.fixed=4,plotwidth="8cm",colgap.forest.left="1cm",colgap.forest.right="1cm",just.forest="right",colgap.left="0.5cm",colgap.right="0.5cm")。

圖1. Meta分析森林圖

森林圖（見圖1）的繪制需要采用forest.meta函數(shù)，將上述步驟處理好的Meta分析處理后的數(shù)據(jù)載入，根據(jù)作者需求，對便簽名稱、文字大小及字體、圖形位置與色彩進(jìn)行設(shè)定即可。其中family為字體設(shè)定函數(shù){“serif”-襯線、“sans”-無襯線、“mono”-等寬}，fontsize為字體大小設(shè)定函數(shù)，lab.e、lab.c為對應(yīng)分組的展示標(biāo)簽，lwd為圖形線段的粗細(xì)。此外，col為顏色定義函數(shù)，其中col.diamond. fixed和col.diamond.random可分別定義固定效應(yīng)模型及隨機(jī)效應(yīng)模型合并統(tǒng)計(jì)量的符號指代，col.square、col.study分別為研究個(gè)體效應(yīng)量（方塊）、可信區(qū)間（線段）的顏色定義函數(shù)。lty為線條虛實(shí)樣式定義函數(shù){“0”-不畫線、“1”-實(shí)線、“2”-虛線、“3”-點(diǎn)線、“4”-點(diǎn)劃線、“5”-長劃線、“6”-點(diǎn)長劃線}，plotwidth為圖形整體寬度的定義函數(shù)。Colgap為森林圖中各個(gè)區(qū)塊之間間隔的定義函數(shù)，此操作較為繁瑣，需多次嘗試方可尋找最佳數(shù)值，其中colgap.forest.left和colgap.forest.right分別為森林圖和左右數(shù)據(jù)信息間隔寬度，colgap.left和colgap.right為各列數(shù)據(jù)信息之間的間隔寬度。just.forest為各列數(shù)據(jù)信息對齊的定義函數(shù){“l(fā)eft”-左對齊、“right”-右對齊、“center”-居中}。

4.討論

隨著循證醫(yī)學(xué)的不斷發(fā)展，其遵循證據(jù)的思想逐漸深入人心，從而憑借薈萃多方數(shù)據(jù)，系統(tǒng)評價(jià)某一問題的Meta分析一度成為研究熱點(diǎn)，統(tǒng)計(jì)計(jì)算相關(guān)數(shù)據(jù)的軟件及方法層出不窮[3]。其中R軟件以開源及自由等優(yōu)勢獲得廣大用戶的歡迎。為避免操作過于復(fù)雜，給初學(xué)者帶來困難，本文僅以Meta程序包進(jìn)行基礎(chǔ)示范。Meta程序包除本文所提及的命令語句之外仍有較多功能未提及，使用者可在R軟件搜索欄中檢索Meta程序包的幫助文件，以獲取高級使用教程，更好地完成Meta分析。