gamma_rate_one_sample {LRTesteR}

Test the rate parameter of a gamma distribution.

Description

Test the rate parameter of a gamma distribution.

Usage

gamma_rate_one_sample(x, rate, alternative = "two.sided", conf.level = 0.95)

Arguments

x

a numeric vector of at least 50 data values.

rate

a number indicating the tested value of the rate parameter.

alternative

a character string specifying the alternative hypothesis, must be one of "two.sided" (default), "greater" or "less".

conf.level

confidence level of the likelihood interval.

Value

An S3 class containing the test statistic, p value, likelihood based confidence interval, and alternative hypothesis.

Source

• https://en.wikipedia.org/wiki/Likelihood-ratio_test

• Yudi Pawitan. In All Likelihood. Oxford University Press.

• Hodd, McKean, and Craig. Introduction to Mathematical Statistics. Pearson.

Examples

library(LRTesteR)

# Null is true
set.seed(1)
x <- rgamma(100, shape = 1, rate = 1)
gamma_rate_one_sample(x, 1, "two.sided")

# Null is false
set.seed(1)
x <- rgamma(100, shape = 1, rate = 2)
gamma_rate_one_sample(x, 1, "greater")