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

Question

我正在运行分位数回归(quantreg包)，并使用texreg创建我的模型的latex输出。

我对bootstrapped s.e很感兴趣。并在摘要的选项中设置“se=”，但是当我使用texreg时，我得到的是“n.i.d”。s.e.

如何更改该选项？

下面是我正在做的事情：

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

如何设置s.e.方法？我从来没有改变过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论坛上。

最好的，菲利普