Question: Merging RandData object with names on the IRanges part
0
gravatar for Ulrike Goebel
10.2 years ago by
Ulrike Goebel110
Ulrike Goebel110 wrote:
Dear list, I would like to do the following: Read an output file of BWA (SAM format) in "chunks" and incrementally build a RangedData object from the chunks (by 'rbind') . Ultimately that should be used to get the number of reads per annotated transcript/region, but this is not the question here. Assume as an example: t1 <- RangedData(IRanges(start=c(7828367, 7828552,4121953), end=c(7828402, 7828587, 4121988)), space=c("Chr1", "Chr1", "Chr3"), mapq=c(1,2,1),flag=c(3,4,5)) I can merge two copies of this by 'rbind(t1,t1)'. But: t2 <- RangedData(IRanges(start=c(7828367, 7828552,4121953), end=c(7828402, 7828587, 4121988), names=c("a", "b", "c")), space=c("Chr1", "Chr1", "Chr3"), mapq=c(1,2,1),flag=c(3,4,5)) (Here, I would like to keep the read names along with their positions in the IRanges object). > rbind(t2,t2) Error in validObject(.Object) : invalid class "RangedData" object: the names of the ranges must equal the rownames Am I doing something completely wrong here ? Or is it confusing two different meanings of 'names' ? BTW, I really like IRanges ! Ulrike > sessionInfo() R version 2.10.0 Under development (unstable) (2009-08-01 r49053) x86_64-unknown-linux-gnu locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8 [5] LC_MONETARY=C LC_MESSAGES=en_US.UTF-8 [7] LC_PAPER=en_US.UTF-8 LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C attached base packages: [1] grid stats graphics grDevices utils datasets methods [8] base other attached packages: [1] ChIPR_1.1.3 MASS_7.3-0 spatstat_1.16-1 [4] deldir_0.0-8 gpclib_1.4-4 mgcv_1.5-5 [7] convert_1.21.1 marray_1.23.0 matchprobes_1.17.0 [10] AnnotationDbi_1.7.11 Biostrings_2.13.29 TeachingDemos_2.4 [13] Ringo_1.9.8 Matrix_0.999375-30 lattice_0.17-25 [16] limma_2.19.2 RColorBrewer_1.0-2 Biobase_2.5.5 [19] IRanges_1.3.56 loaded via a namespace (and not attached): [1] affy_1.23.4 affyio_1.13.3 annotate_1.23.1 [4] DBI_0.2-4 genefilter_1.25.7 nlme_3.1-92 [7] preprocessCore_1.7.4 RSQLite_0.7-1 splines_2.10.0 [10] survival_2.35-4 tools_2.10.0 xtable_1.5-5 -- Dr. Ulrike Goebel Bioinformatics Support Max-Planck Institute for Plant Breeding Research Carl-von-Linne Weg 10 50829 Cologne Germany +49(0) 221 5062 121
iranges • 811 views
ADD COMMENTlink modified 10.2 years ago • written 10.2 years ago by Ulrike Goebel110
Answer: Merging RandData object with names on the IRanges part
0
gravatar for Martin Morgan
10.2 years ago by
Martin Morgan ♦♦ 23k
United States
Martin Morgan ♦♦ 23k wrote:
Ulrike Goebel wrote: > Dear list, > > I would like to do the following: > Read an output file of BWA (SAM format) in "chunks" and incrementally > build a RangedData object from > the chunks (by 'rbind') . Ultimately that should be used to get the > number of reads per annotated transcript/region, but this is not the > question here. Probably you've already worked through this, but, I was wondering whether reading these files in chunks is necessary? It seems like, even in very large files, using scan with it's 'what' (including NULL) argument, you could selectively extract the fairly minimal information required to construct ranges? Martin > Assume as an example: > t1 <- RangedData(IRanges(start=c(7828367, 7828552,4121953), > end=c(7828402, 7828587, 4121988)), space=c("Chr1", "Chr1", "Chr3"), > mapq=c(1,2,1),flag=c(3,4,5)) > > I can merge two copies of this by 'rbind(t1,t1)'. > > But: > t2 <- RangedData(IRanges(start=c(7828367, 7828552,4121953), > end=c(7828402, 7828587, 4121988), names=c("a", "b", "c")), > space=c("Chr1", "Chr1", "Chr3"), mapq=c(1,2,1),flag=c(3,4,5)) > (Here, I would like to keep the read names along with their positions in > the IRanges object). > >> rbind(t2,t2) > Error in validObject(.Object) : > invalid class "RangedData" object: the names of the ranges must equal > the rownames > > Am I doing something completely wrong here ? Or is it confusing two > different meanings of 'names' ? > > > BTW, I really like IRanges ! > > Ulrike >> sessionInfo() > R version 2.10.0 Under development (unstable) (2009-08-01 r49053) > x86_64-unknown-linux-gnu > > locale: > [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C > [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8 > [5] LC_MONETARY=C LC_MESSAGES=en_US.UTF-8 > [7] LC_PAPER=en_US.UTF-8 LC_NAME=C > [9] LC_ADDRESS=C LC_TELEPHONE=C > [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C > > attached base packages: > [1] grid stats graphics grDevices utils datasets methods > [8] base > > other attached packages: > [1] ChIPR_1.1.3 MASS_7.3-0 spatstat_1.16-1 > [4] deldir_0.0-8 gpclib_1.4-4 mgcv_1.5-5 > [7] convert_1.21.1 marray_1.23.0 matchprobes_1.17.0 > [10] AnnotationDbi_1.7.11 Biostrings_2.13.29 TeachingDemos_2.4 > [13] Ringo_1.9.8 Matrix_0.999375-30 lattice_0.17-25 > [16] limma_2.19.2 RColorBrewer_1.0-2 Biobase_2.5.5 > [19] IRanges_1.3.56 > > loaded via a namespace (and not attached): > [1] affy_1.23.4 affyio_1.13.3 annotate_1.23.1 > [4] DBI_0.2-4 genefilter_1.25.7 nlme_3.1-92 > [7] preprocessCore_1.7.4 RSQLite_0.7-1 splines_2.10.0 > [10] survival_2.35-4 tools_2.10.0 xtable_1.5-5 > > > >
ADD COMMENTlink written 10.2 years ago by Martin Morgan ♦♦ 23k
Answer: Merging RandData object with names on the IRanges part
0
gravatar for Patrick Aboyoun
10.2 years ago by
Patrick Aboyoun1.6k
United States
Patrick Aboyoun1.6k wrote:
Ulrike, First of all, I'm glad IRanges is useful for you. Second, thanks for finding a bug in the rbind method for RangedData objects. Because of developer oversight, the duplicate names in the ranges was being handled differently than the duplicate rownames in the values. This has been corrected in a recent svn check-in to IRanges in the BioC 2.5 code line. You can get this updated IRanges package (version 1.3.58) either through svn access or wait 24-48 hours for the updated IRanges package to be placed on bioconductor.org and downloadable via biocLite. > suppressMessages(library(IRanges)) > t1 <- RangedData(IRanges(start=c(7828367, 7828552,4121953), end=c(7828402, 7828587, 4121988)), space=c("Chr1", "Chr1", "Chr3"), mapq=c(1,2,1),flag=c(3,4,5)) > rbind(t1, t1) RangedData: 6 ranges by 2 columns on 2 sequences colnames(2): mapq flag names(2): Chr1 Chr3 > t2 <- RangedData(IRanges(start=c(7828367, 7828552,4121953), end=c(7828402, 7828587, 4121988), names=c("a", "b", "c")), space=c("Chr1", "Chr1", "Chr3"), mapq=c(1,2,1),flag=c(3,4,5)) > rbind(t2, t2) RangedData: 6 ranges by 2 columns on 2 sequences colnames(2): mapq flag names(2): Chr1 Chr3 > sessionInfo() R version 2.10.0 Under development (unstable) (2009-08-05 r49073) i386-apple-darwin9.7.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] stats graphics grDevices utils datasets methods base other attached packages: [1] IRanges_1.3.58 Patrick Ulrike Goebel wrote: > Dear list, > > I would like to do the following: > Read an output file of BWA (SAM format) in "chunks" and incrementally > build a RangedData object from > the chunks (by 'rbind') . Ultimately that should be used to get the > number of reads per annotated transcript/region, but this is not the > question here. > > Assume as an example: > t1 <- RangedData(IRanges(start=c(7828367, 7828552,4121953), > end=c(7828402, 7828587, 4121988)), space=c("Chr1", "Chr1", "Chr3"), > mapq=c(1,2,1),flag=c(3,4,5)) > > I can merge two copies of this by 'rbind(t1,t1)'. > > But: > t2 <- RangedData(IRanges(start=c(7828367, 7828552,4121953), > end=c(7828402, 7828587, 4121988), names=c("a", "b", "c")), > space=c("Chr1", "Chr1", "Chr3"), mapq=c(1,2,1),flag=c(3,4,5)) > (Here, I would like to keep the read names along with their positions > in the IRanges object). > > > rbind(t2,t2) > Error in validObject(.Object) : > invalid class "RangedData" object: the names of the ranges must equal > the rownames > > Am I doing something completely wrong here ? Or is it confusing two > different meanings of 'names' ? > > > BTW, I really like IRanges ! > > Ulrike > > sessionInfo() > R version 2.10.0 Under development (unstable) (2009-08-01 r49053) > x86_64-unknown-linux-gnu > > locale: > [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C > [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8 > [5] LC_MONETARY=C LC_MESSAGES=en_US.UTF-8 > [7] LC_PAPER=en_US.UTF-8 LC_NAME=C > [9] LC_ADDRESS=C LC_TELEPHONE=C > [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C > > attached base packages: > [1] grid stats graphics grDevices utils datasets methods > [8] base > > other attached packages: > [1] ChIPR_1.1.3 MASS_7.3-0 spatstat_1.16-1 > [4] deldir_0.0-8 gpclib_1.4-4 mgcv_1.5-5 > [7] convert_1.21.1 marray_1.23.0 matchprobes_1.17.0 > [10] AnnotationDbi_1.7.11 Biostrings_2.13.29 TeachingDemos_2.4 > [13] Ringo_1.9.8 Matrix_0.999375-30 lattice_0.17-25 > [16] limma_2.19.2 RColorBrewer_1.0-2 Biobase_2.5.5 > [19] IRanges_1.3.56 > > loaded via a namespace (and not attached): > [1] affy_1.23.4 affyio_1.13.3 annotate_1.23.1 > [4] DBI_0.2-4 genefilter_1.25.7 nlme_3.1-92 > [7] preprocessCore_1.7.4 RSQLite_0.7-1 splines_2.10.0 > [10] survival_2.35-4 tools_2.10.0 xtable_1.5-5 > > > >
ADD COMMENTlink written 10.2 years ago by Patrick Aboyoun1.6k
Answer: Merging RandData object with names on the IRanges part
0
gravatar for Michael Dondrup
10.2 years ago by
Michael Dondrup140 wrote:
Hi Ulrike, I think you can use function c(). you then can get the IRanges names preserved in your examples: > ranges(c(t1,t2))[[3]] IRanges object: start end width names [1] 7828367 7828402 36 a [2] 7828552 7828587 36 b Michael Am 20.08.2009 um 17:03 schrieb Ulrike Goebel: > Dear list, > > I would like to do the following: > Read an output file of BWA (SAM format) in "chunks" and > incrementally build a RangedData object from > the chunks (by 'rbind') . Ultimately that should be used to get the > number of reads per annotated transcript/region, but this is not the > question here. > > Assume as an example: > t1 <- RangedData(IRanges(start=c(7828367, 7828552,4121953), > end=c(7828402, 7828587, 4121988)), space=c("Chr1", "Chr1", "Chr3"), > mapq=c(1,2,1),flag=c(3,4,5)) > > I can merge two copies of this by 'rbind(t1,t1)'. > > But: > t2 <- RangedData(IRanges(start=c(7828367, 7828552,4121953), > end=c(7828402, 7828587, 4121988), names=c("a", "b", "c")), > space=c("Chr1", "Chr1", "Chr3"), mapq=c(1,2,1),flag=c(3,4,5)) > (Here, I would like to keep the read names along with their > positions in the IRanges object). > > > rbind(t2,t2) > Error in validObject(.Object) : > invalid class "RangedData" object: the names of the ranges must > equal the rownames > > Am I doing something completely wrong here ? Or is it confusing two > different meanings of 'names' ? > > > BTW, I really like IRanges ! > > Ulrike > > sessionInfo() > R version 2.10.0 Under development (unstable) (2009-08-01 r49053) > x86_64-unknown-linux-gnu > > locale: > [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C > [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8 > [5] LC_MONETARY=C LC_MESSAGES=en_US.UTF-8 > [7] LC_PAPER=en_US.UTF-8 LC_NAME=C > [9] LC_ADDRESS=C LC_TELEPHONE=C > [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C > > attached base packages: > [1] grid stats graphics grDevices utils datasets > methods > [8] base > > other attached packages: > [1] ChIPR_1.1.3 MASS_7.3-0 spatstat_1.16-1 > [4] deldir_0.0-8 gpclib_1.4-4 mgcv_1.5-5 > [7] convert_1.21.1 marray_1.23.0 matchprobes_1.17.0 > [10] AnnotationDbi_1.7.11 Biostrings_2.13.29 TeachingDemos_2.4 > [13] Ringo_1.9.8 Matrix_0.999375-30 lattice_0.17-25 > [16] limma_2.19.2 RColorBrewer_1.0-2 Biobase_2.5.5 > [19] IRanges_1.3.56 > > loaded via a namespace (and not attached): > [1] affy_1.23.4 affyio_1.13.3 annotate_1.23.1 > [4] DBI_0.2-4 genefilter_1.25.7 nlme_3.1-92 > [7] preprocessCore_1.7.4 RSQLite_0.7-1 splines_2.10.0 > [10] survival_2.35-4 tools_2.10.0 xtable_1.5-5 > > > > > -- > Dr. Ulrike Goebel > Bioinformatics Support > Max-Planck Institute for Plant Breeding Research > Carl-von-Linne Weg 10 > 50829 Cologne > Germany > +49(0) 221 5062 121 > > _______________________________________________ > Bioconductor mailing list > Bioconductor at stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor Michael Dondrup, Ph.D. Bergen Center for Computational Science Computational Biology Unit Unifob AS - Thorm?hlensgate 55, N-5008 Bergen, Norway Phone: +47 55584029 Fax: +47 55584295
ADD COMMENTlink written 10.2 years ago by Michael Dondrup140
On Thu, Aug 20, 2009 at 8:33 AM, Michael Dondrup < Michael.Dondrup@bccs.uib.no> wrote: > Hi Ulrike, > > I think you can use function c(). > you then can get the IRanges names preserved in your examples: > > > ranges(c(t1,t2))[[3]] > IRanges object: > start end width names > [1] 7828367 7828402 36 a > [2] 7828552 7828587 36 b > > Yes, this works, but just for the record, c() is not the same as rbind(). The c() concenates chromosome blocks into a single RangedData, without changing any of the blocks, where as the rbind() essentially merges the blocks from the same chromosome across multiple RangedData instances. > > Michael > > Am 20.08.2009 um 17:03 schrieb Ulrike Goebel: > > > Dear list, >> >> I would like to do the following: >> Read an output file of BWA (SAM format) in "chunks" and incrementally >> build a RangedData object from >> the chunks (by 'rbind') . Ultimately that should be used to get the number >> of reads per annotated transcript/region, but this is not the question here. >> >> Assume as an example: >> t1 <- RangedData(IRanges(start=c(7828367, 7828552,4121953), end=c(7828402, >> 7828587, 4121988)), space=c("Chr1", "Chr1", "Chr3"), >> mapq=c(1,2,1),flag=c(3,4,5)) >> >> I can merge two copies of this by 'rbind(t1,t1)'. >> >> But: >> t2 <- RangedData(IRanges(start=c(7828367, 7828552,4121953), end=c(7828402, >> 7828587, 4121988), names=c("a", "b", "c")), space=c("Chr1", "Chr1", "Chr3"), >> mapq=c(1,2,1),flag=c(3,4,5)) >> (Here, I would like to keep the read names along with their positions in >> the IRanges object). >> >> > rbind(t2,t2) >> Error in validObject(.Object) : >> invalid class "RangedData" object: the names of the ranges must equal the >> rownames >> >> Am I doing something completely wrong here ? Or is it confusing two >> different meanings of 'names' ? >> >> >> BTW, I really like IRanges ! >> >> Ulrike >> > sessionInfo() >> R version 2.10.0 Under development (unstable) (2009-08-01 r49053) >> x86_64-unknown-linux-gnu >> >> locale: >> [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C >> [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8 >> [5] LC_MONETARY=C LC_MESSAGES=en_US.UTF-8 >> [7] LC_PAPER=en_US.UTF-8 LC_NAME=C >> [9] LC_ADDRESS=C LC_TELEPHONE=C >> [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C >> >> attached base packages: >> [1] grid stats graphics grDevices utils datasets methods >> [8] base >> >> other attached packages: >> [1] ChIPR_1.1.3 MASS_7.3-0 spatstat_1.16-1 >> [4] deldir_0.0-8 gpclib_1.4-4 mgcv_1.5-5 >> [7] convert_1.21.1 marray_1.23.0 matchprobes_1.17.0 >> [10] AnnotationDbi_1.7.11 Biostrings_2.13.29 TeachingDemos_2.4 >> [13] Ringo_1.9.8 Matrix_0.999375-30 lattice_0.17-25 >> [16] limma_2.19.2 RColorBrewer_1.0-2 Biobase_2.5.5 >> [19] IRanges_1.3.56 >> >> loaded via a namespace (and not attached): >> [1] affy_1.23.4 affyio_1.13.3 annotate_1.23.1 >> [4] DBI_0.2-4 genefilter_1.25.7 nlme_3.1-92 >> [7] preprocessCore_1.7.4 RSQLite_0.7-1 splines_2.10.0 >> [10] survival_2.35-4 tools_2.10.0 xtable_1.5-5 >> >> >> >> >> -- >> Dr. Ulrike Goebel >> Bioinformatics Support >> Max-Planck Institute for Plant Breeding Research >> Carl-von-Linne Weg 10 >> 50829 Cologne >> Germany >> +49(0) 221 5062 121 >> >> _______________________________________________ >> Bioconductor mailing list >> Bioconductor@stat.math.ethz.ch >> https://stat.ethz.ch/mailman/listinfo/bioconductor >> Search the archives: >> http://news.gmane.org/gmane.science.biology.informatics.conductor >> > > Michael Dondrup, Ph.D. > Bergen Center for Computational Science > Computational Biology Unit > Unifob AS - Thormøhlensgate 55, N-5008 Bergen, Norway > Phone: +47 55584029 Fax: +47 55584295 > > > _______________________________________________ > Bioconductor mailing list > Bioconductor@stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: > http://news.gmane.org/gmane.science.biology.informatics.conductor > [[alternative HTML version deleted]]
ADD REPLYlink written 10.2 years ago by Michael Lawrence620
Answer: Merging RandData object with names on the IRanges part
0
gravatar for Ulrike Goebel
10.2 years ago by
Ulrike Goebel110
Ulrike Goebel110 wrote:
Thanks a lot for all your suggestions/help ! I will have a look at Rsamtools (I was already wondering whether something like this exists), and, yes, maybe it is really not necessary to read the file in chunks. Thanks again ! Ulrike Ulrike Goebel wrote: > Dear list, > > I would like to do the following: > Read an output file of BWA (SAM format) in "chunks" and incrementally > build a RangedData object from > the chunks (by 'rbind') . Ultimately that should be used to get the > number of reads per annotated transcript/region, but this is not the > question here. > > Assume as an example: > t1 <- RangedData(IRanges(start=c(7828367, 7828552,4121953), > end=c(7828402, 7828587, 4121988)), space=c("Chr1", "Chr1", "Chr3"), > mapq=c(1,2,1),flag=c(3,4,5)) > > I can merge two copies of this by 'rbind(t1,t1)'. > > But: > t2 <- RangedData(IRanges(start=c(7828367, 7828552,4121953), > end=c(7828402, 7828587, 4121988), names=c("a", "b", "c")), > space=c("Chr1", "Chr1", "Chr3"), mapq=c(1,2,1),flag=c(3,4,5)) > (Here, I would like to keep the read names along with their positions > in the IRanges object). > > > rbind(t2,t2) > Error in validObject(.Object) : > invalid class "RangedData" object: the names of the ranges must equal > the rownames > > Am I doing something completely wrong here ? Or is it confusing two > different meanings of 'names' ? > > > BTW, I really like IRanges ! > > Ulrike > > sessionInfo() > R version 2.10.0 Under development (unstable) (2009-08-01 r49053) > x86_64-unknown-linux-gnu > > locale: > [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C > [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8 > [5] LC_MONETARY=C LC_MESSAGES=en_US.UTF-8 > [7] LC_PAPER=en_US.UTF-8 LC_NAME=C > [9] LC_ADDRESS=C LC_TELEPHONE=C > [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C > > attached base packages: > [1] grid stats graphics grDevices utils datasets methods > [8] base > > other attached packages: > [1] ChIPR_1.1.3 MASS_7.3-0 spatstat_1.16-1 > [4] deldir_0.0-8 gpclib_1.4-4 mgcv_1.5-5 > [7] convert_1.21.1 marray_1.23.0 matchprobes_1.17.0 > [10] AnnotationDbi_1.7.11 Biostrings_2.13.29 TeachingDemos_2.4 > [13] Ringo_1.9.8 Matrix_0.999375-30 lattice_0.17-25 > [16] limma_2.19.2 RColorBrewer_1.0-2 Biobase_2.5.5 > [19] IRanges_1.3.56 > > loaded via a namespace (and not attached): > [1] affy_1.23.4 affyio_1.13.3 annotate_1.23.1 > [4] DBI_0.2-4 genefilter_1.25.7 nlme_3.1-92 > [7] preprocessCore_1.7.4 RSQLite_0.7-1 splines_2.10.0 > [10] survival_2.35-4 tools_2.10.0 xtable_1.5-5 > > > > -- Dr. Ulrike Goebel Bioinformatics Support Max-Planck Institute for Plant Breeding Research Carl-von-Linne Weg 10 50829 Cologne Germany +49(0) 221 5062 121
ADD COMMENTlink written 10.2 years ago by Ulrike Goebel110
Answer: Merging RandData object with names on the IRanges part
0
gravatar for Ulrike Goebel
10.2 years ago by
Ulrike Goebel110
Ulrike Goebel110 wrote:
Sorry, I guess this was a clash with the old version that was still installed. After removing it, and also removing all .Rnw files from the /doc directory, I could build and install it. And now the rbind works ! Have a nice weekend Ulrike Ulrike Goebel wrote: > Hi Patrick, > > just to inform you: I have checked out IRanges 1.3.58 from svn, but > building fails in the vignettes part: > > svn co --username=readonly --password=readonly > https://hedgehog.fhcrc.org/bioconductor/trunk/madman/Rpacks/IRanges > R CMD build IRanges > ... > ** building package indices ... > * DONE (IRanges) > * creating vignettes ... ERROR > > Attaching package: 'IRanges' > > > The following object(s) are masked from package:base : > > cbind, > duplicated, > Map, > mapply, > order, > pmax, > pmax.int, > pmin, > pmin.int, > rank, > rbind, > rep.int, > sapply, > sort, > table, > unique > > > Error: processing vignette 'IRangesOverview.Rnw' failed with diagnostics: > chunk 33 (label=ranges-reduce) > Error in extraArgsAsList(valid_argnames, ...) : > could not find function "anyDuplicated" > Execution halted > > Best, Ulrike -- Dr. Ulrike Goebel Bioinformatics Support Max-Planck Institute for Plant Breeding Research Carl-von-Linne Weg 10 50829 Cologne Germany +49(0) 221 5062 121
ADD COMMENTlink written 10.2 years ago by Ulrike Goebel110
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: 186 users visited in the last hour