Re: [問題] 有關GEE結果的判讀
※ 引述《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
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
09/07 00:52, 4F
→
09/07 00:53, , 5F
09/07 00:53, 5F
推
09/07 02:14, , 6F
09/07 02:14, 6F
→
09/07 02:15, , 7F
09/07 02:15, 7F
→
09/07 02:15, , 8F
09/07 02:15, 8F
推
09/07 22:30, , 9F
09/07 22:30, 9F
→
09/07 22:31, , 10F
09/07 22:31, 10F
→
09/08 00:22, , 11F
09/08 00:22, 11F
→
09/08 00:22, , 12F
09/08 00:22, 12F
→
09/08 00:23, , 13F
09/08 00:23, 13F
討論串 (同標題文章)