R中用bayesQR程序包做贝叶斯分位数回归时遇到如下问题,代码如下
library(bayesQR)
mydata<-read.delim("clipboard") (mydata一共包含11列数据,第一列为injury,是因变量,后面10列是自变量。)
quantvec<-seq(0.05,0.95,0.05)
out<- bayesQR(injury~ 0 + ., data =mydata, quantile =quantvec,nrdaw=50000,keep=10)
pre<-predict(out,X=mydata[,2:11],type='response')
最后一行代码敲完后,出现如下错误:Error in pvec[preds] : 类别为'list'的下标不对。我的目的是得到ROC曲线,但是运行到pre<-predict(out,X=mydata[,2:11],type='response')这一行就出现刚提到的错误了。把代码换成pre<-predict(out,X=mydata[,1:10],type='response')虽然没有错误,但是不是我要的结果。请问这个错误该怎么解决啊?求高手帮忙看看。