[BioC] Question on makeVennDiagram function in ChIPpeakAnno
Noah Dowell
noahd at ucla.edu
Sun Jun 13 03:04:40 CEST 2010
Hello All,
I ran the following function and received the following output but it is unclear to me why I have a negative number for vennCounts (and no p-value). Any ideas? My output for the example in the vignette is correct but running the function with my data I am getting the following:
> makeVennDiagram(RangedDataList(smallAnnoTFA, smallAnnoTFB), NameOfPeaks=c("TFA", "TFB"), maxgap = 0, totalTest= 100, cex = 1, counts.col = "purple")
$p.value
[1] NaN
$vennCounts
TFA TFB Counts
[1,] 0 0 -2667
[2,] 0 1 755
[3,] 1 0 1334
[4,] 1 1 678
attr(,"class")
[1] "VennCounts"
There were 17 warnings (use warnings() to see them)
> warnings()
Warning messages:
1: In .local(query, subject, maxgap, minoverlap, type, select, ... :
argument 'multiple' is deprecated; use 'select'.
2: In .local(query, subject, maxgap, minoverlap, type, select, ... :
argument 'multiple' is deprecated; use 'select'.
3: In .local(query, subject, maxgap, minoverlap, type, select, ... :
argument 'multiple' is deprecated; use 'select'.
Thank you for any help.
Noah
> sessionInfo()
R version 2.11.0 (2010-04-22)
i386-apple-darwin9.8.0
locale:
[1] en_US.UTF-8/en_US.UTF-8/C/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] grid stats graphics grDevices utils datasets methods base
other attached packages:
[1] org.Sc.sgd.db_2.4.1 rtracklayer_1.8.1
[3] RCurl_1.3-1 bitops_1.0-4.1
[5] Starr_1.4.0 affxparser_1.20.0
[7] affy_1.26.0 Ringo_1.12.0
[9] Matrix_0.999375-38 lattice_0.18-5
[11] RColorBrewer_1.0-2 ChIPpeakAnno_1.4.0
[13] limma_3.4.0 org.Hs.eg.db_2.4.1
[15] GO.db_2.4.1 RSQLite_0.8-4
[17] DBI_0.2-5 AnnotationDbi_1.10.0
[19] BSgenome.Ecoli.NCBI.20080805_1.3.16 BSgenome_1.16.0
[21] GenomicRanges_1.0.1 Biostrings_2.16.0
[23] IRanges_1.6.0 multtest_2.4.0
[25] Biobase_2.8.0 biomaRt_2.4.0
loaded via a namespace (and not attached):
[1] affyio_1.16.0 annotate_1.26.0 genefilter_1.30.0 MASS_7.3-5
[5] preprocessCore_1.10.0 pspline_1.0-14 splines_2.11.0 survival_2.35-8
[9] tools_2.11.0 XML_2.8-1 xtable_1.5-6
>
More information about the Bioconductor
mailing list