biomaRt code giving me trouble

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 for me. I hope you could spot any mistakes in my code:

go=c("GO:0051330", "GO:0000080", chrom=c(1,2,"Y"))

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

Error in getBM(attributes = "hgnc_symbol", filters = c("go", "chromosome_name"),  : 
  object 'chrom' not found

getBM(attributes="hgnc_symbol"
+ filters=c("go","chromosome_name")

Error: unexpected symbol in:
"getBM(attributes="hgnc_symbol"
filters"

values=list(go,chrom), mart=ensembl)

Error: unexpected ',' in "values=list(go,chrom),"

• 3,438 views • link

updated 4.6 years ago by Ram

• written 11.9 years ago by chrisclarkson100

The command head(gene$hgnc) returns NULL? What do I have to do to give it a value?

• 0 views • link

updated 4.6 years ago by Ram

• written 11.9 years ago by chrisclarkson100

How did you create it in the first place (or did you)?

• 0 views • link

written 11.9 years ago by Devon Ryan

For what it's worth, section 4.3 of that vignette doesn't seem to work because the example GO IDs in it no longer match to the genes. Presumably this is due to ensembl now using GRCh38 when the vignette was made for hg19 (this is just a guess). If you change the GO ids to GO:0000086 and GO:0051322 and change the "go" filter to "go_id" then you'll get at least partly similar results. If I get a chance, I'll look at the devel branch vignette and see if it's broken too.

• 0 views • link

written 11.9 years ago by Devon Ryan

I've now alerted the package maintainer. The first error that you ran into has already been corrected for the next release of the package, but the GO ids and the filter are still wrong. Hopefully Steffen will have a chance to update things before the next release goes live, which shouldn't be long.

• 0 views • link

written 11.9 years ago by Devon Ryan

1 answer

It looks like you just merged two lines. Instead of this:

go=c("GO:0051330", "GO:0000080", chrom=c(1,2,"Y"))

Try this:

go=c("GO:0051330", "GO:0000080")
chrom=c("1","2","Y")

• 0 views • link

written 11.9 years ago by Devon Ryan

perfect thanks. how do i get the getBM command to work?

• 0 views • link

written 11.9 years ago by chrisclarkson100

By reading help(getBM) and showing us what gene$hgnc contains in the other thread you just started.

• 0 views • link

written 11.9 years ago by Devon Ryan