我正在使用 R 中的 bife 包估计一个固定效应概率模型。我正在尝试将输出提取到可以与观星者或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")
我想
我找到了这个解决方案的替代解决方案,即使可能为时已晚
基本上,首先,我进入了包"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)
我希望它会有所帮助!最好