calcFE: Calculates founder Equivalents

In rmsharp/nprcmanager: Genetic Tools for Colony Management

Description

Part of the Genetic Value Analysis

Usage

calcFE(ped)

Arguments

ped	the pedigree information in datatable format. Pedigree (req. fields: id, sire, dam, gen, population).

Details

It is assumed that the pedigree has no partial parentage

Value

The founder equivalents FE = 1 / sum(p ^ 2), where p is average number of descendants and r is the mean number of founder alleles retained in the gene dropping experiment.

Examples

## Example from Analysis of Founder Representation in Pedigrees: Founder
## Equivalents and Founder Genome Equivalents.
## Zoo Biology 8:111-123, (1989) by Robert C. Lacy
library(nprcgenekeepr)
ped <- data.frame(
id = c("A", "B", "C", "D", "E", "F", "G"),
sire = c(NA, NA, "A", "A", NA, "D", "D"),
dam = c(NA, NA, "B", "B", NA, "E", "E"),
stringsAsFactors = FALSE
)
ped["gen"] <- findGeneration(ped$id, ped$sire, ped$dam)
ped$population <- getGVPopulation(ped, NULL)
pedFactors <- data.frame(
  id = c("A", "B", "C", "D", "E", "F", "G"),
  sire = c(NA, NA, "A", "A", NA, "D", "D"),
  dam = c(NA, NA, "B", "B", NA, "E", "E"),
  stringsAsFactors = TRUE
)
pedFactors["gen"] <- findGeneration(pedFactors$id, pedFactors$sire,
                                    pedFactors$dam)
pedFactors$population <- getGVPopulation(pedFactors, NULL)
fe <- calcFE(ped)
feFactors <- calcFE(pedFactors)

rmsharp/nprcmanager documentation built on April 24, 2021, 3:13 p.m.