Question: GetBM gives controversial results
0
gravatar for lizaveta
6 months ago by
lizaveta0
lizaveta0 wrote:

I want to get chromosome location of SMARCA1 gene.

library("biomaRt")
human = useMart(biomart="ENSEMBL_MART_ENSEMBL", dataset = "hsapiens_gene_ensembl",host="grch37.ensembl.org", path="/biomart/martservice")
getBM(attributes = c("hgnc_symbol","chromosome_name", "start_position", 
                     "end_position","strand"),
     filters = "hgnc_symbol", values = "SMARCA1", mart = human)

I have a table after running this code:

hgnc_symbol chromosome_name start_position end_position strand
1     SMARCA1               X      128580480    128657477     -1

Then, I search for this very gene in the region it should be located:

getBM(attributes = c("hgnc_symbol","chromosome_name", "start_position", "end_position","strand", "transcript_length","ensembl_transcript_id"),
filters = c("chromosome_name","start","end"),
values = list("chrX", 128580480, 128657477),mart = human)

and get nothing:

[1] hgnc_symbol           chromosome_name       start_position        end_position          strand                transcript_length     ensembl_transcript_id
<0 rows> (or 0-length row.names)

How is this possible?

biomart bioconductor R gene • 300 views
ADD COMMENTlink written 6 months ago by lizaveta0

If an answer was helpful you should upvote it, if the answer resolved your question you should mark it as accepted. Upvote|Bookmark|Accept

ADD REPLYlink written 6 months ago by Kevin Blighe39k
7
gravatar for russhh
6 months ago by
russhh4.2k
UK, U. Glasgow
russhh4.2k wrote:

set chromosome_name to X not chrX in your getBM

ADD COMMENTlink written 6 months ago by russhh4.2k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 774 users visited in the last hour