Hi, Currently the build for `recount` is failing in the Linux build machine http://www.bioconductor.org/checkResults/3.4/bioc-LATEST/recount/zin1-buildsrc.html. It's been failing with the same error message for a while. The build message ends with:
Setting options('download.file.method.GEOquery'='auto') Setting options('GEOquery.inmemory.gpl'=FALSE) trying URL 'http://duffel.rail.bio/recount/SRP009615/rse_gene.Rdata' Content type 'application/octet-stream' length 1383643 bytes (1.3 MB) ================================================== downloaded 1.3 MB Quitting from lines 400-406 (recount-quickstart.Rmd) Error: processing vignette 'recount-quickstart.Rmd' failed with diagnostics: UCSC library operation failed Execution halted The lines of code from the recount vignette this refers to are: ## Define expressed regions for study SRP009615, only for chromosome Y regions <- expressed_regions('SRP009615', 'chrY', cutoff = 5L, maxClusterGap = 3000L) ## Briefly explore the resulting regions regions The UCSC error is something I've seen in previous versions of rtracklayer, which have been addressed in rtracklayer 1.33.11. recount version 0.99.25 (which is the one from the build report) had under suggests rtracklayer (>= 1.33.11). I thought that this would mean that the build would fail if rtracklayer below version 1.33.11 was being used. I'm no longer sure, which is why in recount version 0.99.27 I moved rtracklayer to imports to test this. The latest build report shows that rtracklayer 1.33.11 did build correctly in the Linux machine. Everything seems to be working fine in my local tests and in R Travis https://travis-ci.org/leekgroup/recount/builds/150995800. Below is some code showing that everything is working as expected in my local Linux test. So I'm a tad confused and can't reproduce this issue when using rtracklayer version 1.33.11. Which is why I'm guessing that somehow the Linux build machine is not using rtracklayer >= 1.33.11. ## The url redirects as expected to the correct url > rtracklayer:::expandURL('http://duffel.rail.bio/recount/SRP009615/rse_gene.Rdata') [1] "http://idies.jhu.edu/recount/data/SRP009615/rse_gene.Rdata" > library('recount') Loading required package: SummarizedExperiment Loading required package: GenomicRanges Loading required package: stats4 Loading required package: BiocGenerics Loading required package: parallel Attaching package: ‘BiocGenerics’ The following objects are masked from ‘package:parallel’: clusterApply, clusterApplyLB, clusterCall, clusterEvalQ, clusterExport, clusterMap, parApply, parCapply, parLapply, parLapplyLB, parRapply, parSapply, parSapplyLB The following objects are masked from ‘package:stats’: IQR, mad, xtabs The following objects are masked from ‘package:base’: anyDuplicated, append, as.data.frame, cbind, colnames, do.call, duplicated, eval, evalq, Filter, Find, get, grep, grepl, intersect, is.unsorted, lapply, lengths, Map, mapply, match, mget, order, paste, pmax, pmax.int, pmin, pmin.int, Position, rank, rbind, Reduce, rownames, sapply, setdiff, sort, table, tapply, union, unique, unsplit, which, which.max, which.min Loading required package: S4Vectors Attaching package: ‘S4Vectors’ The following objects are masked from ‘package:base’: colMeans, colSums, expand.grid, rowMeans, rowSums Loading required package: IRanges Loading required package: GenomeInfoDb Loading required package: Biobase Welcome to Bioconductor Vignettes contain introductory material; view with 'browseVignettes()'. To cite Bioconductor, see 'citation("Biobase")', and for packages 'citation("pkgname")'. Setting options('download.file.method.GEOquery'='auto') Setting options('GEOquery.inmemory.gpl'=FALSE) > > ## Define expressed regions for study SRP009615, only for chromosome Y > regions <- expressed_regions('SRP009615', 'chrY', cutoff = 5L, + maxClusterGap = 3000L) 2016-08-09 14:16:59 loadCoverage: loading BigWig file http://duffel.rail.bio/recount/SRP009615/bw/mean_SRP009615.bw 2016-08-09 14:17:01 loadCoverage: applying the cutoff to the merged data 2016-08-09 14:17:02 filterData: originally there were 57227415 rows, now there are 57227415 rows. Meaning that 0 percent was filtered. 2016-08-09 14:17:02 findRegions: identifying potential segments 2016-08-09 14:17:02 findRegions: segmenting information 2016-08-09 14:17:02 .getSegmentsRle: segmenting with cutoff(s) 5 2016-08-09 14:17:02 findRegions: identifying candidate regions 2016-08-09 14:17:03 findRegions: identifying region clusters > > ## Briefly explore the resulting regions > regions GRanges object with 808 ranges and 6 metadata columns: seqnames ranges strand | value area <Rle> <IRanges> <Rle> | <numeric> <numeric> 1 chrY [2929794, 2929829] * | 14.7265009482702 530.154034137726 2 chrY [2956678, 2956701] * | 12.8106340567271 307.45521736145 3 chrY [2977203, 2977227] * | 5.34908433914185 133.727108478546 4 chrY [2977957, 2977994] * | 6.46976616508082 245.851114273071 5 chrY [2978850, 2978871] * | 5.7976552139629 127.548414707184 ... ... ... ... . ... ... 804 chrY [26614511, 26614546] * | 7.28189378314548 262.148176193237 805 chrY [26614548, 26614553] * | 5.48768202463786 32.9260921478271 806 chrY [26614779, 26614808] * | 6.64339276949565 199.301783084869 807 chrY [26626808, 26626848] * | 12.6038152648181 516.756425857544 808 chrY [26626971, 26627028] * | 14.1673366941255 821.705528259277 indexStart indexEnd cluster clusterL <integer> <integer> <Rle> <Rle> 1 2929794 2929829 1 36 2 2956678 2956701 2 24 3 2977203 2977227 3 2750 4 2977957 2977994 3 2750 5 2978850 2978871 3 2750 ... ... ... ... ... 804 26614511 26614546 224 298 805 26614548 26614553 224 298 806 26614779 26614808 224 298 807 26626808 26626848 225 221 808 26626971 26627028 225 221 ------- seqinfo: 1 sequence from an unspecified genome > > options(width = 120) > devtools::session_info() Session info ----------------------------------------------------------------------------------------------------------- setting value version R version 3.3.1 Patched (2016-07-11 r70891) system x86_64, linux-gnu ui X11 language (EN) collate en_US.UTF-8 tz <NA> date 2016-08-09 Packages --------------------------------------------------------------------------------------------------------------- package * version date source acepack 1.3-3.3 2014-11-24 CRAN (R 3.3.0) AnnotationDbi 1.35.4 2016-07-22 Bioconductor Biobase * 2.33.0 2016-05-08 Bioconductor BiocGenerics * 0.19.2 2016-07-22 Bioconductor BiocParallel 1.7.6 2016-08-09 cran (@1.7.6) biomaRt 2.29.2 2016-06-07 Bioconductor Biostrings 2.41.4 2016-07-22 Bioconductor bitops 1.0-6 2013-08-17 CRAN (R 3.3.0) BSgenome 1.41.2 2016-07-22 Bioconductor bumphunter 1.13.1 2016-07-22 Bioconductor chron 2.3-47 2015-06-24 CRAN (R 3.3.0) cluster 2.0.4 2016-04-18 CRAN (R 3.3.1) codetools 0.2-14 2015-07-15 CRAN (R 3.3.1) colorspace 1.2-6 2015-03-11 CRAN (R 3.3.0) data.table 1.9.6 2015-09-19 CRAN (R 3.3.1) DBI 0.4-1 2016-05-08 CRAN (R 3.3.0) derfinder 1.7.9 2016-07-22 Bioconductor derfinderHelper 1.7.3 2016-07-22 Bioconductor devtools 1.12.0 2016-06-24 CRAN (R 3.3.1) digest 0.6.9 2016-01-08 CRAN (R 3.3.0) doRNG 1.6 2014-03-07 CRAN (R 3.3.0) downloader 0.4 2015-07-09 cran (@0.4) foreach 1.4.3 2015-10-13 CRAN (R 3.3.0) foreign 0.8-66 2015-08-19 CRAN (R 3.3.1) Formula 1.2-1 2015-04-07 CRAN (R 3.3.0) GenomeInfoDb * 1.9.4 2016-07-22 Bioconductor GenomicAlignments 1.9.6 2016-07-22 Bioconductor GenomicFeatures 1.25.15 2016-07-22 Bioconductor GenomicFiles 1.9.11 2016-07-22 Bioconductor GenomicRanges * 1.25.93 2016-08-09 cran (@1.25.93) GEOquery 2.39.3 2016-06-15 Bioconductor ggplot2 2.1.0 2016-03-01 CRAN (R 3.3.1) gridExtra 2.2.1 2016-02-29 CRAN (R 3.3.0) gtable 0.2.0 2016-02-26 CRAN (R 3.3.0) Hmisc 3.17-4 2016-05-02 CRAN (R 3.3.1) httr 1.2.1 2016-07-03 cran (@1.2.1) IRanges * 2.7.12 2016-08-09 cran (@2.7.12) iterators 1.0.8 2015-10-13 CRAN (R 3.3.0) jsonlite 1.0 2016-07-01 CRAN (R 3.3.1) lattice 0.20-33 2015-07-14 CRAN (R 3.3.1) latticeExtra 0.6-28 2016-02-09 CRAN (R 3.3.0) locfit 1.5-9.1 2013-04-20 CRAN (R 3.3.0) magrittr 1.5 2014-11-22 CRAN (R 3.3.1) Matrix 1.2-6 2016-05-02 CRAN (R 3.3.1) matrixStats 0.50.2 2016-04-24 CRAN (R 3.3.1) memoise 1.0.0 2016-01-29 CRAN (R 3.3.0) munsell 0.4.3 2016-02-13 CRAN (R 3.3.0) nnet 7.3-12 2016-02-02 CRAN (R 3.3.1) pkgmaker 0.22 2014-05-14 CRAN (R 3.3.0) plyr 1.8.4 2016-06-08 CRAN (R 3.3.1) qvalue 2.5.2 2016-07-22 Bioconductor R6 2.1.2 2016-01-26 CRAN (R 3.3.0) RColorBrewer 1.1-2 2014-12-07 CRAN (R 3.3.0) Rcpp 0.12.6 2016-07-19 CRAN (R 3.3.1) RCurl 1.95-4.8 2016-03-01 CRAN (R 3.3.0) recount * 0.99.27 2016-08-09 Github (leekgroup/recount@922e7a6) registry 0.3 2015-07-08 CRAN (R 3.3.0) rentrez 1.0.2 2016-04-21 cran (@1.0.2) reshape2 1.4.1 2014-12-06 CRAN (R 3.3.0) rngtools 1.2.4 2014-03-06 CRAN (R 3.3.0) rpart 4.1-10 2015-06-29 CRAN (R 3.3.1) Rsamtools 1.25.0 2016-07-22 Bioconductor RSQLite 1.0.0 2014-10-25 CRAN (R 3.3.0) rtracklayer 1.33.11 2016-07-22 Github (Bioconductor-mirror/rtracklayer@c12bd86) S4Vectors * 0.11.10 2016-07-22 Bioconductor scales 0.4.0 2016-02-26 CRAN (R 3.3.0) stringi 1.1.1 2016-05-27 CRAN (R 3.3.0) stringr 1.0.0 2015-04-30 CRAN (R 3.3.1) SummarizedExperiment * 1.3.81 2016-08-09 cran (@1.3.81) survival 2.39-5 2016-06-26 CRAN (R 3.3.1) VariantAnnotation 1.19.8 2016-07-22 Bioconductor withr 1.0.2 2016-06-20 CRAN (R 3.3.1) XML 3.98-1.4 2016-03-01 CRAN (R 3.3.0) xtable 1.8-2 2016-02-05 CRAN (R 3.3.0) XVector 0.13.7 2016-08-09 cran (@0.13.7) zlibbioc 1.19.0 2016-05-08 Bioconductor > Note that in my local test, I installed rtracklayer from the Bioconductor-mirror. In my Travis test, rtracklayer gets installed using biocLite(). I know that the latest Windows build failed, but that's because of the recurrent zlibbioc issue. I've seen it build and pass R CMD check before in moscato1. I let this issue linger because we wanted to post our pre-print to bioRxiv before @bioconductor automatically tweeted about the recount package. However, now that the pre-print is up, we'd like to have the recount page http://bioconductor.org/packages/recount/ built (which I know depends on the Linux build working). Best, Leo _______________________________________________ Bioc-devel@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/bioc-devel