Je travaille sur des données imbriquée de ce type.
Code : Tout sélectionner
head(data)
village survey night site location n
1 HKT 1 1 1 I 28
2 HKT 1 1 1 O 32
3 HKT 1 1 2 I 21
4 HKT 1 1 2 O 26
5 HKT 1 1 3 I 12
6 HKT 1 1 3 O 23
'data.frame': 1800 obs. of 6 variables:
$ village : Factor w/ 4 levels "HKT","HPN","KNH",..: 1 1 1 1 1 1 1 1 1 1 ...
$ survey : Factor w/ 9 levels "1","2","3","4",..: 1 1 1 1 1 1 1 1 1 1 ...
$ night : Factor w/ 5 levels "1","2","3","4",..: 1 1 1 1 1 1 1 1 1 1 ...
$ site : Factor w/ 5 levels "1","2","3","4",..: 1 1 2 2 3 3 4 4 5 5 ...
$ location: Factor w/ 2 levels "I","O": 1 2 1 2 1 2 1 2 1 2 ...
$ n : int 28 32 21 26 12 23 20 23 48 67 ...
Je voudrais calculer la moyenne de la colonne n par village et par survey et par nuit. Pour l'instant, j'arrive assez facilement à calculer cette moyenne sur un facteur (ex: village) en utilisant la fonction tapply :
Code : Tout sélectionner
tapply(data$n,data$village,mean)
HKT HPN KNH TOT
8.108966 9.108889 5.615556 13.066667
mais je n'arrive pas à combiner l'analyse sur plusieurs facteurs.
Merci d'avance pour votre aide,
Victor