R/JitterUniqueValues.R

In Mthrun/DataVisualizations: Visualizations of High-Dimensional Data

JitterUniqueValues=function(Data,Npoints=20,min=0.99999,max=1.00001){
  
  vals=unique(Data)
  nu=length(vals)
  n=length(Data)
  Flag=FALSE
  if(is.null(Npoints)){#generate vector of length data specifically
    Npoints=round(n/nu,0)
    Flag=TRUE
  }else{
    if(Npoints<1) Npoints=1
  }
  diff=max-min
  DataJitter=c()
  for(v in vals){
    if(is.finite(v)){
      if(v!=0){
        x <- c(v, v * runif(Npoints-1, min, max))
       
      }else{
        x <- c(v, v + runif(Npoints-1, -diff, diff))
      }
    }else{
      x=c(v,rep(v,Npoints-1))
    }
    DataJitter=c(DataJitter,x)
  }
  if(Flag){#check if rounding had some effect
    if(n!=length(DataJitter)){#
      DataJitter=sample(DataJitter,size = n,replace = T)
    }
  }
  return(DataJitter)
}