2015 年 7 月 20 日更新:
extract.zelig
现在有一个tobit
方法 ( Zelig_4.2-1
)
所以texreg(reg)
现在按预期工作。无论如何,我将把下面的内容留给后代。
在确定了问题的根源后,我更新了extract.zelig
方法并将其传递给了包创建者/维护者 Philip Leifield,他将其合并到了最新的 R-Forge 版本中(可通过 获得install.packages("texreg", repos="http://R-Forge.R-project.org")
)。我不确定它是否在当前的 CRAN 版本中(2015-04-07)...
这是我们需要添加的内容:
else if ("tobit" %in% class(model)) {
coefficient.names <- rownames(s$table)
coefficients <- s$table[, 1]
standard.errors <- s$table[, 2]
significance <- s$table[, 5]
gof <- numeric()
gof.names <- character()
gof.decimal <- logical()
if (include.aic == TRUE) {
aic <- AIC(model)
gof <- c(gof, aic)
gof.names <- c(gof.names, "AIC")
gof.decimal <- c(gof.decimal, TRUE)
}
if (include.bic == TRUE) {
bic <- BIC(model)
gof <- c(gof, bic)
gof.names <- c(gof.names, "BIC")
gof.decimal <- c(gof.decimal, TRUE)
}
if (include.loglik == TRUE) {
lik <- logLik(model)[1]
gof <- c(gof, lik)
gof.names <- c(gof.names, "Log Likelihood")
gof.decimal <- c(gof.decimal, TRUE)
}
if (include.nobs == TRUE) {
n <- nrow(model$data)
gof <- c(gof, n)
gof.names <- c(gof.names, "Num. obs.")
gof.decimal <- c(gof.decimal, FALSE)
}
tr <- createTexreg(coef.names = coefficient.names, coef = coefficients,
se = standard.errors, pvalues = significance, gof.names = gof.names,
gof = gof, gof.decimal = gof.decimal)
return(tr)
}