On Sep 29, 2011, at 10:10 AM, maxbre wrote:
hello all,
this is my reproducible example data frame
test<-structure(list(date = structure(c(1262300400, 1262304000,
1262304000,
1262307600, 1262307600, 1262311200, 1262311200, 1262314800,
1262314800,
1262318400, 1262318400, 1262322000, 1262322000, 1262325600,
1262325600,
1262329200, 1262329200, 1262332800, 1262332800, 1262336400,
1262336400,
1262340000, 1262340000, 1262343600, 1262343600, 1262347200,
1262347200,
1262350800, 1262350800, 1262354400, 1262354400, 1262358000,
1262358000,
1262361600, 1262361600, 1262365200, 1262365200, 1262368800,
1262368800,
1262372400, 1262372400, 1262376000, 1262376000, 1262379600,
1262379600,
1262383200, 1262383200, 1262386800), class = c("POSIXct", "POSIXt"
), tzone = "")), .Names = "date", row.names = c(NA, -48L), class =
"data.frame")
given that, my objective is to create a new variable accounting for
the
julian day (i.e. the number of consecutive days elapsed from a fixed
origin
date)
and these are (among many others) my (most) unsuccessful attempts…
library(chron)
#not working!
test$jday<-julian(test$date, origin=as.POSIXct(test$date))
I thought the error message was reasonably informative:
Error in julian.POSIXt(test$date, origin = as.POSIXct(test$date)) :
'origin' must be of length one
And the obvious fix was to just give it one starting date:
> test$jday<-julian(test$date, origin=as.POSIXct(test$date)[1])
> test$jday
Time differences in days
[1] 0.00000000 0.04166667 0.04166667 0.08333333 0.08333333
0.12500000 0.12500000 0.16666667 0.16666667 0.20833333
[11] 0.20833333 0.25000000 0.25000000 0.29166667 0.29166667 0.33333333
0.33333333 0.37500000 0.37500000 0.41666667
[21] 0.41666667 0.45833333 0.45833333 0.50000000 0.50000000 0.54166667
0.54166667 0.58333333 0.58333333 0.62500000
[31] 0.62500000 0.66666667 0.66666667 0.70833333 0.70833333 0.75000000
0.75000000 0.79166667 0.79166667 0.83333333
[41] 0.83333333 0.87500000 0.87500000 0.91666667 0.91666667 0.95833333
0.95833333 1.00000000
attr(,"tzone")
[1] ""
attr(,"origin")
[1] "2009-12-31 18:00:00 EST"
#not working!
test$jday<-julian(test$date, origin=test$date)
#not working!
test$jday <-julian(as.Date(test$date, "%m"),
as.Date(test$date, "%d"),
as.Date(test$date, "%Y"),
origin=c(month=1, day=0, year=2010))
Exactly the same errors.
any help for this?
Yes:
Read the error messages for meaning and include them in your postings.
--
David Winsemius, MD
West Hartford, CT
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.