[BioC] Rankprod- gene appearing in both up and down regulated lists?

Adam Cornwell [guest] guest at bioconductor.org
Wed Aug 21 19:18:27 CEST 2013


I am using RankProd for differential expression of array data. I've confirmed that all of my input gene names, rownames, etc are all gene-level and contain no duplicates. Yet, when I run RPadvance and look at the output of topGene(rankprodResult, cutoff=0.05, logged = TRUE, logbase=2)$Table1
and
topGene(rankprodResult, cutoff=0.05, logged = TRUE, logbase=2)$Table2
I see that there is one gene which appears in both of those lists- with the same fold change but different p-values for the occurrence in Table1 and Table2. I then went and confirmed that my initial input did not contain duplicates again- which is indeed the case. The foldchange is 1.2744, which seems to indicate that the gene is actually downregulated and should only be in Table2. 
rankprod.table1[rankprod.table1[,"FC:(class1/class2)"] >= 1,]
        gene.index            RP/Rsum FC:(class1/class2)                pfp            P.value 
        19582.0000          2031.3112             1.2744             0.0244             0.0017

What's going on here? I've used this method a few times before, and haven't seen this. I'm using the newest version of RankProd.

(While I'm asking this already- are all the p-values of 0 from underflowing or from the method used to calculate those values? It feels quite strange sometimes to be presenting all these p-values of 0 in my results.)

 -- output of sessionInfo(): 

R version 3.0.1 (2013-05-16)
Platform: x86_64-w64-mingw32/x64 (64-bit)

locale:
[1] LC_COLLATE=English_United States.1252  LC_CTYPE=English_United States.1252    LC_MONETARY=English_United States.1252
[4] LC_NUMERIC=C                           LC_TIME=English_United States.1252    

attached base packages:
[1] grid      parallel  stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
 [1] pd.hugene.2.0.st_3.8.0                heatmap.plus_1.3                      RankProd_2.32.0                      
 [4] hugene20sttranscriptcluster.db_2.12.1 xlsx_0.5.1                            xlsxjars_0.5.0                       
 [7] rJava_0.9-4                           annotate_1.38.0                       org.Hs.eg.db_2.9.0                   
[10] GGally_0.4.4                          reshape_0.8.4                         plyr_1.8                             
[13] ggplot2_0.9.3.1                       biomaRt_2.16.0                        puma_3.2.1                           
[16] mclust_4.2                            VennDiagram_1.6.4                     scatterplot3d_0.3-33                 
[19] annaffy_1.32.0                        KEGG.db_2.9.1                         GO.db_2.9.0                          
[22] RSQLite_0.11.4                        DBI_0.2-7                             AnnotationDbi_1.22.6                 
[25] gplots_2.11.3                         KernSmooth_2.23-10                    caTools_1.14                         
[28] gdata_2.13.2                          gtools_3.0.0                          MASS_7.3-28                          
[31] RColorBrewer_1.0-5                    affyPLM_1.36.0                        preprocessCore_1.22.0                
[34] simpleaffy_2.36.1                     gcrma_2.32.0                          genefilter_1.42.0                    
[37] oligo_1.24.2                          oligoClasses_1.22.0                   affy_1.38.1                          
[40] marray_1.38.0                         limma_3.16.7                          Biobase_2.20.1                       
[43] BiocGenerics_0.6.0                   

loaded via a namespace (and not attached):
 [1] affxparser_1.32.3    affyio_1.28.0        BiocInstaller_1.10.3 Biostrings_2.28.0    bit_1.1-10           bitops_1.0-5        
 [7] codetools_0.2-8      colorspace_1.2-2     dichromat_2.0-0      digest_0.6.3         dynamicTreeCut_1.21  ff_2.2-11           
[13] foreach_1.4.1        GenomicRanges_1.12.4 gtable_0.1.2         impute_1.34.0        IRanges_1.18.3       iterators_1.0.6     
[19] labeling_0.2         munsell_0.4.2        proto_0.3-10         RCurl_1.95-4.1       reshape2_1.2.2       scales_0.2.3        
[25] splines_3.0.1        stats4_3.0.1         stringr_0.6.2        survival_2.37-4      tools_3.0.1          WGCNA_1.27-1        
[31] XML_3.98-1.1         xtable_1.7-1         zlibbioc_1.6.0      

--
Sent via the guest posting facility at bioconductor.org.



More information about the Bioconductor mailing list