Идентификация по группам в data.frame

Question

Я хотел бы определить случай группами, которые просто больше, чем avg plus sd. Например, использование видов как группы и petal.wid как моя переменная в данных диафрагмы.

Каков лучший способ сделать это? создание функции?

Я сделал это, но я не могу установить отношение к оригинальным данным для идентификации дела.

data(iris) 
library(plyr) 
petal.wid.avg <- ddply(iris, .(Species), function(df) 
    return(c(petal.wid.avg=mean(df$Petal.Width), petal.wid.sd=sd(df$Petal.Width))) 
) 
petal.wid.avg$avgsd <- petal.wid.avg$petal.wid.avg + petal.wid.avg$petal.wid.sd 
petal.wid.avg 

Answer 1

Есть много способов сделать это, но функция ave, пожалуй, самый простой.

iris$big <- with(iris, 
    ave(Petal.Width, Species, FUN = function(x) x > mean(x) + sd(x)) 
) 

Вот plyr решение:

iris <- ddply(
    datasets::iris, 
    .(Species), 
    transform, 
    big = Petal.Width > mean(Petal.Width) + sd(Petal.Width) 
) 

---

Baed на комментарии, вот остальная часть решения.

iris <- subset(iris, big) 
iris <- ddply(
    iris, 
    .(Species), 
    transform, 
    smallest = Petal.Width == min(Petal.Width) 
) 
(iris <- subset(iris, smallest)) 

Обратите внимание, что там, где у вас есть связи (как в этом наборе данных), вы не получите уникальную строку «только больше».