tests/testthat/test-circacompare_mixed.R

In circacompare: Analyses of Circadian Data

test_that("circacompare_mixed() works", {
  set.seed(99)
  phi1_in <- 3.15
  mixed_data <- function(n){
    counter <- 1
    for(i in 1:n){
      x <- make_data(k1=0, alpha1=0, phi1=rnorm(1, phi1_in, 0.5), hours=72, noise_sd = 1)
      x$id <- counter
      counter <- counter + 1
      if(i==1){res <- x}else{res <- rbind(res, x)}
    }
    return(res)
  }

  df <- mixed_data(20)

  out <- circacompare_mixed(
    x = df,
    col_time = "time",
    col_group = "group",
    col_outcome = "measure",
    col_id = "id",
    control=list(grouped_params=c("phi"), random_params=c("phi1"))
  )

  phi1_fit <- extract_model_coefs(out$fit)['phi1',]
  phi1_est <- phi1_fit['estimate']
  phi1_ll <- phi1_est - 1.96*phi1_fit['std_error']
  phi1_ul <- phi1_est + 1.96*phi1_fit['std_error']
  expect_true(phi1_in < phi1_ul & phi1_in > phi1_ll)


})