Physical position to gene positions
1
0
Entering edit mode
jebreen20 • 0
@jebreen20-18614
Last seen 5.3 years ago

Dear All, 

I am a mathematician, 

I work in SNPs data, I have the markers (A, B, H, A, ..., -, B, A, H) for p individuals, with its physical positions, I want to map these makers to gene position, how I can use that with R?

 

Regards 

bioconductor molecular biology snps • 1.2k views
ADD COMMENT
0
Entering edit mode

You need to better define what you mean by 'map these markers to gene position', which could mean any number of things. For example

  1. SNPs that overlap exons
  2. SNPs that overlap any part of a gene
  3. Genes that are nearest to the SNP

So what exactly are you after?

ADD REPLY
0
Entering edit mode

Thanks for your repaly, I meant 

  1. Genes that are nearest to the SNP
ADD REPLY
0
Entering edit mode
@james-w-macdonald-5106
Last seen 16 hours ago
United States

Put the positions for your SNPs into a GRanges object, generate a GRanges for all the gene locations using, e.g., the TxDb.Hsapiens.UCSC.hg38.knownGene package) and then use nearest, or distanceToNearest, depending on what you really want. This is somewhat automated in ChipPeakAnno, but isn't particularly difficult to do by hand using the GenomicRanges and GenomicFeatures packages, for which there are vignettes to help you get started.

ADD COMMENT
0
Entering edit mode

This will give me the gen positions in cM?

and what is TxDb.Hsapiens.UCSC.hg38.knownGene ? I do not get that reall.

ADD REPLY
0
Entering edit mode

No, centimorgans aren't unambiguous distances. The distances you (presumably) want are in nucleotides which aren't ambiguous.

As for what the TxDb package I mentioned is, I already told you what vignettes you need to read. You won't get anywhere with R or Bioc if you aren't willing to read documentation and learn for yourself, and the GenomicFeatures vignette is titled 'Making and utilizing TxDb objects', so let me reiterate. Read the vignettes for those two packages.

ADD REPLY
0
Entering edit mode

many Thanks, but not this is what I want, 

maybe my question not clear, I have the genetic map and the physical positions of these markers, I want to find the gene position in cM for these markers. 

ADD REPLY

Login before adding your answer.

Traffic: 762 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6