compare_mixing_2: Graphical comparison of mixing plans based on estimated...

In Mayooran1987/mixingsimulation: Microbiological Risk Assessment of Powder Mixing Process by Simulation Results

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

Description

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.

Usage

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.

See Also

sim_single_pd_stages

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 <- 2000
compare_mixing_2(mulower, muupper, sigma , alpha_in, k, l, r, distribution, UDL, n_sim)

Mayooran1987/mixingsimulation documentation built on May 12, 2023, 12:36 a.m.