Question: 2 peak files in plotAvgProf2 [ChIPseeker]
gravatar for anu014
8 months 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 months ago by zx87547.5k • written 8 months 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 months ago • written 8 months ago by jared.andrews072.5k

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

ADD REPLYlink modified 7 months ago • written 7 months ago by anu014160
gravatar for ZZzzzzhong
8 months ago by
ZZzzzzhong200 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 months ago by ZZzzzzhong200

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 7 months 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: 1439 users visited in the last hour