tests/testthat/test-velocity.R

In dynverse/scvelo: An R wrapper around the great "scvelo" python package

context("testing velocity.R")

test_that("velocity can be calculated", {
  set.seed(1)
  unspliced <- Matrix::rsparsematrix(50, 500, 1, rand.x = runif)
  dimnames(unspliced) <- list(paste0('A', seq_len(nrow(unspliced))), paste0('B', seq_len(ncol(unspliced))))
  spliced <- Matrix::rsparsematrix(50, 500, 1, rand.x = runif)
  dimnames(spliced) <- dimnames(unspliced)

  velocity <- get_velocity(spliced, unspliced)

  expect_is(velocity$scvelo, "anndata._core.anndata.AnnData")
})