19醫(yī)學(xué)應(yīng)用統(tǒng)計(jì)分析-第十九章-多項(xiàng)無(wú)序分類(lèi)變量與數(shù)值變量關(guān)系的分析

第九章

多項(xiàng)有序分類(lèi)變量與數(shù)值變量關(guān)系的分析

主講老師：陳青山數(shù)據(jù)數(shù)據(jù)特征：結(jié)果變量：多項(xiàng)有序分類(lèi)變量；影響變量：數(shù)值變量。統(tǒng)計(jì)方法：有序多分類(lèi)Logistic回歸（本章講解此部分）；Spearman等級(jí)相關(guān)。舉例：根據(jù)鼻咽癌患者的年齡預(yù)測(cè)其手術(shù)的預(yù)后（好轉(zhuǎn)、惡化、死亡）；年氣溫與中暑程度的關(guān)系。實(shí)例例19-1某醫(yī)院為了解就診患者的年齡對(duì)就診滿(mǎn)意度的影響，隨機(jī)抽取176名就診患者進(jìn)行問(wèn)卷調(diào)查，調(diào)查結(jié)果見(jiàn)表19-1。問(wèn)患者的年齡是否影響其就診滿(mǎn)意度？表19-1患者的年齡與就診滿(mǎn)意度情況

注：滿(mǎn)意度中1、2和3分別表示不滿(mǎn)意、態(tài)度中立和滿(mǎn)意。

編號(hào)年齡滿(mǎn)意度155325633241………176252實(shí)例分析數(shù)據(jù)特征：結(jié)果變量：多項(xiàng)有序分類(lèi)變量，病人滿(mǎn)意度（satisfy=1、2、3），有不滿(mǎn)意、態(tài)度中立和滿(mǎn)意三個(gè)等級(jí)；影響變量：數(shù)值變量，年齡（歲）。分析目的：分析就診患者的年齡是否影響就診滿(mǎn)意度。統(tǒng)計(jì)方法：采用有序多分類(lèi)Logistic回歸分析。SAS的計(jì)算程序[操作程序]

例19-1的SAS操作程序SASP19_1：libnamesas"F:\data\sas";/*新建永久邏輯庫(kù)，定義邏輯庫(kù)名（數(shù)據(jù)庫(kù)庫(kù)名）sas，指定保存路徑為F:\data\sas*/

proclogisticdata=sas.d19_1descending;/*調(diào)用logistic過(guò)程,

選擇sas.d19_1數(shù)據(jù)集，設(shè)置反應(yīng)變量的水平排序?yàn)榻敌颉?/modelsatisfy(event='1')=age;/*satisfy作為結(jié)果變量，指定滿(mǎn)意程度以不滿(mǎn)run;意為參照，age作為影響變量*/SAS的計(jì)算結(jié)果反應(yīng)變量的水平排序?yàn)榻敌?，按照satisfy=3的概率擬合模型。本例：x2=0.2291，P=0.6322>0.05。不管模型中結(jié)果變量的分割點(diǎn)在什么位置，模型中的影響變量的系數(shù)

都保持不變，可以進(jìn)行有序多分類(lèi)Logistic回歸。本例：似然比檢驗(yàn)結(jié)果x2=11.0094，P=0.0009<0.05；Score檢驗(yàn)結(jié)果x2=10.8142，P=0.0010<0.05；Wald檢驗(yàn)結(jié)果x2=10.4877，P=0.0012<0.05；P值均小于0.05，故可以認(rèn)為模型成立。年齡的OR值為1.041，其95%可信區(qū)間為（1.016，1.066），有統(tǒng)計(jì)學(xué)意義，表明年齡越大的患者其就診滿(mǎn)意度越高。SPSS的操作步驟[操作步驟]例19-1的SPSS操作程序SPSSP19-1：

打開(kāi)d19-1.sav文件，在SPSS程序中按以下步驟操作：SPSS的計(jì)算結(jié)果應(yīng)用似然比檢驗(yàn)（-2LogLikelihood）對(duì)模型是否成立進(jìn)行檢驗(yàn)，x2=11.099，P=0.001<0.05，故可以認(rèn)為模型成立。通過(guò)公式OR=eβ可求得，年齡的OR值為1.04，其95%可信區(qū)間為（1.016，1.066），有統(tǒng)計(jì)學(xué)意義，表明年齡越大的患者其就診滿(mǎn)意度越高。-2LogLikelihood檢驗(yàn)結(jié)果x2=0.232，P=0.63