MatchGenes in Bumphunter - result list with multiple genes
Entering edit mode
Last seen 3.3 years ago

Is it possible to get a result list with multiple genes per region using matchGenes in bumphunter? If yes, how? Is it also possible to output the start, end and chr in my input file. When I run the code below, I get only one gene per region even though UCSC genome browser shows multiple genes.

My code in R: BiocManager::install("bumphunter") library(bumphunter)

BiocManager::install("TxDb.Hsapiens.UCSC.hg38.knownGene") library(TxDb.Hsapiens.UCSC.hg38.knownGene)

CNACCCMC<- read.csv("FindGeneSymbolTest.csv") CNACCCMC<-makeGRangesFromDataFrame(CNACCC_MC, keep.extra.columns=TRUE, start.field="start", end.field="end", seqnames.field="chr")

genome <- annotateTranscripts(TxDb.Hsapiens.UCSC.hg38.knownGene)

MatchGenesCCCMCTest <- matchGenes(CNACCCMC, genome, type = "any") write.csv(MatchGenesCCCMCTest, "MatchGenesCCCMC_Test.csv", row.names=FALSE)

Thank you in advance.

Kind regards, Hanna Engqvist

cancer • 348 views

Login before adding your answer.

Traffic: 473 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6