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.
