simbc: Simulate a backcross

In kbroman/qtlsim: Software for Testing QTL Analysis via Simulations

Simulate a backcross

Description

Simulates genotype and phenotype data for a backcross experiment with QTLs exactly at marker postions, and acting additively.

Usage

simbc(n.ind=100, n.mar=rep(11,9), mar.sp=rep(10,10*9),
      qtl.chr=c(1,1,2,2,3,4,5), qtl.mar=c(4,8,4,8,6,4,1),
      qtl.dist=rep(0,7), qtl.eff=c(1,1,1,-1,1,1,1)/2,
      sigma=1)

Arguments

n.ind	Number of backcross individuals.
n.mar	Vector indicating the number of markers on each chromosome.
mar.sp	Vector of length sum(n.mar)-n.chr-1, giving the inter-marker spacings in cM.
qtl.chr	Chromosomes on which qtls are sitting.
qtl.mar	Markers to the left of each qtl, numbered 1, 2, ...
qtl.dist	Distance between qtl and marker to left, in cM.
qtl.eff	QTL effects.
sigma	Residual (environmental) standard deviation.

Value

A list with two components: geno, a matrix containing the genotype data, and pheno, a vector of phenotypes.

Author(s)

Karl W Broman, [e-mail]

See Also

convertSim

Examples

dat <- simbc()

kbroman/qtlsim documentation built on May 17, 2023, 11:53 p.m.