用SAS軟件實(shí)現(xiàn)高維列聯(lián)表資料的統(tǒng)計(jì)學(xué)分析_CMH_2檢驗(yàn)_胡良平

1、藥學(xué)服務(wù)與研究； （） 胡良平，等用軟件實(shí)現(xiàn)高維列聯(lián)表資料的統(tǒng)計(jì)學(xué)分析：檢驗(yàn)專題講座用軟件實(shí)現(xiàn)高維列聯(lián)表資料的統(tǒng)計(jì)學(xué)分析：檢驗(yàn)胡良平，沈 寧（軍事醫(yī)學(xué)科學(xué)院研究生部生物醫(yī)學(xué)統(tǒng)計(jì)學(xué)咨詢中心，北京 ）摘 要統(tǒng)計(jì)分析是在 統(tǒng)計(jì)分析方法的基礎(chǔ)上發(fā)展并提出的，用于分層分析，即控制非試驗(yàn)因素后對(duì)保留下來(lái)的由原因變量與結(jié)果變量構(gòu)成的、二維列聯(lián)表資料的統(tǒng)計(jì)學(xué)處理。這種方法可以控制一個(gè)或多個(gè)原因變量，而考察關(guān)注的一個(gè)原因變量對(duì)結(jié)果變量的影響。本文結(jié)合實(shí)例，對(duì)高維列聯(lián)表資料和用軟件進(jìn)行統(tǒng)計(jì)學(xué)處理的方法進(jìn)行了介紹，并給出了結(jié)果解釋。關(guān)鍵詞 統(tǒng)計(jì)學(xué)；檢驗(yàn)；實(shí)現(xiàn)；統(tǒng)計(jì)分析中圖分類號(hào)，文獻(xiàn)標(biāo)志碼文章編號(hào)（）：，（，） （

2、）， ；，（）：在藥學(xué)和臨床研究中，常涉及單變量定性資料的研究，例如分析在晚期非小細(xì)胞肺癌病人中，以不同藥物處理所得獲益病人數(shù)，以此對(duì)比各種藥物的療效。當(dāng)研究結(jié)果變量為定性變量時(shí)，研究者習(xí)慣將資料整理為列聯(lián)表形式：列聯(lián)表、列聯(lián)表和高維列聯(lián)表，后者可以根據(jù)結(jié)果變量的類型分為：結(jié)果為二值變量的高維列聯(lián)表，如考察病人的生存與否（是、否）；結(jié)果為多值有序變量的高維列聯(lián)表，如考察在不同藥物干預(yù)后病人疾病進(jìn)展情況（治愈、穩(wěn)定、進(jìn)展、惡化）；以及結(jié)果為多值名義變量的高維列聯(lián)表，如考察不同抗腫瘤藥物用于不同分型（腸型、彌漫型、混合型 ）病人的人數(shù)構(gòu)成情況。本文將結(jié)合實(shí)例介紹應(yīng)用軟件對(duì)結(jié)果變量為二值變量定性資料

3、實(shí)施檢驗(yàn)。效差異。評(píng)價(jià)指標(biāo)為接受不同方案治療有效的病人人數(shù)，有效定義為完全緩解部分緩解穩(wěn)定，無(wú)效定義為疾病進(jìn)展。在考察過(guò)程中，研究者發(fā)現(xiàn)除藥物外，對(duì)療效可能有影響的因素還有年齡（歲，歲）和病人是否復(fù)治，表 為原作者的調(diào)查表格。表兩組病人化療后有效的人數(shù) 年齡有效人數(shù)治療狀態(tài)（歲）對(duì)照組治療組（）（）初治（）（）（）復(fù)治（）（）（） 呈現(xiàn)實(shí)例以晚期非小細(xì)胞肺癌病人作為研究對(duì)象進(jìn)行的回顧性研究中，分析治療組（多西他賽奧沙利鉑，）與對(duì)照組 （多西他賽 順鉑，）的療作者簡(jiǎn)介胡良平（男），碩士，研究員：本文在參考文獻(xiàn)表的基礎(chǔ)上進(jìn)行了修改，數(shù)據(jù)不代表原作者的結(jié)果 分析思路本研究的分析思路為：?jiǎn)巫兞糠治?，?/p>

