[BioC] No replicates and differential analysis !!

James W. MacDonald jmacdon at med.umich.edu
Wed Jan 25 16:24:31 CET 2006


Nicolas Servant wrote:
> Thanks for your answer,
> But in this case, i have to choose a fold change threshold ! And it is 
> supported that the FC tends to be greater at low expression levels.
> For instance a FC greater than 2 for expression values near 50 is 
> readily seen, but it is low probability to observe FC greater than 2 for 
> expression values near 1000
> So i would like to use a more robust approach.

With only two samples, you are stuck with fold changes. However, you 
might be able to make your results more robust by filtering out those 
genes that you think are too small. I often use kOverA() in the 
genefilter package to do this.

Best,

Jim


> 
> Regards,
> Nicolas S.
> 
> Sean Davis wrote:
> 
> 
>>On 1/25/06 8:34 AM, "Nicolas Servant" <Nicolas.Servant at curie.fr> wrote:
>>
>> 
>>
>>
>>>Hello,
>>>
>>>Does anybody know a R package or function to compare expression level
>>>(affy data) of two groups with no replicates in each group ? In fact,
>>>just compare one array to an other.
>>>The purpose is to find differentially expressed genes.
>>>We cannot used statistical test (not enougth replicates), but we can
>>>used graphical approach based on scatter plot, and outliers detection
>>>approach.
>>>   
>>>
>>
>>Simply take array A and divide it by array B.  Then rank the genes by those
>>ratios.  
>>
>>Sean
>>
>>_______________________________________________
>>Bioconductor mailing list
>>Bioconductor at stat.math.ethz.ch
>>https://stat.ethz.ch/mailman/listinfo/bioconductor
>>
>>
>> 
>>
> 
> 
> 


-- 
James W. MacDonald
Affymetrix and cDNA Microarray Core
University of Michigan Cancer Center
1500 E. Medical Center Drive
7410 CCGC
Ann Arbor MI 48109
734-647-5623



More information about the Bioconductor mailing list