2016-06-24 12 views
0

Betrachtet man dieses Beispiel Datenrahmen:Zusammenfassen von Daten mit zwei Funktionen in dplyr

d <- read.table(text=" 
    trt rep y 
    1 1 30 
    1 1 50 
    1 1 70 
    1 2 0 
    1 2 0 
    1 2 0 
    2 1 10 
    2 1 0 
    2 1 0 
    2 2 5 
    2 2 0 
    2 2 . 
    " 
    , header = TRUE, check.names = F, na.strings = ".") 

Ich versuche, eine Übersichtstabelle mit zwei Operationen des „y“ Variable zu erhalten.

Die erste neue Spalte sollte für jeden rep die einfachen Mittelwert von trt hat:

by_rep1 = d %>% 
    group_by(trt, rep) %>% 
    summarise(sev = mean(na.omit(y))) 

und die zweiten, den Anteil der positiven Werte von trt für jeden rep.

by_rep2 = d %>% 
    group_by(trt, rep) %>% 
    summarise_each(funs(round(mean(.>0, na.rm=TRUE),2)), y) 

ich so lange tun, denn ich bin nicht Ahnung, wie es in einem Schritt zu tun:

inner_join(by_rep1, by_rep2, by = c("trt", "rep")) 

# trt rep mean_y  y 
# (int) (int)  (dbl) (dbl) 
#1  1  1 50.000000 1.00 
#2  1  2 0.000000 0.00 
#3  2  1 3.333333 0.33 
#4  2  2 2.500000 0.50 

hat jemand weiß, wie man in einem einzigen Schritt zu tun, beide Funktionen verbinden?

Antwort

4

Sie können sie in einer einzigen summarize Anweisung setzen:

d %>% group_by(trt, rep) %>% summarise(mean_y = mean(y, na.rm = T), 
             y = round(mean(y > 0, na.rm = T), 2)) 
Source: local data frame [4 x 4] 
Groups: trt [?] 

    trt rep mean_y  y 
    (int) (int)  (dbl) (dbl) 
1  1  1 50.000000 1.00 
2  1  2 0.000000 0.00 
3  2  1 3.333333 0.33 
4  2  2 2.500000 0.50 
1

Wir auch data.table

library(data.table) 
setDT(d)[, .(mean_y = mean(y, na.rm = TRUE), y = round(mean(y > 0, 
       na.rm = TRUE), 2)) , .(trt, rep)] 
# trt rep mean_y y 
#1: 1 1 50.000000 1.00 
#2: 1 2 0.000000 0.00 
#3: 2 1 3.333333 0.33 
#4: 2 2 2.500000 0.50 

Oder dies kann getan werden, indem nur base R

do.call(data.frame, aggregate(y~., d, FUN = function(x) 
    c(mean_y=mean(x, na.rm=TRUE), y=round(mean(x > 0, na.rm=TRUE),2)), na.action=NULL)) 
verwenden können