getSpanningSet: Split Data into Bins
In Sage-Bionetworks/snm: Supervised Normalization of Microarrays

Description Usage Arguments Author(s) Examples

Splits data into bins that span the range of intensities

1	getSpanningSet(snm.obj)

snm.obj

An object of class snm

Brig Mecham <brig.mecham@sagebase.org> and John D. Storey <jstorey@princeton.edu>

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function(snm.obj) {
  m <- rowMeans(snm.obj$M[snm.obj$nulls,])
  steps <- seq(quantile(snm.obj$M[snm.obj$nulls,],0.001),
               quantile(snm.obj$M[snm.obj$nulls,],0.999), length.out=snm.obj$nbins)
  th <- sapply(m, function(x) { 
    sum(x <= steps)
  }) 
  bins <- split(snm.obj$nulls, th)
  bins
  }