beadarray: mRNA QC help
0
0
Entering edit mode
Levi Waldron ▴ 80
@levi-waldron-6357
Last seen 10.2 years ago
Dear Abhi, it looks to me like your chips had low hybridization, based both on the low percentage of probes detected and on the low interquartile range of intensities. Fig. 2 from PMID23136189 (admittedly my publication) shows the range of IQRs for samples within several published studies using Illumina BeadArrays for FFPE tissues. The normal IQR is around 2, which corresponds to around 50-60% present at p<0.05. Your typical IQR is around 1, and I can't directly read percent present but it looks around 20-30%. I would be concerned, and do some additional checks: 1) if you have technical replicates, check them by MA plot and sample-wise correlations, 2) MA plots against the median pseudochip, 3) check for detection and high intensity of the positive control probes (labelled CY3_HYB and HOUSEKEEPING on Illumina chips I've looked at), and 4) run the ffpe::sortedIqrPlot() function to see whether IQR or percent present are related to correlation to median pseudochip in your experiment. Are you using FFPE specimens? Sincerely, Levi On Thu, Jan 23, 2014 at 6:38 PM, Abhishek Pratap <apratap@sagebase.org>wrote: > Hi All > > I am trying to analyze some mRNA illumina bead level data through bead > array. Based on the detection p-values(plot attached): if I read it > right significant # probes > have no expression which concerns me but I could have easily missed > some important step. > > I am attaching the two box plots 1.) intensity and 2.) probe level > detection p-values + code > > Let me know your opinion. > > Thanks! > -Abhi > > beadlevel_data <- readIllumina(dir=dataDir, useImages=F, > illuminaAnnotation="Humanv4") > > #intensity boxplot > boxplot(beadlevel_data, transFun = logGreenChannelTransform, col = "green", > ylab = expression(log[2](intensity)), las=2, outline=FALSE, > main= "Array intensities") > > #summarize : detection p-val boxplot > BSData <- summarize(beadlevel_data) > det = calculateDetection(BSData) > boxplot(det, main='dset1_mRNA_probe_detection_pvals (Null: there is no > expression) ', ylab='p-value(detection score)',xaxt='n', > xlab='individual samples/array' ) > > _______________________________________________ > Bioconductor mailing list > Bioconductor@r-project.org > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: > http://news.gmane.org/gmane.science.biology.informatics.conductor > -- Levi Waldron Assistant Professor of Biostatistics City University of New York School of Public Health, Hunter College 2180 3rd Ave Rm 538 New York NY 10035-4003 phone: 212-396-7747 -- Levi Waldron Assistant Professor of Biostatistics City University of New York School of Public Health, Hunter College 2180 3rd Ave Rm 538 New York NY 10035-4003 phone: 212-396-7747 [[alternative HTML version deleted]]
probe ffpe probe ffpe • 1.0k views
ADD COMMENT

Login before adding your answer.

Traffic: 914 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6