用list 一个循环就可以了
s<-list()
for(i in 1:100)
{
n =100 #产生100个观测值
p =25 #产生25个自变量
beta =rep.int(c(1,0),c(5,p-5)) # 定义自变量系数beta可以改变
sigma =0.5 # 定义残差系数（可以改变）
Sigma =0.5^t(sapply(1:p, function(i, j) abs(i-j), 1:p))#定义产生x矩阵时x之间的协方差阵系数0.5
x=rmvnorm(n, sigma=Sigma) # 生成x矩阵
e=rnorm(n) # 生成残差矩阵
y=c(x %*% beta + sigma * e) # 生成y矩阵
s[[i]]<-y
}