tests/testthat/test_00_ZABI.R

In Stan125/gamlss.dist: Distributions for Generalized Additive Models for Location Scale and Shape

########################################
########       gamlss.dist      ########
########   Distribution tests   ########
###  Zero altered  binomial: ZABI()  ###
########################################

## Context
testthat::context("ZABI")

## Skip this test - currently not important
testthat::skip("Currently not available")

## Seed
set.seed(401)

## Family
fam <- ZABI()

## Random Values
n <- 1e4
mu <- 0.5
sigma <- 0.1
rvec <- rZABI(n, bd = n,  mu, sigma) 

## Empirical Moments
ex_emp <- mean(rvec)
vx_emp <- var(rvec)

## Theoretical moments
ex_theo <- fam$mean(n, mu, sigma)
vx_theo <- fam$variance(n, mu, sigma) 

## Test here if they are about the same
expect_true(abs(ex_emp / ex_theo - 1) < 0.02)  
expect_true(abs(vx_emp / vx_theo - 1) < 0.02)