2016-09-28 191 views
1

我有一套數據,我想做邏輯迴歸建模二元結果變量(治療)的機率,階段作爲序數解釋變量(0,1,2,3 ,4)。 Hba1c是一個連續變量。邏輯迴歸與序數解釋變量

我的課堂陳述是否正確?

我該如何計算序數變量每個級別的比值比?

PROC LOGISTIC data=new; 
class EyeID Therapy (ref ="0") Stage (param = ordinal) Gender (ref="M") Ethnicity (ref="C")/ param = ref; 
model Therapy = Stage Gender age A1c Ethnicity; 
oddsratio Stage; 
run; 

這是輸出:

Odds Ratio Estimates and Wald Confidence Intervals 
Odds Ratio Estimate 95% Confidence Limits 
Stage 1 vs 0 0.873 0.547 1.394 
Stage 2 vs 0 2.434 0.895 6.620 
Stage 3 vs 0 0.915 0.431 1.941 
Stage 4 vs 0 0.356 0.132 0.961 
Stage 2 vs 1 2.788 0.980 7.935 
Stage 3 vs 1 1.048 0.465 2.360 
Stage 4 vs 1 0.408 0.144 1.156 
Stage 3 vs 2 0.376 0.113 1.249 
Stage 4 vs 2 0.146 0.038 0.567 
Stage 4 vs 3 0.389 0.117 1.288 

如果我在報告階段爲序變量,那麼它是正確的,我創建這樣一個表?

Stage 1 vs 0 0.873 0.547 1.394 
Stage 2 vs 1 2.788 0.98 7.935 
Stage 3 vs 2 0.376 0.113 1.249 
Stage 4 vs 3 0.389 0.117 1.288 

我不應該這樣報告,對嗎?這是如果階段是絕對的?

Stage 1 vs 0 0.873 0.547 1.394 
Stage 2 vs 0 2.434 0.895 6.62 
Stage 3 vs 0 0.915 0.431 1.941 
Stage 4 vs 0 0.356 0.132 0.961 
+0

序數模型的累積概率,與參考水平相比並非如此。那是你在找什麼? – Reeza

+0

@Reeza我已更新該帖子。我是否正確解釋輸出? – ybao

回答

1

我不認爲你需要Therapy在類聲明。

沒有樣本數據,我無法測試這個,但我的第一次通過本來就是這樣寫的。

proc logistic data=test; 
class PVDStage (param = ordinal); 
model Therapy(ref = '0') = PVDStage hba1c; 
ODDSRATIO PVDStage; 
run; 

如果您可以提供一些示例數據,我會修改我的答案以確保其正常工作。

+0

謝謝。這很好用! – ybao

+0

ref ='0'應該是event ='0',實際上會導致意想不到的結果。默認值將被使用,SAS不會通過一個錯誤,這看起來很奇怪。 – Reeza

+0

我從來沒有遇到過問題的語法問題。您是否有SAS Note或Doc鏈接來討論它? – DomPazz