simPedMCN: Middle Class Neighborhood pedigree construction

View source: R/simPed.R

simPedMCNR Documentation

Middle Class Neighborhood pedigree construction


Simulates a pedigree for the “middle class neighborhood” mating design (Shabalina, Yampolsky, and Kondrashov 1997).


simPedMCN(pedTemp, g, Nfam = NULL, noff = 2)



A data.frame pedigree of a template pedigree from which the middle class neighborhood design should continue. If NULL, a new pedigree will be created with Nfam families.


Integer number of generations to produce from the middle class neighborhood design


Integer number of families with which to start a new pedigree following the middle class neighborhood design.


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.

See Also

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)

nadiv documentation built on May 29, 2024, 10:40 a.m.