在 R 中使用 texreg 进行分位数回归时如何设置 se=boot？

Question

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

我对 bootstrapped 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 函数，如果这是我需要做的，你能建议一个解释如何做的链接吗？

Answer 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")

Answer 2

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

最好的，菲利普