Nothing

```
rm(list=ls())
library(quadrupen)
## Simulating multivariate Gaussian with blockwise correlation
## and piecewise constant vector of parameters
beta <- rep(c(0,1,0,-1,0), c(25,10,25,10,25))
Soo <- matrix(0.75,25,25) ## bloc correlation between zero variables
Sww <- matrix(0.75,10,10) ## bloc correlation between active variables
Sigma <- bdiag(Soo,Sww,Soo,Sww,Soo) + 0.2
diag(Sigma) <- 1
## This gives a great advantage to the elastic-net
## for support recovery; have a look at image(Sigma)
n <- 100
x <- as.matrix(matrix(rnorm(95*n),n,95) %*% chol(Sigma))
y <- 10 + x %*% beta + rnorm(n,0,10)
lasso <- elastic.net(x,y, lambda2=0)
enet <- elastic.net(x,y, lambda2=10)
breg1 <- bounded.reg(x,y)
breg2 <- bounded.reg(x,y, lambda2=10)
## Plot the Lasso path
plot(lasso, main="Lasso solution path")
## Plot the Elastic-net path
plot(enet, main = "Elastic-net solution path")
## Plot the Elastic-net path (fraction on X-axis, unstandardized coefficient)
plot(enet, standardize=FALSE, xvar="fraction")
## Plot the bounded regression path (fraction on X-axis)
plot(breg1)
## Plot another bounded regression path (fraction on X-axis, unstandardized coefficient)
plot(breg2, standardize=FALSE, xvar="fraction")
```

**Any scripts or data that you put into this service are public.**

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.