Re: [問題] 有關GEE結果的判讀

看板Statistics作者 (Syameroke)時間12年前 (2013/09/07 00:34), 編輯推噓2(2011)
留言13則, 3人參與, 最新討論串2/4 (看更多)
※ 引述《marki (marki)》之銘言: : 各位大大好 : 小弟有一個關於GEE的問題想請教各位大大 : 我有一個data base 是一群慢性腎疾病病人追蹤16年的數據 : 每個病人都是一年追蹤一次 : 記錄的數據是每一年每個病人的年紀、收縮壓、舒張壓、血色素、尿酸、慢性腎疾病分期 : 等等的實驗數據 : 所以算是重複測量 : 小弟想要探討的問題是 : 不同的併發症(例如貧血)在不同的慢性腎疾病分期 (CKD stage)比例是否有不同 : ps. 所以有沒有貧血 應該是類別變項 : 跑出來的結果是 : 變項 B SE p OR OR 95%CI : Lower Upper : Year -0.146 0.038 0.000 0.864 0.802 0.931 : CKDstage 5 2.186 0.646 0.001 8.900 2.507 31.598 : CKDstage 4 2.592 0.624 0.000 13.363 3.930 45.435 : CKDstage 3 0.527 0.299 0.078 1.694 0.942 3.046 : CKDstage 2 -0.248 0.236 0.293 0.781 0.492 1.239 : CKDstage 1 ─ ─ ─ ─ ─ ─ : Year*CKDstage 5 0.119 0.091 0.194 1.126 0.941 1.347 : Year*CKDstage 4 0.009 0.089 0.920 1.009 0.848 1.201 : Year*CKDstage 3 0.163 0.058 0.005 1.177 1.051 1.319 : Year*CKDstage 2 0.111 0.048 0.021 1.118 1.017 1.229 : Year*CKDstage 1 ─ ─ ─ ─ ─ ─ : 我想問的事情是 : 1. 所以根據這樣的分析,我得到的結論是 : 「貧血在CKD stage stage4比CKD stage1容易發生」嗎? : 2. 如果上面的答案是肯定的話,那請問用GEE跑出來的結果和我用每個病人第一次的 : CKD stage與「有/無貧血」跑qui-square跑出來意義有什麼不一樣? 就回答這兩個問題: 1. 不算對也不完全錯,因為我猜你只有看下列這列: CKDstage 4 2.592 0.624 0.000 13.363 3.930 45.435 他告訴你的意思是,當year=0時(也就是baseline),第4期的比第1期的貧血比率高 而隨著時間推移,這個OR值將會越來越小,因為0.009+(-0.146)<0 Year -0.146 0.038 0.000 0.864 0.802 0.931 Year*CKDstage 4 0.009 0.089 0.920 1.009 0.848 1.201 按照這個模型,只要追蹤時間達到19年,那第一期的貧血比率就會高於第四期,不過你的推論 範圍可能只有到16年,所以應該還好,但最後幾年檢定絕對是不會顯著的 2. 這個答案要看你的資料來決定如何回答 我先假設你在model中用的CKD Stage是指"baseline"的Stage,那這樣的話最有可能是你的 "chi-square"檢定整體而言無差異,但其實內部有幾格是有差異的 如果你的model中使用的CKD Stage會隨著時間變換,那就有點麻煩了,因為這個model也許做 了許多過多不太合理的假設(就我所知CKD患者會不可逆的往ESRD走,若貧血也是不可回覆的 ,那這個model提供的數字其實合理性不夠),這可能要拿到原始資料才能回答你 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 61.64.174.131 ※ 編輯: ching0629 來自: 61.64.174.131 (09/07 00:47)

09/07 00:48, , 1F
另外,這個模型非常奇怪,貧血可逆? Year的係數告訴我他可
09/07 00:48, 1F

09/07 00:49, , 2F
逆,我不確定這件事情的真實性,你要多加判斷一下,這個模
09/07 00:49, 2F

09/07 00:49, , 3F
型未必適合你
09/07 00:49, 3F

09/07 00:52, , 4F
當然,搞不好是model沒有錯而是資料錯了,也許統計公司對
09/07 00:52, 4F

09/07 00:53, , 5F
不可逆並不清楚
09/07 00:53, 5F

09/07 02:14, , 6F
year和EPO治療不曉得有沒有相關性...
09/07 02:14, 6F

09/07 02:15, , 7F
如果罹病後期嚴重到一定程度才能健保給付EPO的話...
09/07 02:15, 7F

09/07 02:15, , 8F
那會有干擾作用...
09/07 02:15, 8F

09/07 22:30, , 9F
可以請問一下大大那個19年是怎麼得出來的啊
09/07 22:30, 9F

09/07 22:31, , 10F
另外,大大是怎麼看出來anemia是可逆的啊
09/07 22:31, 10F

09/08 00:22, , 11F
Year的係數為負就代表第一期可逆,若StageN*year的係數比
09/08 00:22, 11F

09/08 00:22, , 12F
year係數的絕對值小就代表第n期是可逆的
09/08 00:22, 12F

09/08 00:23, , 13F
19是這樣來的:(0-2.592)/(-0.146+0.009)
09/08 00:23, 13F
文章代碼(AID): #1IAWG8is (Statistics)
討論串 (同標題文章)
文章代碼(AID): #1IAWG8is (Statistics)