Создать вектор (данные) нормального распределения с выбросами?

Question

В R, как вы генерируете вектор (данные) с выбросами? Отлично, если данные «приемлемы» нормально распределены.

Answer 1

Это действительно зависит от определения «outlier»;

c(rnorm(100), 100, -100) # an egregious example 
    plot(density(c(rnorm(90), rnorm(5, 1)))) # not as egregious 

Answer 2

@DWin прав, что это зависит от того, что вы подразумеваете под «outlier». Для записи я использую то же определение, которое он использует, поэтому я бы использовал (использовал) что-то вроде кода he и @ Ferdinand.kraft, list. Другие иногда означают дату, более экстремальное, чем обычно. Это сложно определить для симуляционного исследования, но общее определение - это точка, более чем в 1,5 раза превышающая межквартильный диапазон за 1-й (3-й) квартиль. Вот простой способ найти, что (я уверен, что будет более эффективные способы):

flag <- 0 
while(flag==0){         
    X <- rnorm(N)          
    bp <- boxplot(X, plot=FALSE) 
    if(length(bp$out)!=0){ 
    flag <- 1 
    } 
}