tests/testthat/test_lqd2dens.R

In fdadensity: Functional Data Analysis for Density Functions by Transformation to a Hilbert Space

library(testthat)

test_that('lqd2dens works on a trivial example, U(0,2)', {
  x = seq(0,1,length.out =512)
  y.lqd = rep(log(2),length.out =512)
  expect_equal( rep(0.5, 512), lqd2dens(lqd=y.lqd, lqdSup=x, dSup = seq(0, 2, length.out = 512) )) 
})

test_that('lqd2dens works on a second trivial example', {
  x =  seq(0,2,length.out = 512)
  y = (seq(0,2,length.out =512) + 0.1)/2.2
  y.lqd <- dens2lqd( dens=y, dSup = x)
  expect_equal( lqd2dens(lqd = y.lqd, lqdSup = seq(0, 1, length.out = length(x)), dSup = x)[c(1,512)], y[c(1,512)], tolerance=1e-3)  
})