similar to: Wavelet reconstruction with coefficients set to zero

Airon, I used R2.4.0 on a Windows XP (SP2) (not Chinese :-)) and it still works: > data = read.csv("u:/airon.csv") > xdata = data$Adj..Close > modwt.la8 = modwt(xdata, "la8", n.level=6) > summary(modwt.la8) Length Class Mode d1 1467 -none- numeric d2 1467 -none- numeric d3 1467 -none- numeric d4 1467 -none- numeric d5 1467 -none- numeric d6 1467

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

Hi: I have encountered problems with imodwt and universal.thresh.modwt and cannot find any reference in R Search. Hope someone can give me some ideas: Starting with modwt.la8 <- modwt(xdata, "la8", n.level=6) <-- this seems to work fine (1) ydata <- imodwt(modwt.la8) will always give ydata as numeric(0) (no values) instead of being a time series data with

Dear R Users, For running wavelet functions using dwt( ), modwt( ), and mra( ), a wavelet filter algorithm is applied. For all these functions, default is "la8" and other possibility is "haar". In related documents, another possibilities like as symlet and coiflet ... are not cited. Besides "la8" and "haar", which wavelet filters can be used? Thank

Hi there, When I run the code below I get the error Error in xy.coords(x, NULL, log = log) :(list) object cannot be coerced to type 'double' Any tips how I can resolve this? > > library("waveslim") > > vols=read.csv(file="C:/Users/ocuk/My Documents/Abs Vol.csv", header=TRUE, > sep=",") > x<-c(vols[,1]) > #x > #data(ibm) >

Hai I'm Dee. I'm trying to write var data from these codes inside excel file. My directory to store the data is *D:\FYP\image* . these are my codes, can you help give an advice or idea with my problem: l*ibrary("biOps") library("waveslim") library("xlsReadWrite") x <- readTiff("D:\\FYP\\image\\SignatureImage\\user186g1.tif") y <-

I take the liberty of pasting the on-line documentation for package "wavelets", functions "dwt" and "wt.filter", pertinent to the wavelet/filter choice. While I apologize for my verbose messages with huge attachments, I encourage whoever loads packages to CRAN to take some minutes to reread the attached documentation. Usage dwt(X, filter="la8", n.levels,

Hai, I'm trying to write these var output data from these codes inside excel file. My directory to store the data is /D:\FYP\image / but receive an error message : /Error in write.xls(mydata, "D:\\FYP\\image.mydata.xls") : object 'mydata' not found/ these are my codes, can you help give an advice or idea with my problem: /library("biOps")

Hello, I get the following message when I try to load a time series of 128 stock returns. Error in modwt(stock, wf = "la8", n.levels = 8, boundary = "periodic") : wavelet transform exceeds sample size in modwt I thought that MODWT had no restrictions on the sample size? Any help? -- View this message in context:

Hi everyone I am using the wavelet and waveslim package in R to find the wave variance.Here is the code return.modwt<-modwt(X, filter="la8", n.levels=5, boundary="periodic", fast=TRUE) return.modwt.var<- wave.variance(return.modwt, type="nongaussian") Where X is a uni-variate time series. I am expecting a matrix with 5 rows(no. of levels) and 3

Hello, everyone I have a question. Assume I have the following zoo object: me.la <- structure(c(1524.75, 1554.5, 1532.25, 1587.5, 1575.25, 1535.5, 1550, 1493.5, 1492.5, 1472.25, 1457.5, 1442.75, 1399, 1535.75, 1565.25, 1543.5, 1598.5, 1586.5, 1547, 1561.5, 1504.75, 1503.75, 1483.75, 1468.75, 1453.75, 1410, 1546.75, 1575.25, 1554, 1609, 1597.5, 1558.5, 1573, 1516.25, 1515.5, 1495, 1480, 1465,

Airon, I'm not sure I have understood your question. Anyway, it seems to me that you'll have to code a little. Take a look at ?modwt. I think you'll have to put zeros each scale a time except the one you want to reconstruct. Then do imodwt. After some small loop, I think you'll get the desired result. HTH, Rogerio. ---------- Cabe?alho original ----------- De: "Airon

Given a time series of length N, I am trying to figure out its maximum DWT decomposition level. With reference to "dwt" function of R package "wavelets", running the provided example I get the following and wonder how the maximum decomposition level (which is not an integer number) is calculated. I tried myself the formula in the on-line documentation but could not get the

I would like to create a vector of pairlist (flag, binary_value) like: (variable ="TrendOff", value = 0) (variable ="MOdwt", value = 1) (variable ? "ZeroPadding", value =1) ................................................ I tried the following syntax but the emcompassing list (that I called "flags") is not made up of pairlists. Instead it

I realized that I may not have answered the question you were asking and that no one else has responded. I can across a similar problem and may have an answer to your question now. If you have both the wavelet coefficients and the scaling coefficients then create a fake sequence of the same length as the original and decompose that sequence using wd form wavethersh with the same wavelet family

Hello, does someone knows how can I get the ordenates of a father wavelet using the package wavethresh? And using other package? Any hints? For example, I'd like knowing what's the value of \phi{j,k}(t), that is, the value of a father wavelet (from some family) at scale j, shift k, point t. Thanks in advance, Rogerio.

hi everyone . I am trying to use some packages but there are some functions have the same name in different package. for example dwt function both in packages wavelets and waveslim dwt(X, filter="la8", n.levels, boundary="periodic") How can I avoid mixing them up ? -- TANG Jie Email: totangjie@gmail.com Tel: 0086-2154896104 Shanghai Typhoon Institute,China

Hi I have written a small introduction to newbies in Danish on ssh and friends. Now some people are questioning my advice and I think they have a point. I am advocating people to use DSA-keys and a config file with this: Protocol 2 ForwardAgent yes ForwardX11 yes Compression yes CompressionLevel 9 and running ssh-agent and ssh-add, and then loggin in without giving keys. One

My wife and I have been working on our own task-listing, time-tracking web application, and we''re really interested in some feedback as we ourselves feel we''ve come a long way. If you have some minutes to kill (or just want to see how the Comet functionality provided by the Juggernaut plugin can work) have a look at http://www.clockingit.com and

> wt.filter("d8") #### HOW COME ???? Error in validObject(.Object) : invalid class "wt.filter" object: invalid object for slot "transform" in class "wt.filter": got class "function", should be or extend class "character" > wt.filter("d10") # OK An object of