[R] Morlet wavelet analysis

[Irina Foss]

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