Sib-pair Command: simulate pedigrees


ClassAnalysis and data manipulation command
Namesimulate pedigrees
Arguments [<nped> [<ngen> [<min_number_of offspring> [<max_number_of_offpsring> ]]]]

Generate a set of nped (default 100) random pedigrees, each of ngen (default 2) generations. The component nuclear families each contain between min_number_of offspring and max_number_of_offpsring offspring (defaulting to a range 0-2). The generated pedigrees are each descended from a single founder couple (with marry-ins).

Example:


set seeds 12234 2323 13143
set loc trait  aff
set loc marker marker
set loc m2     marker
set loc m3     marker
#
# Ten three-generation families with 3-5 children per sibship
#
simulate pedigrees 10 3 3 5
run

sim m3 0.9 0.1
sim marker 0.9 0.1
if (marker=="1/1") then trait=(rand <= 0.15) else trait=(rand <= 0.3)
#
# Censor genotypes in copy of original trait locus
#
m2=marker
if (trait and rand <= 0.6) then m2="x/x" : m3="x/x"
if (not trait and rand <= 0.3) then m2="x/x" : m3="x/x"

set plevel 1
show pedigrees

<< (simulate qtl)Up to index>> (permute)