[BioC] problem with org.Hs.egCHRLOCEND
James W. MacDonald
jmacdon at med.umich.edu
Thu Jun 11 22:43:10 CEST 2009
Hi Kay,
Kay Jaja wrote:
> Does any one know where i can find org.Hs.egCHRLOCEND
> to find the end base pair position for a gene. I have downloaded
> org.Hs.eg.db package and org.Hs.egCHRLOCEND is not there.
That is correct - the gene ends are not part of that package. You can
use biomaRt for this.
> library(biomaRt)
> mart <- useMart("ensembl","hsapiens_gene_ensembl")
Checking attributes ... ok
Checking filters ... ok
> gns <- c("CRIPAK","CAND2","STK25")
> getBM(c("hgnc_symbol","chromosome_name",
"start_position","end_position"), "hgnc_symbol", gns, mart)
hgnc_symbol chromosome_name start_position end_position
1 CRIPAK 4 1378152 1379776
2 CAND2 3 12813171 12851301
3 STK25 2 242083105 242096776
Alternatively, you could use getGene(), which has a simpler interface
but returns things you might not want.
Best,
Jim
> I have tried the following example
>
> Library(org.Hs.eg.db)
>
>> glist
> [1] "CRIPAK" "CAND2" "STK25"
>
>> mget(eglist, org.Hs.egCHR)
> $`285464`
> [1] "4"
>
> $`23066`
> [1] "3"
>
> $`10494`
> [1] "2"
>
> Find the start position:
> > mget(eglist, org.Hs.egCHRLOC)
> $`285464`
> 4
> 1375339
>
> $`23066`
> 3
> 12813170
>
> $`10494`
> 2
> -242083104
>
> And the end positions I get an error
>
>
>> mget(eglist, org.Hs.egCHRLOCEND)
> Error: object "org.Hs.egCHRLOCEND" not found
> Error in mget(eglist, org.Hs.egCHRLOCEND) :
> error in evaluating the argument 'envir' in selecting a method for function 'mget'
>
>
> Thanks for your help
>
>
>
> [[alternative HTML version deleted]]
>
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
--
James W. MacDonald, M.S.
Biostatistician
Douglas Lab
University of Michigan
Department of Human Genetics
5912 Buhl
1241 E. Catherine St.
Ann Arbor MI 48109-5618
734-615-7826
More information about the Bioconductor
mailing list