我尝试使用by
函数来替换数据框中许多变量的异常值,根据组变量进行操作。以下是我的努力,但是我遇到了一个错误。
# R code:
library(tidyverse)
library(dplyr)
# outlier function definition
my_outlier <- function(x){
stats <- boxplot.stats(unlist(x))
outlier <- stats$out
outlier_idx <- which(unlist(x) %in% outlier)
max <- max(x[-outlier_idx]); min <- min(x[-outlier_idx])
x <- ifelse(x>max, max,ifelse(x < min, min, x) )
return(x)
}
# use the above defined func to substitue outliers of 1 variable in a dataframe, according to a Group variable.
group_data <- as_tibble(data.frame(x=c(runif(10),2.5,-2.3,runif(10,1,2),3.5,-1.5), group=c(rep(1,12),rep(2,12)) ) )
View(group_data)
by(group_data$x, group_data$group, my_outlier, simplify=FALSE)
# use the above defined func to substitue outliers of 1+ variable in a dataframe, according to a Group variable.
group_datas <- as_tibble(data.frame(x=c(runif(10),2.5,-2.3,runif(10,1,2),3.5,-1.5),
y=c(runif(10,2,3),4,-1,runif(10,3,4),6,-1),
group=c(rep(1,12),rep(2,12)) ) )
by(group_data[,1:2], group_data$group, my_outlier)
当我使用自定义函数来替换数据框中一个或多个变量的异常值时,根据分组变量进行操作时,出现了错误。
我不知道我的代码哪一部分引起了这个错误。
max <- max(x[-outlier_idx]); min <- min(x[-outlier_idx])
中计算最大值的是哪一列?如果你想包括分组列在内计算整个矩阵,你需要在outlier_idx
后加上逗号,如max <- max(x[-outlier_idx,]); min <- min(x[-outlier_idx,])
。同时请注意,ifelse
只返回标量而不是向量。 - discipulus