Discussion :

[maf]

Bonjour et merci d'avance de votre aide ...

J'ai un fichier dans lequel j'ai deux colonnes, une première V1 dans laquelle il y a des données relatives au temps. (1 = 00h01, 14 = 00h14, 100=01h00, 547=05h47 etc...).

Une deuxième, V2, uniquement avec des chiffres (de 0 à 500 environ)

J'aimerais pouvoir effectuer la moyenne de la colonne V2 par tranche de 15 minutes (relatif à la colonne V1)

J'ai essayé des combinaisons avec breaks, ave, by ... mais je n'y arrive pas.

Merci à celle ou celui qui aurait une solution à cela

[dev_ggy]

Question qui me semble important avant de te données une réponse satisfaisante. Les données de temps (V1) ont elles des écarts constants ou bien il y a des trous dans la chronologie ?

dans le premier cas tu peux

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
data <- rep(1:20,rep(15,20))
Tes_Donnees$data <- data
aggregate(Tes_Donnees, by = list(data) , mean)

Avec Tes_Donnees qui sont triés par ordre croissant.
Dans se cas, tu as la moyenne de tes données pour chaque quart d'heure.

Cordialement.

[maf]

Bonjour,

Il y a effectivement des "trous" car après 59 --> 00h59 on passe directement à 100 --> 01h00.

Et il y a plusieurs fois les mêmes heures car les relevés sont effectués de manière aléatoire. Et certaines heures ne sont pas présentes.

[dev_ggy]

1) dédoublonner:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
 
doublon<-duplicated(as.vector(Tes_Donnees$V1))
data<-data.frame(Tes_Donnees,doublon)
data<-data[data$doublon==FALSE,]

2) Merge

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
rep <- rep(1:20,rep(15,20))
xy<-merge(rep, data, by = V1,  all.x = TRUE )

3) Les sommes

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

aggregate(V1~rep, xy, FUN = mean, na.action = na.omit)

Voila les trois points à développer. Le code est peu être pas tout à fait exacte. A toi de le modifier en conséquence. Mais je pense que les idées sont là.

Cordialement.

[vchouraki]

Bonjour,

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

?tapply

Le problème est plus de convertir des heures (date) en minutes (nombre)... Avez-vous également des dates ou est-que l'ensemble de vos mesures est faite le même jour?

Un exemple :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
# Jeu de data bidon
test <- data.frame(a=rnorm(81), b=c(15:59,500:512,1323:1345))
# test <- data.frame(a=rnorm(240), b=c(0:59,100:159,500:559,1300:1359))
# Formatage de b en date
test$b <- strptime(sprintf("%04d", test$b), format="%H%M")
# Calcul de la difference avec la date "minimum" en minutes
test$b.diff <- as.numeric(difftime(test$b, min(test$b), units="min"))
# Creation d'une variable qui groupe les donnees par tranches de 15 minutes
test$b.cut <- test$b.diff %/% 15
# Calcul de la moyenne par groupe
with(test, tapply(a, b.cut, mean))

HTH

Vincent