simPedMCN | R Documentation |
Simulates a pedigree for the “middle class neighborhood” mating design (Shabalina, Yampolsky, and Kondrashov 1997).
simPedMCN(pedTemp, g, Nfam = NULL, noff = 2)
pedTemp |
A |
g |
Integer number of generations to produce from the middle class neighborhood design |
Nfam |
Integer number of families with which to start a new pedigree following the middle class neighborhood design. |
noff |
Integer number of full-sib offspring produced by each family (must be >=2). |
This creates a pedigree following a breeding design which maintains equal contributions to the next generation by each family in the design. It effectively removes the effect of natural selection which makes it amenable to quantify the contribution of mutations to phenotypic variance over the course of the breeding design.
For a starting pedigree template (pedTemp
), the last generation is used
as parents to begin the breeding design for the next g
generations.
The number of families in the last generation of the template pedigree
(pedTemp
) will be the number of families in each generation.
Alternatively, if no template pedigree is provided (pedTemp=NULL
),
Nfam
number of families will be produced in the first generation from
Nfam
unique sire and Nfam
unique dams.
Either pedTemp
or Nfam
must be NULL
, but not both.
A data.frame
with columns corresponding to: id, dam, sire, sex,
and generation. Sex is M
for males and F
for females. The
first generation produced in the middle class neighborhood scheme is assigned
a value of “1”, with their parents being assigned to generation
0
. If pedTemp
was provided, the generations from this pedigree
will be denoted with negative integers.
Shabalina, S.A, L.Y. Yampolsky, and A.S. Kondrashov. 1997. Rapid decline of fitness in panmictic populations of Drosophila melanogaster maintained under relaxed natural selection. Proc. Natl. Acad. Sci. USA. 94:13034-13039.
simPedHS
, simPedDFC
# No template pedigree provided - start from scrtach
mcn1 <- simPedMCN(pedTemp = NULL, g = 3, Nfam = 4, noff = 2)
# Provide a template pedigree (half-sib design)
hsped <- simPedHS(s = 2, d = 2, n = 4)
mcnHS <- simPedMCN(pedTemp = hsped, g = 3)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.