Hi Everyone!
I have a problem with using rDGIdb looking for possible drug candidates. I want to use the option of VCF input and I am following the code described in https://bioconductor.org/packages/release/bioc/vignettes/rDGIdb/inst/doc/vignette.pdf. However it always gives me an error in the end:
library(VariantAnnotation)
library(TxDb.Hsapiens.UCSC.hg19.knownGene)
library(org.Hs.eg.db)
txdb <- TxDb.Hsapiens.UCSC.hg19.knownGene
vcf <- readVcf("somatic-wgs.vcf.gz", "hg19")
seqlevels(vcf) <- paste("chr", seqlevels(vcf), sep = "")
rd <- rowRanges(vcf)
loc <- locateVariants(rd, txdb, CodingVariants())
And the error is: Error in compatibleSeqnames(rep(seqnames(x), elementNROWS(y)), seqnames(y@unlistData)) : Level set of 'x' must be subset of that of 'y', or vice versa
I checked the seqlevels of my vcf file and TxDB as well:
seqlevels(vcf))
[1] "chr1" "chr2" "chr3" "chr4" "chr5"
[6] "chr6" "chr7" "chr8" "chr9" "chr10"
[11] "chr11" "chr12" "chr13" "chr14" "chr15"
[16] "chr16" "chr17" "chr18" "chr19" "chr20"
[21] "chr21" "chr22" "chrX" "chrY" "chrMT"
[26] "chrGL000207.1" "chrGL000226.1" "chrGL000229.1" "chrGL000231.1" "chrGL000210.1"
[31] "chrGL000239.1" "chrGL000235.1" "chrGL000201.1" "chrGL000247.1" "chrGL000245.1"
[36] "chrGL000197.1" "chrGL000203.1" "chrGL000246.1" "chrGL000249.1" "chrGL000196.1"
[41] "chrGL000248.1" "chrGL000244.1" "chrGL000238.1" "chrGL000202.1" "chrGL000234.1"
[46] "chrGL000232.1" "chrGL000206.1" "chrGL000240.1" "chrGL000236.1" "chrGL000241.1"
[51] "chrGL000243.1" "chrGL000242.1" "chrGL000230.1" "chrGL000237.1" "chrGL000233.1"
[56] "chrGL000204.1" "chrGL000198.1" "chrGL000208.1" "chrGL000191.1" "chrGL000227.1"
[61] "chrGL000228.1" "chrGL000214.1" "chrGL000221.1" "chrGL000209.1" "chrGL000218.1"
[66] "chrGL000220.1" "chrGL000213.1" "chrGL000211.1" "chrGL000199.1" "chrGL000217.1"
[71] "chrGL000216.1" "chrGL000215.1" "chrGL000205.1" "chrGL000219.1" "chrGL000224.1"
[76] "chrGL000223.1" "chrGL000195.1" "chrGL000212.1" "chrGL000222.1" "chrGL000200.1"
[81] "chrGL000193.1" "chrGL000194.1" "chrGL000225.1" "chrGL000192.1"
(seqlevels(txdb))
[1] "chr1" "chr2" "chr3"
[4] "chr4" "chr5" "chr6"
[7] "chr7" "chr8" "chr9"
[10] "chr10" "chr11" "chr12"
[13] "chr13" "chr14" "chr15"
[16] "chr16" "chr17" "chr18"
[19] "chr19" "chr20" "chr21"
[22] "chr22" "chrX" "chrY"
[25] "chrM" "chr1_gl000191_random" "chr1_gl000192_random"
[28] "chr4_ctg9_hap1" "chr4_gl000193_random" "chr4_gl000194_random"
[31] "chr6_apd_hap1" "chr6_cox_hap2" "chr6_dbb_hap3"
[34] "chr6_mann_hap4" "chr6_mcf_hap5" "chr6_qbl_hap6"
[37] "chr6_ssto_hap7" "chr7_gl000195_random" "chr8_gl000196_random"
[40] "chr8_gl000197_random" "chr9_gl000198_random" "chr9_gl000199_random"
[43] "chr9_gl000200_random" "chr9_gl000201_random" "chr11_gl000202_random"
[46] "chr17_ctg5_hap1" "chr17_gl000203_random" "chr17_gl000204_random"
[49] "chr17_gl000205_random" "chr17_gl000206_random" "chr18_gl000207_random"
[52] "chr19_gl000208_random" "chr19_gl000209_random" "chr21_gl000210_random"
[55] "chrUn_gl000211" "chrUn_gl000212" "chrUn_gl000213"
[58] "chrUn_gl000214" "chrUn_gl000215" "chrUn_gl000216"
[61] "chrUn_gl000217" "chrUn_gl000218" "chrUn_gl000219"
[64] "chrUn_gl000220" "chrUn_gl000221" "chrUn_gl000222"
[67] "chrUn_gl000223" "chrUn_gl000224" "chrUn_gl000225"
[70] "chrUn_gl000226" "chrUn_gl000227" "chrUn_gl000228"
[73] "chrUn_gl000229" "chrUn_gl000230" "chrUn_gl000231"
[76] "chrUn_gl000232" "chrUn_gl000233" "chrUn_gl000234"
[79] "chrUn_gl000235" "chrUn_gl000236" "chrUn_gl000237"
[82] "chrUn_gl000238" "chrUn_gl000239" "chrUn_gl000240"
[85] "chrUn_gl000241" "chrUn_gl000242" "chrUn_gl000243"
[88] "chrUn_gl000244" "chrUn_gl000245" "chrUn_gl000246"
[91] "chrUn_gl000247" "chrUn_gl000248" "chrUn_gl000249"
Any suggestions? Thank you!