[R] Morlet wavelet analysis

Irina Foss Irina.Foss at thurso.uhi.ac.uk
Thu Jun 4 15:49:45 CEST 2009


Dear,

I am using "cwt "function from Rwave package to perform Morlet wavelet analysis.

d1<-c(1.31673209591515, -0.171333455797346, -1.67618546079420, 
-0.931604651010362, -0.954614183224057, -1.19821794547758, 0.516096363353144, 
-0.0432811977638559, 0.737764943619919, 0.438046629673177, -0.208607167908743, 
-0.3091308821262, -1.42473112931594, 0.234125312118165, -0.307047554490597, 
0.686201101422384, 1.40104418156371, -0.722985280679127, -1.47295298153946, 
-1.49070808616222, -0.294899345331949, -0.623882988251866, 0.95580872220543, 
0.159226599747532, 0.258592921036826, -0.893605825950617, -1.48090318423471, 
-0.599454337019121, -0.315308442503905, 1.48221661623468, 1.07313149452230, 
0.0330503877044541, 1.07332650281700, -1.37121080562242, -0.181568897416963, 
1.55169357402125, -0.502404116831855, 2.19007041840759, 0.177374817139349, 
1.37270541339308, 0.179612250662539, 1.33781757389541)

d1<-ts(d1deltat=1,start=c(1960,1))
 library(Rwave)
D1-cwt(d1,noctave=8,nvoice=1)

After I did Morlet wavelet analysis by using function "wavelet" from dplR package.

 library(dplR)

D2<-wavelet(d1,Dt=1)

But the wavelets from two tests are very different. What is the reason??

Thank you,

Irina




Irina Foss
Environmental Research Institute
North Highland College
UHI Millennium Institute
Castle Street
Thurso, Caithness
Scotland
KW14 7JD
United Kingdom

Tel:  +44 (0) 1847 889 587
Fax: +44 (0) 1847 890 014
http://www.eri.ac.uk/sub/staff.php?id=11




More information about the R-help mailing list