MatchGenes in Bumphunter - result list with multiple genes
0
0
Entering edit mode
@hanna_engqvist-21246
Last seen 5.4 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 • 941 views
ADD COMMENT

Login before adding your answer.

Traffic: 735 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6