4、涉及兩個(gè)變量的回歸分析。結(jié)果變量為治療是否有效，屬于定性資料，也稱為結(jié)果變量為二值（是、否）變量的三維列聯(lián)表資料。設(shè)計(jì)類型為類似于兩因素析因設(shè)藥學(xué)服務(wù)與研究； （） 胡良平，等用軟件實(shí)現(xiàn)高維列聯(lián)表資料的統(tǒng)計(jì)學(xué)分析：檢驗(yàn)計(jì)，因素包括藥物干預(yù)、年齡、既往接受治療的狀態(tài)。據(jù)前述分析可以選擇的方法有：加權(quán) 檢驗(yàn)（限三維列聯(lián)表）、檢驗(yàn)、回歸分析和對(duì)數(shù)線性模型分析。本文僅介紹 檢驗(yàn)的表 軟件程序編寫(xiě)及說(shuō)明 軟件輸入過(guò)程說(shuō)明實(shí)現(xiàn)和結(jié)果解釋。 原始數(shù)據(jù) 表向 準(zhǔn)型表格的 根據(jù)資料類型以及要選擇的統(tǒng)計(jì)學(xué)方法，將研究者的原始表轉(zhuǎn)化為用于統(tǒng)計(jì)學(xué)分析的標(biāo)準(zhǔn)型（見(jiàn)表）。在轉(zhuǎn)化表格的過(guò)程中，應(yīng)該把因素不同水平組；數(shù)據(jù)步開(kāi)

5、始：建立數(shù)據(jù)集并命名為；、和分別代表年齡、復(fù)治情況、治療方案和療效；變量為頻數(shù)，即指示輸入的數(shù)據(jù)為頻數(shù)；合下二值結(jié)果變量各水平下的頻數(shù)全部呈現(xiàn)出來(lái)。表結(jié)果變量為二值變量的四維列聯(lián)表資料的標(biāo)準(zhǔn)型 年齡復(fù)治治療方案病人人數(shù)（例）分層與否（藥物配伍）有效無(wú)效歲 是 多西他賽 奧沙利鉑（治療組）； ；開(kāi)始輸入數(shù)據(jù)；按照表中數(shù)據(jù)出現(xiàn)的次序，輸入數(shù)據(jù)；調(diào)用過(guò)程；第二個(gè) 語(yǔ)句后面列出個(gè)變量，分別為個(gè)原因變量（年齡和復(fù)治情況）和個(gè)結(jié)果變量（治療效果）；程序結(jié)束。多西他賽順鉑（對(duì)照組）否 多西他賽奧沙利鉑（治療組）多西他賽順鉑（對(duì)照組）歲 是 多西他賽奧沙利鉑（治療組）多西他賽順鉑（對(duì)照組）否 多西他賽奧沙利鉑

6、（治療組）多西他賽順鉑（對(duì)照組） 分析方法的 本研究涉及了 個(gè)原因變量：年齡、復(fù)治與否、治療方案；結(jié)果變量為個(gè)，即療效（以獲益人數(shù)表示），所以該資料屬于結(jié)果變量為二值變量的四維列聯(lián)表資料。研究者的主要目的是研究不同治療方案對(duì)療效有何影響，所以以上個(gè)原因變量中的年齡和復(fù)治情況，就成為重要非試驗(yàn)因素。研究者最關(guān)心的是“治療方案”對(duì)療效的影響。當(dāng)只想消除或減弱非試驗(yàn)因素對(duì)結(jié)果的影響，重點(diǎn)考察試驗(yàn)因素的試驗(yàn)效應(yīng)時(shí)，可以選擇加權(quán)檢驗(yàn)和 檢驗(yàn)，排除其他非試驗(yàn)因素的影響，觀察不同治療方案的療效之間的差別是否具有統(tǒng)計(jì)學(xué)意義。 用 件分析表中的 料 程序編寫(xiě) 程序編寫(xiě)由數(shù)據(jù)步和過(guò)程步組成。在數(shù)據(jù)步中， 表示 歲

