在 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 論壇。

最好的，菲利普