x <- c(rnorm(200, 0.4, 0.05), rnorm(200, 0.1, 0.05))
x[round(runif(20)*400)] <- rnorm(20, 0.2, 0.05)
x[round(runif(20)*400)] <- rnorm(20, 0.3, 0.05)
y <- denoise(x)
plot(x)
points(y, pch=20, col="red")
Run the code above in your browser using DataLab