sim.haplotype {rres}

Simulate artificial haplotypes.

Description

sim.haplotype returns haplotypes of the specified number of SNPs simulated under linkage equilibrium.

Usage

sim.haplotype(freq, nhaplo)

Arguments

freq

vector of values between 0 and 1.

nhaplo

positive integer.

Details

freq are reference allele frequencies. nhaplo haplotypes are simulated independently.

Value

A matrix of nhaplo rows and length(freq) columns. Reference alleles are coded 1, alternate alleles are coded 2.

Examples

nsnp = 7 # number of SNPs
freq = runif(nsnp, 0.05, 0.95)
nhaplo = 4 # number of founder haplotypes
sim.haplotype(freq, nhaplo)