Compare the bias of estimated sigma and average model size of two methods
1 2 3 4 5 6 7 8 9 10 11 12 | ## Not run:
n<-200;p<-2000
data<-datageneration(n,p,0)
X<-data$X
Y<-data$Y
lam0<-sqrt(2*log(p)/n)
mcp<-scalescad_mcp(X,Y,lam0,method="mcp")
scad<-scalescad_mcp(X,Y,lam0,method="scad")
r1<-resultC(mcp$hsigma,mcp$coefficients)
r2<-resultC(scad$hsigma,scad$coefficients)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.