R/est.sd.ctrl.R

In warnes/exp.ssize: Estimate Effectivenes of Microarry Sample Size Estimation via Simulation

"est.sd.ctrl" <-
function(sample.ctrl, nrep)
{

ngenes <- dim(sample.ctrl)[1]
nrep <- nrep[nrep <= dim(sample.ctrl)[2]]

est.sd <- matrix(0, ngenes, length(nrep));

  for ( j in 1:length(nrep))
   {
   est.sd[,j] <- sd(t(sample.ctrl[,1:nrep[j]]))
   colnames(est.sd) <- as.character(nrep)
   est.sd
   } # end of for j

}