1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
|
\name{TOthreshda2}
\alias{TOthreshda2}
\title{Data analytic wavelet thresholding routine}
\usage{
TOthreshda2(ywd, alpha = 0.05, verbose = FALSE, return.threshold = FALSE)
}
\arguments{
\item{ywd}{The \code{\link{wd.object}} that you wish to threshold.}
\item{alpha}{The smoothing parameter which is a p-value }
\item{verbose}{Whether messages get printed}
\item{return.threshold}{If TRUE then the threshold value gets returned
rather than the actual thresholded object}
}
\description{
This function might be better called using the regular
\code{\link{threshold}} function using the \code{op2} policy.
Corresponds to the wavelet thresholding routine developed by
Ogden and Parzen (1994) Data dependent wavelet thresholding in nonparametric
regression with change-point applications. \emph{Tech Rep 176},
University of South Carolina, Department of Statistics.
}
\details{
The TOthreshda2 method operates in a similar fashion to
\code{\link{TOthreshda1}} except that it takes the cumulative sum
of squared coefficients, creating a sample "Brownian bridge" process,
and then using the standard Kolmogorov-Smirnov statistic in testing.
In this situation, the level of the hypothesis tests, alpha, has default
value 0.05. Note that the choice of alpha controls the smoothness of
the resulting wavelet estimator -- in general, a relatively large alpha
makes it easier to include coefficients, resulting in a more wiggly
estimate; a smaller alpha will make it more difficult to include
coefficients, yielding smoother estimates.
}
\value{
Returns the threshold value if \code{return.threshold==TRUE} otherwise
returns the shrunk set of wavelet coefficients.
}
\seealso{\code{\link{threshold}},\code{\link{TOthreshda1}}, \code{\link{wd}}}
\author{Todd Ogden}
\keyword{smooth}
|
