Question: Rankprod- gene appearing in both up and down regulated lists?
0
gravatar for Guest User
6.2 years ago by
Guest User12k
Guest User12k wrote:
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 • 688 views
ADD COMMENTlink written 6.2 years ago by Guest User12k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 253 users visited in the last hour