x<-read.csv("d1.csv");
attach(x);
library(quantreg);
"lprq"<-function(GPR,CRB,h,m=10,tau=0.2){
  xx<-seq(min(GPR),max(GPR),length=m)
  fv<-xx
  dv<-xx
  for(i in 1:length(xx)){
    z<-GPR-xx[i]
    wx<-dnorm(z/h)
    r<-rq(CRB~z,weights=wx,tau=tau,ci=FALSE)
    fv[i]<-r$coef[1]
    dv[i]<-r$coef[2]}
  list(xx=xx,fv=fv,dv=dv)}
plot(GPR,CRB,xlab="GPR",ylab="CRB",col=5,font=2,font.lab=2)
title(main="tau=0.2",cex.main=1)
hs<-c(1,2,3,4)
for(i in hs){
  h=hs[i]
  fit<-lprq(GPR,CRB,h=h,tau=0.2)
  lines(fit$xx,fit$fv,lty=i,lwd=3,col=i)}
legend(45,-70,c("h=0.05","h=0.1","h=0.15","h=0.2"),col=1:length(hs),lty=1:length(hs))
大家好，这是我的R语言程序，想利用数据进行非参数分位数回归，在运行到倒数第二步，也就是  fit<-lprq(GPR,CRB,h=h,tau=0.2)
  lines(fit$xx,fit$fv,lty=i,lwd=3,col=i)}的时候，出错，显示Error in rq.fit.br(wx, wy, tau = tau, ...) : Singular design matrix。
小白一枚，请大家指教，谢谢！
数据请见附件