context("meanR")

test_that('statistics', {

  expect_identical(length(reps), length(R))
  expect_identical(length(R), nrow(rs$CP))

  ## Try computation and check that all values are coherent
  expect_silent(MeanR <- rs$meanR)

  expect_true(inherits(MeanR, "meanR"))
  expect_true(MeanR$FStat > 0)
  expect_true(MeanR$p.value >= 0)
  expect_true(MeanR$p.value <= 1)

  ## Try computation with random data
  reps <- sample(1:10, length(reps), replace = TRUE)
  R <- rnorm(length(reps))
  CP <- rWishart(1, length(reps), diag(1, length(reps)))[,,1]

  expect_silent(MeanR <- fitSurface(data, fit, transforms, "loewe",
                               R, CP, reps, statistic = "meanR")$meanR)
  expect_true(MeanR$FStat > 0)
  expect_true(MeanR$p.value >= 0)
  expect_true(MeanR$p.value <= 1)

  expect_silent(summary(MeanR))

})