library(knitr)library(MAVE)n<-400; p<-10 sigma<-diag(p)x<-rmvnorm(n,rep(0,p),sigma)beta1<-c(rep(1,2),rep(0,p-2))beta2<-c(rep(0,2),rep(1,2),rep(0,p-4))beta3<-c(rep(0,4),rep(-1,2),rep(0,p-6))beta4<-c(rep(0,6),1,-1