在Rmarkdown中将bife对象的结果输出到Latex？答案

【问题标题】：Outputting the results from bife object to Latex in Rmarkdown?在Rmarkdown中将bife对象的结果输出到Latex？
【发布时间】：2019-06-20 03:12:21
【问题描述】：

我正在使用 R 中的 bife 包估计一个固定效应概率模型。我正在尝试将输出提取到我可以与 stargazer 或 texreg 一起使用的东西中，这样我就可以使用 Rmarkdown 将它们输出到论文中来创建一个 LaTeX 对象。我知道我可以手动提取系数和标准误差等，但我想知道是否没有更有效的方法可以将此对象强制转换为适用于任一包的对象。

这是一个可重现的例子：

install.packages("bife")
library(bife)
data("iris")
iris$big <- ifelse(iris$Sepal.Length > median(iris$Sepal.Length),1,0)
output <- bife(big ~ Sepal.Width + Petal.Length | Species, data=iris, "logit")

【问题讨论】：

你有没有找到一个好的解决方案？
很遗憾没有。
我知道这不是一个直接的答案，因为您需要使用bife:: 和一个概率，但我最终使用survival::clogit() 作为我的固定效应 logit，并且在 Stargazer 中有效。

标签： r-markdown stargazer texreg

【解决方案1】：

我想我找到了一个替代解决方案，即使可能为时已晚

基本上，首先，我去了包“texreg”的存储库，发现了这个函数：

extract.bife <- function(model,
                         include.loglik = TRUE,
                         include.deviance = TRUE,
                         include.nobs = TRUE,
                         ...) {
  s <- summary(model)
  coefficient.names <- rownames(s$cm)
  co <- s$cm[, 1]
  se <- s$cm[, 2]
  pval <- s$cm[, 4]

  gof <- numeric()
  gof.names <- character()
  gof.decimal <- logical()
  if (include.loglik == TRUE) {
    lik <- logLik(model)
    gof <- c(gof, lik)
    gof.names <- c(gof.names, "Log Likelihood")
    gof.decimal <- c(gof.decimal, TRUE)
  }
  if (include.deviance == TRUE) {
    gof <- c(gof, deviance(model))
    gof.names <- c(gof.names, "Deviance")
    gof.decimal <- c(gof.decimal, TRUE)
  }
  if (include.nobs == TRUE) {
    n <- s$nobs["nobs"]
    gof <- c(gof, n)
    gof.names <- c(gof.names, "Num. obs.")
    gof.decimal <- c(gof.decimal, FALSE)
  }

  tr <- createTexreg(
    coef.names = coefficient.names,
    coef = co,
    se = se,
    pvalues = pval,
    gof.names = gof.names,
    gof = gof,
    gof.decimal = gof.decimal
  )
  return(tr)
}

因此，对于您的示例，只需将其应用于您的模型并使用函数 texreg，您可能会得到一个 Latex-"like" 输出

tr <- extract.bife(output)
texreg(tr)

我希望它会有所帮助！最好的

【讨论】：