dat <- unlist(read.table(textConnection( gsub("\n", " ", "0 0 0 0 0 -0.270340178359072 -0.367828308500142 0.666980581124872 1.87088147328446 2.14548907684624 -0.480859313143032 -0.72911654245842 -0.490308602315934 -0.66152028906509 -0.221049033806403 0.367003418871239 0.631073992586373 0.0487728723414486 0.762655178750436 0.78574757843331 0.338239686422963 0.784206454089066 -2.14265084073625 2.11325193044223 0.186018356196443 0 0 0 0 0 0 0 0 0 0 0.519132472499234 -2.604783141655 -0.244519550114012 -1.6570790528784 3.34184602886343 2.10361226260999 1.9796808733979 -0.822247322003058 1.06850578033292 -0.678811824405992 0.804225748913681 0.57363964388698 0.437113583759113 0.437208643628268 0.989892093383503 1.76545983424176 0.119483882364649 -0.222311941138971 -0.74669456611669 -0.0663660879732063 0 0 0 0 0") ), as.is = T)) test_that("test subseries", { expect_equal(subseries(dat, 0, 10), as.vector(dat[1:10])) expect_equal(subseries(dat, 17, 21), as.vector(dat[18:21])) expect_error(subseries(dat, -1, 21)) expect_error(subseries(dat, 1, 12221)) expect_error(subseries(dat, -1, 12221)) })