在R中获得多个模型系数的函数的置信区间

Question

我已经在R（coxme程序包）中运行了混合效果Cox模型。在我的模型中，我具有性别（男性/女性）和父母身份（有孩子/没有）的变量，两者之间存在相互作用。

model <- coxme(Surv(spell.length, event) ~ sex + parent + sex:parent + (1 | id),...)

我想要四组具有置信区间的样板模型系数：男性/没有孩子；女性/没有孩子；男性/儿童女性/儿童。

我想知道如何做，因为最后一组是两个模型变量的函数。有人知道您可以在R中执行此操作吗？

谢谢。

更新：

对不起，应该给出一个数据示例。

library(coxme)
set.seed(1000)
df <- data.frame(
    id = 1:1000,
    spell.length = runif(1000, min = 0, max = 50),
    event = ifelse(runif(1000)<=0.8,0,1),
    sex = factor(x = rep(c("male","female"),500)),
    parent = factor(x=ifelse(runif(1000)<=0.5,"child","no child"))
    ) 
coxme(Surv(spell.length,event) ~ sex + parent + sex:parent + (1|id), df)

Answer 1

在线阅读后，我想答案是要记住...

并且标准误差只是估计量方差的平方根。

在R中，您可以使用vcov函数提取方差和协方差。例如：

data(mtcars)
model <- lm(mpg ~ disp + wt, data=mtcars)
var <- vcov(model)[2,2]+vcov(model)[3,3]+2*vcov[2,3]
se <- sqrt(var)