average_fragment_length: Define the average fragment length

In federicogiorgi/vulcan: VirtUaL ChIP-Seq data Analysis using Networks

Description

A function to get average fragment length from a ChIP-Seq experiment

Usage

average_fragment_length(bam.files, plot = TRUE, max.dist = 550)

Arguments

bam.files	a vector of BAM files locations
plot	logical. Should a plot be generated?
max.dist	numeric. Maximum fragment length accepted. Default=550

Value

nothing

Examples

library(vulcandata)
sheetfile<-'deleteme.csv'
vulcandata::vulcansheet(sheetfile)
a<-read.csv(sheetfile,as.is=TRUE)
bams<-a$bamReads
unlink(sheetfile)
average_fragment_length(bams,plot=TRUE)

federicogiorgi/vulcan documentation built on May 27, 2019, 7:25 a.m.