Graphical comparison of mixing plans based on estimated probability of detection at the end of the mixing with different revolutions

Source: R/compare_mixing_2.R

This function provides a graphical display to compare mixing plans using a different number of revolutions based on the estimated probability of detection at the end of mixing.

compare_mixing_2(
  mulower,
  muupper,
  sigma,
  alpha_in,
  k,
  l,
  r,
  distribution,
  UDL,
  n_sim
)

Arguments

mulower

the lower value of the mean concentration (\(\mu\)) for use in the graphical display's x-axis.

muupper

the upper value of the mean concentration (\(\mu\)) for use in the graphical display's x-axis.

sigma

the standard deviation of the colony-forming units (CFUs) in the mixed sample on the logarithmic scale (default value 0.8)

alpha_in

concentration parameter at the initial stage

k

number of small portions / primary samples

l

number of revolutions / stages

r

the rate of the concentration parameter changes at each mixing stage

distribution

what suitable distribution type we have employed for simulation such as "Poisson-Type A" or "Poisson-Type B" or "Lognormal-Type A" or "Lognormal-Type B" or "Poisson lognormal-Type A" or "Poisson lognormal-Type B"

UDL

the upper decision limit, which depends on the type of microorganisms and testing regulations.

n_sim

number of simulations

Value

graphical display compares mixing plans using a different number of revolutions based on the estimated probability of detection at the end of the mixing.

References

  • Nauta, M.J., 2005. Microbiological risk assessment models for partitioning and mixing during food handling. International Journal of Food Microbiology 100, 311-322.

Examples

mulower <- 0.1
muupper <- 200
sigma <- 0.8
alpha_in <- 0.01
k <- 30
l <- c(500,5000)
r <- 0.01
distribution <-  "Poisson lognormal-Type B"
UDL <- 0
n_sim <- 20000
compare_mixing_2(mulower, muupper, sigma , alpha_in, k, l, r, distribution, UDL, n_sim)