2017-10-26 5 views
0

Ich zähle Z-Scores für mehrere Variablen von zwei Gruppen.Zähle mehrere Z-Scores mit Dplyr Mutate Dpply/Lapply

Hier ist ein Beispiel:

data = mtcars

Die Variablen ich die Z-Werte erhalten möchten:

vars <- c("mpg", "disp", "hp", "drat", "wt", "qsec")

Z-Score für eine variable Zählung (in Betrieb):

mtcars %>% group_by(am, vs) %>% mutate(z_mpg = (mpg - mean(mpg))/sd(mpg))

Das Problem ist, dass ich nicht mit dem vorherigen Code arbeiten oder arbeiten kann, um alle "vars" -Variablen auszuführen, damit ich alle Z-Scores auf einmal bekomme.

Wenn Sie eine Idee haben, wie Sie dies tun, indem Sie Daten normalisieren (Mittelwert 0, SD 1), während Sie die Gruppen berücksichtigen, anstatt z-Scoring, würde mir das auch helfen.

Danke!

Antwort

0

Sie würden mutate_at verwenden und funs verwenden, um Ihre Z-Score-Funktion zu definieren. In diesem Fall wird . verwendet, um die Spalte anzugeben, die Sie mutieren.

mtcars %>% 
    group_by(am, vs) %>% 
    mutate_at(.cols = vars, funs(z = (. - mean(.))/sd(.))) 

Source: local data frame [32 x 17] 
Groups: am, vs [4] 

    mpg cyl disp hp drat wt qsec vs am gear carb  mpg_z_ disp_z_  hp_z_ drat_z_  wt_z_  qsec_z_ 
    <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>  <dbl>  <dbl>  <dbl>  <dbl>  <dbl>  <dbl> 
1 21.0  6 160.0 110 3.90 2.620 16.46  0  1  4  4 0.3118089 -0.4852978 -0.7168218 -0.1024154 -0.48795905 0.60787578 
2 21.0  6 160.0 110 3.90 2.875 17.02  0  1  4  4 0.3118089 -0.4852978 -0.7168218 -0.1024154 0.03595488 1.12105734 
3 22.8  4 108.0 93 3.85 2.320 18.61  1  1  4  1 -1.1710339 0.9679756 0.5147599 -0.7890520 0.66286051 -0.09519147 
4 21.4  6 258.0 110 3.08 3.215 19.44  1  0  3  1 0.2659345 1.6870444 0.3753676 -1.0547191 0.05956492 -0.36000354 
5 18.7  8 360.0 175 3.15 3.440 17.02  0  0  3  2 1.3156017 0.0331832 -0.5745432 0.1266602 -0.86434641 -0.15281032 
6 18.1  6 225.0 105 2.76 3.460 20.22  1  0  3  1 -1.0695190 1.0153670 0.1364973 -1.7435153 0.76407410 0.17268463 
7 14.3  8 360.0 245 3.21 3.570 15.84  0  0  3  4 -0.2703291 0.0331832 1.5237884 0.3872184 -0.69514319 -1.62477907 
8 24.4  4 146.7 62 3.69 3.190 20.00  1  0  4  2 1.4799831 -0.5783405 -1.9177872 0.2582986 -0..02243925 
9 22.8  4 140.8 95 3.92 3.150 22.90  1  0  4  2 0.8324905 -0.6984282 -0.3412433 0.7533708 -0.12734568 2.00294653 
10 19.2  6 167.6 123 3.92 3.440 18.30  1  0  4  4 -0.6243679 -0.1529447 0.9964304 0.7533708 0.70656315 -1.13854778 
# ... with 22 more rows 
+0

Danke! @ Jake Kaupp, das funktioniert! Kann ich so Daten normalisieren? 'mtcars%>% group_by (am, vs)%>% mutate_at (.cols = vars, Spaß (n = (. - min (.))/(Max (.) - min (.)))) ' – J5H

+0

Wenn Sie sich die Zeit genommen haben, den Code zu schreiben, sollten Sie wissen, ob es funktioniert hat! –

Verwandte Themen