| Title: | A plyranges and tximeta workflow |
|---|---|
| Description: | An extended workflow using the plyranges and tximeta packages for fluent genomic data analysis. Use tximeta to correctly import RNA-seq transcript quantifications and summarize them to gene counts for downstream analysis. Use plyranges for clearly expressing operations over genomic coordinates and to combine results from differential expression and differential accessibility analyses. |
| Authors: | Stuart Lee [aut, cre] (ORCID: <https://orcid.org/0000-0003-1179-8436>), Michael Love [aut, ctb] |
| Maintainer: | Stuart Lee <[email protected]> |
| License: | MIT + file LICENSE |
| Version: | 1.24.0 |
| Built: | 2026-10-04 17:15:25 UTC |
| Source: | https://github.com/bioc/fluentGenomics |
Download, retrieve and cache macrophage ATAC seq results
cache_atac_se(verbose = FALSE)cache_atac_se(verbose = FALSE)
verbose |
Provide progress bar for downloading and updating the cache? (Default is FALSE) |
This code generates the SummarizedExperiment object used for ATAC seq peaks data analysis as discussed in the vignette.
the path to the .rds file containing the SummarizedExperiment
https://zenodo.org/record/1188300/files/
if (interactive()) { path_to_se <- cache_atac_se() atac <- readRDS(path_to_se) }if (interactive()) { path_to_se <- cache_atac_se() atac <- readRDS(path_to_se) }
A GRanges object containing approximately 300,000 differential accessiblity peaks estimated via limma. A description of how this object was generated can be found in the vignette.
peakspeaks
An object of class GRanges of length 296220.
#' @format A GRanges with 296220 rows and 3 metadata columns:
the ATAC peak identifer
the estimated log2 fold change
the estimated differential accessibility FDR adjusted p-value