Promedio por grupos en R


Los ejemplos siguientes muestran como obtener estadísticos por grupo en R usando la función vectorizada tapply y la función aggregate. Se puede verificar se pueden usar para obtener un promedio por grupos o por contingencia de grupos.

> head(iris)
  Sepal.Length Sepal.Width Petal.Length Petal.Width Species
1          5.1         3.5          1.4         0.2  setosa
2          4.9         3.0          1.4         0.2  setosa
3          4.7         3.2          1.3         0.2  setosa
4          4.6         3.1          1.5         0.2  setosa
5          5.0         3.6          1.4         0.2  setosa
6          5.4         3.9          1.7         0.4  setosa

> attach(iris)

> tapply(Sepal.Length, Species, mean)
    setosa versicolor  virginica
     5.006      5.936      6.588

> aggregate(Sepal.Length~Species, data=iris, mean)
     Species Sepal.Length
1     setosa        5.006
2 versicolor        5.936
3  virginica        6.588

> tapply(Sepal.Length, list(Species,Petal.Width>1), mean)
              FALSE    TRUE
setosa     5.006000      NA
versicolor 5.414286 6.02093
virginica        NA 6.58800

> aggregate(Sepal.Length~Species+(Petal.Width>1), data=iris, mean)
     Species Petal.Width > 1 Sepal.Length
1     setosa           FALSE     5.006000
2 versicolor           FALSE     5.414286
3 versicolor            TRUE     6.020930
4  virginica            TRUE     6.588000


Entradas populares de este blog

Valores perdidos (NA) en R: identificación y tratamiento (I)

Ejemplos de tablas en Stata

Métodos y técnicas de Investigación Económica