7、，表示 歲；表 示 復(fù) 治， 表示初治； 表示多西他 賽奧 沙利鉑，表示多西他賽順鉑；表示治療有效，表示治療無(wú)效，見(jiàn)表。 結(jié)果解釋 第部分結(jié)果（圖、表 、表 ）是不考慮年齡和復(fù)治情況的基礎(chǔ)上，考察藥物方案與療效之間的關(guān)系。由圖結(jié)果可見(jiàn)，多西他賽奧沙利鉑組的有效率（）較多西他賽順鉑組的有效率（）略高，但種檢驗(yàn)方法的結(jié)果均（見(jiàn)表）。表為不控制年齡和復(fù)治情況下藥物方案與療效間的關(guān)系（相對(duì)風(fēng)險(xiǎn)的估計(jì)值）。由表 結(jié)果可見(jiàn)上述四格表資料的優(yōu)勢(shì)比（）、相對(duì)危險(xiǎn)度（）及其 置信區(qū)間。本例為病例對(duì)照研究，應(yīng)該選擇第一行結(jié)果即案例對(duì)照（優(yōu)比），即 值為，置信區(qū)間的下限和上限分別為和。過(guò)程表一合計(jì)頻數(shù)百分比行百分比

8、列百分比合計(jì)圖不控制年齡和復(fù)治情況下藥物方案與療效間的關(guān)系：兩組頻數(shù)及百分比： 藥學(xué)服務(wù)與研究 ；（）表 “”表的統(tǒng)計(jì)量（不控制年齡和復(fù)治情況下藥物方案與療效間的關(guān)系） “”（）統(tǒng)計(jì)量自由度值概率卡方似然比卡方連續(xù)校正卡方卡方系數(shù)列聯(lián)系數(shù)統(tǒng)計(jì)量表 相對(duì)風(fēng)險(xiǎn)的估計(jì)值（控制年齡和復(fù)治情況下藥物方案與療效間的關(guān)系） （）普通相對(duì)風(fēng)險(xiǎn)的估計(jì)值（行行）研究類型方法值置信限案例對(duì)照 （優(yōu)比） （第列風(fēng)險(xiǎn)） （第列風(fēng)險(xiǎn)） 表 相對(duì)風(fēng)險(xiǎn)的估計(jì)值（不控制年齡和復(fù)治情況下藥物方案與療效間的關(guān)系） （）研究類型值置信限案例對(duì)照（優(yōu)比）（第列風(fēng)險(xiǎn)）（第列風(fēng)險(xiǎn)）第部分結(jié)果顯示在消除年齡和復(fù)治與否兩個(gè)因素后，治療方案對(duì)療

9、效的影響。 檢驗(yàn)的種具體方法，即非零相關(guān)統(tǒng)計(jì)量（），行均值得分差值統(tǒng)計(jì)量（）和一般聯(lián)系統(tǒng)計(jì)量（）得到 值均為，說(shuō)明兩種治療方案療效之間的差異有統(tǒng)計(jì)學(xué)意義。表為控制年齡和復(fù)治情況下藥物方案與療效間的關(guān)系。由表可見(jiàn)，在消除年齡和復(fù)治與否兩個(gè)因素后，用 法和法分別對(duì)兩種設(shè)計(jì)類型，即病例對(duì)照研究（）和隊(duì)列研究（）得到、 及置信區(qū)間結(jié)果。本研究為病例對(duì)照研究，選擇案例對(duì)照（優(yōu)比），得出結(jié)果，置信區(qū)間的上限和下限分別為和。在控制年齡和復(fù)治情況下，對(duì)各層優(yōu)勢(shì)比是否滿足 齊 性 的 檢 驗(yàn) 結(jié) 果 顯 示， ，即說(shuō)明各控制因素分層 （年齡和復(fù)治）之間效應(yīng)方向上不一致，即證明年齡因素與療法有交互作用，此時(shí)運(yùn)用檢

