

Study with the several resources on Docsity
Earn points by helping other students or get them with a premium plan
Prepare for your exams
Study with the several resources on Docsity
Earn points to download
Earn points by helping other students or get them with a premium plan
Material Type: Assignment; Class: LINEAR STATISTICL MODELS; Subject: Statistics; University: University of South Carolina - Columbia; Term: Fall 2008;
Typology: Assignments
1 / 2
This page cannot be seen from the preview
Don't miss anything!


STAT 714, FALL 2008 HOMEWORK 5
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 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.
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