Как удалить выбросы из нескольких столбцов в r

Часто вам может потребоваться удалить выбросы из нескольких столбцов одновременно в R.

Распространенный способ определить наблюдение как выброс — это если оно в 1,5 раза превышает межквартильный размах выше третьего квартиля (Q3) или в 1,5 раза превышает межквартильный размах ниже первого квартиля (Q1).

Используя это определение, мы можем использовать следующие шаги для создания простой функции для выявления выбросов, а затем применить эту функцию к нескольким столбцам в кадре данных R.

Шаг 1: Создайте фрейм данных.

Начнем с создания фрейма данных в R:

 df <- data.frame(index=c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10),
                 var1=c(4, 4, 5, 4, 3, 2, 8, 9, 4, 5),
                 var2=c(1, 2, 4, 4, 6, 9, 7, 8, 5, 29),
                 var3=c(9, 9, 9, 5, 5, 3, 4, 5, 11, 34))

Шаг 2: Определите функцию выброса.

Далее давайте определим функцию, которая может идентифицировать выбросы, и функцию, которая затем может удалять выбросы:

 outliers <- function (x) {

  Q1 <- quantile(x, probs =.25)
  Q3 <- quantile(x, probs =.75)
  iqr = Q3-Q1

 upper_limit = Q3 + (iqr*1.5)
 lower_limit = Q1 - (iqr*1.5)

 x > upper_limit | x < lower_limit
}

remove_outliers <- function (df, cols = names (df)) {
  for (col in cols) {
    df <- df[!outliers(df[[col]]),]
  }
  df
}

Шаг 3. Примените функцию выброса к фрейму данных.

Наконец, давайте применим эту функцию к нескольким столбцам фрейма данных, чтобы удалить выбросы:

 remove_outliers(df, c('var1', 'var2', 'var3'))

  index var1 var2 var3
1 1 4 1 9
2 2 4 2 9
3 3 5 4 9
4 4 4 4 5
5 5 3 6 5
9 9 4 5 11

Дополнительные руководства по R можно найти здесь .