r - 估计零通胀模型时如何用 R 获取系数

Question

可能很简单，但我想知道，如何在使用 Zeroinfl 命令时获取系数？

treatment <- factor(rep(c(1, 2), c(43, 41)), 
    levels = c(1, 2),labels = c("placebo", "treated"))
improved <- factor(rep(c(1, 2, 3, 1, 2, 3), c(29, 7, 7, 13, 7, 21)),
    levels = c(1, 2, 3),labels = c("none", "some", "marked"))

numberofdrugs <- rpois(84, 2)
healthvalue <- rpois(84,0.5)
y <- data.frame(healthvalue,numberofdrugs, treatment, improved)

require(pscl)
ZIP<-zeroinfl(healthvalue~numberofdrugs+treatment+improved, y)

summary(ZIP)

我通常使用 ZIP$coef[1] 来获取系数，但不幸的是在这里你获取了一大堆。那么如何从 ZIP 模型中获取单个系数呢？

Answer 1

使用coef提取函数列出一个长向量中的所有系数，然后您可以使用单索引表示法来选择它们:

coef(ZIP)[1]

count_(Intercept) 
        0.1128742 

或者，您需要首先选择要从中获取系数的模型:

ZIP$coef$count[1]

(Intercept) 
  0.1128742

ZIP$coef[[1]][1]

(Intercept) 
  0.1128742