AOQL_grab_B: Construction of AOQ curve and calculate AOQL value based on average microbial counts

Description

AOQL_grab_B provides the AOQ curve and calculates AOQL value based on average microbial counts.

Usage

AOQL_grab_B(c, r, t, distribution,llim, K, m, sd)

Arguments

acceptance number

number of primary increments in a grab sample or grab sample size

number of grab samples

distribution

what suitable microbiological distribution we have used such as 'Poisson gamma' or 'Lognormal'or 'Poisson lognormal'

llim

the upper limit for graphing the arithmetic mean of cell count

dispersion parameter of the Poisson gamma distribution (default value 0.25)

microbiological limit with default value zero, generally expressed as number of microorganisms in specific sample weight

standard deviation of the lognormal and Poisson-lognormal distributions on the log10 scale (default value 0.8)

Value

AOQ curve and AOQL value based on average microbial counts

Details

Since $P_a$ is the probability of acceptance, $\lambda$ is the arithmetic mean of cell count and the outgoing contaminated arithmetic mean of cell count of primary increments is given by $AOQ$ as the product $\lambda P_a$. The quantity $AOQL$ is defined as the maximum proportion of outgoing contaminated primary increments and is given by $$AOQL ={\max_{\lambda \geq 0}}{\lambda P_a}$$

Examples

Run this code

# NOT RUN {
  c <-  0
  r <-  25
  t <-  30
  distribution <- 'Poisson lognormal'
  llim <- 0.20
  AOQL_grab_B(c, r, t, distribution, llim)
# }