inst/scripts/create_example_mbs.R
In CuppenResearch/MutationalPatterns: Comprehensive genome-wide analysis of mutational processes
library(tidyverse)
library(VariantAnnotation)
ref_genome <- "BSgenome.Hsapiens.UCSC.hg19"
library(ref_genome, character.only = TRUE)
# Get grl
grl <- readRDS("inst/states/blood_grl.rds")
# Get mbs
grl_mbs <- get_mut_type(grl, "mbs")
saveRDS(grl_mbs, "inst/states/blood_grl_mbs.rds")
# Count contexts
mbs_counts <- count_mbs_contexts(grl_mbs)
saveRDS(mbs_counts, "inst/states/blood_mbs_counts.rds")
CuppenResearch/MutationalPatterns documentation built on Nov. 23, 2022, 4:13 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
library(tidyverse)
library(VariantAnnotation)
ref_genome <- "BSgenome.Hsapiens.UCSC.hg19"
library(ref_genome, character.only = TRUE)
# Get grl
grl <- readRDS("inst/states/blood_grl.rds")
# Get mbs
grl_mbs <- get_mut_type(grl, "mbs")
saveRDS(grl_mbs, "inst/states/blood_grl_mbs.rds")
# Count contexts
mbs_counts <- count_mbs_contexts(grl_mbs)
saveRDS(mbs_counts, "inst/states/blood_mbs_counts.rds")
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.