【问题标题】:effect modification from coxph in R从R中的coxph进行效果修改
【发布时间】:2018-04-18 20:37:48
【问题描述】:

我正在研究 R 中时变 cox 模型上下文中的效果修改。我的模型如下:

fit1.1<- with (data, coxph(Surv(startdays, enddays, event) ~ age*strata(Region) + BMI.time.var + cluster(id), ties="breslow"))

我得到以下输出:

我真正感兴趣的是从这些输出中计算每个区域中年龄的效果估计值和 95% CI,以便在效果修改表中呈现。为此,我希望能够估计:

对年龄的影响估计:东北部、南部、西部、中西部。

我无法识别允许这样做的代码。

【问题讨论】:

    标签: r interaction cox-regression


    【解决方案1】:

    你可以试试 confint(fit1.1)。

    另外,如果您需要从 fit1.1 中提取任何其他信息,请尝试使用 str(fit1.1) 或 str(summary(fit1.1)) 查看拟合对象的结构。

    【讨论】:

    • 嗨 Len,我没有完全描述我在寻找什么。我很抱歉。我真正感兴趣的是从这些输出中计算每个区域的效果估计值和 95% 的年龄 CI,以便在效果修改表中呈现。为此,我希望能够估计:对以下地区年龄的影响估计:东北部、南部、西部、中西部。
    • 没问题。如果没有数据样本,这有点难以解决,但我认为“age:strate(Region)Northeast”(以及“age:strate(Region)South”和“age:strate(Region)West”)给出 age*region 相对于 "age:strate(Region)Midwest" 的交互项.如果您只需要每个区域的年龄估计值和 95% CI 置信区间,为什么不将数据框拆分为每个区域的数据框,然后在该步骤之后进行 cox 分析?
    猜你喜欢
    • 1970-01-01
    • 2015-09-11
    • 1970-01-01
    • 2022-10-25
    • 1970-01-01
    • 2019-07-06
    • 1970-01-01
    • 1970-01-01
    • 2023-04-10
    相关资源
    最近更新 更多