cat(cli::col_yellow("test of inverse Gamma:\n"))

data("wafers")
HLig <- HLfit( y ~X1+X2+X1*X3+X2*X3+I(X2^2)+(1|batch),
       family=Gamma(log),rand.family=inverse.Gamma(log),
       resid.model= ~ X3+I(X3^2) ,data=wafers)
testthat::expect_equal(HLig$APHLs$p_v,-1157.523,tolerance=1e-3)

# late addition: only test of Gamma(identity), now that lambda range is automatically controlled
HLgi <- HLfit( y ~X1+X2+X1*X3+X2*X3+I(X2^2)+(1|batch),
               family=Gamma(log),rand.family=Gamma(identity),
               resid.model= ~ X3+I(X3^2) ,data=wafers)
testthat::expect_equal(logLik(HLgi),c(P_bv=-1175.070),tolerance=1e-3)