以《多分类有序反应变量Logistic 回归及其应用》(作者高歌,同济大学学报,第31 卷第10 期 2003 年10 月)的例子为例:
作者随机抽取 2002 年国家执业医师临床实践技能考试 30 岁以下本科学历考生 5 000 人.确定 9 个自变量如下: ①性别 x 1 (男, x 1 = 0 ;女, x1 = 1) ; ②工作单位级别 x 2 (二级及以下医疗机构, x2 = 0 ;三级医院, x 2 = 1) ; ③考试评分方法 x3 (考官独立评分, x3 = 0 ;考官集体评分, x 3 = 1) ; ④考官回避制度 x4 (考官未回避本单位考生, x4 = 0 ;考官已回避本单位考生, x 4 = 1) ; ⑤考题重复次数 x5 (重复 3 次及以上, x5 =0 ;重复3 次以下, x 5 = 1) ; ⑥合格率控制措施 x 6 (有控制, x6 = 0 ;无控制, x6 = 1) ; ⑦考试基地级别 x 7 (二级甲等医院, x7 = 0 ;三级医院, x 7 = 1) ; ⑧模拟人质量 x8 (质量差, x 8 = 0 ;质量好, x 8 = 1) ; ⑨ “病人” 体征x 9 (无阳性体征, x 9 = 0 ;有阳性体征, x9 = 1) .以考生的考试成绩 y (优, y = 1 ;中, y = 2 ;不合格, y = 3)为三分类有序反应变量.
根据文章的信息,我试图写出作者的SAS程序,也不知道对不对。
Proc
logistic
;
Class
x1 x2 x3 x4 x5 x6 x7 x8 x9/param=ref des;
model y= x1 x2 x3 x4 x5 x6 x7 x8 x9/selection=stepwise
sls=0.1
sle=0.1;
Run;
课本上说:因变量为多值变量时,logistic回归分析要求多个回归方程中自变量的系数是相等的。因此对于多值变量的logistic回归模型,要求做平行性检验,test for the proportional odds assumption。
而且用逐步回归删选自变量结果并不会出现score test for the proportional odds assumption。我想问的是,在这样的情况下还需要做平行性检验吗?该如何做?程序怎么写?
是不是把筛选后的自变量模型重新运行,验证score test for the proportional odds ssumption。