## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

## ----part1--------------------------------------------------------------------
library(LRTesteR)
library(statmod)

set.seed(1)
x <- rnorm(50)

gaussian_mu_test(x = x, mu = 0, alternative = "two.sided")

## ----part2--------------------------------------------------------------------
empirical_mu_test(x = x, mu = 0, alternative = "two.sided")

## ----part3--------------------------------------------------------------------
set.seed(1)
x <- rinvgauss(50)

inverse_gaussian_mu_test(x = x, mu = 1, alternative = "two.sided")

## ----part4--------------------------------------------------------------------
empirical_mu_test(x = x, mu = 1, alternative = "two.sided")

## ----part5--------------------------------------------------------------------
set.seed(1)
x <- rnorm(50)

gaussian_variance_test(x = x, sigma.squared = 1, alternative = "two.sided")

## ----part6--------------------------------------------------------------------
empirical_variance_test(x = x, sigma.squared = 1, alternative = "two.sided")

## ----part7--------------------------------------------------------------------
set.seed(1)
x <- rnorm(100)
fctr <- c(rep(1, 50), rep(2, 50))
fctr <- factor(fctr, levels = c("1", "2"))

gaussian_mu_one_way_test(x = x, fctr = fctr, conf.level = .95)

## ----part8--------------------------------------------------------------------
empirical_mu_one_way_test(x = x, fctr = fctr, conf.level = .95)

## ----part9--------------------------------------------------------------------
set.seed(1)
x <- rinvgauss(100)
fctr <- c(rep(1, 50), rep(2, 50))
fctr <- factor(fctr, levels = c("1", "2"))

inverse_gaussian_mu_one_way_test(x = x, fctr = fctr, conf.level = .95)

## ----part10-------------------------------------------------------------------
empirical_mu_one_way_test(x = x, fctr = fctr, conf.level = .95)

## ----part11-------------------------------------------------------------------
set.seed(1)
x <- rnorm(100)
fctr <- c(rep(1, 50), rep(2, 50))
fctr <- factor(fctr, levels = c("1", "2"))

gaussian_variance_one_way_test(x = x, fctr = fctr, conf.level = .95)

## ----part12-------------------------------------------------------------------
empirical_variance_one_way_test(x = x, fctr = fctr, conf.level = .95)