inst/unitTests/test_c_utils.R

library(genoset)
library(RUnit)
#library(inline)

test_leftBound <- function() {
  #.C("leftBound", 1:10, 4, 10, 2)
}

#test_widthToStartEnd <- function() {
#  w = rep(3, 5)
#  s = numeric(5)
#  e = numeric(5)
#  myend = cumsum(w)
#  mystart = myend - 2
#  f <- cfunction(c(width="integer", start="numeric", end="numeric"), "
#                 widthToStartEnd(INTEGER(width), REAL(start), REAL(end), LENGTH(end));
#                 return(R_NilValue);
#", convention=".Call", includes="/gne/home/phaverty/Projects/R/genoset-package/genoset/src/genoset.h")
#  f(w, s, e, length(s))
#  checkEquals(s, mystart)
#  checkEquals(e, myend)
#
#}

test_RleViews_viewMeans <- function(){

  rle = Rle(1:5, rep(2, 5))
  ir = IRanges(start=c(3, 4), end=c(3, 9), names=c("GENE1", "GENE2"))
  myview = Views(rle, ir)
  s = start(ir)
  e = end(ir)
  one = .Call("RleViews_viewMeans", myview, TRUE, PACKAGE = "IRanges")
  two = .Call("rangeMeans_rle", s, e, as.numeric(runValue(rle)), runLength(rle), TRUE, PACKAGE = "genoset")
  checkEquals(one, two, checkNames=FALSE)

  rle2 = Rle(c(1, NA, 3, NA, 5), rep(2,5))
  myview2 = Views(rle2, ir)
  one2 = .Call("RleViews_viewMeans", myview2, TRUE, PACKAGE = "IRanges")
  two2 = .Call("rangeMeans_rle", s, e, as.numeric(runValue(rle2)), runLength(rle2), TRUE, PACKAGE = "genoset")
  checkEquals(one2,two2, checkNames=FALSE)

  one3 = .Call("RleViews_viewMeans", myview2, FALSE, PACKAGE = "IRanges")
  two3 = .Call("rangeMeans_rle", s, e, as.numeric(runValue(rle2)), runLength(rle2), FALSE, PACKAGE = "genoset")
  checkEquals(one3,two3, checkNames=FALSE)
}

test_numCallable_rle <- function(){
  rle = Rle(c(3L, 4L, 1L, 5L, 9L), c(4, 2, 3, 4, 5))
  ir = IRanges(start=c(3, 4, 7, 9), end=c(4, 8, 8, 17), names=c("GENE1", "GENE2", "GENE3", "GENE4"))
  s = start(ir)
  e = end(ir)
  one = c(2L, 3L, 0L, 8L)
  two = .Call("numCallable_rle", s, e, runValue(rle), runLength(rle), 3L, PACKAGE = "genoset")
  checkEquals(one, two, checkNames=FALSE)
  checkEquals(one, numCallable(rle, ir, 3L), checkNames=FALSE)
  checkEquals(one, numCallable(rle, cbind(start(ir), end(ir)), 3L), checkNames=FALSE)
  real.mat = cbind(start(ir), end(ir))
  storage.mode(real.mat) = "double"
  checkEquals(one, numCallable(rle, real.mat, 3L), checkNames=FALSE)
}

Try the genoset package in your browser

Any scripts or data that you put into this service are public.

genoset documentation built on Nov. 17, 2017, 9:22 a.m.