来自zelig package(R)的turit-ing Tobit输出

时间:2015-04-02 00:16:25

标签: r texreg r-zelig

这是一个奇怪的问题,但这里有:

我正在尝试将我的模型结果输出到带有texreg的TeX表中。

reg <- zelig(Y ~ X, model = "tobit", below = 0, above = Inf)

但是,我从texreg收到错误:

texreg(reg)
  

.local(model, ...)中的错误:         目前仅支持以下Zelig模型:              logit,ls,mlogit,ologit,probit,relogit。

我的问题基本上是:这是来自Zelig还是来自texreg的错误?

1 个答案:

答案 0 :(得分:3)

更新2015-07-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)
    }