R软件中有些包设计到LASSO回归的，然而这些包仅仅包含了一个响应变量，即y是nX1维，而我面临的是多个响应变量，即y是nXp维，求相应的程序。数据可以任意生成如下：x=matrix(rnorm(100,30),100,30)y=x[,1:4]%*%matrix(1:8,4)