Entering edit mode
Tim Smith
★
1.1k
@tim-smith-1532
Last seen 10.2 years ago
Hi All,
Sorry for the naive question! I was trying to retrieve some
coordinates (start and end positions) from biomaRt and I'm not sure if
I'm doing things rightâ¦
Problem definition: For a gene, retrieve the 5'UTR and exon
coordinates for the most common isoform of the gene.
***************
library(biomaRt)
ensembl = useMart("ensembl",dataset="hsapiens_gene_ensembl")
getAtt <- c('chromosome_name', 'start_position', 'end_position',
'strand','exon_chrom_start','exon_chrom_end',
      '5_utr_start','5_utr_end','3_utr_start','3_utr_end')
elocs <- getBM(attributes=getAtt,filters="hgnc_symbol",value="ZMYM4",m
art=ensembl)
print(elocs)
**************
However, this would give me the coordinates for all isoforms and it
would be difficult to get the coordinates for the most common isoform.
How can I identify the most common isoform?
many thanks!
[[alternative HTML version deleted]]