LNNB-snippet.R
From Organic Design wiki
p <- 1000 n <- 10 mu <- defaults["mu"] sigma <- defaults["sigma"] tau <- defaults["tau"]
means <- c(rep(rnorm(2*m1, mu, tau), each=slides),rep(rnorm(m0 , mu, tau), each=p))
X <- 2^(rnorm(n*p, means, sigma)) dim(X) <- c(p, n) X <- t(X)