rm(list=ls())

source("wavelet_function_r.txt")

data=matrix(scan("LF.txt"), ncol=2, byrow=T) #data to be analyzed should be 2 cols - year in 1st and values in 2nd


lvl=.95
lvl2=.9

XL=as.vector(data[,2])


fin=wavelet(data[,2], 1, 1, .25, 6)

source("ken_sig.txt") #significance function

zz=noise(lvl2, XL,"w",  fin, length(data[,2]))
xx=noise(lvl,XL, "w",  fin, length(data[,2]))

source("wave_plot.txt")