library(ffpeExampleData)
data(lumibatch.GSE17565)
QC <- sampleQC(lumibatch.GSE17565,xaxis="index",cor.to="pseudochip",QCmeasure="IQR")
##sort samples
QCvsRNA <- data.frame(inputRNA.ng=lumibatch.GSE17565$inputRNA.ng,rejectQC=QC$rejectQC)
QCvsRNA <- QCvsRNA[order(QCvsRNA$rejectQC,-QCvsRNA$inputRNA.ng),]
##QC rejects samples with lowest input RNA concentration\n
par(mgp=c(4,2,0))
dotchart(log10(QCvsRNA$inputRNA.ng),
QCvsRNA$rejectQC,
xlab="log10(RNA conc. in ng)",
ylab="rejected?",
col=ifelse(QCvsRNA$rejectQC,"red","black"))
Run the code above in your browser using DataLab