Homework Assignment 1 - Statistical Learning | STAT 542, Assignments of Statistics

Material Type: Assignment; Professor: Liang; Class: Statistical Learning; Subject: Statistics; University: University of Illinois - Urbana-Champaign; Term: Fall 2008;

Typology: Assignments

Pre 2010

Uploaded on 03/11/2009

koofers-user-l3j
koofers-user-l3j 🇺🇸

9 documents

1 / 2

Toggle sidebar

This page cannot be seen from the preview

Don't miss anything!

bg1
STAT 542 Fall 2008
Homework Assignment 1
Due Tuesday, September 16
1. Assume a linear regression model with Normal errors. Take the error variance σ2known.
Show that the AIC criterion and Mallow’s Cpare equivalent. Recall that
AIC = log likelihood evaluated at MLE+dim of the model
Cp= residual sum of squares + 2σ2dim of the model
2. In a simple linear regression problem, we observe npairs of observations (y1,x1),...,(yn, xn),
where yiRand xiRp. The least squares estimate of the p-dimensional regression
coefficient is given by
ˆ
β= arg minβRp
n
X
j=1
(yjxt
jβ)2.
A basic tool for examining the fit is the leave-one-out residual: consider fitting the model
omitting the ith observation with the corresponding LS estimate defined as
ˆ
β[i]= arg min
n
X
j:j6=i
(yjxt
jβ)2,
get a prediction for the omitted observation, ˆy[i]=xt
iˆ
β[i], then the leave-one-out residual
for the ith observation is (yiˆy[i]).
We will show that
yiˆy[i]=yixt
iˆ
β
1Hii
,
where Hii denotes the (i, i)th entry of the projection matrix H. That is, the leave-one-
out residual is a re-scale of the original residual, hence it is not necessary to re-fit the
model each time an observation is omitted.
We will prove this result via the following steps. Throughout, we assume all the design
matrices are of ful l rank, and all the LS estimates exist and are unique.
(a) Consider a new data set of nobservations where y
1= ˆy[1] =xt
1ˆ
β[1] and the others
are the same as the original data. Let ˆ
βdenote the corresponding LS estimate
based on the new data. Show that
ˆ
β=ˆ
β[1].
Hint: Show that ˆ
β[1] minimizes (y
1xt
1β)2+Pn
j=2(yjxt
jβ)2.
(b) Show that
ˆy[1] =Pn
j=2 H1jyj
1H11
.
(c) Show that
y1ˆy[1] =y1xt
1ˆ
β
1H11
.
pf2

Partial preview of the text

Download Homework Assignment 1 - Statistical Learning | STAT 542 and more Assignments Statistics in PDF only on Docsity!

STAT 542 Fall 2008

Homework Assignment 1

Due Tuesday, September 16

  1. Assume a linear regression model with Normal errors. Take the error variance σ

2 known.

Show that the AIC criterion and Mallow’s Cp are equivalent. Recall that

AIC = −

log likelihood evaluated at MLE

dim of the model

C

p

= residual sum of squares + 2σ

2

dim of the model

  1. In a simple linear regression problem, we observe n pairs of observations (y 1

, x 1

),... , (y n

, x n

where yi ∈ R and xi ∈ R

p

. The least squares estimate of the p-dimensional regression

coefficient is given by

β = arg min β∈R

p

n ∑

j=

(y j

− x

t

j

β)

2 .

A basic tool for examining the fit is the leave-one-out residual: consider fitting the model

omitting the ith observation with the corresponding LS estimate defined as

β [i]

= arg min

n ∑

j:j 6 =i

(y j

− x

t

j

β)

2 ,

get a prediction for the omitted observation, ˆy [i]

= x

t

i

β [i]

, then the leave-one-out residual

for the ith observation is (y i

− ˆy [i]

We will show that

y i

− yˆ [i]

y i

− x

t

i

β

1 − H

ii

where H ii

denotes the (i, i)th entry of the projection matrix H. That is, the leave-one-

out residual is a re-scale of the original residual, hence it is not necessary to re-fit the

model each time an observation is omitted.

We will prove this result via the following steps. Throughout, we assume all the design

matrices are of full rank, and all the LS estimates exist and are unique.

(a) Consider a new data set of n observations where y

1

= ˆy [1]

= x

t

1

β [1]

and the others

are the same as the original data. Let

β

∗ denote the corresponding LS estimate

based on the new data. Show that

β

=

β [1]

Hint: Show that

β [1]

minimizes (y

1

− x

t

1

β)

2

n

j=

(y j

− x

t

j

β)

2 .

(b) Show that

yˆ [1]

n

j=

H

1 j

y j

1 − H 11

(c) Show that

y 1

− ˆy [1]

y 1 − x

t

1

β

1 − H 11

  1. Get the Boston Housing data from the course website.

(a) Load the data into R, Splus or MATLAB. Do the following transformations on some

of the variables.

X

1

= log(X 1 ), X

2

= X 2 / 10 , X

3

= log X 3 , X

5

= log(X 5 ), X

6

= log(X 6 ),

X

7

= X

  1. 5

7

4 , X

8

= log(X 8

), X

9

= log(X 9

), X

10

= log(X 10

X

11

= e

  1. 4 X 11 / 10

3 , X

12

= X

12

/ 100 , X

13

X

13

, X

14

= log(X 14

(b) Fit a multiple linear regression model to predict the last variable (log of the median

value of owner-occupied homes in $1000’s) from the other variables. Summarize

your analysis.

(c) Search through all the models and select the best sub-model using C p

, AIC and

BIC. Compare the results.

(d) Use the Zheng-Loh model selection method and compare to (c).

(e) Report the 10-fold cross-validation errors for the full model, and the stepwise search

procedure using AIC and BIC.

In addition to your answers, please hand in a printout of necessary code.