Description Usage Arguments Details Value Author(s) Examples
Simulate recombination from Poisson process
1 2 | simuPoisson(parentsGeno,chr,cM,N)
|
parentsGeno |
genotype matrix for parents on all chromosomes: a data.frame or matrix with two rows, each row represent the marker gentoypes for one parent. Genotypes must be coded additively: the heterozygotes genotype must be coded as half the two homozygotes gentoypes |
chr |
a vector of size p for the chromosome numbers of all markers |
cM |
a vector of size p for the centiMorgan map positions for each marker |
N |
total number of lines to simulate |
mapchr |
map for a single chromosome |
simulation of genotypes is done for each individual chromosome. minimum and maximum number of recombinations was set up based on the 95% quantile of Poisson distribution minrecmb=qpois(0.025,chrlength/100,lower.tail=T) maxrecmb=qpois(0.975,chrlength/100,lower.tail=T)
simuPoisson |
returns a matrix with N rows, each row is the simulated genotype for an individual |
Lian Lian
1 2 3 | data(parentsGeno)
data(map)
dat=simuPoisson(parentsGeno,map$chr,map$cM,10)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.