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

Question

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

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

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

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

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

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

Answer 1

После файла справки ?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)