setwd("C:/Users/Sérgio/Documents/FisicaAmbiental/Wavelet")

require(WaveletComp)


x <- read.csv('Tsolo.txt',header = TRUE)

print(nrow(x))


my.data <- data.frame(x = x)
my.w <- analyze.wavelet(my.data,
	loess.span = 0,
	dt = 1, dj = 1/250,
	lowerPeriod = 1,
	upperPeriod = 256,
	make.pval = TRUE, n.sim = 10)
wt.image(my.w, n.levels = 250,
	legend.params = list(lab = "wavelet power levels") )