[BioC] Very Different group sizes in DE of two conditions

Gordon K Smyth smyth at wehi.EDU.AU
Thu Jul 5 02:22:32 CEST 2012


> Date: Tue,  3 Jul 2012 03:08:25 -0700 (PDT)
> From: "Moriah [guest]" <guest at bioconductor.org>
> To: bioconductor at r-project.org, moriahcohen at gmail.com
> Subject: [BioC] Very Different group sizes in  DE of two conditions
>
> Hi,
>
> I would like to conduct DE analysis of miRNA and mRNA using DESeq and 
> EdgeR on two conditions (Breast cancer subtypes) the problem is that for 
> one condition I have 122 samples and for the second only 22.  I would 
> like to know if it is O.K to run DE on such different size groups ?

Yes, it is no problem whatsoever.

Gordon

> What I tried to do is to run 100 different permutations of 22 samples 
> out of 122 and to try to run DE and search for differentially expressed 
> miRNA in all 100 different permutations against the 22 from the second 
> condition.
>
> I got only one miRNA that is differentially expressed in DESeq in all 
> permutaions but when I run 122 samples against 22 samples I got no such 
> miRNA. There is a difference between EdgeR results as well.
>
> Is my way of creating samller groups is applicable? Or is there any 
> other way to deal with such different group sizes?
>
> Thanks!
>
> -- output of sessionInfo():
>
> R version 2.15.1 (2012-06-22)
> Platform: x86_64-pc-linux-gnu (64-bit)
>
> locale:
> [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C
> [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8
> [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8
> [7] LC_PAPER=C                 LC_NAME=C
> [9] LC_ADDRESS=C               LC_TELEPHONE=C
> [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
>
> attached base packages:
> [1] stats     graphics  grDevices utils     datasets  methods   base
>
> other attached packages:
> [1] edgeR_2.6.7    limma_3.12.1   DESeq_1.8.3    locfit_1.5-8   Biobase_2.14.0
>
> loaded via a namespace (and not attached):
> [1] annotate_1.34.0      AnnotationDbi_1.18.1 BiocGenerics_0.2.0
> [4] DBI_0.2-5            genefilter_1.38.0    geneplotter_1.34.0
> [7] grid_2.15.1          IRanges_1.14.3       lattice_0.20-6
> [10] RColorBrewer_1.0-5   RSQLite_0.11.1       splines_2.15.1
> [13] stats4_2.15.1        survival_2.36-14     xtable_1.7-0
>>

______________________________________________________________________
The information in this email is confidential and intend...{{dropped:4}}



More information about the Bioconductor mailing list