Introduction to dasper
1 November 2022
Package
dasper 1.9.0
1 Basics
1.1 Install dasper
R is an open-source statistical environment which can be easily modified to enhance its functionality via packages. dasper is a R package available via the Bioconductor repository for packages. R can be installed on any operating system from CRAN after which you can install dasper by using the following commands in your R session:
if (!requireNamespace("BiocManager", quietly = TRUE)) {
install.packages("BiocManager")
}
BiocManager::install("dzhang32/dasper")
## Check that you have a valid Bioconductor installation
BiocManager::valid()1.2 Required knowledge
The expected input of dasper are junctions reads (e.g. directly outputted from an aligner such as STAR or extracted from a BAM file (e.g. using megadepth) and coverage in the form of BigWig files (which can be generated from BAM files using megadepth or RSeQC). dasper is based on many other packages and in particular in those that have implemented the infrastructure needed for dealing with RNA-sequencing data. The packages SummarizedExperiment and GenomicRanges are used throughout, therefore familiarity with these packages will greatly help in interpreting the output of dasper.
If you are asking yourself the question “Where do I start using Bioconductor?” you might be interested in this blog post. Or if you find the structure of a SummarizedExperiment unclear, you might consider checking out this manual.
1.3 Asking for help
As package developers, we try to explain clearly how to use our packages and in which order to use the functions. But R and Bioconductor have a steep learning curve so it is critical to learn where to ask for help. The blog post quoted above mentions some but we would like to highlight the Bioconductor support site as the main resource for getting help: remember to use the dasper tag and check the older posts. Other alternatives are available such as creating GitHub issues and tweeting. However, please note that if you want to receive help you should adhere to the posting guidelines. It is particularly critical that you provide a small reproducible example and your session information so package developers can track down the source of the error.
1.4 Citing dasper
We hope that dasper will be useful for your research. Please use the following information to cite the package and the overall approach. Thank you!
## Citation info
citation("dasper")
#>
#> To cite package 'dasper' in publications use:
#>
#> dzhang32 (2022). _Detecting abberant splicing events from
#> RNA-sequencing data_. doi:10.18129/B9.bioc.dasper
#> <https://doi.org/10.18129/B9.bioc.dasper>,
#> https://github.com/dzhang32/dasper - R package version 1.9.0,
#> <http://www.bioconductor.org/packages/dasper>.
#>
#> dzhang32 (2020). "Detecting abberant splicing events from
#> RNA-sequencing data." _bioRxiv_. doi:10.1101/TODO
#> <https://doi.org/10.1101/TODO>,
#> <https://www.biorxiv.org/content/10.1101/TODO>.
#>
#> To see these entries in BibTeX format, use 'print(<citation>,
#> bibtex=TRUE)', 'toBibtex(.)', or set
#> 'options(citation.bibtex.max=999)'.2 Quick start guide
2.1 Workflow
The above workflow diagram gives a top-level overview of the available functions within dasper and describes the order in which they are intended to be run. This are broadly split into 4 categories:
- Process junctions functions are prefixed with a junction_. They will load in your junctions into an
RangedSummarizedExperimentobject, annotate your junctions using reference annotation, filter out junctions for those that likely originate from technical noise and normalize your junction counts to allow for comparison between samples. - Process coverage functions are prefixed with a coverage_. They annotate your junctions with coverage from associated regions then normalize this coverage to allow for comparison between samples.
- Outlier detection functions are prefixed with a outlier_. They will use a outlier detection algorithm (isolation forest) to detect the most outlier/abnormal/aberrant junctions in terms of their counts and associated coverage.
- Sashimi plots can be generated using the function
plot_sashimi. This enables you to select a gene, transcript or region of interest and plot the splicing and coverage across that region in the form of sashimi plot.
2.2 Example
dasper includes wrapper functions for the 3 major analysis steps in the workflow - processing junction data, processing coverage data, then performing outlier detection. If you are familiar with dasper or uninterested with the intermediates, you can go from start to finish using these wrappers.
2.2.1 Setup
First, we need to install the system dependency megadepth, which is required for the loading of coverage from BigWig files. The easiest way to do this is to run install_megadepth from the megadepth
megadepth::install_megadepth()
#> The latest megadepth version is 1.2.0
#> This is not an interactive session, therefore megadepth has been installed temporarily to
#> /tmp/Rtmp2Hddaj/megadepthdasper requires reference annotation that can either be inputted as a GTF path or pre-loaded into a TxDb object as shown below.
# use GenomicState to load txdb (GENCODE v31)
ref <- GenomicState::GenomicStateHub(version = "31", genome = "hg38", filetype = "TxDb")[[1]]
#> snapshotDate(): 2022-10-26
#> downloading 1 resources
#> retrieving 1 resource
#> loading from cache
#> Loading required package: GenomicFeatures
#> Loading required package: BiocGenerics
#>
#> Attaching package: 'BiocGenerics'
#> The following objects are masked from 'package:stats':
#>
#> IQR, mad, sd, var, xtabs
#> The following objects are masked from 'package:base':
#>
#> Filter, Find, Map, Position, Reduce, anyDuplicated, aperm, append,
#> as.data.frame, basename, cbind, colnames, dirname, do.call,
#> duplicated, eval, evalq, get, grep, grepl, intersect, is.unsorted,
#> lapply, mapply, match, mget, order, paste, pmax, pmax.int, pmin,
#> pmin.int, rank, rbind, rownames, sapply, setdiff, sort, table,
#> tapply, union, unique, unsplit, which.max, which.min
#> Loading required package: S4Vectors
#> Loading required package: stats4
#>
#> Attaching package: 'S4Vectors'
#> The following objects are masked from 'package:base':
#>
#> I, expand.grid, unname
#> Loading required package: IRanges
#> Loading required package: GenomeInfoDb
#> Loading required package: GenomicRanges
#> Loading required package: AnnotationDbi
#> 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")'.RNA-seq data in the form of junctions and BigWig files are the required raw input to dasper. For illustration, in this vignette we will use an BigWig file from the GTEx v6 data publicly hosted by recount2. In reality, users should use the BigWig files that correspond to the same samples as those your junction data originate from.
# Obtain the urls to the remotely hosted GTEx BigWig files
url <- recount::download_study(
project = "SRP012682",
type = "samples",
download = FALSE
)
#> Setting options('download.file.method.GEOquery'='auto')
#> Setting options('GEOquery.inmemory.gpl'=FALSE)
# cache the file using BiocFileCache for faster retrieval
bw_path <- dasper:::.file_cache(url[1])A small set of example junctions are included within the package, which can be used for testing the functionality of dasper.
library(dasper)
junctions_example
#> class: RangedSummarizedExperiment
#> dim: 97671 5
#> metadata(0):
#> assays(1): raw
#> rownames: NULL
#> rowData names(0):
#> colnames(5): count_1 count_2 gtex_54020 gtex_51495 gtex_56930
#> colData names(2): samp_id case_control2.2.2 Running dasper
Here, we run the 3 wrapper functions (junction_process, coverage_process and outlier_process) on junctions_example. These functions are designed to be compatible with usage of the %>% pipe from tidyverse. The time-limiting steps of this analysis can be run parallel using BiocParallel as shown.
library(magrittr)
outlier_scores <-
junction_process(junctions_example, ref) %>%
coverage_process(ref,
coverage_paths_case = rep(bw_path, 2),
coverage_paths_control = rep(bw_path, 3),
bp_param = BiocParallel::MulticoreParam(5)
) %>%
outlier_process(bp_param = BiocParallel::MulticoreParam(5))3 Comprehensive usage guide
In this section, we will run each of the 1-9 functions shown in workflow. This can be helpful for users that want to understand or modify the intermediates of the dasper pipeline or are only interested in a executing a specific step (e.g. annotating junctions using junction_annot). If you want to follow along and run this code in your R session, make sure you have followed the instructions in the section setup.
3.1 Junction processing
The first step is to load in patient (and control) junction data using junction_load. By default, this expects junction files in the format outputted by the STAR aligner (SJ.out). If the user’s junctions are in a different format you can accommodate for this by adjusting the load_func argument. The output of this will be junctions stored in a RangedSummarizedExperiment object required for the downstream dasper functions.
Additionally, users can choose to add a set of control junctions from GTEx v6 (publicly available through recount2) using the argument controls. This may be useful if the user only has RNA-seq data from only a single or small number of patients, which is often the case in a diagnostic setting. In this example, we arbitrarily chose to use the GTEx junctions from the fibroblast tissue. However, we would recommend matching the control tissue to the tissue from which your patient RNA-seq data was derived. The current available tissues that be selected via controls are “fibroblasts”, “lymphocytes”, “skeletal_muscle” and “whole_blood”.
junctions_example_1_path <-
system.file(
"extdata",
"junctions_example_1.txt",
package = "dasper",
mustWork = TRUE
)
junctions_example_2_path <-
system.file(
"extdata",
"junctions_example_2.txt",
package = "dasper",
mustWork = TRUE
)
# only keep chromosomes 21 + 22 for speed
junctions <-
junction_load(
junction_paths = c(
junctions_example_1_path,
junctions_example_2_path
),
controls = "fibroblasts",
chrs = c("chr21", "chr22")
)
#> [1] "2022-11-01 17:40:23.594079 - Loading junctions for sample 1/2..."
#> [1] "2022-11-01 17:40:23.651343 - Loading junctions for sample 2/2..."
#> [1] "2022-11-01 17:40:23.684362 - Adding control junctions..."
#> [1] "2022-11-01 17:40:23.690613 - Downloading and importing fibroblasts junction data..."
#> adding rname 'https://www.dropbox.com/s/kh9fvknltai89fo/GTEx_v6_junctions_fibroblast.rda?dl=1'
#> [1] "2022-11-01 17:41:28.765661 - Tidying and storing junction data as a RangedSummarizedExperiment..."
#> [1] "2022-11-01 17:41:30.780076 - done!"
junctions
#> class: RangedSummarizedExperiment
#> dim: 97671 286
#> metadata(0):
#> assays(1): raw
#> rownames: NULL
#> rowData names(0):
#> colnames(286): count_1 count_2 ... gtex_56472 gtex_57475
#> colData names(2): samp_id case_controlNext, we will annotate our junctions using junction_annot. Using reference annotation inputted via ref, this will classify junctions into the categories “annotated”, “novel_acceptor”, “novel_donor”, “novel_combo”, “novel_exon_skip”, “ambig_gene” and “unannotated”. Additionally, the genes/transcripts/exons that each junction overlaps will be stored in the SummarizedExperiment::rowData. This information is necessary for the downstream dasper functions but can also be useful independently if user’s are interested in for example, only retrieving the annotated/unannotated junctions.
# take only the first 5 samples (2 cases, 3 controls)
# to increase the speed of the vignette
junctions <- junctions[, 1:5]
junctions <- junction_annot(junctions, ref)
#> [1] "2022-11-01 17:41:30.821883 - Obtaining co-ordinates of annotated exons and junctions..."
#> [1] "2022-11-01 17:41:44.475088 - Getting junction annotation using overlapping exons..."
#> [1] "2022-11-01 17:41:46.863388 - Tidying junction annotation..."
#> [1] "2022-11-01 17:41:47.310188 - Deriving junction categories..."
#> [1] "2022-11-01 17:41:51.745606 - done!"
head(SummarizedExperiment::rowData(junctions))
#> DataFrame with 6 rows and 14 columns
#> in_ref gene_id_start
#> <logical> <CharacterList>
#> 1 FALSE
#> 2 TRUE ENSG00000277117.4
#> 3 TRUE ENSG00000277117.4
#> 4 FALSE
#> 5 FALSE
#> 6 FALSE
#> tx_name_start exon_id_start
#> <CharacterList> <IntegerList>
#> 1
#> 2 ENST00000612610.4,ENST00000620481.4,ENST00000623960.3,... 695012,695013
#> 3 ENST00000612610.4,ENST00000620481.4,ENST00000623960.3,... 695014,695015
#> 4
#> 5
#> 6
#> strand_start exon_width_start gene_id_end
#> <CharacterList> <IntegerList> <CharacterList>
#> 1
#> 2 + 291,291 ENSG00000277117.4
#> 3 + 165,165 ENSG00000277117.4
#> 4
#> 5
#> 6
#> tx_name_end exon_id_end
#> <CharacterList> <IntegerList>
#> 1
#> 2 ENST00000612610.4,ENST00000620481.4,ENST00000623960.3,... 695014,695015
#> 3 ENST00000612610.4,ENST00000620481.4,ENST00000623960.3,... 695016,695017
#> 4
#> 5
#> 6
#> strand_end exon_width_end gene_id_junction strand_junction type
#> <CharacterList> <IntegerList> <CharacterList> <CharacterList> <factor>
#> 1 * unannotated
#> 2 + 165,165 ENSG00000277117.4 + annotated
#> 3 + 36,36 ENSG00000277117.4 + annotated
#> 4 * unannotated
#> 5 * unannotated
#> 6 * unannotatedWe will then filter our junctions using junction_filter. You have the option of choosing to filter by count, width, type (annotation category) and whether junctions overlap a particular genomic region. The default settings shown below filter only by count and require a junction to have at least 5 reads in at least 1 sample to be kept in.
junctions <- junction_filter(junctions,
count_thresh = c(raw = 5),
n_samp = c(raw = 1)
)
#> [1] "2022-11-01 17:41:51.8268 - Filtering junctions..."
#> [1] "2022-11-01 17:41:51.827031 - by count..."
#> [1] "2022-11-01 17:41:52.203893 - done!"
junctions
#> class: RangedSummarizedExperiment
#> dim: 4873 5
#> metadata(0):
#> assays(1): raw
#> rownames: NULL
#> rowData names(14): in_ref gene_id_start ... strand_junction type
#> colnames(5): count_1 count_2 gtex_54020 gtex_51495 gtex_56930
#> colData names(2): samp_id case_controlWe can then normalize our raw junctions counts into a proportion-spliced-in using junction_norm. This will first cluster each junction by finding all other junctions that share an acceptor or donor site with it. Then, calculate the normalized counts by dividing the number of reads supporting each junction with the total number of reads in it’s associated cluster.
junctions <- junction_norm(junctions)
#> [1] "2022-11-01 17:41:52.271308 - Clustering junctions..."
#> [1] "2022-11-01 17:41:52.403194 - Normalising junction counts..."
#> [1] "2022-11-01 17:41:52.533041 - done!"
# save a separate object for plotting
junctions_normed <- junctions
# show the raw counts
tail(SummarizedExperiment::assays(junctions)[["raw"]])
#> count_1 count_2 gtex_54020 gtex_51495 gtex_56930
#> [4868,] 0 0 2 4 6
#> [4869,] 0 0 39 20 54
#> [4870,] 0 0 5 2 1
#> [4871,] 0 0 82 57 125
#> [4872,] 0 9 0 0 0
#> [4873,] 0 12 0 0 0
# and those after normalisation
tail(SummarizedExperiment::assays(junctions)[["norm"]])
#> count_1 count_2 gtex_54020 gtex_51495 gtex_56930
#> [4868,] 0 0.0000000 1.00000000 1.00000000 1.000000000
#> [4869,] 0 0.0000000 1.00000000 1.00000000 1.000000000
#> [4870,] 0 0.0000000 0.05747126 0.03389831 0.007936508
#> [4871,] 0 0.0000000 0.94252874 0.96610169 0.992063492
#> [4872,] 0 0.4285714 0.00000000 0.00000000 0.000000000
#> [4873,] 0 0.5714286 0.00000000 0.00000000 0.000000000Finally, we need to score each patient junction by how much it’s (normalized) counts deviate from the count distribution of the same junction in the control samples using junction_score. By default, this uses a z-score measure, however this can be modified to a user-inputted function by adjusting the argument score_func. These junction scores with be stored in the SummarizedExperiment::assays “score”.
junctions <- junction_score(junctions)
#> [1] "2022-11-01 17:41:52.569125 - Calculating the direction of change of junctions..."
#> [1] "2022-11-01 17:41:52.627648 - Generating junction abnormality score..."
#> [1] "2022-11-01 17:41:52.780416 - done!"
names(SummarizedExperiment::assays(junctions))
#> [1] "raw" "norm" "direction" "score"3.2 Coverage processing
We then move on to loading and normalising the coverage across regions associated with each junction using coverage_norm. Namely, these are the 2 exonic regions flanking each junction and the intron inbetween. Given that there are 3 regions for each junction and we need to obtain coverage for every junction from every sample, this step can be very computationally intense to run. Here, dasper allows parrallelisation across samples using BiocParallel and uses the tool megadepth developed by Chris Wilks, which is significantly faster than other tools (rtracklayer and pyBigWig) for loading coverage from BigWig files (see runtime comparison).
coverage <- coverage_norm(
junctions,
ref,
coverage_paths_case = rep(bw_path, 2),
coverage_paths_control = rep(bw_path, 3),
bp_param = BiocParallel::SerialParam()
)
#> [1] "2022-11-01 17:41:52.803602 - Obtaining exonic and intronic regions to load coverage from..."
#> [1] "2022-11-01 17:41:53.086782 - Obtaining regions to use to normalise coverage..."
#> [1] "2022-11-01 17:41:57.609373 - Loading coverage..."
#> [1] "2022-11-01 17:42:14.451366 - Normalising coverage..."
#> [1] "2022-11-01 17:42:14.452229 - done!"coverage_score can then be used to compare the coverage associated with each junction to the coverage distribution corresponding to the same region in controls. After which, junctions should have the SummarizedExperiment::assays “coverage_score”.
junctions <- coverage_score(junctions, coverage)
#> [1] "2022-11-01 17:42:14.471277 - Generating coverage abnormality score..."
#> [1] "2022-11-01 17:42:14.832069 - Obtaining regions with greatest coverage dysruption..."
#> [1] "2022-11-01 17:42:15.083153 - done!"
names(SummarizedExperiment::assays(junctions))
#> [1] "raw" "norm" "direction" "score"
#> [5] "coverage_region" "coverage_score"3.3 Outlier processing
The last step in the dasper pipeline is to use the “scores” and “coverage scores” as input into an outlier detection model. This involves using outlier_detect to score each junction in each sample by how aberrant it looks. The outputted scores will be stored in the SummarizedExperiment::assays “outlier_score”. This step can be can be parallelised using BiocParallel. The fit of the isolation forest may fluctuate depending on it’s intial random_state. For reproducibility you can avoid this by setting the parameter random_state to an fixed integer (for more details on the isolation forest model see https://scikit-learn.org/stable/modules/generated/sklearn.ensemble.IsolationForest.html).
junctions <- outlier_detect(
junctions,
bp_param = BiocParallel::SerialParam(),
random_state = 32L
)
#> [1] "2022-11-01 17:42:15.142973 - generating outlier scores for sample 1/2"
#> + '/home/biocbuild/.cache/R/basilisk/1.11.0/0/bin/conda' 'create' '--yes' '--prefix' '/home/biocbuild/.cache/R/basilisk/1.11.0/dasper/1.9.0/env_sklearn' 'python=3.8.13' '--quiet' '-c' 'conda-forge'
#> + '/home/biocbuild/.cache/R/basilisk/1.11.0/0/bin/conda' 'install' '--yes' '--prefix' '/home/biocbuild/.cache/R/basilisk/1.11.0/dasper/1.9.0/env_sklearn' 'python=3.8.13'
#> + '/home/biocbuild/.cache/R/basilisk/1.11.0/0/bin/conda' 'install' '--yes' '--prefix' '/home/biocbuild/.cache/R/basilisk/1.11.0/dasper/1.9.0/env_sklearn' '-c' 'conda-forge' 'python=3.8.13' 'scikit-learn=0.23.2' 'pandas=0.25.3'
#> [1] "2022-11-01 17:42:59.085052 - fitting outlier detection model with parameters: behaviour=deprecated, bootstrap=FALSE, contamination=auto, max_features=1, max_samples=auto, n_estimators=100, n_jobs=NULL, random_state=32, verbose=0, warm_start=FALSE"
#> [1] "2022-11-01 17:42:59.529579 - fitting outlier detection model with parameters: behaviour=deprecated, bootstrap=FALSE, contamination=auto, max_features=1, max_samples=auto, n_estimators=100, n_jobs=NULL, random_state=32, verbose=0, warm_start=FALSE"
#> [1] "2022-11-01 17:42:59.534888 - generating outlier scores for sample 2/2"
#> [1] "2022-11-01 17:42:59.877022 - fitting outlier detection model with parameters: behaviour=deprecated, bootstrap=FALSE, contamination=auto, max_features=1, max_samples=auto, n_estimators=100, n_jobs=NULL, random_state=32, verbose=0, warm_start=FALSE"
#> [1] "2022-11-01 17:43:00.321362 - fitting outlier detection model with parameters: behaviour=deprecated, bootstrap=FALSE, contamination=auto, max_features=1, max_samples=auto, n_estimators=100, n_jobs=NULL, random_state=32, verbose=0, warm_start=FALSE"
#> [1] "2022-11-01 17:43:00.352384 - done!"
names(SummarizedExperiment::assays(junctions))
#> [1] "raw" "norm" "direction" "score"
#> [5] "coverage_region" "coverage_score" "outlier_score"Finally, outlier_aggregate will aggregate the junction-level “outlier_scores” to a cluster level. The returned outlier_scores object contains a DataFrame with each row detailing a cluster in terms of how aberrant it looks. The gene_id_cluster column describes the gene linked to each cluster, derived from the gene annotation returned by junction_annot. Splicing events are ranked in the rank column with 1 referring to the most aberrant splicing event in that patient. Details of the specific junctions that look the most aberrant for each cluster can be found in the junctions column.
outlier_scores <- outlier_aggregate(junctions, samp_id_col = "samp_id", )
#> [1] "2022-11-01 17:43:00.458193 - Aggregating outlier scores to cluster level... "
#> [1] "2022-11-01 17:43:01.195034 - Annotating clusters with gene details..."
#> [1] "2022-11-01 17:43:01.250908 - done!"
head(outlier_scores)
#> DataFrame with 6 rows and 6 columns
#> samp_id cluster_index mean_outlier_score rank gene_id_cluster
#> <character> <character> <numeric> <numeric> <CharacterList>
#> 1 samp_1 4123 -0.200174 1 ENSG00000100243.21
#> 2 samp_1 1286 -0.199448 2 ENSG00000142156.14
#> 3 samp_1 3993 -0.189962 3 ENSG00000100418.8
#> 4 samp_1 1000 -0.158565 4 ENSG00000160213.7
#> 5 samp_1 3456 -0.156321 5 ENSG00000100129.18
#> 6 samp_1 3462 -0.155961 6 ENSG00000100129.18
#> junctions
#> <list>
#> 1 4123: 1:-0.219912:...,4124:-1:-0.180436:...,...
#> 2 1287:-1:-0.225970:...,1288: 1:-0.172927:...,...
#> 3 3993: 1:-0.134843:...,3994:-1:-0.245081:...,...
#> 4 1000: 1:-0.119976:...,1001:-1:-0.197153:...,...
#> 5 3456: 1:-0.121972:...,3457:-1:-0.190670:...,...
#> 6 3462: 1:-0.105657:...,3463:-1:-0.206265:...,...3.4 Sashimi plots
In order to help further inspection of the candidate genes with abberrant splicing returned by dasper, we also include functions to visualise splicing across genes/transcripts/regions of interest. For example, you could visualise the gene with the most aberrant cluster in samp_1.
# take gene_id of the cluster ranked 1
gene_id <- unlist(outlier_scores[["gene_id_cluster"]][1])
plot_sashimi(
junctions_normed,
ref,
case_id = list(samp_id = "samp_1"),
gene_tx_id = gene_id,
gene_tx_col = "gene_id",
count_label = FALSE
)
Often regions of splicing can be very complex within a gene/transcript. In which case, you may consider zooming in on a specific region of interest using the argument region. Furthermore, the thickness of the lines represents the normalized counts of each junction however in this way, it can be difficult to accurately differentiate junctions that have similar counts. Users may want to add a label representing the count of each junction using the argument count_label.
plot_sashimi(junctions_normed,
ref,
case_id = list(samp_id = "samp_1"),
gene_tx_id = gene_id,
gene_tx_col = "gene_id",
region = GenomicRanges::GRanges("chr22:42620000-42630000"),
count_label = TRUE
)
4 Reproducibility
The dasper package (dzhang32, 2022) was made possible thanks to:
- R (R Core Team, 2022)
- BiocStyle (Oleś, 2022)
- RefManageR (McLean, 2017)
- knitr (Xie, 2022)
- rmarkdown (Allaire, Xie, McPherson, Luraschi, Ushey, Atkins, Wickham, Cheng, Chang, and Iannone, 2022)
- sessioninfo (Wickham, Chang, Flight, Müller, and Hester, 2021)
- testthat (Wickham, 2011)
This package was developed using biocthis.
Code for creating the vignette
## Create the vignette
library("rmarkdown")
system.time(render("dasper.Rmd", "BiocStyle::html_document"))
## Extract the R code
library("knitr")
knit("dasper.Rmd", tangle = TRUE)Date the vignette was generated.
#> [1] "2022-11-01 17:43:06 EDT"Wallclock time spent generating the vignette.
#> Time difference of 3.357 minsR session information.
#> ─ Session info ───────────────────────────────────────────────────────────────────────────────────────────────────────
#> setting value
#> version R Under development (unstable) (2022-10-25 r83175)
#> os Ubuntu 22.04.1 LTS
#> system x86_64, linux-gnu
#> ui X11
#> language (EN)
#> collate C
#> ctype en_US.UTF-8
#> tz America/New_York
#> date 2022-11-01
#> pandoc 2.9.2.1 @ /usr/bin/ (via rmarkdown)
#>
#> ─ Packages ───────────────────────────────────────────────────────────────────────────────────────────────────────────
#> package * version date (UTC) lib source
#> abind 1.4-5 2016-07-21 [2] CRAN (R 4.3.0)
#> AnnotationDbi * 1.61.0 2022-11-01 [2] Bioconductor
#> AnnotationHub 3.7.0 2022-11-01 [2] Bioconductor
#> assertthat 0.2.1 2019-03-21 [2] CRAN (R 4.3.0)
#> backports 1.4.1 2021-12-13 [2] CRAN (R 4.3.0)
#> base64enc 0.1-3 2015-07-28 [2] CRAN (R 4.3.0)
#> basilisk 1.11.0 2022-11-01 [2] Bioconductor
#> basilisk.utils 1.11.0 2022-11-01 [2] Bioconductor
#> bibtex 0.5.0 2022-09-25 [2] CRAN (R 4.3.0)
#> Biobase * 2.59.0 2022-11-01 [2] Bioconductor
#> BiocFileCache 2.7.0 2022-11-01 [2] Bioconductor
#> BiocGenerics * 0.45.0 2022-11-01 [2] Bioconductor
#> BiocIO 1.9.0 2022-11-01 [2] Bioconductor
#> BiocManager 1.30.19 2022-10-25 [2] CRAN (R 4.3.0)
#> BiocParallel 1.33.0 2022-11-01 [2] Bioconductor
#> BiocStyle * 2.27.0 2022-11-01 [2] Bioconductor
#> BiocVersion 3.17.0 2022-11-01 [2] Bioconductor
#> biomaRt 2.55.0 2022-11-01 [2] Bioconductor
#> Biostrings 2.67.0 2022-11-01 [2] Bioconductor
#> bit 4.0.4 2020-08-04 [2] CRAN (R 4.3.0)
#> bit64 4.0.5 2020-08-30 [2] CRAN (R 4.3.0)
#> bitops 1.0-7 2021-04-24 [2] CRAN (R 4.3.0)
#> blob 1.2.3 2022-04-10 [2] CRAN (R 4.3.0)
#> bookdown 0.29 2022-09-12 [2] CRAN (R 4.3.0)
#> brio 1.1.3 2021-11-30 [2] CRAN (R 4.3.0)
#> broom 1.0.1 2022-08-29 [2] CRAN (R 4.3.0)
#> BSgenome 1.67.0 2022-11-01 [2] Bioconductor
#> bslib 0.4.0 2022-07-16 [2] CRAN (R 4.3.0)
#> bumphunter 1.41.0 2022-11-01 [2] Bioconductor
#> cachem 1.0.6 2021-08-19 [2] CRAN (R 4.3.0)
#> car 3.1-1 2022-10-19 [2] CRAN (R 4.3.0)
#> carData 3.0-5 2022-01-06 [2] CRAN (R 4.3.0)
#> checkmate 2.1.0 2022-04-21 [2] CRAN (R 4.3.0)
#> cli 3.4.1 2022-09-23 [2] CRAN (R 4.3.0)
#> cluster 2.1.4 2022-08-22 [2] CRAN (R 4.3.0)
#> cmdfun 1.0.2 2020-10-10 [2] CRAN (R 4.3.0)
#> codetools 0.2-18 2020-11-04 [2] CRAN (R 4.3.0)
#> colorspace 2.0-3 2022-02-21 [2] CRAN (R 4.3.0)
#> cowplot 1.1.1 2020-12-30 [2] CRAN (R 4.3.0)
#> crayon 1.5.2 2022-09-29 [2] CRAN (R 4.3.0)
#> curl 4.3.3 2022-10-06 [2] CRAN (R 4.3.0)
#> dasper * 1.9.0 2022-11-01 [1] Bioconductor
#> data.table 1.14.4 2022-10-17 [2] CRAN (R 4.3.0)
#> DBI 1.1.3 2022-06-18 [2] CRAN (R 4.3.0)
#> dbplyr 2.2.1 2022-06-27 [2] CRAN (R 4.3.0)
#> DelayedArray 0.25.0 2022-11-01 [2] Bioconductor
#> deldir 1.0-6 2021-10-23 [2] CRAN (R 4.3.0)
#> derfinder 1.33.0 2022-11-01 [2] Bioconductor
#> derfinderHelper 1.33.0 2022-11-01 [2] Bioconductor
#> desc 1.4.2 2022-09-08 [2] CRAN (R 4.3.0)
#> digest 0.6.30 2022-10-18 [2] CRAN (R 4.3.0)
#> dir.expiry 1.7.0 2022-11-01 [2] Bioconductor
#> doRNG 1.8.2 2020-01-27 [2] CRAN (R 4.3.0)
#> downloader 0.4 2015-07-09 [2] CRAN (R 4.3.0)
#> dplyr 1.0.10 2022-09-01 [2] CRAN (R 4.3.0)
#> ellipsis 0.3.2 2021-04-29 [2] CRAN (R 4.3.0)
#> evaluate 0.17 2022-10-07 [2] CRAN (R 4.3.0)
#> fansi 1.0.3 2022-03-24 [2] CRAN (R 4.3.0)
#> farver 2.1.1 2022-07-06 [2] CRAN (R 4.3.0)
#> fastmap 1.1.0 2021-01-25 [2] CRAN (R 4.3.0)
#> filelock 1.0.2 2018-10-05 [2] CRAN (R 4.3.0)
#> foreach 1.5.2 2022-02-02 [2] CRAN (R 4.3.0)
#> foreign 0.8-83 2022-09-28 [2] CRAN (R 4.3.0)
#> Formula 1.2-4 2020-10-16 [2] CRAN (R 4.3.0)
#> fs 1.5.2 2021-12-08 [2] CRAN (R 4.3.0)
#> generics 0.1.3 2022-07-05 [2] CRAN (R 4.3.0)
#> GenomeInfoDb * 1.35.0 2022-11-01 [2] Bioconductor
#> GenomeInfoDbData 1.2.9 2022-10-29 [2] Bioconductor
#> GenomicAlignments 1.35.0 2022-11-01 [2] Bioconductor
#> GenomicFeatures * 1.51.0 2022-11-01 [2] Bioconductor
#> GenomicFiles 1.35.0 2022-11-01 [2] Bioconductor
#> GenomicRanges * 1.51.0 2022-11-01 [2] Bioconductor
#> GenomicState 0.99.15 2021-10-13 [2] Bioconductor
#> GEOquery 2.67.0 2022-11-01 [2] Bioconductor
#> ggplot2 3.3.6 2022-05-03 [2] CRAN (R 4.3.0)
#> ggpubr 0.4.0 2020-06-27 [2] CRAN (R 4.3.0)
#> ggrepel 0.9.1 2021-01-15 [2] CRAN (R 4.3.0)
#> ggsci 2.9 2018-05-14 [2] CRAN (R 4.3.0)
#> ggsignif 0.6.4 2022-10-13 [2] CRAN (R 4.3.0)
#> glue 1.6.2 2022-02-24 [2] CRAN (R 4.3.0)
#> gridExtra 2.3 2017-09-09 [2] CRAN (R 4.3.0)
#> gtable 0.3.1 2022-09-01 [2] CRAN (R 4.3.0)
#> here 1.0.1 2020-12-13 [2] CRAN (R 4.3.0)
#> highr 0.9 2021-04-16 [2] CRAN (R 4.3.0)
#> Hmisc 4.7-1 2022-08-15 [2] CRAN (R 4.3.0)
#> hms 1.1.2 2022-08-19 [2] CRAN (R 4.3.0)
#> htmlTable 2.4.1 2022-07-07 [2] CRAN (R 4.3.0)
#> htmltools 0.5.3 2022-07-18 [2] CRAN (R 4.3.0)
#> htmlwidgets 1.5.4 2021-09-08 [2] CRAN (R 4.3.0)
#> httpuv 1.6.6 2022-09-08 [2] CRAN (R 4.3.0)
#> httr 1.4.4 2022-08-17 [2] CRAN (R 4.3.0)
#> interactiveDisplayBase 1.37.0 2022-11-01 [2] Bioconductor
#> interp 1.1-3 2022-07-13 [2] CRAN (R 4.3.0)
#> IRanges * 2.33.0 2022-11-01 [2] Bioconductor
#> iterators 1.0.14 2022-02-05 [2] CRAN (R 4.3.0)
#> jpeg 0.1-9 2021-07-24 [2] CRAN (R 4.3.0)
#> jquerylib 0.1.4 2021-04-26 [2] CRAN (R 4.3.0)
#> jsonlite 1.8.3 2022-10-21 [2] CRAN (R 4.3.0)
#> KEGGREST 1.39.0 2022-11-01 [2] Bioconductor
#> knitr 1.40 2022-08-24 [2] CRAN (R 4.3.0)
#> labeling 0.4.2 2020-10-20 [2] CRAN (R 4.3.0)
#> later 1.3.0 2021-08-18 [2] CRAN (R 4.3.0)
#> lattice 0.20-45 2021-09-22 [2] CRAN (R 4.3.0)
#> latticeExtra 0.6-30 2022-07-04 [2] CRAN (R 4.3.0)
#> lifecycle 1.0.3 2022-10-07 [2] CRAN (R 4.3.0)
#> limma 3.55.0 2022-11-01 [2] Bioconductor
#> locfit 1.5-9.6 2022-07-11 [2] CRAN (R 4.3.0)
#> lubridate 1.8.0 2021-10-07 [2] CRAN (R 4.3.0)
#> magrittr 2.0.3 2022-03-30 [2] CRAN (R 4.3.0)
#> Matrix 1.5-1 2022-09-13 [2] CRAN (R 4.3.0)
#> MatrixGenerics 1.11.0 2022-11-01 [2] Bioconductor
#> matrixStats 0.62.0 2022-04-19 [2] CRAN (R 4.3.0)
#> megadepth 1.9.0 2022-11-01 [2] Bioconductor
#> memoise 2.0.1 2021-11-26 [2] CRAN (R 4.3.0)
#> mime 0.12 2021-09-28 [2] CRAN (R 4.3.0)
#> munsell 0.5.0 2018-06-12 [2] CRAN (R 4.3.0)
#> nnet 7.3-18 2022-09-28 [2] CRAN (R 4.3.0)
#> pillar 1.8.1 2022-08-19 [2] CRAN (R 4.3.0)
#> pkgconfig 2.0.3 2019-09-22 [2] CRAN (R 4.3.0)
#> pkgload 1.3.1 2022-10-28 [2] CRAN (R 4.3.0)
#> plyr 1.8.7 2022-03-24 [2] CRAN (R 4.3.0)
#> plyranges 1.19.0 2022-11-01 [2] Bioconductor
#> png 0.1-7 2013-12-03 [2] CRAN (R 4.3.0)
#> prettyunits 1.1.1 2020-01-24 [2] CRAN (R 4.3.0)
#> progress 1.2.2 2019-05-16 [2] CRAN (R 4.3.0)
#> promises 1.2.0.1 2021-02-11 [2] CRAN (R 4.3.0)
#> purrr 0.3.5 2022-10-06 [2] CRAN (R 4.3.0)
#> qvalue 2.31.0 2022-11-01 [2] Bioconductor
#> R6 2.5.1 2021-08-19 [2] CRAN (R 4.3.0)
#> rappdirs 0.3.3 2021-01-31 [2] CRAN (R 4.3.0)
#> RColorBrewer 1.1-3 2022-04-03 [2] CRAN (R 4.3.0)
#> Rcpp 1.0.9 2022-07-08 [2] CRAN (R 4.3.0)
#> RCurl 1.98-1.9 2022-10-03 [2] CRAN (R 4.3.0)
#> readr 2.1.3 2022-10-01 [2] CRAN (R 4.3.0)
#> recount 1.25.0 2022-11-01 [2] Bioconductor
#> RefManageR * 1.4.0 2022-09-30 [2] CRAN (R 4.3.0)
#> rentrez 1.2.3 2020-11-10 [2] CRAN (R 4.3.0)
#> reshape2 1.4.4 2020-04-09 [2] CRAN (R 4.3.0)
#> restfulr 0.0.15 2022-06-16 [2] CRAN (R 4.3.0)
#> reticulate 1.26 2022-08-31 [2] CRAN (R 4.3.0)
#> rjson 0.2.21 2022-01-09 [2] CRAN (R 4.3.0)
#> rlang 1.0.6 2022-09-24 [2] CRAN (R 4.3.0)
#> rmarkdown 2.17 2022-10-07 [2] CRAN (R 4.3.0)
#> rngtools 1.5.2 2021-09-20 [2] CRAN (R 4.3.0)
#> rpart 4.1.19 2022-10-21 [2] CRAN (R 4.3.0)
#> rprojroot 2.0.3 2022-04-02 [2] CRAN (R 4.3.0)
#> Rsamtools 2.15.0 2022-11-01 [2] Bioconductor
#> RSQLite 2.2.18 2022-10-04 [2] CRAN (R 4.3.0)
#> rstatix 0.7.0 2021-02-13 [2] CRAN (R 4.3.0)
#> rstudioapi 0.14 2022-08-22 [2] CRAN (R 4.3.0)
#> rtracklayer 1.59.0 2022-11-01 [2] Bioconductor
#> S4Vectors * 0.37.0 2022-11-01 [2] Bioconductor
#> sass 0.4.2 2022-07-16 [2] CRAN (R 4.3.0)
#> scales 1.2.1 2022-08-20 [2] CRAN (R 4.3.0)
#> sessioninfo * 1.2.2 2021-12-06 [2] CRAN (R 4.3.0)
#> shiny 1.7.3 2022-10-25 [2] CRAN (R 4.3.0)
#> stringi 1.7.8 2022-07-11 [2] CRAN (R 4.3.0)
#> stringr 1.4.1 2022-08-20 [2] CRAN (R 4.3.0)
#> SummarizedExperiment 1.29.0 2022-11-01 [2] Bioconductor
#> survival 3.4-0 2022-08-09 [2] CRAN (R 4.3.0)
#> testthat 3.1.5 2022-10-08 [2] CRAN (R 4.3.0)
#> tibble 3.1.8 2022-07-22 [2] CRAN (R 4.3.0)
#> tidyr 1.2.1 2022-09-08 [2] CRAN (R 4.3.0)
#> tidyselect 1.2.0 2022-10-10 [2] CRAN (R 4.3.0)
#> tzdb 0.3.0 2022-03-28 [2] CRAN (R 4.3.0)
#> utf8 1.2.2 2021-07-24 [2] CRAN (R 4.3.0)
#> VariantAnnotation 1.45.0 2022-11-01 [2] Bioconductor
#> vctrs 0.5.0 2022-10-22 [2] CRAN (R 4.3.0)
#> vroom 1.6.0 2022-09-30 [2] CRAN (R 4.3.0)
#> waldo 0.4.0 2022-03-16 [2] CRAN (R 4.3.0)
#> withr 2.5.0 2022-03-03 [2] CRAN (R 4.3.0)
#> xfun 0.34 2022-10-18 [2] CRAN (R 4.3.0)
#> XML 3.99-0.12 2022-10-28 [2] CRAN (R 4.3.0)
#> xml2 1.3.3 2021-11-30 [2] CRAN (R 4.3.0)
#> xtable 1.8-4 2019-04-21 [2] CRAN (R 4.3.0)
#> XVector 0.39.0 2022-11-01 [2] Bioconductor
#> yaml 2.3.6 2022-10-18 [2] CRAN (R 4.3.0)
#> zlibbioc 1.45.0 2022-11-01 [2] Bioconductor
#>
#> [1] /tmp/RtmpWA0OdR/Rinst22b50923468d06
#> [2] /home/biocbuild/bbs-3.17-bioc/R/library
#>
#> ─ Python configuration ───────────────────────────────────────────────────────────────────────────────────────────────
#> python: /home/biocbuild/.cache/R/basilisk/1.11.0/dasper/1.9.0/env_sklearn/bin/python
#> libpython: /home/biocbuild/.cache/R/basilisk/1.11.0/dasper/1.9.0/env_sklearn/lib/libpython3.8.so
#> pythonhome: /home/biocbuild/.cache/R/basilisk/1.11.0/dasper/1.9.0/env_sklearn:/home/biocbuild/.cache/R/basilisk/1.11.0/dasper/1.9.0/env_sklearn
#> version: 3.8.13 | packaged by conda-forge | (default, Mar 25 2022, 06:04:18) [GCC 10.3.0]
#> numpy: /home/biocbuild/.cache/R/basilisk/1.11.0/dasper/1.9.0/env_sklearn/lib/python3.8/site-packages/numpy
#> numpy_version: 1.23.4
#>
#> NOTE: Python version was forced by use_python function
#>
#> ──────────────────────────────────────────────────────────────────────────────────────────────────────────────────────5 Bibliography
This vignette was generated using BiocStyle (Oleś, 2022) with knitr (Xie, 2022) and rmarkdown (Allaire, Xie, McPherson et al., 2022) running behind the scenes.
Citations made with RefManageR (McLean, 2017).
[1] J. Allaire, Y. Xie, J. McPherson, et al. rmarkdown: Dynamic Documents for R. R package version 2.17. 2022. URL: https://github.com/rstudio/rmarkdown.
[2] M. W. McLean. “RefManageR: Import and Manage BibTeX and BibLaTeX References in R”. In: The Journal of Open Source Software (2017). DOI: 10.21105/joss.00338.
[3] A. Oleś. BiocStyle: Standard styles for vignettes and other Bioconductor documents. R package version 2.27.0. 2022. URL: https://github.com/Bioconductor/BiocStyle.
[4] R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. Vienna, Austria, 2022. URL: https://www.R-project.org/.
[5] H. Wickham. “testthat: Get Started with Testing”. In: The R Journal 3 (2011), pp. 5–10. URL: https://journal.r-project.org/archive/2011-1/RJournal_2011-1_Wickham.pdf.
[6] H. Wickham, W. Chang, R. Flight, et al. sessioninfo: R Session Information. R package version 1.2.2. 2021. URL: https://CRAN.R-project.org/package=sessioninfo.
[7] Y. Xie. knitr: A General-Purpose Package for Dynamic Report Generation in R. R package version 1.40. 2022. URL: https://yihui.org/knitr/.
[8] dzhang32. Detecting abberant splicing events from RNA-sequencing data. https://github.com/dzhang32/dasper - R package version 1.9.0. 2022. DOI: 10.18129/B9.bioc.dasper. URL: http://www.bioconductor.org/packages/dasper.