R/pedigree.trim.R

In kinship2: Pedigree Functions

# Automatically generated from all.nw using noweb

pedigree.trim <- function(removeID, ped){
## trim subjects from a pedigree who match the removeID 
## trim relation matrix as well

  if(!("pedigree" %in% class(ped)))
    stop("Must be a pegigree object.\n")

  rmidx <- match(removeID, ped$id)
  if(length(rmidx)>0) {
    pedtrimmed <- ped[-rmidx]
    return(pedtrimmed)
  } else {
  return(ped)
  }
}