Description Usage Arguments Details Value Note Author(s) References See Also Examples

Starting with a noise-contaminated signal, we decompose it using a 'nondecimated' lifting algorithm (i.e. by applying the modified lifting transform following several random paths), shrink all the obtained detail coefficients and invert each transform to give an estimated signal. The average of all these estimates is the final proposal for estimating the true (unknown) signal.

1 2 |

`x` |
Vector of any length (possibly irregularly spaced) that gives the grid locations where the signal is observed. |

`f` |
Vector of the same length as |

`J` |
Number of trajectories to be used by the nondecimated lifting algorithm. |

`Pred` |
The type of regression to be used in the prediction step of the modified lifting algorithm. Choices are linear, quadratic or cubic (respectively, |

`neighbours` |
Number of neighbours to be used for defining the neighbourhood of each point that has to be removed. If ( |

`closest` |
If ( |

`intercept` |
Specifies whether ( |

`nkeep` |
Number of scaling points we want at the end of the application of the transform. The usual choice is |

`trule` |
The type of Bayesian shrinkage technique, with possible choices posterior median ( |

`verbose` |
A boolean indicating whether extra information should be printed. |

`do.orig` |
A boolean indicating whether the original |

`returnall` |
A boolean indicating whether the function returns useful variables or just the denoised datapoints. |

Essentially, this function applies `J`

times the modified lifting algorithm (that can be found in `fwtnpperm`

), and removes the noise from all sets of detail coefficients by using empirical Bayes shrinkage (of package 'EbayesThresh'). Inverting (by means of the function `invtnp`

of the package 'adlift') each transform consequently results in `J`

estimates of the (unknown) true signal. The average of these estimators is our proposed estimator. The functions that appear as choices for `Pred`

can be found in the package 'adlift'.

`vec` |
Matrix whose rows give the trajectories to be used by the nondecimated lifting algorithm. |

`ghatnat` |
Vector that gives the estimated true signal given by denoising using the lifting scheme that establishes its own order for removing the points (but with the same specification for prediction stage and neighbourhood as the modified algorithm), rather than a randomly generated order. |

`aveghat` |
Estimated signal, obtained as the average of the individual estimates from the random trajectory runs. |

Use this function together with the "adlift" and "EbayesThresh" packages available from CRAN.

Marina Knight (marina.knight@bristol.ac.uk)

See the paper 'A "nondecimated" lifting transform.' by Knight, M.I. and Nason, G.P. (2009) for further details.

`denoiseperm`

, `fwtnpperm`

, `fwtnpperm`

, and also `invtnp`

of package 'adlift'

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 | ```
# construct the grid
x<-runif(256)
# construct the true, normally unknown, signal
g<-make.signal2("blocks",x=x)
# generate noise with mean 0 and signal-to-noise ratio 5
noise<-rnorm(256,mean=0,sd=sqrt(var(g))/5)
# generate a noisy version of g
f<-g+noise
# decide on a number of random trajectories to be used (below J=100, in paper J=20,30), and apply
# the nondecimated lifting transform to the noisy signal (x,f)
#
# below we apply the modified lifting transform J times, each time following a different path,
# and using adaptive prediction with neighbourhoods of size 2 in closest configuration;
# all details are then thresholded using posterior medians and the algorithms inverted
# the aggregate estimator of g proposed by our method is found in out$aveghat
out<-nlt(x,f,J=10,Pred=AdaptPred,neighbours=2,closest=TRUE,intercept=TRUE,nkeep=2,trule="median")
``` |

```
Loading required package: EbayesThresh
Loading required package: adlift
**********************************************
adlift: a package to perform wavelet lifting schemes
--- Written by Matt Nunes and Marina Knight ---
Current package version: 1.3-3 ( 2017-09-13 )
-+ packaged by MAN +-
**********************************************
adlift 1.3-3 loaded
Attaching package: 'adlift'
The following object is masked from 'package:EbayesThresh':
postmean.cauchy
**********************************************
nlt: a package to perform nondecimated wavelet lifting
--- Written by Marina Knight and Matt Nunes ---
Current package version: 2.1-3 ( 11/11/2012 )
-+ packaged by MAN +-
**********************************************
nlt 2.1-3 loaded
1 ...
2 ...
3 ...
4 ...
5 ...
6 ...
7 ...
8 ...
9 ...
10 ...
```

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.