[程式] SPSS GLM算incidence rate ratio

看板Statistics作者 (Mutex)時間9年前 (2015/03/19 21:23), 編輯推噓0(007)
留言7則, 2人參與, 最新討論串1/1
目的: 希望使用SPSS 分別算出組別的incidence rate 然後再用Incidence Rate Ratio(IRR)比較兩者的差別 參考如右圖 http://imgur.com/6jFBVkj
這是我測試的data http://ppt.cc/jZaJ 目前組別的Event, PY, rate已使用aggregate分別求出 group event_sum follow_sum rate=[(event_sum/follow_sum)*1000] 1 4 279.8 14.3 2 5 805.6 6.2 照說IIR應該會約落在14.3/6.2=2.3左右(reference:2) 或 6.2/14.3=0.4(reference:1) 但在參考下面這網址使用GLM的poisson算IRR時,算出來的IRR卻只有1.708 http://www.ats.ucla.edu/stat/spss/dae/poissonreg.htm 我使用的語法和結果 http://imgur.com/rP6sSS3
研究了好一陣子都弄不出來 不知道是哪裡出錯了 懇請大師達人們賜教 感激不盡 -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 36.226.185.112 ※ 文章網址: https://www.ptt.cc/bbs/Statistics/M.1426771394.A.37F.html

03/19 22:55, , 1F
剛嘗試用stata算的 http://imgur.com/0AYtq5N
03/19 22:55, 1F

03/19 23:23, , 2F
你所謂 rate=[(event_sum/follow_sum)*1000] 和你的
03/19 23:23, 2F

03/19 23:23, , 3F
迴歸模型,我看不出二者的關係。
03/19 23:23, 3F

03/19 23:24, , 4F
1.708 是從 exp(beta_group) = exp(0.535) 來的。
03/19 23:24, 4F

03/19 23:25, , 5F
我在猜,你是不是想讓follow設定成 offset 而不是 IV
03/19 23:25, 5F

03/19 23:27, , 6F
是的話,IRR就會是14.3/6.2=2.3 了
03/19 23:27, 6F

03/19 23:29, , 7F
所以我猜,follow按你的本意應該是暴露量而不是共變數。
03/19 23:29, 7F
文章代碼(AID): #1L2it2D_ (Statistics)