Your code is nearly unreadable, probably because you posted in HTML instead of
plain text. Plain text is the standard expected on R-help.
The only apparent date/time conversion in your code is the use of as.POSIXct().
You are missing a closing quote mark after the date.
If I fix that, then on my system as.POSIXct() succeeds:
> seq(as.POSIXct("2002-01-01"), by = "hour", length.out = 5)
[1] "2002-01-01 00:00:00 PST" "2002-01-01 01:00:00 PST"
[3] "2002-01-01 02:00:00 PST" "2002-01-01 03:00:00 PST"
[5] "2002-01-01 04:00:00 PST"
Your warning (not error) messages refer to an unknown timezone. So check your
timezone, for example:
> Sys.getenv('TZ')
[1] "US/Pacific"
I hope this helps.
-Don
On 12/13/16, 9:33 AM, "R-help on behalf of Muhammad Anees"
<r-help-boun...@r-project.org on behalf of a...@aneconomist.com> wrote:
Dear All,
I tried to find an answer everywhere including stackoverflows and manuals
but there is no specific solution that might have helped. Hence, I am sure
you will bear with me for this first query.
I am trying to replicate production of Wavelet Coherency using the
following codes: I dont know why the two errors/warnings appear and my
resultant figures dont show up the time as I use the format of Date as
2000-01-01
Code used is:
> library(WaveletComp)> x<- API> y<- TFS> my.data = data.frame(x = x, y =
y)> my.wc = analyze.coherency(my.data, my.pair = c("x","y"),+
loess.span = 0,+ dt = 1/24, dj = 1/100,+
lowerPeriod = 1/2,+ make.pval =
T, n.sim = 10)> wc.image(my.wc, n.levels = 250,+ siglvl.contour = 0.1,
siglvl.arrow = 0.05,+ legend.params = list(lab = "cross-wavelet power
levels"))> wc.image(my.wc, n.levels = 250, color.key = "interval",+
siglvl.contour = 0.1, siglvl.arrow = 0.05, which.arrow.sig = "wt",+
legend.params = list(lab = "cross-wavelet power levels"))> wc.image(my.wc,
which.image = "wc", n.levels = 250,+ siglvl.contour = 0.1,
siglvl.arrow = 0.05,+ legend.params = list(lab = "wavelet coherence
levels"))> my.date = seq(as.POSIXct("2002-01-01),+ by = "hour",
length.out = 24*96)########################################
########################################
Warning messages:1: In strptime(xx, f <- "%Y-%m-%d %H:%M:%OS", tz = tz) :
unknown timezone '%F %T'2: In as.POSIXct.POSIXlt(x) : unknown
timezone '%F %T'3: In strptime(x, f, tz = tz) : unknown timezone '%F
%T'4: In as.POSIXct.POSIXlt(as.POSIXlt(x, tz, ...), tz, ...) :
unknown timezone '%F %T'
########################################
########################################
I would be happier to see a suggestion for removing this warning and
get the results as the time axis showing Date as in the data on the
figures of my output from above.
Regards
------------------------------
Muhammad Anees
Assistant Professor,
The Econometrician &
Statistical Consultant @
www.aneconomist.com
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
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.
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
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.
