R/characterizeBreeds.R

"characterizeBreeds" <- function(phen){
  BreedNames <- unique(phen$Breed)
  breed <- setNames(vector("list", length(BreedNames)), BreedNames)

  for(b in BreedNames){
    breed[[b]] <- list(Groups        = if(any(is.na(phen$Sex[phen$Breed==b]))){c("cohort")}else{c("male","female")},
                       hasCandidates = any(phen$isCandidate[phen$Breed==b]))
  }
  breed
}

Try the optiSel package in your browser

Any scripts or data that you put into this service are public.

optiSel documentation built on Jan. 16, 2021, 5:30 p.m.