|
|
1.1 root 1: All the smoothers we shall describe can easily be {\em robustified} by replacing the averaging or least-squares operation by a more robust procedure.
2: known as {\em Hanning}, {em splitting} and {\em twicing}, in various combinations.
3: called {\em regression smoothers}.
4: Many of the smoothers that we discuss are {\em linear} (section~3.4.2), and this facilitates
5: an approximate assessment of their {\em degrees of freedom} (section~3.5).
6: Thus the weight in the fit at $x_0$ which is associated with the point $x_j$ is $S_{0j}$, and this sequence of weights is known as the {\em equivalent kernel} at $x_0$.
7: The {\em loess} smooth on the other hand has a strictly local neighbourhood yet the weights die down smoothly to zero.
8: itself; hence the name {\sl smoother}.
9: An important property of a smoother is its {\sl nonparametric} nature; that is,
10: We call the estimate produced by a smoother a {\sl smooth}.
11: setting, usually referred to as {\sl scatterplot
12: The simplest smoother occurs in the case of a {\sl categorical}
13: While the reader might not normally think of this as {\sl smoothing}, this simple
14: {\sl local averaging}, that is,
15: The averaging is done in {\sl neighbourhoods} around the target value.
16: question of which {\sl brand} of smoother to use, because smoothers
17: adjustable {\sl smoothing parameter}.
18: Thus there is a {\sl fundamental tradeoff between bias and variance},
19: In a sense the regression line is an {\sl infinitely smooth} function, and
20: are {\sl close} to $x_i$?
21: This is called a {\sl symmetric nearest neighbourhood} and
22: the {\sl running mean}
23: of which side they are on; this is called a {\sl nearest neighbourhood}.
24: also called a {\sl moving average}, and is popular for evenly-spaced time-series data.
25: in practice it does not work very well. It tends to be so wiggly that it hardly deserves the name {\sl smoother.}
26: The {\sl running
27: the use of a {\sl weighted} least-squares fit in each
28: ^{Cleveland's (1979)} implementation of a locally-weighted running-lines smoother, {\sl loess},
29: A kernel smoother uses an explicitly defined set of local weights, defined by the {\sl kernel}, to
30: function only of its {\sl metric} distance from $x_0$, while the weights used by the nearest-neighbour smoothers are typically a function of both {\sl metric} and {\sl rank} distance.
31: Their {\sl equivalent kernels} are one way to compare
32: All the smoothers studied in this chapter are {\sl linear} in $Y$, which means that the fit at a point $x_0$ can be written as $\gsmooth(x_0)=\sum_{j=1}^n S_{0j}y_j$, and the $S_{0j}$ depend on all the $x_i$ and on the smoothing parameter $\lambda$.
33: We do this using the {\sl equivalent degrees of freedom}, which we describe in the next chapter.
34: Regression splines offer a compromise by representing the fit as a {\sl piecewise}
35: of {\sl knots} or breakpoints, $\xi_1,\ldots,\xi_K$.
36: A variant of polynomial splines are the natural splines; although they are defined for all piecewise polynomials of odd degree, we discuss the natural {\sl cubic} splines.
37: A very simple approach (referred to as cardinal splines) requires a single parameter, the {\sl number} of interior knots.
38: In summary, regression splines are attractive because of their computational neatness, {\sl when the knots are given}.
39: that minimizer is a {\sl natural cubic spline} with knots at the unique values of $x_i$
40: This would result in $n+2$ parameters, although the constraints on each end bring it down to $n$. We'll see however that the coefficients are estimated in a constrained way as well, and this can bring the {\sl effective} dimension down dramatically.
41: Since the columns of $\bB$ are the evaluated $B$-splines, in order from left to right and evaluated at the {\sl sorted} values of $X$, and the cubic $B$-splines have local support, $\bB$ is lower 4-banded.
42: Let $\bN$ be an $n\times n$ nonsingular {\sl natural-spline} basis matrix for representing the solution (Exercise~2.5).
43: to use {\sl local averaging}.
44: One might say, then, that a cubic smoothing spline is approximately a {\sl kernel}
45: Here we describe in more detail the locally-weighted smoother of ^{Cleveland (1979)}, currently called {\sl loess} in the S statistical-computing language.
46: \item{(iii)}Weights $w_i$ are assigned to each point in $\NN(x_0)$, using the {\sl tri-cube} weight function:
47: as a percentage or {\sl span} of the data points, is the smoothing parameter.
48: Nearest neighbourhoods work satisfactorily with {\sl loess} at the endpoints, however,
49: The first two smoothers require a definition of {\sl nearest neighbours}
50: {\sl Nearest} is determined by a distance measure and for this
51: dimensions: the so-called {\sl thin-plate spline} is one such
52: Another generalization is known as multivariate {\sl tensor product} splines.
53: In addition, the emphasis in the book is not on smoothers {\sl per~se}, but
54: \exercise {\sl Updating formula for running-line smooth.}
55: \exercise {\sl Basis for natural splines.} Suppose $B$ is an $n\times (K+4)$ matrix containing the evaluations of the cubic $B$-spline basis functions with $K$ interior knots evaluated at the $n$ values
56: \exercise {\sl Derivation of smoothing splines; ^{Reinsch (1967)}.} Consider the following optimization problem: minimize
57: \exercise {\sl Semi-parametric regression; ^^{Green, P.J.}^^{Jennison, C. }^^{Seheult, A.} Green, Jennison, and Seheult (1985).} Suppose we have a set of $n$ observations of $p$ predictors arranged in
58: Construct an appropriate penalized residual sum of squares, and show that the minimizers must satisfy the following pair of {\sl estimating} equations:
59: \exercise {\sl Efficient kernel smoothing; ^{Silverman (1982)}, ^{H\"ardle (1986)}}
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.