J'ai besoin d'un coup de pouce pour pouvoir calculer le nombre d'événements A par mois pour un capteur donné.
Mon jeu de données se présente comme suit:
Code : Tout sélectionner
a=seq(as.Date("2000/1/1"), by = "day", length.out = 1000)
n=20
date=sample(a,n)
tab=structure(list(
capteur1 = structure(sample(1:2,size =n,replace = T), .Label = c('A','B'), class = "factor"),
capteur2 = structure(sample(1:2,size =n,replace = T), .Label = c('A','B'), class = "factor"),
capteur3 = structure(sample(1:2,size =n,replace = T), .Label = c('A','B'), class = "factor")),
.Names = c("capteur1","capteur2","capteur3"), class = "data.frame", row.names = c(NA,-20L))
tab=cbind(date,tab)
tab
Pour le capteur 3 par exemple, j'aimerai avoir un résultat du genre:
Mois Annee Nbr A
XX YYYY n1
TT ZZZZ n2
................
J'ai essayé de le faire avec la fonction 'aggregate':
Code : Tout sélectionner
annee=substr(tab$date,1,4)
mois=substr(tab$date,6,7)
jour=substr(tab$date,9,10)
aggregate(tab$capteur3,list(mois,annee),FUN)
# FUN étant la fonction que je recherche
Merci !