Hello, dear friend, you can consult us at any time if you have any questions, add WeChat: daixieit

496 Homework 2

Due Tuesday, Feb. 14

Guidelines to Turn the HW in:

(b) The solution to each problem must contain a summary of the main conclusions or final answers, handwritten or typed using either word or R Markdown (or a combination of these).

(b) If you use R to obtain the solutions, explicitly refer to the R output. You also need to include the R code as a separate image or included in the word or Markdown le.

Problems.

1. Let X, X1 , . . . , Xn  f and dene

γ := ╱尸 f (x)fˆn (x)dx= f (x) fˆn (x)dx.

Since  fˆn (x)= (K(x  X)) =: fsm (x), we have γ =    f (x)fsm (x)dx.  In

class, we claimed that

  i1  fˆ(i)(Xi )\ = γ,

where fˆ(i)(x) is the leave-one-out kernel density estimator. Show the validity of this claim. Hint: First, compute  fˆ(i)(Xi )|Xi.

2.  Get the data on fragments of glass collected in forensic work from canvas.

(a) Estimate the density of the rst variable (refractive index) using a his- togram and a kernel density estimator. Experiment with 2 or 3 different binwidths and bandwidths.

(b) Now, use cross-validation to choose the amount of smoothing.  Add the

condence intervals fˆn (x) z (x)  and generate a graph sim-

ilar to the rst graph in Section 2.5 of Garc´ıa-Portgu´es’ notes.

3. Exercise 2.13 in A Short Course on Nonparametric Curve Estimation” by Garc´ıa-Portugu´es.

Remark:  Only consider the bimodal density  (see Section 2.4 in the notes), and the bandwidth selectors RT (Rule of Thumb), UCV (least-squares cross- validation method as seen in class), and DPI (the Sheather-Jones method).

4. In class, we computed the coverage probabilities of the CIs fˆn (x)z (x) when used to estimate fsm (x) := 匝K(x  X). This was done when the band-

width h is set to a fixed value. Carry out the same analysis when the bandwidth h is chosen according to (a) the rule of thumb method, and according to (b) cross-validation.  Generate the same graph as in class for these two cases (see graph above Exercise 2.3 in Section 2.5 of Garc´ıa-Portugu´es notes).

5. In class, we computed the coverage probabilities of the CI

fˆn (x) zα/2  fˆn (x)   K2 (u)du

when applied to estimate both fsm (x) = 匝(Kh (x  X) and f (x).  We also proposed the following unbiased” estimator:

f˜n (x; h, ζ) = fˆn (x; h)          (fˆn (x; ζh) n (x; h))2             

for some xed ζ > 1.

(a)  Compute the coverage probabilities of the following estimator

f˜n (x; h, ζ) zα/2  fˆn (x)   K2 (u)du

when applied to estimate f (x), and generate a graph similar to that in Section 2.5 of Garc´ıa-Portugu´es (see graph above Exercise 2.3 therein). Choose h according to the RT method and experiment with a few values of ζ and comment on the results.

(b) Repeat (a) with the CI

f˜n (x; h, ζ) zα/2  n (x; h, ζ)   K2 (u)du