原文链接:
http://tecdat.cn/?p=10997tecdat.cn模拟假数据集
set.seed(0820)
n <- 50
p <- 25
beta <- rep(0,p)
beta[1:5] <- 1:5/5
X <- matrix(rnorm(n*p),n,p)
X <- scale(X)
Xb <- X%*%beta
Y <- X%*%beta+rnorm(n)
Y <- Y-mean(Y)
plot(cor(X,Y),xlab="j",ylab="Cor(Y,X_j)",main="Sample correlations",cex=2)
![18541ccee6e5de39662ddab0138ce091.png](https://i-blog.csdnimg.cn/blog_migrate/4bee960d81e46dd25e83ba43f5d75e33.png)
![5138c3a59fe700e19fc0d596ad1d5fbc.png](https://i-blog.csdnimg.cn/blog_migrate/3385ccab4a7ec56111dbf82477d43444.png)
标准最小二乘法
summary(ols)
##
## Call:
##
## Resid