[程式] SAS卡方問題

看板Statistics作者 (希)時間9年前 (2016/06/23 21:04), 編輯推噓1(108)
留言9則, 3人參與, 最新討論串1/1
大家晚安~~有個問提想請問大家~~ 我用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://www.ptt.cc/bbs/Statistics/M.1466687096.A.010.html

06/23 21:41, , 1F
你用proc logistic去跑
06/23 21:41, 1F

06/23 22:20, , 2F
有~跑出來了~想再請教logistic跑出來的or & 95CI為
06/23 22:20, 2F

06/23 22:23, , 3F
OR:1.416 95CI:1.383-1.45 但這樣的解釋是2組比1組的得
06/23 22:23, 3F

06/23 22:24, , 4F
病為1.416倍,3組比2組的得病也是1.416倍嗎?
06/23 22:24, 4F

06/23 22:25, , 5F
真的很感謝你的回答!!!
06/23 22:25, 5F

06/23 23:58, , 6F
用dummy variable可以檢視非線性的可能性
06/23 23:58, 6F

06/24 22:37, , 7F
M大謝謝您的回答,雖然不太懂您的意思(不好意思統計太弱)
06/24 22:37, 7F

06/24 22:38, , 8F
但已經私下請教1樓I大,問題已經解決了,但真的都很感謝
06/24 22:38, 8F

06/24 22:38, , 9F
你們的回覆!!!謝謝你們!!!
06/24 22:38, 9F
文章代碼(AID): #1NQzvu0G (Statistics)