Entering edit mode
Guest User ★ 12k@guest-user-4897
Last seen 7.4 years ago
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:  LC_COLLATE=English_United States.1252 LC_CTYPE=English_United States.1252 LC_MONETARY=English_United States.1252  LC_NUMERIC=C LC_TIME=English_United States.1252 attached base packages:  grid parallel stats graphics grDevices utils datasets methods base other attached packages:  pd.hugene.2.0.st_3.8.0 heatmap.plus_1.3 RankProd_2.32.0  hugene20sttranscriptcluster.db_2.12.1 xlsx_0.5.1 xlsxjars_0.5.0  rJava_0.9-4 annotate_1.38.0 org.Hs.eg.db_2.9.0  GGally_0.4.4 reshape_0.8.4 plyr_1.8  ggplot2_0.9.3.1 biomaRt_2.16.0 puma_3.2.1  mclust_4.2 VennDiagram_1.6.4 scatterplot3d_0.3-33  annaffy_1.32.0 KEGG.db_2.9.1 GO.db_2.9.0  RSQLite_0.11.4 DBI_0.2-7 AnnotationDbi_1.22.6  gplots_2.11.3 KernSmooth_2.23-10 caTools_1.14  gdata_2.13.2 gtools_3.0.0 MASS_7.3-28  RColorBrewer_1.0-5 affyPLM_1.36.0 preprocessCore_1.22.0  simpleaffy_2.36.1 gcrma_2.32.0 genefilter_1.42.0  oligo_1.24.2 oligoClasses_1.22.0 affy_1.38.1  marray_1.38.0 limma_3.16.7 Biobase_2.20.1  BiocGenerics_0.6.0 loaded via a namespace (and not attached):  affxparser_1.32.3 affyio_1.28.0 BiocInstaller_1.10.3 Biostrings_2.28.0 bit_1.1-10 bitops_1.0-5  codetools_0.2-8 colorspace_1.2-2 dichromat_2.0-0 digest_0.6.3 dynamicTreeCut_1.21 ff_2.2-11  foreach_1.4.1 GenomicRanges_1.12.4 gtable_0.1.2 impute_1.34.0 IRanges_1.18.3 iterators_1.0.6  labeling_0.2 munsell_0.4.2 proto_0.3-10 RCurl_1.95-4.1 reshape2_1.2.2 scales_0.2.3  splines_3.0.1 stats4_3.0.1 stringr_0.6.2 survival_2.37-4 tools_3.0.1 WGCNA_1.27-1  XML_3.98-1.1 xtable_1.7-1 zlibbioc_1.6.0 -- Sent via the guest posting facility at bioconductor.org.