library(LRTesteR)
library(statmod)
# Null is true
set.seed(1)
x <- rinvgauss(n = 100, mean = 1, shape = 2)
inverse_gaussian_shape_one_sample(x, 2, "two.sided")
# Null is false
set.seed(1)
x <- rinvgauss(n = 100, mean = 1, shape = 2)
inverse_gaussian_shape_one_sample(x, 1, "greater")
Run the code above in your browser using DataLab