test_that("Poisson residual no error", {
  a <- c(0, 1, 2, 3, 4, 5, 6)
  b <- c(1, 2, 3, 4, 5, 6, 7)
  c <- c(0, 1, 0, 0, 0, 1, 0)
  d <- c(3, 4, 5, 6, 7, 8, 9)
  e <- c(0, 0, 0, 1, 1, 1, 1)
  df <- data.table("a" = a, "b" = b, "c" = c, "d" = d, "e" = e)
  time1 <- "a"
  time2 <- "b"
  event <- "c"
  pyr <- "a"
  names <- c("d")
  term_n <- c(0)
  tform <- c("loglin")
  keep_constant <- c(0)
  a_n <- c(-0.1)
  modelform <- "M"

  # tfile <- file(paste(tempfile(), ".txt",sep="" ),open = "wt")
  control <- list("ncores" = 2, "lr" = 0.95, "maxiter" = -1, "halfmax" = 1, "epsilon" = 1e-9, "deriv_epsilon" = 1e-9, "step_max" = 1.0, "change_all" = TRUE, "thres_step_max" = 1.0, "verbose" = 0)
  poisres <- PoisRun(Poisson(a, c) ~ loglinear(d, 0), df, control = control, a_n = a_n)
  #  model_control <- list("pearson" = FALSE, "deviance" = FALSE)
  expect_no_error(Residual(poisres, df, pearson = FALSE, deviance = FALSE))
  expect_no_error(Residual(poisres, df, pearson = TRUE, deviance = FALSE))
  expect_no_error(Residual(poisres, df, pearson = FALSE, deviance = TRUE))
  expect_error(Residual(poisres, df, pearson = TRUE, deviance = TRUE))
  #  expect_no_error(RunPoissonRegression_Residual(df, pyr, event, names, term_n, tform, keep_constant, a_n, modelform, control, strat_col = "rand", model_control))
  #  model_control <- list("pearson" = TRUE, "deviance" = FALSE)
  #  expect_no_error(RunPoissonRegression_Residual(df, pyr, event, names, term_n, tform, keep_constant, a_n, modelform, control, strat_col = "rand", model_control))
  #  model_control <- list("pearson" = FALSE, "deviance" = TRUE)
  #  expect_no_error(RunPoissonRegression_Residual(df, pyr, event, names, term_n, tform, keep_constant, a_n, modelform, control, strat_col = "rand", model_control))
  #  #
  #  model_control <- list("pearson" = FALSE, "deviance" = FALSE, "strata" = TRUE)
  #  expect_no_error(RunPoissonRegression_Residual(df, pyr, event, names, term_n, tform, keep_constant, a_n, modelform, control, strat_col = "e", model_control))
  c <- c(0, 0, 0, 0, 0, 1, 0)
  e <- c(0, 0, 0, 1, 1, 1, 1)
  df <- data.table("a" = a, "b" = b, "c" = c, "d" = d, "e" = e)
  poisres_s <- PoisRun(Poisson_Strata(a, c, e) ~ loglinear(d, 0), df, control = control, a_n = a_n)
  expect_no_error(Residual(poisres_s, df, pearson = FALSE, deviance = FALSE))
  #  expect_no_error(RunPoissonRegression_Residual(df, pyr, event, names, term_n, tform, keep_constant, a_n, modelform, control, strat_col = "e", model_control))
  #  #
  #  keep_constant <- c(1)
  #  expect_error(RunPoissonRegression_Residual(df, pyr, event, names, term_n, tform, keep_constant, a_n, modelform, control, strat_col = "e", model_control))
  #  keep_constant <- c(0)
  #  names <- c("CONST")
  #  expect_no_error(RunPoissonRegression_Residual(df, pyr, event, names, term_n, tform, keep_constant, a_n, modelform, control, strat_col = "e", model_control))
  #  #
  a <- c(0, 1, 2, 3, 4, 5, 6)
  b <- c(1, 2, 3, 4, 5, 6, 7)
  c <- c(0, 0, 0, 0, 0, 0, 0)
  d <- c(3, 4, 5, 6, 7, 8, 9)
  e <- c(0, 0, 0, 1, 1, 1, 1)
  df <- data.table("a" = a, "b" = b, "c" = c, "d" = d, "e" = e)
  expect_error(Residual(poisres_s, df, pearson = FALSE, deviance = FALSE))
  #  expect_error(RunPoissonRegression_Residual(df, pyr, event, names, term_n, tform, keep_constant, a_n, modelform, control, strat_col = "e", model_control))
})