如何解決醫(yī)藥分析中的重復(fù)測量問題

添加日期：2010年11月27日 閱讀：1046

    重復(fù)測量(Repeated Measures)是醫(yī)學(xué)領(lǐng)域十分常見的一種情形，具體是指對同一研究對象的同一個測量指標(biāo)或同一組測量指標(biāo)，在不同時間點上進(jìn)行多次重復(fù)測量，以用于分析測量指標(biāo)隨時間變化而產(chǎn)生的變化特點。這類測量資料在臨床和流行病學(xué)研究中比較常見。例如，為研究各種新藥對高血壓病人的治療效果，需要在服藥后一小時、二小時、四小時、六小時等的時點分別測量受試者的血壓，以分析其血壓的變動情況。醫(yī)學(xué)研究是應(yīng)用統(tǒng)計分析方法*多的領(lǐng)域之一，統(tǒng)計分析方法也能幫助很好地解決上述問題。 
    粗看一下，很多人可能覺得運(yùn)用t檢驗或方差分析作各測量時間點上的差別檢驗就能較好地解決上述問題。但實際上這樣的處理方法是錯誤的，很可能會得出錯誤的結(jié)論。原因是應(yīng)用傳統(tǒng)的t檢驗或方差分析有一個前提條件，即數(shù)據(jù)的獨立性，而醫(yī)學(xué)領(lǐng)域的重復(fù)測量數(shù)據(jù)往往存在高度的相關(guān)性。例如張三服藥二小時后的舒張壓是90，那么四小時后他的舒張壓也很可能在90附近，不會差得太遠(yuǎn)，這就是數(shù)據(jù)之間的相關(guān)性。如果此時去測李四，他的舒張壓可能是70到110之間的任意值，因為這個數(shù)值與張三二小時后的舒張壓值沒有任何關(guān)系。因此，當(dāng)t檢驗或方差分析的前提條件被破壞時，我們必須使用一套新的統(tǒng)計分析方法——重復(fù)測量的方差分析來校正原先的結(jié)果。這時，我們必須借助高級統(tǒng)計分析軟件來實現(xiàn)，下文的分析中，我們講借助全球高端六西格瑪質(zhì)量管理軟件JMP進(jìn)行。順便提一下，JMP軟件是全球?qū)ｉT面向六西格瑪和質(zhì)量管理的高端統(tǒng)計分析軟件，是FDA(美國食品藥品管理局)**指定的專業(yè)統(tǒng)計軟件之一，已為全球絕大多數(shù)知名藥企用于質(zhì)量改進(jìn)和藥品研發(fā)，客戶包括拜爾、羅氏、禮來、諾和諾德、阿斯利康、輝瑞、葛蘭素史克、衛(wèi)材制藥、武田制藥、賽諾菲安萬特等等。
一般來說，重復(fù)測量的方差分析可分為三個主要步驟： 
    第*， 球形性檢驗 
    這里的球形性指的是協(xié)方差矩陣的球形性質(zhì)，即矩陣的主對角線元素（方差）相等、非主對角線元素（協(xié)方差）為零。常用Mauchly氏法檢驗球形性。當(dāng)球形性沒有得到滿足時，必須對與時間有關(guān)的F統(tǒng)計量的分子、分母自由度進(jìn)行調(diào)整，以便減少犯第*類錯誤的概率。 
    第二， 方差分解及F檢驗 
    方差分解及F檢驗的概念相信大家都有所了解，但與傳統(tǒng)方差分析不同的是，重復(fù)測量數(shù)據(jù)的誤差結(jié)構(gòu)分成兩個層次：一層誤差來自不同時間點的測量, 另一層誤差來自不同的觀察對象。而每一層的誤差可能與若干個解釋變量有關(guān), 從而構(gòu)成了相對復(fù)雜的誤差結(jié)構(gòu)。一般用來自不同觀察對的誤差檢驗組間效應(yīng)（如不同藥物），用來自不同時間點的誤差分別檢驗組內(nèi)效應(yīng)（如不同時間）和組間效應(yīng)與組內(nèi)效應(yīng)之間的交互作用。 
    第三， 繪制輪廓圖 
    除了定量分析之外，還可以用可視化的方式（如輪廓圖）來繪制測量指標(biāo)在組間變量的不同水平上、不同測量時間點時的變化曲線（具體形式參見下文），以此反映組間效應(yīng)的顯著程度和指標(biāo)隨時間變化的趨勢。測量時間間隔可以是等距的，也可以是不等距的。 
    下面用一個典型的生物醫(yī)學(xué)實例來具體詮釋一下“重復(fù)測量的方差分析”的整個過程，在涉及具體的統(tǒng)計計算時，一律采用JMP進(jìn)行。采用JMP的原因，一方面是確保分析的**性；另一方面是因為下面使用的功能都是其他同類軟件所不具備的；另外JMP在分析結(jié)果的可視化、分析能力等方面較同類軟件相比都有很大的優(yōu)勢。 
    案例分析 一項藥物代謝動力學(xué)研究，目的是對比某種藥物的不同劑型在體內(nèi)的代謝速度。劑型分膠囊型和片劑型。將16名受試對象隨機(jī)分為兩組，每組8名。一組給予膠囊，另一組給予片劑，分別在服藥后1小時、2小時、4小時、6小時及8小時測定血中的藥物濃度。測定結(jié)果如以下圖一所示，試選擇合適的數(shù)理統(tǒng)計方法加以分析。 
        責(zé)任編輯：小季    m.atm-sprinta.com    2010-11-27 15:24:23

文章來源：