[問題]Hazard ratio與發生密度

看板Statistics作者 (努力賺錢)時間13年前 (2012/10/29 20:54), 編輯推噓0(0推 0噓 12→)
留言12則, 3人參與, 最新討論串1/2 (看更多)
我想要請問，當觀察組的發生密度(新個案/總追蹤人年*1000)數值小於參考組時， 為何觀察組的Hazard ratio會大於1,雖然不顯著。相反的，我的分析結果，也有 觀察組發生密度較參考組高，但Hazard ratio的數值卻小於1的情況？ 我所使用的軟體為SAS9.2版，但因為X變相資料是利用多段不等機率抽樣的營養調查， 故使用Sudaan 9 來校正抽樣誤差。 但是即便我不使用Sudaan語法，這個結果仍然存在。我看過歷屆學長姊的論文，發現 也有想同的問題存在，但在這種情況下大多是不顯著的。 還有這只是curde model，並未校正其他變項。 想請問這樣的情況要怎麼解釋，還是我可能在哪個步驟出錯?? 這個問題不處理好，我一整個沒心思跑其他資料了!!!T^T 請高手相救，感恩!! -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 111.250.63.230
→
10/29 21:03, , 1F
hazard ratio指的是什麼?發生密度又是怎麼算的?
10/29 21:03, 1F
→
10/29 21:08, , 2F
風險比(HR);發生密度(新發個案數/總人年數*1000)!!
10/29 21:08, 2F
→
10/29 21:49, , 3F
文中的 "發生密度" 我猜是類似 p.d.f., 以 f(x) 表示之, 則
10/29 21:49, 3F
→
10/29 21:50, , 4F
對應的 hazard rate 是 f(x)/(1-F(x)). 並不一定 f(x)<g(x)
10/29 21:50, 4F
→
10/29 21:50, , 5F
就保證 f(x)/(1-F(x))<g(x)/(1-G(x)).
10/29 21:50, 5F
→
10/29 22:17, , 6F
如果是y大說的ratio of hazard functions,那麼不太會是個常數
10/29 22:17, 6F
→
10/29 22:19, , 7F
所以想先搞清楚原po在問什麼
10/29 22:19, 7F
→
10/29 23:06, , 8F
我是想要比較吃素的人是否有較低的慢性病罹患風險!因此我先
10/29 23:06, 8F
→
10/29 23:07, , 9F
出有吃素及沒吃素其罹病的個數，及兩組的總追蹤人年數。
10/29 23:07, 9F
→
10/29 23:12, , 10F
再以SAS proc survive的語法分析HR
10/29 23:12, 10F
→
10/30 10:13, , 11F
PROC Lifetest吧? 有點懷疑你知不知道自己在做什麼分析
10/30 10:13, 11F
→
10/30 10:14, , 12F
你的data長什麼樣子, SAS code怎麼寫的?
10/30 10:14, 12F
‣ 更多 ofine 的文章
文章代碼(AID): #1GZdoAKo (Statistics)
‣ 更多 ofine 的文章
文章代碼(AID): #1GZdoAKo (Statistics)