例解回归分析第五版10.4题目R语言代码
// An highlighted block
load("RABE5.RData") # 载入数据
A <- P290 # 载入数据
# (a)
f1 <- function(X){
(X - mean(X))/sqrt((sum((X - mean(X)) ^ 2)))
} # 标准化变换
S <- apply(A, 2, f1) # 计算标准化后的X和Y
S <- as.data.frame(S) #转换成数据框
lm1 <- lm(S$Y ~ S$X1 + S$X2 + S$X3 + S$X4 + S$X5 + S$X6 +0 )
summary(lm1) # R2很高,但很多系数不显著,数据存在共线性
# (b)
XiShu <- coefficients(lm1)
f2 <- function(X){
sqrt(