Comment créer des périodes ordonnées mensuellement à partir de datevector sans boucler?

Question

Alors que la classe difftime est assez simple à utiliser lorsque l'on veut obtenir des différences annuelles ou quotidiennes entre les dates. Les différences mensuelles n'étaient pas aussi simples pour moi. Étant donné que l'arithmétique norme R permet des objets soustractions date Boucler dans un vecteur serait possible, mais here's une solution qui évite mise en boucle:

# generate reproducible example 
x<-seq(as.Date("2010-11-15"),as.Date("2011-05-20"),"months") 
y<-seq(as.Date("2010-04-15"),as.Date("2012-05-20"),"months") 
z<-seq(as.Date("2012-08-15"),as.Date("2013-05-20"),"months") 

d <- c(z,x,y) 

# my function – suggestions are welcome! 
getPeriods <- function(datevector){ 
x <- floor((as.numeric(datevector)-as.numeric(min(datevector)))/30) +1 
return(x) 
} 

# This returns a vector of monthly ordered periods. 

Y at-il une meilleure façon plus native de le faire? Par exemple, puis-je utiliser seq() en combinaison avec longueur() en quelque sorte? Je n'ai pas pu le faire parce que l'argument "to" de seq ne permet pas les vecteurs. Tout commentaire apprécié, d'ailleurs j'espère que cela peut aider quelqu'un d'autre ...

Answer 1

La solution standard consiste à convertir en POSIXlt et à partir de là.

getPeriods2 <- function(datevector){ 
    tmp <- as.POSIXlt(datevector) 
    tmp <- tmp$year*12 + tmp$mon 
    tmp - min(tmp) + 1 
} 

Ce ne tient pas compte des différences de jours, ce qui signifie que la différence entre le 31 mars et avril 1 est 1 dans ma fonction, mais est 0 dans le vôtre.

Maintenant que les mois ne sont pas de longueurs égales, la définition d'une "période mensuelle" est de toute façon délicate. Surtout quand vous courez sur quelques années, où vous pourriez avoir le problème de l'année bissextile. Une différence mensuelle n'est pas objective: dans votre définition, une différence mensuelle est de 30 jours, ce qui peut être discuté. Certaines personnes préfèrent 365/12 comme différence, de sorte qu'il reste correct lorsque les périodes s'étendent sur plusieurs années.

Voir par exemple votre fonction:

> tt <- as.Date(c("1998-01-01","2012-01-01")) 
> tt 
[1] "1998-01-01" "2012-01-01" 

> getPeriods(tt) 
[1] 1 171 

> diff(getPeriods(tt))/12 
[1] 14.16667 

> diff(getPeriods2(tt))/12 
[1] 14 

Ainsi, selon votre fonction, il y a 2 mois plus de 14 ans que il est en réalité. Mon idée personnelle: utiliser des mois comme des mois, pas comme des périodes de 30 jours. Je ne dis pas que votre fonction est mauvaise, c'est une approche qui est plus souvent utilisée. Mais en ce qui me concerne, vous allez sûrement avoir des ennuis à un moment donné.

Answer 2

Le paquet zoo possède de nombreuses fonctions utiles pour les séries temporelles, telles que les fonctions yearmon et as.yearmon. Les dates internes sont converties à l'année avec des mois décimaux (month.index/12), donc si vous voulez une sortie entière, vous devez probablement multiplier par 12 et arrondir. Avec votre vecteur de date, d:

require(zoo) 
d.ym <- as.yearmon(d) 
#Edit: 
#to get the number of months from the first record: 
round(as.numeric(d.ym -d.ym[1])*12) 

# calculate successive differences as integer months 
diffdates <- round(diff(as.numeric(d.ym))*12) 
diffdates 

L'effort impliqué dans l'apprentissage du zoo en vaut la peine. C'est un paquet mature. Il y a un grand nombre de questions sur R-help que Grothendieck répond par rapport aux fonctions du zoo.

Answer 3

Vous pouvez également le faire avec l'ensemble lubrifiant. Nous devrions d'abord les transformer en objets POSIXt. L'heure d'été créera des surprises pour nous (c'est-à-dire que les dates ne seront pas exactement séparées d'un mois), définissons donc le fuseau horaire sur Temps universel coordonné, qui ne participe pas à l'heure d'été.

d <- as.POSIXlt(d, tz = "UTC") 

Maintenant, nous soustrayons pour obtenir les intervalles de temps impliqués, et divisons par un mois.

(d - min(d))/months(1) 

Ceci nous donne le nombre de mois après la première date à laquelle chaque date se produit. (Je pense que c'est ce que vous recherchez, mais je m'excuse si j'ai mal interprété vos intentions). Ce type de division est dans la dernière version de lubridate, que j'ai soumis à cran ce matin. Cela peut donc prendre quelques jours pour être mis à jour. Je m'excuse.