poisson_lambda_one_way {LRTesteR} | R Documentation |
Test the equality of lambda parameters of poisson distributions.
Description
Test the equality of lambda parameters of poisson distributions.
Usage
poisson_lambda_one_way(x, fctr, conf.level = 0.95)
Arguments
x |
a numeric vector of at least 50 data values per group. |
fctr |
a factor vector indicating groups. |
conf.level |
overall confidence level of the likelihood intervals. Uses Bonferroni correction. |
Details
All lambdas are equal. (lambda_1 = lambda_2 ... lambda_k).
Alternative: At least one lambda is not equal.
Value
An S3 class containing the test statistic, p value, list of likelihood based confidence intervals, overall confidence level, individual confidence level of each interval and alternative hypothesis.
Source
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 <- rpois(150, 1)
fctr <- c(rep(1, 50), rep(2, 50), rep(3, 50))
fctr <- factor(fctr, levels = c("1", "2", "3"))
poisson_lambda_one_way(x, fctr, .95)
# Null is false
set.seed(1)
x <- c(rpois(50, 1), rpois(50, 2), rpois(50, 3))
fctr <- c(rep(1, 50), rep(2, 50), rep(3, 50))
fctr <- factor(fctr, levels = c("1", "2", "3"))
poisson_lambda_one_way(x, fctr, .95)
[Package LRTesteR version 1.1.1 Index]