Convertir une date vecteur en jour Julien dans la R

J'ai une colonne de dates dans le format:
16Jun10
et je voudrais extraire le jour Julien.
J'ai plusieurs années.

J'ai essayé les fonctions de julian et mdy.date et il ne semble pas fonctionner.

OriginalL'auteur user3166363 | 2014-01-28

21

Essayez ce qui suit pour convertir de la classe character(c'est à dire le texte) à la classe POSIXlt, puis l'extrait de jour Julien (yday):
```
tmp <- as.POSIXlt("16Jun10", format = "%d%b%y")
tmp$yday
# [1] 166
```
Pour plus de détails sur les paramètres de la fonction:
```
?POSIXlt
?DateTimeClasses
```
Une autre option est d'utiliser un Date classe, et ensuite utiliser format pour extraire un jour julien (notez que cette classe de définir les jours de julien entre 1:366, tandis que POSIXlt est 0:365):
```
tmp <- as.Date("16Jun10", format = "%d%b%y")
format(tmp, "%j")
# [1] "167"
```
Cette réponse sommet parle de Jour de l'Année/Ordinale date et pas Julian Date. Les dates juliennes sont jour compte de 4713 avant jc (!!). Pourquoi est-ce répondre à de nombreuses reprises upvoted?.
Voir Terminologie paragraphe - "Le terme Julian date peut aussi se référer, en dehors de l'astronomie, de la journée de l'exercice nombre (plus correctement, l'ordinal date) dans le calendrier Grégorien, en particulier dans la programmation d'ordinateur...". Je suppose que les programmeurs ont voté plus que les astronomes ont diminué voté.
Une bonne :). Note Andrei réponse ci-dessous. Il fait calcule liés à l'astronomie les dates Juliennes. Notez également que la balise julian se réfère à la liés à l'astronomie sens. Et maintenant, je me rends..

OriginalL'auteur Marc in the box
17

De la même façon:
```
require(lubridate)
x = as.Date('2010-06-10')
yday(x)

[1] 161
```
Également noter que l'utilisation de lubridate:
```
> dmy('16Jun10')
[1] "2010-06-16 UTC"
```
Pourriez-vous m'aider à obtenir une date Julienne comme 112161 pour 2012/06/09 (aaaa/mm/jj)?

OriginalL'auteur John
5

Vous pouvez utiliser R insol paquet qui a un JD(x, inverse=FALSE) fonction qui convertit POSIXct Julian Nombre de Jours (JDN).

insol paquet a également JDymd(year,month,day,hour=12,minute=0,sec=0) personnalisé pour les dates.

Pour afficher la totalité de la Date Julienne (JD) vous devez définir options(digits=16).

OriginalL'auteur Andrei

my.data = read.table(text = "
     OBS  MONTH1  DAY1  YEAR1
       1       3     1   2012
       2       3    31   2012
       3       4     1   2012
       4       4    30   2012
       5       5     1   2012
       6       5    31   2012
       7       6     1   2012
       8       6    30   2012
       9       7     1   2012
      10       7    31   2012    
", header = TRUE, stringsAsFactors = FALSE)

my.data$MY.DATE1 <- do.call(paste, list(my.data$MONTH1, my.data$DAY1, my.data$YEAR1))
my.data$MY.DATE1 <- as.Date(my.data$MY.DATE1, format=c("%m %d %Y"))

my.data$my.julian.date <- as.numeric(format(my.data$MY.DATE1, "%j"))
my.data

Retourne, qui, techniquement, est incorrecte car les dates Juliennes de ne pas retourner à 1 le premier jour de chaque mois de janvier:

http://en.wikipedia.org/wiki/Julian_day

Ci-dessous les dates sont Ordinale dates:

   OBS MONTH1 DAY1 YEAR1   MY.DATE1 my.julian.date
1    1      3    1  2012 2012-03-01             61
2    2      3   31  2012 2012-03-31             91
3    3      4    1  2012 2012-04-01             92
4    4      4   30  2012 2012-04-30            121
5    5      5    1  2012 2012-05-01            122
6    6      5   31  2012 2012-05-31            152
7    7      6    1  2012 2012-06-01            153
8    8      6   30  2012 2012-06-30            182
9    9      7    1  2012 2012-07-01            183
10  10      7   31  2012 2012-07-31            213

OriginalL'auteur Mark Miller

Voici mes R versions de code écrit à l'origine en APL et converti à J. Nous appelons cette pseudo-Julien, car il est conçu uniquement pour les dates après le 15 octobre 1582, qui est quand la réforme du calendrier, dans certaines parties du monde Occidental, arbitrairement changé la date.

#* toJulian: convert 3-element c(Y,M,D) timestamp into pseudo-Julian day number.
toJulian<- function(TS3)
{   mm<- TS3[2]
    xx<- 0
    if( mm<=2) {xx<- 1}
    mm<- (12*xx)+mm
    yy<- TS3[1]-xx
    nc<- floor(0.01*yy)
    jd<- floor(365.25*yy)+floor(30.6001*(1+mm))+TS3[3]+1720995+(2-(nc-floor(0.25*nc)))
    return(jd)
#EG toJulian c(1959,5,24) -> 2436713
#EG toJulian c(1992,12,16) -> 2448973
}

Voici la fonction inverse:

#* toGregorian: convert pseudo-Julian day number to timestamp in form c(Y,M,D)
# (>15 Oct 1582).  Adapted from "Numerical Recipes in C" by Press,
# Teukolsky, et al.
toGregorian<- function(jdn)
{   igreg<- 2299161       # Gregorian calendar conversion day c(1582,10,15).
    ja<- floor(jdn)
    xx<- 0
    if(igreg<=ja){xx<- 1}
    jalpha<- floor((floor((xx*ja)-1867216)-0.25)/36524.25)
    ja<- ((1-xx)*ja) + ((xx*ja)+1+jalpha-floor(0.25*jalpha))
    jb<- ja+1524
    jc<- floor(6680+((jb-2439870)-122.1)/365.25)
    jd<- floor(365.25*jc)
    je<- floor((jb-jd)/30.6001)
    id<- floor((jb-jd)-floor(30.6001*je))
    mm<- floor(je-1)
    if(12<mm){mm<- mm-12}
    iyyy<- floor(jc-4715)
    if(mm>2){iyyy<- iyyy-1}
    if(0>iyyy){iyyy<- iyyy-1}
    gd<- c(iyyy, mm, id)
    return(gd)
#EG toGregorian 2436713 -> c(1959,5,24)
#EG toGregorian 2448973 -> c(1992,12,16)
}

OriginalL'auteur DevonMcC

Vous devez vous connecter pour publier un commentaire.