R/mcirktau.R

"mcirktau" <- function(data, posorder, ws=NULL){
  
  # DATE WRITTEN: 22 Ene 2013          LAST REVISED:  22 Ene 2013
  # AUTHOR: Sandra Barragan.
  # DESCRIPTION: Computes the mean Circular Kendall's tau between
  # a circular order and a data set.
  # REFERENCE:  Fisher (1993) 
  
  ntaus <- apply(data, 1, cirKendall, posorder)
  if(is.null(ws)){ws <- rep(1, nrow(data))}
  return(list(mtau=weighted.mean(ntaus, w=ws, na.rm=TRUE), ntaus))
}

Try the isocir package in your browser

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

isocir documentation built on Aug. 17, 2023, 9:07 a.m.