Доверительный интервал для квантильной регрессии с использованием бутстрапа, категориальных переменных

Я пытаюсь получить доверительные интервалы начальной загрузки для квантильной регрессии. Мои факторы реакции регрессии - это в основном категориальные переменные. Когда я запускаю код boot.rq, я получаю следующую ошибку:

NAs, введенные coercionError в boot.rq.xy (x, y, U, tau): NA / NaN / Inf в вызове внешней функции (arg 7)

Ниже приведен пример, в котором я хотел бы вычислить доверительные интервалы для типа;

library(car)
library(quantreg)
library(boot)
newdata<-Prestige
mod2 = rq(income ~ education+ prestige + women + type, data=newdata)
dat2 <- newdata[,c("education","women","prestige","census","type")]
mod3=boot.rq(cbind(1,dat2),newdata$income,tau=0.5, R=10000)

Иллюстрация из предыдущего обсуждения на этом форуме;

Доверительный интервал для квантильной регрессии с использованием начальной загрузки


person ALEXANDER MACHARIA    schedule 13.11.2019    source источник