punc.R

punc <- function(x){
  len = length(x)
  pp <- which(x > 1.5*IQR(x) + quantile(x,3/4))
  pplen <- length(pp)
  group <- rep(1,len)
  group[pp] <- 1
  if(pplen==0){
    return(group)
    break
  }
  if(pplen==1){
    group[pp] <- 2
    group[(pp+1):len] <- 3
    return(group)
    break
  }
  for(i in 1:(pplen-1)){
    group[pp[i]] <- i + 1
    if(pp[i+1]- pp[i] !=1){
      group[(pp[i]+1):(pp[i+1]-1)] <- i+2 
    }
  }
  group[pp[pplen]] <- i+3
  group[(pp[pplen]+1):(len)] <- i+4
  return(as.numeric(as.factor(group)))
}