#### 7.8 a<-read.table("C:/Documents and Settings/Mandy Jin/Desktop/TA 4315/618.txt") y<-a$V1 x1<-a$V2 x2<-a$V3 x3<-a$V4 x4<-a$V5 f1<-lm(y~x1+x2+x3+x4) e1<-f1$resi SSE1<-sum(e1^2) f2<-lm(y~x1+x4) e2<-f2$resi SSE2<-sum(e2^2) SSR<-SSE2-SSE1 ########### 7.16 b<-read.table("C:/Documents and Settings/Mandy Jin/Desktop/TA 4315/605.txt") y<-b$V1 x1<-b$V2 x2<-b$V3 n<-length(y) y.star<-1/sqrt(n-1)*(y-mean(y))/sd(y) x1.star<-1/sqrt(n-1)*(x1-mean(x1))/sd(x1) x2.star<-1/sqrt(n-1)*(x2-mean(x2))/sd(x2) lm(y.star~x1.star+x2.star) ##########7.27 a<-read.table("C:/Documents and Settings/Mandy Jin/Desktop/TA 4315/618.txt") y<-a$V1 x1<-a$V2 x2<-a$V3 x3<-a$V4 x4<-a$V5 f3<-lm(y~x1+x4) ### see f3 #### c f4<-lm(y~x4) SSR4<-anova(f4)$Sum[1] #### this is SSR(X4) f5<-lm(y~x3+x4) anova(f5) SSR43<-anova(f5)$Sum[2] #### SSR(X4|X3) f6<-lm(y~x1) SSR1<-anova(f6)$Sum[1] #### SSR(X1) f7<-lm(y~x3+x1) SSR13<-anova(f7)$Sum[2] ### SSR(X1|X3)