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

Question

让我们从Survival Package

中获取条件logisitc回归示例

使用以下命令

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]