diff --git a/vignettes/recount-workflow.Rmd b/vignettes/recount-workflow.Rmd index fe16c87..5602f75 100644 --- a/vignettes/recount-workflow.Rmd +++ b/vignettes/recount-workflow.Rmd @@ -785,9 +785,6 @@ rse_er_scaled <- scale_counts(rse_er) ## To highlight that we scaled the counts rm(rse_er) - -## Remove files we no longer need -unlink("SRP045638/bw", recursive = TRUE) ``` Now that we have a scaled count matrix for the expressed regions, we can proceed with the DE analysis just like we did at the gene and exon feature levels (Figures \@ref(fig:erdeanalysis1), \@ref(fig:erdeanalysis2), \@ref(fig:erdeanalysis3), and \@ref(fig:erdeanalysis4)). @@ -861,6 +858,10 @@ bws <- paste0( colData(rse_er_scaled)$bigwig_file ) +## Workaround to https://github.com/lawremi/rtracklayer/issues/83: use the local +## files we already downloaded +bws <- gsub("http://duffel.rail.bio/recount/", "", bws) + ## Note that they are also present in the recount_url data.frame bws <- recount_url$url[match( colData(rse_er_scaled)$bigwig_file, @@ -1000,6 +1001,7 @@ session_info() ```{r clean_up, echo = FALSE} ## Delete big files unlink(dir("SRP045638", "rse_", full.names = TRUE)) +unlink("SRP045638/bw", recursive = TRUE) ```