average_fragment_length {vulcan} | R Documentation |
A function to get average fragment length from a ChIP-Seq experiment
average_fragment_length(bam.files, plot = TRUE, max.dist = 550)
bam.files |
a vector of BAM files locations |
plot |
logical. Should a plot be generated? |
max.dist |
numeric. Maximum fragment length accepted. Default=550 |
nothing
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)