Nothing
gammas <- function(individuals, pedigree) {
stopifnot(is.character(individuals) || is.numeric(individuals))
if (is.numeric(individuals))
storage.mode(individuals) <- "integer"
if (is.integer(individuals) && any(individuals <= 0))
stop("individuals, if integer-valued, must be positive-valued")
stopifnot(is.matrix(pedigree))
stopifnot(ncol(pedigree) == 3)
stopifnot(is.character(pedigree) || is.numeric(pedigree))
if (is.numeric(pedigree))
storage.mode(pedigree) <- "integer"
if (is.integer(pedigree) && any(pedigree <= 0))
stop("pedigree, if integer-valued, must be positive-valued")
stopifnot(typeof(individuals) == typeof(pedigree))
stopifnot(individuals %in% pedigree)
indall <- sort(unique(pedigree))
founders <- setdiff(indall, pedigree[ , 1])
result <- matrix(NA_real_, length(founders), length(individuals))
geneset <- as.integer(2)
for (i in seq(along = founders))
for (j in seq(along = individuals)) {
names(geneset) <- founders[i]
result[i, j] <- descent(individuals[j], pedigree, geneset)
}
rownames(result) <- founders
colnames(result) <- individuals
result
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.