simrapd: Simulate RAPD data

In kbroman/fingers: Identifying Clusters of Related Individuals

Simulate RAPD data

Description

Simulates RAPD data for a set of sibling families.

Usage

simrapd(n.sib = rep(15,10), p = c(rep(0.125,8),rep(0.175,5),rep(0.225,5),
           rep(0.275,8),rep(0.325,3),rep(0.375,4),
           rep(0.475,4),rep(0.575,3)))

Arguments

n.sib	A vector giving the number of siblings per family (length is the number of families).
p	A vector of frequencies of the band allele at each marker (length is the number of markers).

Details

The RAPDs are assumed to be in Hardy-Weinberg equilibrium.

Value

A matrix of dimension (n.ind x n.mar), giving the RAPD phenotypes for each individual at each marker, with 1 indicating a band and 0 indicating no band.

Author(s)

Karl W Broman [e-mail]

References

BL Apostol, WC Black IV, BR Miller, P Reiter, BJ Beaty (1993) Estimation of the number of full sibling families at an oviposition site using RAPD-PCR markers: applications to the mosquito Aedes aegypti. Theor Appl Genet 86:991-1000.

See Also

simulfams

Examples

data <- simrapd(rep(20,5), p=runif(40, 0.1, 0.6))

kbroman/fingers documentation built on May 17, 2023, 11:50 p.m.