Learn R Programming

hsrecombi (version 1.0.1)

bestmapfun: Best fitting genetic-map function

Description

Approximation of mixing parameter of system of map functions

Usage

bestmapfun(theta, dist_M)

Value

list (LEN 2)

mixing

mixing parameter of system of genetic mapping functions

mse

minimum value of target function (theta - dist_M)^2

Arguments

theta

vector of recombination rates

dist_M

vector of genetic positions

Details

The genetic mapping function that fits best to the genetic data (recombination rate and genetic distances) is obtained from Rao's system of genetic-map functions. The corresponding mixing parameter is estimated via 1-dimensional constrained optimisation. See vignette for its application to estimated data.

References

Rao, D.C., Morton, N.E., Lindsten, J., Hulten, M. & Yee, S (1977) A mapping function for man. Human Heredity 27: 99-104. tools:::Rd_expr_doi("10.1159/000152856")

Examples

Run this code
  theta <- seq(0, 0.5, 0.01)
  gendist <- -log(1 - 2 * theta) / 2
  bestmapfun(theta, gendist)

Run the code above in your browser using DataLab