作者marki (marki)
看板Statistics
标题[问题] 想请问交互作用
时间Thu Sep 3 15:15:05 2020
抱歉最近问题有点多，想请问交互作用的问题
我有想要探讨在两个族群:不同的GA_G (0 or 1)其随着Year_G增加
AKI的机会是不是都会下降
我的logistic regression模型是
AKI ~ Year_G + factor(GA_G) + 一堆调整变数 + Year_G*factor(GA_G)
用R跑出来的结果是
Estimate Std. Error z value Pr(>|z|)
(Intercept) -2.454804 0.793533 -3.094 0.00198 **
Year_G -0.163630 0.073466 -2.227 0.02593 *
factor(GA_G)1 -1.551497 0.566584 -2.738 0.00618 **
Year_G:factor(GA_G)1 0.173219 0.132713 1.305 0.19182
....(调整变数的结果我就不列了)
想请问的是
1. 虽然交互作用项不显着，可是当GA_G从0变成1的时候，是不是AKI本来随着
Year_G会逐渐下降会变成会逐渐上升呢?
2. 如果真的是这样，那我还可以说AKI的机会在两组中都会显着地随着Year_G的增加
而下降吗?
3. 如果不行的话，要怎麽在paper上描述这样的结果呢?我有办法计算出在GA_G=1的这
个族群的p value或是判断他有没有显着吗?
4. 如果我想要画出一个Y轴为Probabilitiy，X轴为Year_G的GA_G=1的图，请问下述的程式码有
没有问题?
curve((exp(intercept + Year_G_coef*x + GA_G_coe *1 + 调整变项 + inter_coe2
*x*1)/(1+exp(intercept + Year_G_coe*x + GA_G_coe *1 + 调整变项 + inter_coe *x*1)))
以上
非常感谢
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.116.253.195 (台湾)
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Statistics/M.1599117307.A.9EB.html
※ 编辑: marki (140.116.253.195 台湾), 09/03/2020 15:51:58
1F：推 locka: 你的程式码看起来只是linear regression 不是 logistic reg 09/03 16:24
2F：→ locka: ression 耶? 09/03 16:24
3F：→ marki: 我算logistic是glm(AKI~..., family="binominal")所以应该 09/03 16:42
4F：→ marki: 是logistic拉 只是我不是用predict函数 09/03 16:43
5F：推 raiderho: 能否请你简略说明 AKI, GA_G, Year_G 的意思. 09/05 18:19