[程式] SAS code請教
各位版上的學長學姐們:
小妹，目前念護理所一年級，在使用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://www.ptt.cc/bbs/Statistics/M.1402883061.A.E5C.html
→
06/18 11:19, , 1F
06/18 11:19, 1F
→
06/18 11:21, , 2F
06/18 11:21, 2F
→
06/18 11:22, , 3F
06/18 11:22, 3F
→
06/18 11:22, , 4F
06/18 11:22, 4F