Call main Fortran part of sequoia, and convert its output to a list with dataframes.
SeqParSib(
ParSib = "par",
Specs = NULL,
ErrM = NULL,
GenoM = NULL,
LhIN = NULL,
AgePriors = NULL,
Parents = NULL,
quiet = FALSE
)
either "par" to call parentage assignment, or "sib" to call the rest of the algorithm.
a named vector with parameter values, as generated by
SeqPrep
.
3x3 matrix with genotyping error probabilities; rows=actual, columns=observed
matrix with genotype data, size nInd x nSnp
life history data: ID - sex - birth year
matrix with agepriors, size Specs["nAgeClasses"] by 8.
matrix with rownumbers of assigned parents, size nInd by 2
suppress messages
A list with
the pedigree
Non-assigned likely relatives
Non-assigned parent-parent-offspring trios
Info on dummies
Total log-likelihood per iteration
For a detailed description of the output see sequoia