rGTs {KinMixLite} | R Documentation |
Simulate random genotype profiles and DNA samples for related individuals
Description
Simulate random genotype profiles and DNA samples for arbitrarily related individuals
Usage
rGTs(nreps,IBD,db,DNA,sex=rep(0,ncontr),nU=0)
Arguments
nreps |
Integer, number of replicates |
IBD |
Specification of relationships, as in |
db |
Data frame, database of alleles and their frequencies, for each marker; for format, see |
DNA |
Integer vector, numbers of DNA cells for the respective individuals, can be NA |
sex |
Integer vector, sex of the respective contributors: 1=male, 2=female, 0=unspecified |
nU |
Integer, include also this number of unrelated individuals |
Details
Genotype profiles are generated randomly, using the allele frequency database db
, under the relationships specified by the IBD
argument. In accordance with the underlying biology, allele values for the AMEL marker (if this is one of the markers included) are not influenced by relationships with other individuals; however they are influenced by the sex of the individuals, where this is known. Information on sex can be specified by the optional argument sex
: a male is given the profile X-Y, a female X-X, and an individual with unspecified sex X-X or X-Y with equal probabilities.
Value
Data frame with variables Sim
, Sample.name
, Marker
, Allele
, and DNA
, suitable for input to simExtraction
, etc. See package pcrsim
.
Author(s)
Peter Green (P.J.Green@bristol.ac.uk)
Examples
data(test2data)
data(NGMDyes)