作者ching0629 (Syameroke)
看板Statistics
標題Re: [問題] 有關GEE結果的判讀
時間Sat Sep 7 00:34:13 2013
※ 引述《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)
1F：→ ching0629:另外,這個模型非常奇怪,貧血可逆? Year的係數告訴我他可 09/07 00:48
2F：→ ching0629:逆,我不確定這件事情的真實性,你要多加判斷一下,這個模 09/07 00:49
3F：→ ching0629:型未必適合你 09/07 00:49
4F：→ ching0629:當然,搞不好是model沒有錯而是資料錯了,也許統計公司對 09/07 00:52
5F：→ ching0629:不可逆並不清楚 09/07 00:53
6F：推 anovachen:year和EPO治療不曉得有沒有相關性... 09/07 02:14
7F：→ anovachen:如果罹病後期嚴重到一定程度才能健保給付EPO的話... 09/07 02:15
8F：→ anovachen:那會有干擾作用... 09/07 02:15
9F：推 marki:可以請問一下大大那個19年是怎麼得出來的啊 09/07 22:30
10F：→ marki:另外，大大是怎麼看出來anemia是可逆的啊 09/07 22:31
11F：→ ching0629:Year的係數為負就代表第一期可逆,若StageN*year的係數比 09/08 00:22
12F：→ ching0629:year係數的絕對值小就代表第n期是可逆的 09/08 00:22
13F：→ ching0629:19是這樣來的:(0-2.592)/(-0.146+0.009) 09/08 00:23