inst/tinytest/test_brquasiFitControl.R

In ikosmidis/brquasi: Bias reduction in quasi likelihood estimation

## Test that alias works
brquasi1 <- brquasi_control()
brquasi2 <- brquasiControl()
expect_identical(brquasi1, brquasi2)


## Test if type is passed correctly in brquasiFitControl
types <- c("M", "iRBM", "eRBM", "MPQL_trace")
for (type in types) {
    brquasi_c <- brquasi_control(type = type)
    expect_equal(brquasi_c$type, type)
}
expect_error(brquasiControl(type = "magic_bias_reduction_method"), pattern = "'arg' should be one of")


## Test that error is produced if epsilon is <= 0
expect_error(brquasiControl(epsilon = -1), pattern = "'epsilon' must be > 0")
expect_error(brquasiControl(epsilon = 0), pattern = "'epsilon' must be > 0")
expect_error(brquasiControl(lambda = -1), pattern = "'lambda' must be >= 0")

## Test disp_factor
expect_error(brquasiControl(disp_factor = "n-p-1000"), pattern = "'arg' should be one of")
expect_identical(brquasiControl(disp_factor = "n-p")$disp_factor, "n-p")
expect_identical(brquasiControl(disp_factor = "n")$disp_factor, "n")