Entering edit mode
Hi,
I'm trying to run bioconductor packages with Cricetulus griseus(Chinese hamster) genome annotation. As there is no compiled annotation for Cricetulus griseus, I'm trying to create a custom annotation.
But there's an error that I don't know how to fix it.
I attached codes and the results below:
## Load libraries
library(ChIPQC)
library(GenomicRanges)
## load CriGri-PICR annotation
library(GenomicFeatures)
TxDb.Crigri.PICR <- makeTxDbFromGFF(file = "Cricetulus_griseus_picr.CriGri-PICR.104.gtf.gz", format = "gtf")
## Load sample data
samples <- read.csv('meta/CHO_chipseq.csv')
View(samples)
## Create ChIPQC object
chipObj <- ChIPQC(samples, annotation=TxDb.Crigri.PICR)
The output error messages were like this:
> ## Load sample data
> samples <- read.csv('meta/CHO_chipseq.csv')
> View(samples)
> ## Create ChIPQC object
> chipObj <- ChIPQC(samples, annotation=TxDb.Crigri.PICR)
Tp5-H3K4me1 H3K4me1 1 macs
Tp5-H3K4me3 H3K4me3 1 narrow
Tp5-H3K27ac H3K27ac 1 narrow
Checking chromosomes:
[1] "RAZU01000001.1"
Compiling annotation...
Error in GeneAnnotation == "hg19" :
comparison (1) is possible only for atomic and list types
How can I fix this error?
Thanks ```
Hi Respected Team,
I am also encoutering the same error given below:
Bam file has 24 contigs Calculating coverage histogram for chr1
Calculating SSD for chr1
Calculating unique positions per strand for chr1
Calculating shift for chr1
300 / 300
Error in names(res) <-
*vtmp*
: 'names' attribute [10] must be the same length as the vector [2] In addition: Warning message: stop worker failed: wrong args for environment sub assignmentthis is my command in R for chipqc:
Below is information of environment initialize
Please help me to resolve this issue.
Best regards Gopal