多變量重復測量方差分析統計報表輸出的SAS宏實現

甘秀敏 蔡 強 聶紹發△

重復測量是指對同一觀察對象的同一觀察指標在不同時間點上進行多次測量，用于分析觀察指標在不同時點上的變化規律。重復測量的數據在醫學研究中非常常見，而重復測量方差分析是用以處理重復測量連續性數據主要的統計分析方法，為了更準確、方便地得到該統計方法所得結果的統計報表，本文通過以GLM廣義線性模型為基礎，編寫SAS宏程序實現多變量重復測量方差分析統計報表的輸出。

方法與步驟

1.SAS宏編寫的主要方法及步驟

通過SAS軟件中的GLM過程步對重復測量的數據進行分析，采用結果輸出系統(output delivery system，ODS)技術，并將所有統計分析結果根據其類型儲存為相應的SAS數據庫，包含進行球形性檢驗的數據庫(sphericity)和進行組間、組內比較結果的數據庫(model ANOVA)。依據球形性檢驗的結果在model ANOVA中選出所需要的統計量及是否校正過的P值〔1〕，并對挑選后的結果進行標準格式的編輯轉換，進而得到輸出統計報表所需的標準化SAS數據庫，最后采用ODS技術及Report過程步將最終的數據庫規范化地輸出到Word文檔中。

2.SAS宏程序參數設置及注釋

將統計分析的關鍵參數轉化為SAS宏的參數，在宏中按照相應的規則進行SAS程序的編寫，得到所需要的最終宏程序。相關宏參數及其意義解釋見表1，SAS程序及相關解釋見文后附錄。

表1 宏參數名稱及注釋

結 果

SAS宏程序的應用實例:欲比較兩種治療口腔粘膜炎導致口腔粘膜疼痛藥物的療效情況，將200例病人隨機分成兩組，分別接受兩種治療，在5個醫院同時進行多中心藥物臨床試驗，并分別在用藥前、用藥后每兩天測定一次，共測定5次。具體宏參數根據實際的數據情況分別賦值。

對上述宏變量賦值后運行宏程序，即在路徑為“C:Documents and Settings”的文件夾下自動生成rtf格式的REPORT文件，統計結果報表也被規范化輸出到REPORT文件中(表2)。

表2 重復測量方差分析結果(1)

在藥物試驗過程中也經常遇到一些單中心的臨床試驗，這時只需將宏變量“modelright”賦值改為“%let modelright=%str(group center);”，即可靈活地得到需要的統計表格(見表3)。

表3 重復測量方差分析結果(2)

討 論

在使用SAS軟件進行數據統計分析時，SAS軟件會給出大量的統計分析結果，而根據實際需要，研究者往往僅需其中的一小部分統計分析結果。為了正確挑選出統計結果，并輸出成固定格式的統計報表，需要研究者有深厚的統計學基礎及細心的操作，否則極易在重復多次的統計工作中出現錯漏。

SAS宏程序自動產生統計報表可使統計工作變得相對方便、快捷、準確，目前已得到越來越多的應用，特別是在藥物臨床試驗的數據統計中能夠減少不必要的重復工作進而極大的提高統計工作的效率〔2－4〕。本文采用GLM過程步對重復測量的數據進行了多變量的方差分析，所編寫的SAS宏程序能夠靈活地應對分析模型中自變量及應變量個數的動態變化，并使用ODS技術對統計結果進行了處理和自動輸出，可以在很大程度上減少因人為重復操作而導致的統計報告錯漏的發生，保證統計輸出結果的準確性，使重復多次分析的工作量得到大大的降低，可為廣大衛生工作者提供一個有力的統計分析工具。

1．Glenn A．Walker．Common Statistical Methods for Clinical Research with SAS Example．Second Edition．Cary，NC:SAS Institute Inc．，2002．

2．童新元，張高魁，姚晨．定性指標的SAS統計分析報表——SAS軟件在新藥臨床試驗統計分析中的應用(一)．中國衛生統計，2003，20(1):50-51．

3．田曉燕，金丕煥，張文彤．應用SAS軟件的動態數據交換技術自動生成統計報告．中華預防醫學雜志，2005，39(1):54-55．

4．薛富波，于莉莉，李嬋娟，等．新藥臨床試驗中log-rank檢驗的SAS宏實現．中國臨床藥理學與治療學，2006，11(2):227-231．