"sample.genes.ind.2" <-
function( n, sd.ctrl, var.ratio, delta, distrn = "normal")
{ ## start of the fn sample.genes.ind
## note: here, delta is a number, can change it to a vector
ngenes.ind <- length(sd.ctrl)
sample <- matrix(0, nr = ngenes.ind, nc = 2*n)
sd.trt <- sqrt(var.ratio) * sd.ctrl
sample[, 1:n] <- matrix(rnorm(n*ngenes.ind , rep(0,ngenes.ind), sd.ctrl), nr = ngenes.ind, nc = n)
# get data for the the control group
if(length(delta)== 1) delta <- rep(delta, ngenes.ind)
sample[, (n+1):(2*n)]<- matrix(rnorm(n*ngenes.ind , delta, sd.trt), nr = ngenes.ind, nc = n)
# get data for the the treatment group, the diff from null group being the delta != 0
sample
} ## the end of the function sample.genes.ind
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.