Question: 2 peak files in plotAvgProf2 [ChIPseeker]
0
gravatar for anu014
5 days ago by
anu014150
India
anu014150 wrote:

Hello Biostars!,

Can someone tell me how to plot 2 peak files in a single plot using plotAvgProf2. I was using below command but it gave me just one sample:

plotAvgProf2(c(peak2,peak3), TxDb=txdb, upstream=3000, downstream=3000,weightCol="fold_enrichment", xlab="Genomic Region (5'->3')", ylab = "Read Count Frequency"

Please help me out..

Thanks :)

ADD COMMENTlink modified 4 days ago by zx87545.3k • written 5 days ago by anu014150

I think the peak files need to be passed in as a list.

Also, be sure to recognize that all of ChIPSeeker's average peak profiles just show the number/frequency of peaks, not the number of reads. The vignette is slightly confusing in that regarding since all of the y-axis labels are "Read Count Frequency", though that's not really what's being shown.

ADD REPLYlink modified 5 days ago • written 5 days ago by jared.andrews071.3k
2
gravatar for ZZzzzzhong
4 days ago by
ZZzzzzhong90
ZZzzzzhong90 wrote:
files <- list(peak2,peak3)
promoter <- getPromoters(TxDb=txdb, upstream=3000, downstream=3000)
tagMatrixList <- lapply(files, getTagMatrix, windows=promoter)
plotAvgProf(tagMatrixList, xlim=c(-3000, 3000))
ADD COMMENTlink written 4 days ago by ZZzzzzhong90
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1888 users visited in the last hour