5 Problems on Linear Statistical Models - Homework 5 | STAT 714, Assignments of Statistics

Material Type: Assignment; Class: LINEAR STATISTICL MODELS; Subject: Statistics; University: University of South Carolina - Columbia; Term: Fall 2008;

Typology: Assignments

Pre 2010

Uploaded on 09/02/2009

koofers-user-3rv
koofers-user-3rv 🇺🇸

9 documents

1 / 2

Toggle sidebar

This page cannot be seen from the preview

Don't miss anything!

bg1
STAT 714, FALL 2008 HOMEWORK 5
1. Consider our general linear model y=Xb +e, where E(e) = 0and cov(e) = σ2I.
Let PXdenote the perpendicular projection matrix onto C(X). Let b
ydenote the vector
of least squares fitted values and b
edenote the vector of least squares residuals. Compute
each of the following.
(a) E(b
y)
(b) cov(b
y)
(c) E(b
e)
(d) cov(b
e)
(e) cov(b
y,b
e).
2. The observed tension, y, in a nonextensible string required to maintain a body of
unknown weight, w, in equilibrium on a smooth inclined plane of angle θ, 0 < θ < π/2,
is a random variable with mean E(y) = wsin θ. For Nknown values θ1, θ2, ..., θN, set by
the experimenter and a given body, the observed data are y1, y2, ..., yN.
(a). Find b
w, the least squares estimator of w, the weight of this body.
(b) Compute E(b
w) and var( b
w). You may assume that y1, y2, ..., yNare independent.
(c) Let b
y1,b
y2, ..., b
yNdenote the least squares fitted values. Is it necessarily true that
PN
i=1(yib
yi) = 0? Explain.
3. Define the matrices
X=
1100
1100
1010
1001
and W=
1 1
1 1
1 0
1 0
(a) Show that C(W) C(X).
(b) Take y= (1,0,1,2)0. Express yas the sum of three vectors: one in C(W), one in
C(W)
C(X), and one in another space. In what space does the 3rd vector lie?
4. (a) For the simple linear regression model in Example 2.1 (notes, pp 57-59), show
algebraically that PXand PWare both equal to
1
N+(x1x)2
Pi(xix)2
1
N+(x1x)(x2x)
Pi(xix)2· · · 1
N+(x1x)(xNx)
Pi(xix)2
1
N+(x1x)(x2x)
Pi(xix)2
1
N+(x2x)2
Pi(xix)2· · · 1
N+(x2x)(xNx)
Pi(xix)2
.
.
..
.
.....
.
.
1
N+(x1x)(xNx)
Pi(xix)2
1
N+(x2x)(xNx)
Pi(xix)2· · · 1
N+(xNx)2
Pi(xix)2.
N×N
.
In regression analysis, this matrix is called the hat matrix.
(b) Compute the trace of this matrix.
(c) Use the ceramic data from HW 4 and fit both the centred and uncentred simple linear
regression models. Report least squares estimates for both models. Also, show that the
PAGE 1
pf2

Partial preview of the text

Download 5 Problems on Linear Statistical Models - Homework 5 | STAT 714 and more Assignments Statistics in PDF only on Docsity!

STAT 714, FALL 2008 HOMEWORK 5

  1. Consider our general linear model y = Xb + e, where E(e) = 0 and cov(e) = σ^2 I. Let PX denote the perpendicular projection matrix onto C(X). Let ŷ denote the vector of least squares fitted values and ̂e denote the vector of least squares residuals. Compute each of the following. (a) E(̂y) (b) cov( ŷ ) (c) E(̂e) (d) cov(̂ e) (e) cov( ŷ , ̂e).
  2. The observed tension, y, in a nonextensible string required to maintain a body of unknown weight, w, in equilibrium on a smooth inclined plane of angle θ, 0 < θ < π/2, is a random variable with mean E(y) = w sin θ. For N known values θ 1 , θ 2 , ..., θN , set by the experimenter and a given body, the observed data are y 1 , y 2 , ..., yN. (a). Find ŵ , the least squares estimator of w, the weight of this body. (b) Compute E(ŵ ) and var(ŵ ). You may assume that y 1 , y 2 , ..., yN are independent. (c) Let ∑ ŷ 1 , ŷ 2 , ..., ŷN denote the least squares fitted values. Is it necessarily true that N i=1(yi^ −^ ŷi) = 0? Explain.
  3. Define the matrices

X =

  

   and W =

  

  

(a) Show that C(W) ⊂ C(X). (b) Take y = (1, 0 , 1 , 2)′. Express y as the sum of three vectors: one in C(W), one in C(W)⊥C(X), and one in another space. In what space does the 3rd vector lie?

  1. (a) For the simple linear regression model in Example 2.1 (notes, pp 57-59), show algebraically that PX and PW are both equal to

    

1 N +^ ∑(x^1 −x)^2 i(xi−x)^2

1 N +^

(x∑ 1 −x)(x 2 −x) i(xi−x)^2

· · · (^) N^1 + (x∑^1 −x)(xN^ −x) i(xi−x)^2 1 N +^

(x∑ 1 −x)(x 2 −x) i(xi−x)^2

1 N +^ ∑(x^2 −x)^2 i(xi−x)^2

· · · (^) N^1 + (x∑^2 −x)(xN^ −x) i(xi−x)^2 .. .

1 N +^

(x∑ 1 −x)(xN −x) i(xi−x)^2

1 N +^

(x∑ 2 −x)(xN −x) i(xi−x)^2

· · · (^) N^1 + (xN^ −x) ∑^2 i(xi−x)^2

    

N ×N

In regression analysis, this matrix is called the hat matrix. (b) Compute the trace of this matrix. (c) Use the ceramic data from HW 4 and fit both the centred and uncentred simple linear regression models. Report least squares estimates for both models. Also, show that the

PAGE 1

STAT 714, FALL 2008 HOMEWORK 5

fitted values and residuals are the same for both fits. Attach all relevant output, suitably commented, if you use SAS, R, etc.

  1. Consider the one-way ANOVA model yij = μ + αi + eij , for i = 1, 2 , ..., a and j = 1 , 2 , ..., ni, so that the design matrix is

XN ×p =

   

(^1) n 1 1 n 1 0 n 1 · · · (^0) n 1 (^1) n 2 0 n 2 1 n 2 · · · (^0) n 2 .. .

(^1) na 0 na 0 na · · · (^1) na

   

where p = a + 1 and N =

∑ i ni. (a) Show that the perpendicular projection matrix onto C(X) is given by the N × N matrix PX = Blk Diag(n− i 1 Jni×ni ),

where Jni×ni is the ni × ni matrix of ones and “Blk Diag” stands for “block diagonal.” (b) Show that

y′(PX − P 1 )y =

∑^ a

i=

ni(yi+ − y++)^2.

This is the called the corrected treatment (model) sum of squares. (c) Take a = 3 and n 1 = n 2 = n 3 = 2. Break up C(PX − P 1 ) into a − 1 = 2 orthogonal subspaces following the discussion on pp 64 (notes). With your orthogonal subspaces (and their associated ppms, say, M 1 and M 2 ), verify that

y′(PX − P 1 )y = y′M 1 y + y′M 2 y,

using the observed data y = (1, 0 , 2 , 1 , 3 , 4)′.

PAGE 2