This is a test version of Biostars. For the public version, visit https://www.biostars.org.
CHIPseeker error now
tagMatrix_1 <- getTagMatrix(peak_1, weightCol =NULL, windows = promoter_1)

Error in split.default(1:length(windows), as.factor(seqnames(windows))) : 
  group length is 0 but data length > 0
In addition: Warning message:
In .Seqinfo.mergexy(x, y) :
  The 2 combined objects have no sequence levels in common. (Use
  suppressWarnings() to suppress this warning.)
r
• 3,096 views • link
updated by zx8754 zx8754 • written by lzhizhao9179 • lzhizhao9179

It is a good idea to Google your error/warning messages, if you did the search then you would see that there is an issue at GitHub with a solution.

Seems like the problem was the chromosome column being labelled (1 instead of chr1)

Fixed with:

newStyle <- mapSeqlevels(seqlevels(peak), "UCSC")
peakchr <- renameSeqlevels(peak, newStyle)
• 1 views • link
written by zx8754 zx8754

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • summarizeOverlaps warning A. thaliana differential expression
    written by eva.lucarelli • eva.lucarelli
    Hello, I am currently working on a differential expression tutorial regarding gene expression in Arabidopsis thaliana. At some point in the exercise, we use the …
  • CNVRs visualization with karyoploteR
    written by rthapa • rthapa
    Hi, I am using karyoploteR to visualize CNV. I am getting following error. But I do have the genomic region in a bed file. Could …
  • Error in Gviz BiomartGeneRegionTrack
    written by codezy • codezy
    Dear friends, I would like to run the "BiomartGeneRegionTrack" in package Gviz to get the annotation file. But it seems I can not use the …
  • Error while annotating with Chipseeker
    written by hamaor • hamaor
    I followd the Diffbind pipeline for a ChipSeq data: peaksets2 = dba(sampleSheet="SampleSheet.csv",config = data.frame(RunParallel=TRUE, AnalysisMethod=DBA_DESEQ2, bCorPlot=FALSE, bUsePval=TRUE, minQCth=10),minOverlap = 2) peaksets2 <- dba.count(peaksets2) peaksets2 <- …
  • Error using Tabix
    written by AS • AS
    Hi, I am trying to use vcf files downloaded from public database (NCBI SRA Bioproject:PRJNA486011). I am using only R to do that because I …
  • Something is wrong with getTagMatrix (A command in ChIPseeker)?
    written by 1067648804 • 1067648804
    I am using ChIPseeker for data nalysis. When I am running the getTagMatrix command, errors occured as follows: > promoter <- getPromoters(TxDb=GhTM_1, upstream=3000, downstream=3000) > …
  • RNA Seq Workflow
    written by jadepinket • jadepinket
    I am using the RNA-seq workflow: gene-level exploratory analysis and differential expression. I have summerizeoverlap output 0 so <- summarizeOverlaps(ebg, "Epigenomics02-96.alignments.bam", mode="Union", + singleEnd=FALSE, + …
  • Annotate DMR information
    written by yassimaryam • yassimaryam
    We have used “GRCh37/hg19” for reference genome in our project. We convert output file to GenomicRange object that has been 3 columns included that: (Seqnames, …
  • intansv - problem at the annotation step
    written by rajini28m • rajini28m
    Hi Everyone, I am doing structural variant detection using different programs such as DELLY, Softsearch, Svseq. After detection I am integrating the SV calls from …
  • problem with methylKit of the annotation
    written by hua.peng1314 hua.peng1314
    I am dealing with RRBS data follow [the manual][1]. I met some problem when doing annotation.The problem showed below. ``` > read.transcript.features("mm9.refseq.genes.bed",remove.unsual=TRUE,up.flank=1500,down.flank=500,unique.prom=TRUE) > diffAnn <- …