### 2.Example.imputation.R
set.seed(1)
cv.index <- cvtype(n=1024, cv.bsize=2, cv.kfold=8, cv.random=T)$cv.index 

snr <- 5
testdata <- heav(1024)
x <- testdata$x
y <- testdata$meanf + rnorm(n=1024, 0, testdata$sdf / snr)

yimpute <- cvimpute.by.wavelet(y=y, impute.index=cv.index)$yimpute
par(mar=0.1+c(4,4,2,1))
plot(y, yimpute, xlab="Observations", ylab="Imputed Values", 
                         main="Heavisine", cex=0.5);abline(0,1)