To get started analyzing my m6A sequencing, I used the Guitar R package to generate a feature distribution map, but I got an error report Error in stop_if_wrong_length("'seqnames'", ans_len) : 'seqnames' must have the length of the object to construct (1) or length 1
here is my code
library(BayesPeak) rm(list=ls()) options(stringsAsFactors = F) library(Guitar) package.version("Guitar")
a <- read.table("NC_summits.bed")
stBedFiles <- list("NC_summits.bed") txdb <- makeTxDbFromGFF(file = "Homo_sapiens.GRCh38.105.gtf.gz", format="gtf", dataSource="Ensembl", organism="Homo sapiens") GuitarPlot(txTxdb = txdb, stBedFiles = stBedFiles, headOrtail = TRUE, enableCI = FALSE, mapFilterTranscript = TRUE, pltTxType =c("mrna"))
Thank you for any advice,
Kind regards