I am trying to plot the cross-wavelet power spectrum of these two time series :
date <- lubridate::ymd(yrar, truncated = 2L) my.data <- data.frame(date,