pdf("rpd4_10r.pdf") alldata <- matrix( c(44,89.47,44.609,11.37,62,178,182, 40,75.07,45.313,10.07,62,185,185, 44,85.84,54.297,8.65,45,156,184, 42,68.15,59.571,8.17,40,166,172, 38,89.02,49.874,9.22,55,178,180, 47,77.45,44.811,11.63,58,176,176, 40,75.98,45.681,11.95,70,176,180, 43,81.19,49.091,10.85,64,162,170, 44,81.42,39.442,13.08,63,174,176, 38,81.87,60.055,8.63,48,170,186, 44,73.03,50.541,10.13,45,168,168, 45,87.66,37.388,14.03,56,186,192, 45,66.45,44.754,11.12,51,176,176, 47,79.15,47.273,10.6,47,162,164, 54,83.12,51.855,10.33,50,166,170, 49,81.42,49.156,8.95,44,180,185, 51,69.63,40.836,10.95,57,168,172, 51,77.91,46.672,10,48,162,168, 48,91.63,46.774,10.25,48,162,164, 49,73.37,50.388,10.08,67,168,168, 57,73.37,39.407,12.63,58,174,176, 54,79.38,46.08,11.17,62,156,176, 52,76.32,45.441,9.63,48,164,166, 50,70.87,54.625,8.92,48,146,186, 51,67.25,45.118,11.08,48,172,172, 54,91.63,39.203,12.88,44,168,172, 51,73.71,45.79,10.47,59,186,188, 57,59.08,50.545,9.93,49,148,160, 49,76.32,48.673,9.4,56,186,188, 48,61.24,47.92,11.5,52,170,176, 52,82.78,47.467,10.5,53,170,172),byrow=T,ncol=7) Y <- alldata[,3] X <- alldata[,4:7] ix <- diag(1,ncol=nrow(X),nrow=nrow(X)) X0 <- matrix(rep(1,nrow(X)),ncol=1) X1 <- X[,1] X2 <- X[,2] X3 <- X[,3] X4 <- X[,4] R0 <- t(Y) %*% X0 %*% solve(t(X0) %*% X0) %*% t(X0) %*% Y X01 <- cbind(X0,X1) R01 <- t(Y) %*% X01 %*% solve(t(X01) %*% X01) %*% t(X01) %*% Y X013 <- cbind(X0,X1,X3) R013 <- t(Y) %*% X013 %*% solve(t(X013) %*% X013) %*% t(X013) %*% Y X0132 <- cbind(X0,X1,X3,X2) R0132 <- t(Y) %*% X0132 %*% solve(t(X0132) %*% X0132) %*% t(X0132) %*% Y X01324 <- cbind(X0,X1,X3,X2,X4) R01324 <- t(Y) %*% X01324 %*% solve(t(X01324) %*% X01324) %*% t(X01324) %*% Y R1324_0 <- R01324-R0 R1_0 <- R01-R0 R3_01 <- R013-R01 R2_013 <- R0132-R013 R4_0132 <- R01324-R0132 SSRESID <- t(Y) %*% (ix - X01324 %*% solve(t(X01324)%*%X01324) %*% t(X01324)) %*% Y DFRESID <- nrow(X)-ncol(X01324) MSRESID <- SSRESID[1,1]/DFRESID R1324_0 R1_0 R3_01 R2_013 R4_0132 SSRESID DFRESID MSRESID # Compute Sequential Sums of Squares DF1324_0 <- ncol(X01324)-ncol(X0) MS1324_0 <- R1324_0/DF1324_0 F1324_0 <- MS1324_0/MSRESID R324_01 <- R01324-R01 DF324_01 <- ncol(X01324)-ncol(X01) MS324_01 <- R324_01/DF324_01 F324_01 <- MS324_01/MSRESID R24_013 <- R01324-R013 DF24_013 <- ncol(X01324)-ncol(X013) MS24_013 <- R24_013/DF24_013 F24_013 <- MS24_013/MSRESID R4_0132 <- R01324-R0132 DF4_0132 <- ncol(X01324)-ncol(X0132) MS4_0132 <- R4_0132/DF4_0132 F4_0132 <- MS4_0132/MSRESID R1324_0 DF1324_0 F1324_0 R324_01 DF324_01 F324_01 R24_013 DF24_013 F24_013 R4_0132 DF4_0132 F4_0132 # Compute Partial Sums of Squares X0234 <- cbind(X0,X2,X3,X4) R0234 <- t(Y) %*% X0234 %*% solve(t(X0234) %*% X0234) %*% t(X0234) %*% Y R1_0234 <- R01324-R0234 F1_0234 <- R1_0234/MSRESID X0124 <- cbind(X0,X1,X2,X4) R0124 <- t(Y) %*% X0124 %*% solve(t(X0124) %*% X0124) %*% t(X0124) %*% Y R3_0124 <- R01324-R0124 F3_0124 <- R3_0124/MSRESID X0134 <- cbind(X0,X1,X3,X4) R0134 <- t(Y) %*% X0134 %*% solve(t(X0134) %*% X0134) %*% t(X0134) %*% Y R2_0134 <- R01324-R0134 F2_0134 <- R2_0134/MSRESID R1_0234 F1_0234 R3_0124 F3_0124 R2_0134 F2_0134 R4_0132 F4_0132 dev.off()