Familias2linkdat {paramlink} | R Documentation |
Convert 'Familias' output to linkdat objects
Description
Familias is a widely used program for computations in forensic genetics. The
function documented here facilitates the use of paramlink
for
specialized computations which are not implemented in Familias, e.g.
conditional simulations.
Usage
Familias2linkdat(familiasped, datamatrix, loci)
readFamiliasLoci(loci)
connectedComponents(ID, FID, MID)
Arguments
familiasped |
A |
datamatrix |
A data frame with two columns per marker (one for each allele) and one row per individual. |
loci |
A |
ID |
An integer vector: Individual ID. |
FID |
An integer vector: ID of father. |
MID |
An integer vector: ID of mother. |
Details
The Familias program represents pedigrees and marker data in a way that
differs from paramlink
in several ways, mostly because of
paramlink
's stricter definition of a 'pedigree'. In paramlink
,
a pedigree must be connected, have numerical IDs, and each member must have
either 0 or 2 parents present in the pedigree. None of this is required by
FamiliasPedigree
objects. The conversion function
Familias2linkdat
takes care of all of these potential differences: It
converts each FamiliasPedigree
into a list of connected linkdat objects,
additional parents are added where needed, and non-numerical ID labels are
stored in the plot.labels
slot of the linkdat object(s).
Value
A linkdat
object, or a list of such.
Author(s)
Magnus Dehli Vigeland, Thore Egeland
References
Windows Familias is freely available from https://familias.name.
Examples
x = nuclearPed(1)