Rankprod- gene appearing in both up and down regulated lists?
0
0
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: [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.
RankProd RankProd • 871 views
ADD COMMENT

Login before adding your answer.

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