Группировка по факторным переменным и нахождение выбросов в каждой группе

is_outlier.py

is_outlier = lambda group: (np.abs(group - group.mean()) > 2 * group.std()).astype(int)
factors = df.select_dtypes(exclude=[np.number]).columns.tolist()
df['is_outlier'] = df.groupby(factors).transform(is_outlier)

is_outlier.r

tr1 <- split(test[, sapply(test, is.numeric)], lapply(test[, sapply(test, is.factor)], factor))    
tr2 <- lapply(tr1, function(x) ifelse(abs(x - mean(x)) < 2 * sd(x), 0, 1))    
test$is_outlier <- unsplit(tr2,  lapply(test[,sapply(test, is.factor)], factor))