simuPoisson: simulate genotypes from Poisson process

In lian0090/simuPoisson: Simulate F2 populations

Description

Simulate recombination from Poisson process

Usage

	simuPoisson(parentsGeno,chr,cM,N)
	

Arguments

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

Details

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)

Value

simuPoisson

returns a matrix with N rows, each row is the simulated genotype for an individual

Author(s)

Lian Lian

Examples

	data(parentsGeno)
	data(map)
dat=simuPoisson(parentsGeno,map$chr,map$cM,10)

lian0090/simuPoisson documentation built on May 21, 2019, 6:11 a.m.