2016-09-28 189 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