plot of chunk
headband

Cette fonction permet de calculer une statistique par niveau d’un facteur. Par exemple, considérons les données accessibles ici:

## Warning: attempt to set 'sep' ignored

##    site taille
## 1 site1  3.125
## 2 site1  1.162
## 3 site1  4.222
## 4 site1  2.748
## 5 site1  3.625
## 6 site1  3.892

Pour connaître la moyenne (par exemple) de la taille en fonction du site, on peut procéder comme suit:

tapply(taille,site,"mean")

## site1 site2 site3 site4 site5 site6 
## 3.327 4.770 3.035 2.510 3.637 2.189

Il est possible de calculer toutes sortes de métriques, comme la moyenne, la médiane, la variance, le minimum ou le maximum…

tapply(taille,site,"var")

##  site1  site2  site3  site4  site5  site6 
## 1.1134 2.5378 6.4559 0.8254 4.6700 2.8530

Si vous souhaitez préciser un argument supplémentaire pour la fonction qui vous intéresse (par exemple la fonction quantile) voici comment procéder:

## pour le quantile d'ordre 10% :
quantile(taille,probs=0.1)

##    10% 
## 0.9214

tapply(taille,site,"quantile",probs=0.1)

##  site1  site2  site3  site4  site5  site6 
## 2.0285 3.1482 0.7346 1.2759 1.4641 0.6844