truncationAssignment {m6Aboost} | R Documentation |
An function for assigning the truncation read counts from bigWig to the GRanges peaks.
truncationAssignment(object, bw_positive, bw_negative, sampleName = "")
object |
A GRanges object which should contains all the single nucleotide peaks of miCLIP2 experiment. |
bw_positive |
A path to the bigWig file of truncation read counts at the positive strand that output from the preprocess in the m6Aboost pipeline. |
bw_negative |
A path to the bigWig file of truncation read counts at the negative strand that output from the preprocess in the m6Aboost pipeline. |
sampleName |
The column name that user would like to use for indicating the name of the sample. |
A GRanges object with the truncation read counts.
You Zhou
if (.Platform$OS.type != "windows") { testpath <- system.file("extdata", package = "m6Aboost") test <- readRDS(file.path(testpath, "test.rds")) truncationBw_p <- file.path(testpath, "truncation_positive.bw") truncationBw_n <- file.path(testpath, "truncation_negative.bw") test <- truncationAssignment(test, bw_positive=truncationBw_p, bw_negative=truncationBw_n, sampleName = "WT1") }