10、驗(yàn)得到的結(jié)果并不可靠，應(yīng)該慎重作出結(jié)論，最好選擇其他統(tǒng)計(jì)學(xué)方法（如回歸）進(jìn)行驗(yàn)證，此處的結(jié)果和結(jié)論僅供參考。統(tǒng)計(jì)學(xué)和專業(yè)結(jié)論 調(diào)整了年齡和復(fù)治情況之后，兩種治療方法療效之間的差異無(wú)統(tǒng)計(jì)學(xué)意義。 討 論 加 的 用條件和 不同的研究目的決定同樣的資料，可以選擇不同的統(tǒng)計(jì)學(xué)方法進(jìn)行分析。對(duì)于結(jié)果變量為二值變量的高維列聯(lián)表資料可以選擇的方法有加權(quán) 檢驗(yàn)、檢驗(yàn)、回歸分析和對(duì)數(shù)線性模型分析。加權(quán) 檢驗(yàn)和檢驗(yàn)的應(yīng)用條件是當(dāng)列聯(lián)表的多個(gè)因素不適合采用簡(jiǎn)單的 “合并 （壓縮 ）”的方法進(jìn)行處理，但又不想采用復(fù)雜的對(duì)數(shù)線性模型或 回歸模型進(jìn)行分析時(shí)，可以考慮這兩種方法。但是要注意這兩種方法不能考察被控制的原因變

11、量對(duì)結(jié)果變量的影響，并且加權(quán) 檢驗(yàn)只能應(yīng)用于三維列聯(lián)表，。值得一提的是，從現(xiàn)在軟件的輸出結(jié)果中可看出， 檢驗(yàn)基本上涵蓋了通常所說(shuō)的加權(quán) 檢驗(yàn)。 檢驗(yàn) 檢驗(yàn)法可以計(jì)算 種統(tǒng)計(jì)量，即非零相關(guān)統(tǒng)計(jì)量、行均值得分差值統(tǒng)計(jì)量和一般聯(lián)系統(tǒng)計(jì)量。它們之間的區(qū)別在于：非零相關(guān)統(tǒng)計(jì)量要求行變量（原因變量）與列變量（結(jié)果變量）均為有序變量，如行變量為藥物濃度（高、中、低），列變量為療效（好、中、差）。行均值得分差值統(tǒng)計(jì)量是當(dāng)行變量為無(wú)序變量（如民族可分為漢族、回族、滿族），而列變量為有序變量（如血紅蛋白濃度可分為高、中、低），計(jì)算該統(tǒng)計(jì)量可以獲得校正的（即消除某個(gè)或某些分層變量的影響，當(dāng)僅有一個(gè)定性原因變量時(shí)，就

12、是普通的秩和檢驗(yàn)結(jié)果） 秩和檢驗(yàn)結(jié)果。一般聯(lián)系統(tǒng)計(jì)量是當(dāng)行變量與列變量均為無(wú)序變量時(shí)計(jì)算該統(tǒng)計(jì)量，希望回答兩名義變量之間藥學(xué)服務(wù)與研究； （） 胡良平，等用軟件實(shí)現(xiàn)高維列聯(lián)表資料的統(tǒng)計(jì)學(xué)分析：檢驗(yàn)是否存在關(guān)聯(lián)性，用通俗的語(yǔ)言表達(dá)就是考察各行結(jié)合臨床，結(jié)合統(tǒng)計(jì)和專業(yè)知識(shí)，對(duì)結(jié)果給出合理和上的頻數(shù)分布規(guī)律是否相同。需要注意的是，當(dāng)結(jié)謹(jǐn)慎的解釋，作出令人信服且經(jīng)得起考驗(yàn)的結(jié)論。果變量為二值變量時(shí)，得出的以上種統(tǒng)計(jì)量的數(shù)【參考文獻(xiàn)】值結(jié)果一致。檢驗(yàn)在結(jié)果部分的解釋中， 裘光賢，周兆春，王 婕 多西他賽聯(lián)合奧沙利鉑治療晚期非小細(xì)胞肺癌的療效和安全性藥學(xué)服務(wù)與研究，檢驗(yàn)旨在考察不同分層中，治療組與， （）

