hypersampleplan (version 0.1.1)

hypersampleplan.CL: Compute upper and lower confidence limits for the number of defective in an attribute sampling plan

Description

Attribute sampling is an important field of statistical quality control. When a lot is submitted for inspection of quality control, a sampling plan must specify both the number of samples to be drawn from the lot as well as the acceptance number which is the maximum number of defective items found in the sample that would still make the lot acceptable. Since the sample is not free of defective, there are probabilities of accepting a lot which is actually not acceptable and of rejecting one which is acceptable. The hypergeometric distribution using Chebyshev Polynomials forms the basis for calculating those exact probabilities.

This algorithm is proposed to compute upper and lower confidence limits for the number of defective in a lot given (n, x,N)

Usage

hypersampleplan.CL(n, x, N)

Arguments

n

sample size

x

the observed number of defective items in a random sample

N

lot size

Value

a matrix that provides one-sided confidence limits. For a lower confidence limit the confidence is CL.Lower; for an upper confidence limit the confidence is CL.Upper; for a two-sided limit it is CL.Upper - CL.Lower.

The value of k is interpreted as: For lower confidence limits, one can be (1 - Conf.Limits.1)% sure that the number of defective units in the lot is at least k1. For upper confidence limits, one can be (Conf.Limits.2)% sure that the number of defective units in the lot is at most k2.

References

Alvo, M., & Cabilio, P. (2000). Calculation of hypergeometric probabilities using Chebyshev polynomials. The American Statistician, 54(2), 141-144.

Odeh, R. (1983). Attribute sampling plans, tables of tests and confidence limits for proportions (Vol. 49). CRC Press.

Examples

# When one observes 12 defective items in a sample of size 160 taken from a lot size 1000,
# it will give the results for the one-sided 97.5% upper and lower confidence limits for
# the number of defective items in the lot.
hypersampleplan.CL(160, 12, 1000)