老师，这个是我所有的代码，就是到了画boxplot就开始出现如上报错，我不知道咋解决
可能是你的数据有问题吧，用下这段代码试试:
################################################### #从芯片的原始数据获得表达数据 ################################################### GSE <-"GSE100935" files = list.files("D:/data", pattern = "GSM") files AffyBatch=ReadAffy(filenames=files) n.Sample = length(AffyBatch$sample) SampleNames=sampleNames(AffyBatch) SampleNames #获得样品名字，去掉不必要的信息 SampleNames<-substr(SampleNames,1,10) SampleNames ################第一个样本的灰度图############## pdf(file="image.1.pdf",width = 9,height = 24) image(AffyBatch[,1]) dev.off() ##################原始数据荧光强度box分图############## ###基于affy包的boxplot,函数默认不显示异常值 ###绘制箱线图 ##################################### #获取样品数量，并设置图片颜色 cols = rainbow(n.Sample) #绘图 pdf(file="RawData_expressionDistribution.pdf", width=24, height=18) par(cex = 0.7) if(n.Sample>40) par(cex = 0.5) boxplot(AffyBatch,col = cols,xaxt="n") axis(side=1,at=1:n.Sample,SampleNames,las=2) dev.off() #################RNAdeg 完整性查看######################### RNAdeg=AffyRNAdeg(AffyBatch) #查看 summaryAffyRNAdeg(RNAdeg) #绘图 pdf(file="RNAdeg.pdf",width = 12,height = 9) plotAffyRNAdeg(RNAdeg,cols= cols) legend("topleft",legend=SampleNames,col= cols,lty=1,box.lty=0,bg=NA) box() dev.off()
老师，还是会出现报错，这是为什么呀w(ﾟДﾟ)w