EBSeq/R/crit_fun.R

   1 crit_fun<-function (PPEE, thre)
   2 {
   3     y <- cumsum(sort(PPEE))/(1:length(PPEE))
   4     mm <- y < thre
   5     index <- sum(mm)
   6     if (index > 0) {
   7         out <- 1 - sort(PPEE)[index]
   8             }
   9     if (index == 0) {
  10                         out <- 1
  11                                     }
  12     names(out) <- NULL
  13     return(out)
  14 }
  15