[R] ??: Re:??: Re: Need help in waveslim package: imodwt and universal.thresh.modwt

rdporto1 rdporto1 at terra.com.br
Wed Nov 15 21:33:42 CET 2006


Airon, 

I don't think you have to find an English computer 'cause the
following must work in your Chinese one :-)

Let me explain. First of all, change your lines to

xdata <- ckhdat$Adj..Close[1:1447]
#names(ckhdwt.la8) <- c("w1", "w2", "w3", "w4", "w5","w6", "v6")

note the # sign, i.e., DO NOT change the names before
the function imodwt.

This is because the function imodwt looks for the names
created by the modwt function. If you need to change names,
do it AFTER the reconstruction.

I hope that it helps you.

Rogerio

---------- Cabeçalho original -----------

De: "Airon Yiu" aironyiu at yahoo.com.hk
Para: "rdporto1" rdporto1 at terra.com.br
Cópia: "r-help" r-help at stat.math.ethz.ch
Data: Wed, 15 Nov 2006 23:49:57 +0800 (CST)
Assunto: ??: Re:??: Re:[R] Need help in waveslim package: imodwt and  universal.thresh.modwt

> Hi Rogerio:
>    
>   Here is what I tried.  I have used only data points 1 to 1447 but the same problems appear if I used all the data in the file I sent .
>    
>   >library(waveslim)
> >infile <-"C:\\airon.csv"
> >ckhdat <- read.csv(infile,header=TRUE, quote="")
>   ...
> 1464       23-Dec-05,80.8,81,80.7,80.95,1538304,80.95
> 1465        28-Dec-05,80.85,81,80.3,80.7,3728116,80.7
> 1466      29-Dec-05,80.8,80.95,80.3,80.4,3145493,80.4
> 1467    30-Dec-05,80.4,80.2,78.85,79.65,7508611,79.65
> >xdata <- ckhdata$adjcls[1:1447]
>   > ckhdwt.la8 <- modwt(xdata, "la8",  n.levels = 6)
> > names(ckhdwt.la8) <- c("w1", "w2", "w3", "w4", "w5","w6", "v6")
> > 
> > ydata <- imodwt(ckhdwt.la8)
> > ydata
> numeric(0)
> > thre_wc_univ <- universal.thresh.modwt(ckhdwt.la8, max.level = 4, hard = FALSE)
> ???abs(wc.fine) : ??????????
> > 
> 
>   Note that the error message for universal.thresh.modwt  is in Chinese.  It roughly means "Error at abs(wc.find) : variable cannot have non-numeric value" 
>    
>   I need to find an WinXP (English) machine to see if the same issue appers.
>    
>   Thks



More information about the R-help mailing list