tests/testthat/test-gof.R

In AndrMenezes/unitModalReg: Unit Modal Parametric Regression Models

test_that("multiplication works", {
  set.seed(1212)
  n     <- 200
  betas <- c(1, 2)
  X     <- cbind(1, x1=runif(n))
  eta   <- c(X %*% betas)
  mu    <- exp(eta) / (1 + exp(eta))
  phi   <- exp(6.0)
  y     <- rbeta(n, shape1=mu * phi + 1, shape2=(1-mu) * phi + 1)
  data  <- data.frame(y=y, x1=X[,2])

  m1 <- unitModalReg(y ~ x1, data = data, family = "betaMode")
  m2 <- unitModalReg(y ~ x1, data = data, family = "unitGompertz")
  m3 <- unitModalReg(y ~ x1, data = data, family = "unitGamma")
  m4 <- unitModalReg(y ~ x1, data = data, family = "Kumaraswamy")

  gof(m1, m2, m3, m4)


})