my.date <- seq(as.POSIXct("2014-10-14 00:00:00", format="%F %T"),
by="hour",
length.out=6*128*24)
library('WaveletComp')
x <- periodic.series(start.period=1*24, length=6*128*24)
my.data <- data.frame(date=my.date, x=x)
my.wt <- analyze.wavelet(my.data, "x",
loess.span=0,
dt=1/24, dj=1/20,
lowerPeriod=1/4,
make.pval=TRUE, n.sim=10)
layaout2 <- layout(matrix(c(1, 2, 3), nrow=3))
wt.image(my.wt, main="wavelet power spectrum",
legend.params=list(lab="wavelet power levels (quantiles)",
lab.line=3.5,
label.digits=2),
periodlab="period (days)")
p1 <- recordPlot()
wt.image(my.wt, main="wavelet power spectrum",
legend.params=list(lab="wavelet power levels (quantiles)",
lab.line=3.5,
label.digits=2),
periodlab="period (days)")
p2 <- recordPlot()
plot_grid(p1, p2,
labels='AUTO',
hjust=0, vjust=1)
我弄错了
Error in if (labels) { : missing value where TRUE / FALSE is required