[BioC] PICS : makeRangedDataOutput
Zhu, Lihua (Julie)
Julie.Zhu at umassmed.edu
Fri Oct 21 23:52:41 CEST 2011
Arnaud,
The makeRangedDataOutput outputs null when it should output a RangedData
object with peaks having score >=2. Here is the example you ran in Bioc2011
for ER dataset. Did I miss something? Thanks!
Best regards,
Julie
RD <- as(pics, "RangedData")
RD
RangedData with 304 rows and 1 value column across 2 spaces
space ranges | score
<factor> <IRanges> | <numeric>
pics1 chr21 [42669615, 42669741] | 98.27149
pics2 chr21 [42659913, 42660054] | 88.79954
pics3 chr21 [41898880, 41899010] | 85.83956
pics4 chr21 [45574708, 45574826] | 81.39958
pics5 chr21 [45631173, 45631311] | 79.91959
pics8 chr21 [15504544, 15504699] | 59.93969
pics11 chr21 [42659615, 42659758] | 46.61976
pics12 chr21 [15299644, 15299796] | 46.12643
pics13 chr21 [39208579, 39208729] | 43.65977
... ... ... ... ...
pics292 chr22 [33122494, 33122758] | 1.874657
pics293 chr22 [43400868, 43401220] | 1.849990
pics294 chr22 [38733327, 38733569] | 1.553992
pics298 chr22 [14702383, 14702743] | 1.479992
pics299 chr22 [17040892, 17041186] | 1.479992
pics300 chr22 [22827324, 22827654] | 1.479992
pics302 chr22 [15229177, 15229465] | 1.331993
pics303 chr22 [15233278, 15233637] | 1.257994
pics304 chr22 [15231777, 15232120] | 1.109994
myFilter <- list(score=c(2, Inf), delta=c(50, 300), se=c(0,50))
RD <-makeRangedDataOutput(pics, type="bed", filter=myFilter)
RD
null
R version 2.14.0 beta (2011-10-17 r57293)
Platform: i386-apple-darwin9.8.0/i386 (32-bit)
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/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] ChipSeq_1.0.0 MotIV_1.7.0 rGADEM_2.1.0
seqLogo_1.19.0 BiocInstaller_1.1.28
[6] snowfall_1.84 snow_0.3-7 PICS_1.7.0
BSgenome_1.21.7 Biostrings_2.21.11
[11] GenomicRanges_1.5.51 IRanges_1.11.31
loaded via a namespace (and not attached):
[1] lattice_0.19-33 tools_2.14.0
More information about the Bioconductor
mailing list