我正在尝试获取分位数回归的自举置信区间。我的回归反应因素主要是分类变量。当我运行代码boot.rq时，出现以下错误；
由boot.rq.xy（x，y，U，tau）中的coercionError引入的NA： 外部函数调用中的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)
该插图来自该论坛先前的讨论；