2

我正在运行分位数回归（包quantreg）并texreg用于创建我的模型的乳胶输出。

我对引导 se 感兴趣并在摘要选项中设置 se="boot" 但是当我使用 texreg 时，我得到 "nid" se

如何更改该选项？

这是我正在做的事情：

tm3 <- rq(nback ~ cara + mat + dut + e_brown + e_green + e_blue + ins_no + ins_1
          + ins_5 + lit + sci + lan + gender_f + partner,  
          tau = 0.75, data=tru_all, model=TRUE)
summary(tm3, se = "boot")
texreg(tm3)

我尝试查看 texreg （在控制台上键入 texreg）并在第 35-38 行发现

for (i in 1:length(models)) {
cf <- models[[i]]@coef
se <- models[[i]]@se
pv <- models[[i]]@pvalues

如何设置 se 方法？我从未更改过 R 函数，如果这是我需要做的，你能建议一个链接来解释如何做吗？

4

2 回答 2

1

我已经更新了 texreg 源代码。您对“...”参数的建议应在下一个版本中实施。如果您有更多建议，您可能希望将它们发布在 texreg 论坛中。

最好的，菲利普

于 2012-11-28T09:06:44.210 回答
1

部分成功。我不太擅长 S4 类，但我能够通过下载包源并编辑前几行来获得一些工作extract.rq：

extract.rq <- function(model, include.nobs=TRUE, include.percentile=TRUE,...) {
  co <- summary(model, cov=TRUE,...)$coef[,1]
  names <- rownames(summary(model, cov=TRUE,...)$coef)
  se <- summary(model, cov=TRUE,...)$coef[,2]
  pval <- summary(model, cov=TRUE,...)$coef[,4]

  n <- length(summary(model,...)$resid)
  tau <- summary(model,...)$tau

然后从修改后的源构建和安装包似乎可以工作，如下所示：

library(quantreg)
library(texreg)
data(stackloss)
m <- rq(stack.loss ~ stack.x,.5)

texreg(m,se = "boot")
texreg(m,se = "iid")
于 2012-11-21T16:35:43.747 回答