从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")

但 none 有效

Answer 1

另一种方法是使用 summary 函数。可以看到summary把模型的系数当成了矩阵

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

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

summary(B)$coefficients[1,1]

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

obtain individual coefficent value from clogit in R

r

logistic-regression

survival-analysis

coefficients