Entering edit mode
Hello ChIPseeker, I'm wondering what the argument weightcol='V5' is doing in the getTagMatrix() function? Is it doing a weighting algorithm for the values of column 'V5'? Thanks! Best, YJ
region <- getBioRegion(TxDb=txdb, upstream=rel(0.2), downstream=rel(0.2), by="gene", type="body")
tagMatrix <- getTagMatrix(peak=peak, upstream=rel(0.2), downstream= rel(0.2), windows=region, type="body", by="gene", TxDb=txdb, weightCol = "V5", nbin = 800, ignore_strand = FALSE)
plotPeakProf(tagMatrix=tagMatrix, conf=0.95, xlab = "Genomic Region (5'->3')", ylab = "Peak Count Frequency", facet = "none", free_y = TRUE)