test_that("Test Genetic Interaction score calculations", {
  testthat::skip_on_cran()
  gimap_dataset <- get_example_data("gimap") %>%
    gimap_filter() %>%
    gimap_annotate(cell_line = "HELA") %>%
    gimap_normalize(
      timepoints = "day"
    ) %>%
    calc_gi()

  results <- data.frame(
    rep = c("Day05_RepA_early", "Day22_RepA_late", "Day22_RepB_late", "Day22_RepC_late"),
    intercept = as.numeric(round(c(-0.5110, 0.00210, 0.0105, 0.00501), 3)),
    slope = as.numeric(round(c(0.483, 0.659, 0.646, 0.650), 3))
  )

  gimap_dataset$overall_results$intercept <- round(gimap_dataset$overall_results$intercept, 3)
  gimap_dataset$overall_results$slope <- round(gimap_dataset$overall_results$slope, 3)

  testthat::expect_true(all.equal(gimap_dataset$overall_results$rep, results$rep))
  testthat::expect_identical(
    round(gimap_dataset$overall_results$intercept[1], 3),
    round(results$intercept[1], 3)
  )
  testthat::expect_identical(gimap_dataset$overall_results$rep[1], results$rep[1])

  testthat::expect_identical(
    round(gimap_dataset$gi_scores$mean_expected_cs[1], 3),
    round(0.360, 3)
  )
  testthat::expect_identical(
    round(gimap_dataset$gi_scores$mean_observed_cs[1], 3),
    round(-4.411, 3)
  )
  testthat::expect_identical(
    round(gimap_dataset$gi_scores$mean_gi_score[1], 3),
    round(-4.0740, 3)
  )
  testthat::expect_identical(
    round(gimap_dataset$gi_scores$p_val[1], 3),
    round(0.001203125, 3)
  )
})