r - 从R中的clogit获得个体系数值

Question

让我们以生存包中的条件逻辑回归示例为例

并使用以下命令

library(survival)

data(logan)

resp <- levels(logan$occupation)
n <- nrow(logan)
indx <- rep(1:n, length(resp))
logan2 <- data.frame(logan[indx,],
                     id = indx,
                     tocc = factor(rep(resp, each=n)))
logan2$case <- (logan2$occupation == logan2$tocc)
B <- clogit(case ~ tocc + tocc:education + strata(id), logan2)

现在我们可以生成回归参数，但假设我们特别想要toccfarm-1.896 的值。

我们如何仅输出此内容或将其另存为x另存为。

当我们使用

B$coefficients

我们得到所有的回归系数。

我尝试过类似的东西

B$coefficients[1,]
B$coefficients(term=1)
B$coefficients("toccfarm")

但没有一个有效

Answer 1

另一种方法可以是使用该summary函数。您可以看到，通过汇总，模型的系数被视为矩阵。

> is.matrix(summary(B)$coefficients)
[1] TRUE

此时，您可以存储summary(B)$coefficients在一个对象中，然后根据需要对其进行子集化。

summary(B)$coefficients[1,1]