r语言回归分析_基于R语言实现LASSO回归分析

原文链接：

http://tecdat.cn/?p=10997​tecdat.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)

​

标准最小二乘法

 summary(ols)

## 
## Call:

## 
## Resid