Question: 2 peak files in plotAvgProf2 [ChIPseeker]
gravatar for anu014
8 weeks ago by
anu014160 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 8 weeks ago by zx87546.0k • written 8 weeks ago by anu014160

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 8 weeks ago • written 8 weeks ago by jared.andrews071.6k

Thanks a lot for this information! ACTUALLY, I was thinking that it's read count only..

ADD REPLYlink modified 26 days ago • written 26 days ago by anu014160
gravatar for ZZzzzzhong
8 weeks ago by
ZZzzzzhong120 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 8 weeks ago by ZZzzzzhong120

Thanks a lot ! It was really helpful. Can you tell me how to change the colors of the samples in the plot? It's always giving me blue n brown color as default..

ADD REPLYlink written 26 days ago by anu014160
Please log in to add an answer.


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