作者cubicin (希)
看板Statistics
標題[問題] SAS卡方問題
時間Thu Jun 23 21:04:54 2016
大家晚安~~有個問提想請問大家~~
我用SAS跑卡方的時候
2*2table 沒有問題,但大於三組*2 時,OR&95CI就跑不出來了
例如我想看
有沒有糖尿病的人得到高血壓的有無差別
所以我用卡方去檢定也想看Crude OR 95ci
程式如下:
proc freq data =m.analysis ORDER=data;
tables DM*HT/ CHISQ relrisk;
run;
可以跑出來我要的卡方值以及OR & 95CI
但換成4組時(例如年齡我分為4組 65以下 65-74 75-85 86以上)
proc freq data =m.analysis ORDER=data;
tables age_gp*HT/ CHISQ relrisk;
run;
就只跑出卡方值 OR和95CI就跑不出來了
請問是多組別不能這樣跑嗎?還是有別的撰寫方法?
例如我是想知道
65歲以下(第一組)與65-74(第二組)得到HT的勝算比
65歲以下(第一組)與75-85(第三組)得到HT的勝算比
==>有人這樣看的嗎?
還是我想錯了?
或是像這總情況是怎麼處理?就直接檢定沒有在看OR的嗎?
感謝大家的回答了!!!!謝謝!!
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 223.136.60.163
※ 文章網址: https://webptt.com/m.aspx?n=bbs/Statistics/M.1466687096.A.010.html
1F:推 imaltar: 你用proc logistic去跑 06/23 21:41
2F:→ cubicin: 有~跑出來了~想再請教logistic跑出來的or & 95CI為 06/23 22:20
3F:→ cubicin: OR:1.416 95CI:1.383-1.45 但這樣的解釋是2組比1組的得 06/23 22:23
4F:→ cubicin: 病為1.416倍,3組比2組的得病也是1.416倍嗎? 06/23 22:24
5F:→ cubicin: 真的很感謝你的回答!!! 06/23 22:25
6F:→ mosenay: 用dummy variable可以檢視非線性的可能性 06/23 23:58
7F:→ cubicin: M大謝謝您的回答,雖然不太懂您的意思(不好意思統計太弱) 06/24 22:37
8F:→ cubicin: 但已經私下請教1樓I大,問題已經解決了,但真的都很感謝 06/24 22:38
9F:→ cubicin: 你們的回覆!!!謝謝你們!!! 06/24 22:38