if (requireNamespace("logspline", quietly = TRUE) &&
    requireNamespace("KernSmooth", quietly = TRUE) &&
    requireNamespace("mclust", quietly = TRUE)) {
  context("estimate_density")

  test_that("estimate_density", {
    library(logspline)
    library(KernSmooth)
    library(mclust)

    set.seed(333)
    x <- distribution_normal(500, 1)

    # Methods
    density_kernel <- estimate_density(x, method = "kernel")
    density_logspline <- estimate_density(x, method = "logspline")
    density_KernSmooth <- estimate_density(x, method = "KernSmooth")
    density_mixture <- estimate_density(x, method = "mixture")

    testthat::expect_equal(mean(density_kernel$y - density_logspline$y), 0, tol = 0.1)
    testthat::expect_equal(mean(density_kernel$y - density_KernSmooth$y), 0, tol = 0.1)
    testthat::expect_equal(mean(density_kernel$y - density_mixture$y), 0, tol = 0.1)
  })
}