trying to retrieve data on go term from ensembl through biomaRt

Trying to retrieve data on genes pertaining to the GO term mentioned.... it won't work...

ensembl = useMart("ensembl",dataset="hsapiens_gene_ensembl")
go=c("GO:0006306")
getBM(attributes="hgnc_symbol", "hgnc_id"
       filters=c("go", "chromosome_name"),
       values=list(go, chrom), mart=ensembl)

Returns:

getBM(attributes="hgnc_symbol", "hgnc_id"
+        filters=c("go", "chromosome_name"),
Error: unexpected symbol in:
"getBM(attributes="hgnc_symbol", "hgnc_id"
       filters"
>        values=list(go, chrom), mart=ensembl)
Error: unexpected ',' in "       values=list(go, chrom),"
> getBM(attributes="hgnc_symbol", "hgnc_id"
+        filters=c("go", "chromosome_name"),
Error: unexpected symbol in:
"getBM(attributes="hgnc_symbol", "hgnc_id"
       filters"
>        values=list(go, chrom))

any tips?

Thanks

• 3,164 views • link

updated 4.6 years ago by Ram

• written 11.9 years ago by chrisclarkson100

1 answer

You forgot an enclosure.

getBM(attributes="hgnc_symbol", "hgnc_id"

should be

getBM(attributes=c("hgnc_symbol", "hgnc_id"),

• 0 views • link

written 11.9 years ago by Devon Ryan

Also,

filters=c("go", "chromosome_name"),

should be:

filters=c("go_id", "chromosome_name"),

• 0 views • link

updated 4.6 years ago by Ram

• written 11.9 years ago by Devon Ryan

Hey Devon Ryan,

So my new code is:

getBM(attributes="hgnc_symbol", "hgnc_id")
       filters=c("go", "chromosome_name"),
       values=list(go, chrom), mart=ensembl)

and it returns:

getBM(attributes="hgnc_symbol", "hgnc_id")
Error in martCheck(mart) : 
  You must provide a valid Mart object. To create a Mart object use the function: useMart.  Check ?useMart for more information.
>        filters=c("go", "chromosome_name"),
Error: unexpected ',' in "       filters=c("go", "chromosome_name"),"
>        values=list(go, chrom),
Error: unexpected ',' in "       values=list(go, chrom),"
>        mart=ensembl)
Error: unexpected ')' in "       mart=ensembl)"

• 0 views • link

updated 4.6 years ago by Ram

• written 11.9 years ago by chrisclarkson100

You forgot a comma on the first line :)

• 0 views • link

updated 4.6 years ago by Ram

• written 11.9 years ago by Devon Ryan

More posts like this

Error: unexpected symbol using BiomaRt in R
written 5.1 years ago by andreadcass •

I'm trying to use biomaRt in R to make a table containing HGNC symbols and UnitProt/SwissProt IDs. Here is my code: ensembl = useMart("ENSEMBL_MART_ENSEMBL", dataset="hsapiens_gene_ensembl") …
gene ensemble ID- no symbol
written 5.4 years ago by Rob •

Hello friends, I used this code to convert ensemble gene ID to gene symbol: df <- read.csv("gene_Ens_ID.csv") library('biomaRt') mart <- useDataset("hsapiens_gene_ensembl", useMart("ensembl")) attributeNames <- c('ensembl_gene_id','external_gene_name','description','gene_biotype', …
Finding subcellular localization of genes using biomaRt - filter not found
written 6.1 years ago by andyincolorado •

I'm trying to find the subcellular localizations of a list of genes using biomaRt but I keep getting the error: `Query ERROR: caught BioMart::Exception::Usage: Filter …
annotation - biomaRt - getBM - multiple entrez ID
written 7.6 years ago by Learner

Hi I am trying to follow what is said in here , however, without any success https://support.bioconductor.org/p/52407/ when I do it as described here, I …
Converting Ensembl gene id/gene id version to hgnc symbol using Biomart r package
written 7.9 years ago by sahar850 •

Hi, I need to convert data from TCGA in the form of ensembl gene id version to hgnc symbol using Biomat r package. After creating …
GetBM gives controversial results
written 7.9 years ago by lizaveta •

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", …
Biomart error while convertint transcript ID.
written 9.3 years ago by always_learning

library("biomaRt") ensembl<- useMart("ensembl", dataset="hsapiens_gene_ensembl") value =read.csv("Refseq.txt", strip.white=TRUE) colnames(value) = c("name") values <- as.vector(value$name) getBM(attributes=c("refseq_mrna", "ensembl_gene_id", "hgnc_symbol"), filters = "refseq_mrna", values = values, mart= ensembl) and …
Empty gene names from ensembl query using biomart
written 11.7 years ago by kindlychung •

Here is the code: ``` require(biomaRt) ensembl = useMart("ensembl", ,dataset="hsapiens_gene_ensembl") # filters = listFilters(ensembl) # filters[1:5, ] # attributes = listAttributes(ensembl) # ensembl_attr = cbind(attributes$name, …
struggling with an R command
written 11.9 years ago by chrisclarkson100

I have established ensembl as my database but I can not get this code to work.... Any help? ``` results <- getBM(attributes = c("ensembl_gene_id", "hgnc_symbol"), …
biomaRt code giving me trouble
written 11.9 years ago by chrisclarkson100

Trying to get task number 4.3 from biomaRt vignette to work for me. I have established ensembl as the database but from there nothing works …