R/testIDATs.R

In ttriche/miser: the poor man's methylation suite (cf. `sesamizer`)

#' make sure annotations and so forth are in place before doing a big run...
#' 
#' @param   samps   two or more rows of the `samps` data.frame with `Basename`
#' @param   ...     other arguments, currently unused 
#'
#' @return          the status of an attempt to annotate those IDATs 
#' 
#' @import  minfi
#' 
#' @export
testIDATs <- function(samps, ...) { 
  rgSet <- read.metharray.exp(base=".", targets=samps[1:2, , drop=FALSE])
  colnames(rgSet) <- rgSet$subject
  res <- sesamize(rgSet)
  is(res, "GenomicRatioSet")
}