simulfams: 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

simulfams(n.sib=sample(5:20,size=sample(5:20,size=1),replace=TRUE),
          p=runif(sample(5:15,size=1),min=0.1,max=0.6))

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)

Laura Plantinga and Karl 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

simrapd

Examples

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

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