Aggregate TSSs

Aggregate overlapping TSSs to generate positional scores.

tss_aggregate(experiment)

Arguments

experiment	TSRexploreR object.

Value

TSRexploreR object with GRanges of aggregated TSSs.

Details

'import_bams' generates a GRanges object of non-aggregated TSSs. This function aggregates overlapping TSSs into a sum total score per genomic position.

Examples

bam_file <- system.file("extdata", "S288C.bam", package="TSRexploreR")
assembly <- system.file("extdata", "S288C_Assembly.fasta", package="TSRexploreR")
samples <- data.frame(sample_name="S288C", file_1=bam_file, file_2=NA)

exp <- tsr_explorer(sample_sheet=samples, genome_assembly=assembly)
exp <- exp %>% 
  import_bams(paired=TRUE) %>%
  tss_aggregate
#> Warning: NAs introduced by coercion