пакет вейвлета Inverse DWT не может восстановить серию?

Я использую пакет вейвлетов и заметил, что когда я пытаться

library("wavelets")
x <- rnorm(100)
y <- idwt(dwt(x))
plot(x, y)

реконструкция y, по-видимому, не равна оригиналу x.

Этого следует ожидать?

Для некоторого контекста я пытаюсь выполнить (регуляризованную) логистическую регрессию, используя вейвлет-преобразования набора рядов. Затем я хочу отобразить коэффициент регрессии обратно в исходное пространство временных рядов, чтобы увидеть, какие временные точки использовались при различении.

Но я даже не могу восстановить оригинальную серию. Я может быть совершенно неправильно понимаю вещи, может ли кто-нибудь пролить свет на это?


r wavelet
person Krisrs1128    schedule 03.09.2015    source источник
comment
Это работает, например, с y <- idwt(dwt(x, n.levels=3, boundary="reflection", fast=FALSE)) (см. ?dwt).   -  person    schedule 03.09.2015

Ответы (1)


arrow_upward
4
arrow_downward

Следуя файлу справки ?dwt, вы можете изменить свой сценарий, например:

library(wavelets)
set.seed(42)
x <- rnorm(100)
y <- idwt(dwt(x, n.levels=3, boundary="reflection", fast=FALSE))
plot(x, y)
abline(0,1)

введите здесь описание изображения

person Community    schedule 03.09.2015