tests/testthat/test-drle.R

In kuwisdelu/matter: Out-of-core statistical computing and signal processing

require(testthat)
require(matter)

context("drle-class")

test_that("drle indexing - logical", {

	x <- c(rep(TRUE,10),rep(FALSE,10),rep(c(TRUE,FALSE),5))
	y <- drle(x)

	expect_equal(x, y[])
	expect_equal(x[1], y[1])
	expect_equal(x[1:30], y[1:30])
	expect_equal(x[30:1], y[30:1])

})

test_that("drle indexing - integer", {

	x <- c(rep(1L,10),10:1,1:10)
	y <- drle(x)

	expect_equal(x, y[])
	expect_equal(x[1], y[1])
	expect_equal(x[1:30], y[1:30])
	expect_equal(x[30:1], y[30:1])

})

test_that("drle indexing - double", {

	x <- c(1,1,1,1,1,6,7,8,9,10,21,32,33,34,15)
	y <- drle(x)

	expect_equal(x, y[])
	expect_equal(x[1], y[1])
	expect_equal(x[1:15], y[1:15])
	expect_equal(x[15:1], y[15:1])

})

test_that("drle indexing - NAs", {

	x <- c(rep(1L,6),NA,6:10,rep(NA,3),101:105)
	y <- drle(x)

	expect_equal(x, y[])
	expect_equal(x[1], y[1])
	expect_equal(x[1:20], y[1:20])
	expect_equal(x[20:1], y[20:1])
	expect_equal(x[c(16,NA,17)], y[c(16,NA,17)])

	x <- as.double(x)
	y <- drle(x)

	expect_equal(x, y[])
	expect_equal(x[1], y[1])
	expect_equal(x[1:20], y[1:20])
	expect_equal(x[20:1], y[20:1])
	expect_equal(x[c(16,NA,17)], y[c(16,NA,17)])

})

test_that("drle endomorphic subsetting", {

	x <- c(rep(1L, 10), 10:1, 1:10)
	y <- drle(x)

	expect_equal(y, y[1:30,drop=NA])
	expect_equal(x[1:30], y[1:30,drop=NA][])
	expect_equal(x[30:1], y[30:1,drop=NA][])
	expect_equal(x[c(7,9,11,13)], y[c(7,9,11,13),drop=NA][])
	expect_equal(x[c(7:13,17:23)], y[c(7:13,17:23),drop=NA][])

	x <- c(1,1,1,1,1,6,7,8,9,10,21,32,33,34,15)
	y <- drle(x)

	expect_equal(y, y[1:15,drop=NA])
	expect_equal(x[1:15], y[1:15,drop=NA][])
	expect_equal(x[15:1], y[15:1,drop=NA][])
	expect_equal(x[c(1,3,5,7)], y[c(1,3,5,7),drop=NA][])
	expect_equal(x[c(3:7,11:14)], y[c(3:7,11:14),drop=NA][])

})

test_that("drle combining", {

	x1 <- c(1,1,1,1,1,6,7,8,9,10,21,32,33,34,15)
	x2 <- c(rep(1L, 10), 10:1, 1:10)
	y1 <- drle(x1)
	y2 <- drle(x2)

	x3 <- c(x1, x2)
	y3 <- c(y1, y2)

	expect_equal(x3, y3[])
	expect_equal(x3[1], y3[1])
	expect_equal(x3[1:45], y3[1:45])
	expect_equal(x3[45:1], y3[45:1])

})

test_that("drle factor", {

	x <- factor(c("a","a","a","b","b","b","c","b","b","a"))
	y <- drle(x)

	expect_equal(x, y[])
	expect_equal(x[1], y[1])
	expect_equal(x[1:10], y[1:10])
	expect_equal(x[10:1], y[10:1])

	x2 <- rev(x)
	y2 <- drle(x2)

	expect_equal(c(x, x2), c(y, y2)[])
	expect_equal(c(x, x2, x), c(y, y2, y)[])

	y3 <- y[1:6,drop=NULL]
	y4 <- droplevels(y3)

	expect_is(y3, "drle")
	expect_equal(y3[], x[1:6])
	expect_equal(as.character(y4[]), as.character(x[1:6]))

	expect_equal(levels(y3), levels(x))
	expect_equal(levels(y4), levels(droplevels(x[1:6])))

})