output from getBM() function in biomaRt package
0
0
Entering edit mode
k.askland ▴ 20
@kaskland-6850
Last seen 9.1 years ago
United States

I have been able to use the method suggested in response to my last post (biomaRt SNP (rs ID) -to-HGNC symbol mapping question) as an alternative means to obtaining rsID-HGNC symbol maps.  However, I would ideally like to compare several methods for obtaining gene mapping information for a large set of SNPs. One method I'd like to use is the getBM function from the biomaRt package by accessing the ensembl mart. However, while I can input a list of SNP rsIDs with their chromosomal locations (i.e., chromosome, start and end BP positions) with the script below, there does not appear to be a way to retain that information in the output file. The attributes list available for biomaRt does not appear to include rsID or marker location. I would like to know if there's a way to add the marker information to the set of attributes?

ChromLocMap <- getBM(c('hgnc_symbol','chromosome_name','start_position','end_position', 'band'), 
        filters = c('chromosome_name', 'start', 'end'),
        values=as.list(rsIDlist[,c(as.character("chr"),"startBP","endBP")]), mart=ensembl)

biomaRt • 1.3k views
ADD COMMENT

Login before adding your answer.

Traffic: 574 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