[R] Random number from density()

[Stephen Weigand]

On 4/2/07, Arianne ALBERT <Arianne.ALBERT at cefe.cnrs.fr> wrote:
> Hello,
>
> I'm writing some genetic simulations in R where I would like to place
> genes along a chromosome proportional to the density of markers in a
> given region. For example, a chromosome can be presented as a list of
> marker locations:
>
> Chr1<-c(0, 6.5, 17.5, 26.2, 30.5, 36.4, 44.8, 45.7, 47.8, 48.7, 49.2,
> 50.9, 52.9, 54.5, 56.5, 58.9, 61.2, 64.1)
>
> Where the numbers refer to the locations of markers along the
> chromosome. As you can see, there are a lot of markers around 50, but
> they are less dense elsewhere. I would like to take this information to
> randomly select a location on the chromosome to put a gene, where we are
> more likely to pick a location with high marker density (instead of
> using a uniform distribution between 0 and 64.1).
>
> So far I've used density(Chr1) to generate the probability density, but
> can this also be used to generate a random number? All the help I can
> find suggests getting the quantile function (the reciprocal of the
> integral of the pdf), however, it doesn't seem as though density()
> returns a pdf that can be manipulated in this way. Any suggestions?
>
> Thanks in advance,
>
> Arianne
>
This thread may be helpful:

http://tolstoy.newcastle.edu.au/R/e2/help/07/01/9139.html