# NOT RUN {
## normal data
x <- rnorm(100)
mean(x)
median(x)
skippedMean(x)
sd(x)
mad(x)
skippedSD(x)
## Tukey's gross error model
## (1-eps)*Norm(mean, sd = sigma) + eps*Norm(mean, sd = 3*sigma)
ind <- rbinom(100, size = 1, prob = 0.1)
x.err <- (1-ind)*x + ind*rnorm(100, sd = 3)
mean(x.err)
median(x.err)
skippedMean(x.err)
sd(x.err)
mad(x.err)
skippedSD(x.err)
# }
Run the code above in your browser using DataLab