13、：對(duì)照 組 療 效 之 差 是 否 一 致。 在 本 研 究 中，如，在歲和歲的兩個(gè)層中，考察各層中是否均有治療組療效優(yōu)于劣于對(duì)照組，或在歲與歲層中 治療組與對(duì)照組中療效，（）： ，優(yōu)劣的方向是否一致。如果方向不一致，則證明年 胡良平 統(tǒng)計(jì)學(xué)三型理論在實(shí)驗(yàn)設(shè)計(jì)中的應(yīng)用 北京：人齡因素與療法有交互作用，這時(shí)運(yùn)用檢驗(yàn) 民軍醫(yī)出版社，：得到的結(jié)果是不可靠的，應(yīng)該考慮選擇其他統(tǒng)計(jì)學(xué)方法。所以，該檢驗(yàn)其實(shí)是對(duì)檢驗(yàn)結(jié)果的 ：可靠性進(jìn)行補(bǔ)充說(shuō)明 。，：綜上，在進(jìn)行單變量分析時(shí)，對(duì)變量性質(zhì)進(jìn)行準(zhǔn) 胡良平醫(yī)學(xué)統(tǒng)計(jì)學(xué) 運(yùn)用三型理論分析定量與定性資料 北京：人民軍醫(yī)出版社，：確判斷，對(duì)試驗(yàn)設(shè)計(jì)類型清晰了解，根據(jù)試驗(yàn)

14、目的選 擇合適的統(tǒng)計(jì)學(xué)方法是必要的。對(duì)于初學(xué)者，將原 ：始的試驗(yàn)記錄表格，清晰地整理為資料的標(biāo)準(zhǔn)型，對(duì) ，：檢驗(yàn)分析中，應(yīng)收稿日期修回日期于統(tǒng)計(jì)學(xué)分析尤其重要。 掌握數(shù)據(jù)資料的類型，應(yīng)用軟件操作時(shí)要注重本文編輯 陽(yáng)凌燕檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪（頁(yè)）定、有無(wú)副作用發(fā)生，此后如繼續(xù)服用，每個(gè)月就診一次。上接第口服胺碘酮也可導(dǎo)致肝炎或肝硬化，其發(fā)生率。胺碘該例患者未出現(xiàn)相關(guān)臨床癥狀，僅在實(shí)驗(yàn)室檢查中發(fā)現(xiàn)異酮及其代謝產(chǎn)物去乙基胺碘酮可以明顯抑制，胺碘常，因此建議臨床醫(yī)師在使用胺碘酮初期應(yīng)提高監(jiān)測(cè)頻率。酮含有的苯并呋喃環(huán)可破壞線粒

15、體氧化和電子傳遞過(guò)程 是胺碘酮及其代謝產(chǎn)物可以明顯抑制，而血脂康是通，胺碘酮引起肝臟毒性的主要結(jié)構(gòu)。其肝損害的機(jī)制可能過(guò)代 謝，建 議 老 年 患 者 在 同 時(shí) 服 用 這 兩 種 經(jīng)為：（）胺碘酮可引起肝細(xì)胞線粒體毒性，主要通過(guò)影響氧代謝的藥物時(shí)，需考慮藥物的相互作用，不宜大劑化及氧化磷酸化的解偶聯(lián)，從而抑制電子傳遞鏈上的復(fù)合量給藥。酶。（）大量活性氧的產(chǎn)生加速細(xì)胞凋亡。胺碘酮在抑制電【參考文獻(xiàn)】子傳遞鏈的同時(shí)可產(chǎn)生大量的活性氧，使得線粒體通透性轉(zhuǎn)換孔道開(kāi)放，繼而引起線粒體膜上的細(xì)胞色素 釋放，觸發(fā)黃興福，楊艷敏，朱 俊 靜脈胺碘酮與嚴(yán)重肝功能損害中凋亡過(guò)程 從而加快細(xì)胞凋亡國(guó)心臟起搏與心電生理雜志， （）：。，血脂康膠囊主要成分為紅曲，化學(xué)成分較為復(fù)雜，主要，為他汀類、不飽和脂肪酸、多種人體氨基酸、類固醇及小量黃， （）：酮類物質(zhì)，有調(diào)血脂、保護(hù)血管內(nèi)皮、抑制過(guò)氧化損傷、阻滯馬慧珍，郭美姿，張 偉 血脂康引起藥物性肝炎藥物不良血管平滑肌細(xì)胞增殖和遷移等作用。紅曲含有洛伐他汀，在反應(yīng)雜志，（）：血脂康膠囊臨床應(yīng)用中國(guó)專家共識(shí)中標(biāo)示出每粒膠囊，）中洛伐他汀含量為，后者經(jīng)代謝，