inst/unitTests/test_findSpectrum.R

In CellTrails: Reconstruction, visualization and analysis of branching trajectories

test_findSpectrum <- function(){
  data(exSCE)
  dat <- SingleCellExperiment(assay=list(logcounts=logcounts(exSCE)))
  se <- embedSamples(dat)

  #TEST
  d <- findSpectrum(se$eigenvalues, frac=30)
  RUnit::checkEquals(d, seq_len(10))
  latentSpace(dat) <- se$components[, d]
  RUnit::checkException(latentSpace(dat) <- se$components[seq_len(2), d])
}