Hello,

taus<-seq(from = .05, to = .95, by = 0.05)

QR.coef <- coef(rq(S~B,tau=taus))

print(QR.coef)

lprq <- function(x, y, h, m=19 , tau=.5)

{

xx <- seq(min(x),max(x),length=m)

fv <- xx

dv <- xx

for(i in 1:length(xx)) {

z <- x - xx[i]

wx <- dnorm(z/h)

r <- rq(y~z, weights=wx, tau=tau, ci=FALSE)

fv[i] <- r$coef[1.]

dv[i] <- r$coef[2.]

}

data.frame(dv = dv)

}

QQR.coef <- as.data.frame(matrix(0, ncol = 19, nrow = 19))

#Note: 0.05 in below loop is the bandwidth that can be adjusted

for (i in 1:19){

x<-lprq(B,S,0.05,tau=taus[i])

QQR.coef[,i]<-x

}

the rq in the lprq function of mine is throwing the below error. Am I doing something wrong?

**Error:** Error in rq.fit.br(wx, wy, tau = tau, ...) : Singular design matrix Error