sslasso在处理高维数据表现不错，sslasso 的模拟该如何运行呢？

rm(list())
set.seed(0)
library(MASS)
library(Matrix)
n=100;p=1000;rho=0.9;bdim=50;k=20
beta0=rep(0,p)
indb0=c(1,51,101,151,201,251)
beta0[indb0]=c(-2.5,-2,-1,-1.5,2,2.5)/sqrt(3)
ex=rep(0,p)
sig=rho^matrix(1,bdim,bdim);diag(sig)=1
Sigx=sig
for(i in 1:k-1){
 Sigx=bdiag(Sigx,sig)
}
X=mvrnorm(n,ex,Sigx)

y=as.vector(X%*%beta()+rnorm(n))

你好，我是有问必答小助手，非常抱歉，本次您提出的有问必答问题，技术专家团超时未为您做出解答


本次提问扣除的有问必答次数，将会以问答VIP体验卡（1次有问必答机会、商城购买实体图书享受95折优惠）的形式为您补发到账户。


因为有问必答VIP体验卡有效期仅有1天，您在需要使用的时候【私信】联系我，我会为您补发。