Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- N <- 1000
- P <- 100
- Sigma = diag(sample(1:10, P, P))
- X <- mvrnorm(N, runif(P, -10, 10), Sigma)
- p <- rbinom(N, 1, 0.1)
- beta <- p * rnorm(P, 5, sqrt(5)) + (1-p) * rnorm(P, 0, sqrt(.1))
- e <- rnorm(N, 0, sqrt(10))
- y <- X%*%beta + e
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement