Working with low abundance probesets !!
2
0
Entering edit mode
Sharon Anbu ▴ 480
@sharon-anbu-1524
Last seen 10.2 years ago
Hi All, I am working with the Affy data set for which the max fold change is around 1.7. Most of the probesets have a fold change between 1.2 and 1.5. I have done the standard data analysis procedure reccommanded by many people in Bioconductor (starting with RMA extraction, difference between 2-groups, t-test, p-val & corrected p-val etc). After doing this, I have got only 20 significant genes, which has no biological significance with the experiment. Now, I am thinking of trying other measures instead of difference between groups. What kinds of measures I should try? Is there any method addressing the issue of detecting low abundance probesets in Bioconductor? Thanks in Advance. Regards, Sharon
affy affy • 908 views
ADD COMMENT
0
Entering edit mode
@sean-davis-490
Last seen 3 months ago
United States
On 12/1/05 8:50 AM, "Sharon Anbu" <sharonanandhi at="" gmail.com=""> wrote: > Hi All, > > I am working with the Affy data set for which the max fold change is > around 1.7. Most of the probesets have a fold change between 1.2 and > 1.5. I have done the standard data analysis procedure reccommanded by > many people in Bioconductor (starting with RMA extraction, difference > between 2-groups, t-test, p-val & corrected p-val etc). After doing > this, I have got only 20 significant genes, which has no biological > significance with the experiment. How many arrays have you run? If you know that you are looking for fold changes of 1.2-1.7 (I'm not sure how you know which genes are true and showing these changes, but...), you will likely need a stable system (little biological variability within groups) or a "relatively" large n. You could do a simple power calculation to determine what your power is to detect a difference between your two groups given a "typical" variance and with different numbers of slides. If you have good power with the number of slides that you are using, then you may be seeing true negative results. If not, then increasing the number of arrays might be a good bet. All that said, if you have a set of genes that you are interested in, you could go directly from "fold-change" on the array (ignore statistics such as p-value altogether) and move directly to validation with PCR-based methods. Sean
ADD COMMENT
0
Entering edit mode
An embedded and charset-unspecified text was scrubbed... Name: not available Url: https://stat.ethz.ch/pipermail/bioconductor/attachments/20051202/ 344e25df/attachment.pl
ADD REPLY
0
Entering edit mode
Dear List I am trying to find out if the mean expression values of Q-PCR differ for a group of 192 genes in three different sample groups. I am mainly interested in the F value, not any particular contrasts. I have created an exprSet out of my QPCR data using the following >myEset Expression Set (exprSet) with 192 genes 40 samples phenoData object with 1 variables and 40 cases varLabels cov1: read from file How can I proceed now? I thought of using linear models and then obtain the F value from fit3$F.p.value for the fit, but I got the following errors: >treatments = factor(c (0,1,1,1,1,1,0,0,1,1,0,0,2,0,0,1,1,0,2,2,1,0,0,0,2,2,0,0,1,0,1,0,0,1,0 , 0,0,0,0,0),labels = c("N0","N1","N2")) >design = model.matrix(~0+treatments) >fit <- lmFit(exprs(myEset),design) >contrast.matrix = makeContrasts(N1-N0,levels = design) >fit2 = contrasts.fit(fit,contrast.matrix) >fit3 = eBayes(fit2) Error in ebayes(fit = fit, proportion = proportion, stdev.coef.lim = stdev.coef.lim) : No residual degrees of freedom in linear model fits Best, David
ADD REPLY
0
Entering edit mode
Sharon Anbu ▴ 480
@sharon-anbu-1524
Last seen 10.2 years ago
An embedded and charset-unspecified text was scrubbed... Name: not available Url: https://stat.ethz.ch/pipermail/bioconductor/attachments/20051201/ 31b73920/attachment.pl
ADD COMMENT

Login before adding your answer.

Traffic: 922 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