作者vm212 (啦啦啦)
看板Statistics
标题[问题] SAS code请教
时间Mon Jun 16 09:44:18 2014
各位版上的学长学姐们:
小妹，目前念护理所一年级，在使用SAS软体跑统计遇到了一些问题，希望能请学长学姐们，好心的协助我解决。
我想要利用SAS logistic analysis 去看X1, X2,X3因子与Y疾病的关系。所有的资料都是biary资料，我也作了dummy variable的处理了。
由於我作的是回朔性研究，是翻病历建立病患资料，有许多以前的病历资料都已经不存在，所以有许多遗漏值必须处理，而这些遗漏值我查了一下遗漏值的定义，是属於MNAR (missing not at random)，可以用MCMC method进行插补。
我现在的想法是先处理遗漏值，在进行logistic analysis。
我所写出的sas code 如下: (SAS 9.0)
DATA Logistic;
INPUT Y X1 X2 X3;
CARDS;
0 . 0 1
0 1 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 1 . 0
0 0 0 0
0 0 0 0
0 1 0 0
0 1 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 1 0 0
0 0 0 0
0 1 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 1 0 0
0 0 0 0
0 0 0 0
0 . 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 .
0 1 0 0
0 . 0 0
0 0 0 0
0 1 0 0
0 0 0 0
0 0 0 0
0 1 0 0
0 0 0 0
0 . 0 0
0 1 0 0
0 1 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 1 0
0 0 0 0
1 1 1 0
1 . 0 0
1 . 0 .
1 . . 0
1 1 1 0
1 1 0 1
1 1 0 0
1 . 0 1
1 . 0 1
1 1 0 1
1 0 1 0
1 1 0 1
1 . 0 1
1 1 0 1
1 . 0 1
1 . . 1
1 1 1 0
1 1 0 0
1 . 1 .
1 1 1 0
1 1 0 1
1 . 0 0
1 1 . 1
1 1 0 1
1 . 0 0
1 1 0 0
1 1 0 0
1 . 0 0
1 0 0 0
;
proc mi data=Logistic nimpute=2;
var X1 X2 X3;
run;
proc logistic data=logistic;
class Y(ref='0')/param=ref;
model X1(event='1')=Y;
run;
问题 1:
X1、X2、X3 missing data，总共有24个，请问为什麽只有进行2个地方的插补呢?
http://ppt.cc/xy0-
http://ppt.cc/baD7
问题 2:
进行logistic分析时，为什麽总data只剩29个?似乎没有计算到插补的地方?
http://ppt.cc/QilC
不知道是不是code设计哪边出了问题? 不知道能不能请学长姐提供一些改进的方法?谢谢~
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 42.71.234.6
※ 文章网址: http://webptt.com/cn.aspx?n=bbs/Statistics/M.1402883061.A.E5C.html
1F：→ andrew43:插补的解读错了。结果是「仅X2缺失共有2比」…… 06/18 11:19
2F：→ andrew43:至於第二个问题，我猜是你没有引入插补後的资料。 06/18 11:21
3F：→ andrew43:大概和 out=xxx 或 data=xxx 这种引数有关。 06/18 11:22
4F：→ andrew43:所以你的回归只采用了插补前的data。 06/18 11:22