作者marki (marki)
看板Statistics
标题[问题] GEE的margins的问题
时间Sat May 28 14:37:13 2022
各位前辈好
小弟目前用手上的一个dataset 利用R 的glmtoolbox::glmgee跑了一个logistic模型
程式码大致如下
cfit <- glmtoolbox::glmgee(DV~ IV1 + ...,
id=id, corstr="exchangeable", family=binomial(),
data=dat)
然後用这个模型预测我新建的数据，想看在我建立的数据下，可能发生事件的Probability
newdata <- data.frame(IV1 = c(1,0,1), ....)
predict(cfit, newdata = newdata, se.fit = T, type ="response")
然後出现下面的表格
IV1 ......................... fit se.fit
IV1 1 ......................... 0.09 0.02
.
.
.
问题来了。
如果我想要用这个Table (也就是有SE的状况) 计算95% CI 要怎麽计算阿?
好像不是很单纯用0.09 +/- 1.96* 0.02
因为如果照下面这个网页
https://reurl.cc/k1o85K
最後一个表来说
https://imgur.com/a/AFSnOUI
他的Margin 也都不是Std. Err. * 1.96 (而且每一行相除也不是常数)
会是 Std. Err * 1.96/ sqrt (n)吗?
想请问在统计学上要怎麽处理?
我尝试找过文章或R package都找不太到
谢谢
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.116.253.195 (台湾)
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Statistics/M.1653719835.A.B4F.html
1F：→ andrew43: predict用type="link"然後从odd转换回probability 05/29 00:21
2F：→ andrew43: 在odd之下，算好+/-1.96se的区间，再转换回机率 05/29 00:23