tests/testthat/test-ICA.R

In kusumikakd/EEG: A collection of utilities for EEG analysis

context("ICA analyses")

demo_epochs <- electrode_locations(demo_epochs,
                                   montage = "biosemi64alpha",
                                   overwrite = TRUE)

demo_SOBI <- run_ICA(demo_epochs, pca = 10)
demo_fastic <- run_ICA(demo_epochs, method = "fastica", pca = 10)
demo_fica <- run_ICA(demo_epochs, method = "fica", pca = 10)

test_that("topoplots for ICA work", {

  skip_on_appveyor()
  skip_on_travis()
  vdiffr::expect_doppelganger("topographical plot for SOBI",
                              topoplot(demo_SOBI,
                                       component = "Comp001"))
  vdiffr::expect_doppelganger("topographical plot for fastica",
                              topoplot(demo_fastic,
                                       component = "Comp001"))
  vdiffr::expect_doppelganger("topographical plot for fica",
                              topoplot(demo_fica,
                                       component = "Comp001"))
})

test_that("ICA timecourses work", {
  skip_on_appveyor()
  skip_on_travis()
  vdiffr::expect_doppelganger("timecourse over one component",
                              plot_timecourse(demo_SOBI,
                                              component = "Comp001"))
  vdiffr::expect_doppelganger("fastica timecourse",
                              plot_timecourse(demo_fastic,
                                              component = "Comp002"))
  vdiffr::expect_doppelganger("fastica psd",
                              plot_psd(demo_fastic))
})

test_that("component removal works", {
  skip_on_appveyor()
  skip_on_travis()
  demo_rm <- apply_ica(demo_epochs, demo_SOBI, 1)
  demo_re <- apply_ica(demo_SOBI)
  vdiffr::expect_doppelganger("removed one component",
                              plot_butterfly(demo_rm))
  vdiffr::expect_doppelganger("reconstruct all",
                              plot_butterfly(demo_re))
})

test_that("artefact detect works", {
  skip_on_appveyor()
  skip_on_travis()
  expect_equal(ar_acf(demo_SOBI, plot = FALSE), character(0))
  expect_equal(ar_chanfoc(demo_fastic, plot = FALSE), "Comp006")
  expect_equal(ar_trialfoc(demo_fica, plot = FALSE), "Comp004")

})