R/preProcessX.r

In modehunt: Multiscale Analysis for Density Functions

preProcessX <- function(X.raw, lower = -Inf, upper = Inf){

n.raw <- length(X.raw)

## case 1: (a,b) = (0,inf)
if (lower==0 & upper==Inf){
    n <- n.raw+1
    X <- 1:n*0
    X[2:n] <- sort(X.raw)}

## case 2: (a,b) = (0,1)
if (lower==0 & upper==1){
    n <- n.raw+2
    X <- 1:n*0
    X[2:(n-1)] <- sort(X.raw)
    X[n] <- 1}

## case 3: (a,b) = (-inf,0)
if (lower==-Inf & upper==0){
    n <- n.raw+1
    X <- 1:n*0
    X[1:(n-1)] <- sort(X.raw)}

## case 4: (a,b) = (a,b)
if (lower > -Inf & upper < Inf){
    n <- n.raw + (lower > -Inf) + (upper < Inf)
    X <- 1:n*0
    if (lower > -Inf){
        X[1] <- lower
        X[2:(n.raw+1)] <- sort(X.raw)} else {X[1:n.raw] <- sort(X.raw)}
    if (upper < Inf){X[n] <- upper}}

## case 5: (a,b) = (-inf,inf)
if (lower==-Inf & upper==Inf){
    n <- n.raw
    X <- sort(X.raw)}

## remove ties via perturbing a little bit
d <- diff(X)
if (min(d) == 0){
    X <- X + 10^(-3) * mean(abs(d)) * rnorm(n)
    X <- sort(X)
    cat('Warning: Data contain ties and have been perturbed a litte!')
}


X <- matrix(sort(X), ncol = 1)
return(X)
}