Как установить se=boot при использовании texreg для квантильной регрессии в R?

Question

Я использую квантильную регрессию (пакет quantreg) и используя texreg создать выход латекса из моих моделей.

Я заинтересован в начальной загрузке se и set 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

Я обновил исходный код texreg. Ваше предложение с аргументом "..." должно быть реализовано в следующем выпуске. Если у вас есть другие предложения, вы можете опубликовать их на форуме texreg.

Лучший, Филипп

Answer 2

Частичный успех. Я не очень разбираюсь в классах 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")