R/Multi.Jack.R

#'
#' @title Jack-knife indices in n topologies m times.
#'
#' @description
#' The function calculates the indices values for a MultiData list  m (=replicates) times 
#'
#' @param  MultiData is the list of Trees and distributions to evaluate, a list object.
#' 
#' @param times in the number of times to repeat the process, an integer.
#'
#' @param jtip is the proportion of terminals to delete, real (range 0-1).
#'
#' @param jtopol is the proportion of topologies to delete, real (range 0-1).
#'
#' @returns the rankings
#'
#' @examples
#'  ## get the library
#'  library(jrich)
#'  
#'  ## load the data
#'  data(Multitaxon1) 
#' 
#' Multi.Jack(Multitaxon1, jtip=0.25)
#'
#'@author Miranda-Esquivel Daniel R.
#'
#'




Multi.Jack <- 
function (MultiData = MultiData, times = 100, jtip = 0, jtopol = 0) {
  
  results <- list() 
  
  for(i in 1:times){
    results[[i]] <- Rank.Indices(Multi.Index.Calc(MultiData = MultiData, 
                                 jtip = jtip, jtopol = jtopol))
  }
  
  return(results)
}

Try the jrich package in your browser

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

jrich documentation built on May 2, 2019, 3:45 a.m.