R/dconcordant.R

In MultNonParam: Multivariate Nonparametric Methods

#' Calculate the probability atom of the count of concordant pairs among indpendent pairs of random variables.
#' @param ss Integer number of pairs
#' @param nn number of pairs
#' @return real probability
#' @export
dconcordant<-function(ss,nn){
   if(nn<=20){
      if(ss<0) dc<-0
      if(ss>(nn*(nn-1)/2)) dc<-0
      if((ss>=0)&(ss<=(nn*(nn-1)/2))) dc<-.Fortran("dconcordant", as.integer(ss), as.integer(nn), dc=as.double(0.0), PACKAGE="MultNonParam")$dc
   }else{
     dc<-NA
   }
   return(dc)
}