0

我正在使用 R 中的 bife 包估计一个固定效应 logit 模型。我需要估计估计的边际效应并将它们导出到乳胶。但是，支持 bife 估计的包 texreg 不允许导出“bifeAPEs”类的对象。有什么解决办法吗？

这是一个可重现的示例：

#Require packages
packages<-c("bife", "texreg")
lapply(packages, require, character.only = TRUE)

#Dataset
data("iris")
iris$big <- ifelse(iris$Sepal.Length > median(iris$Sepal.Length),1,0)

#Model
output <- bife(big ~ Sepal.Width + Petal.Length | Species, data=iris, "logit")

#Output
apes_stat <- get_APEs(output)
class(apes_stat)
extract(output)
extract(apes_stat)
4

1 回答 1

0

也许您可以做的是创建一个data.frame平均部分效果，然后data.frame使用 将其转换为乳胶表xtable，如下所示：-

ape_data <- data.frame(APES = attr(apes_stat$delta, "names"),
                       delta = unname(apes_stat$delta))
xtab_ape <- xtable::xtable(ape_data)
xtab_ape

如果您不想rownames在 Latex 中打印，则可以添加以下命令：-

print(xtab_ape, include.rownames = FALSE)
于 2021-08-06T08:53:13.533 回答