These functions provides true concentration level in the original sample when diluted samples collected from a heterogeneous batch.

true_concentration_heterogeneous(mu, sd, a, b, f, u, USL, n_sim)

true_concentration_heterogeneous_multiple(mu, sd, a, b, f, u, USL, n_sim)

true_concentration_curves_heterogeneous(mu_low, mu_high, sd, a, b, f, u, USL, n_sim)

Arguments

mu

the mean concentration (on the log scale).

sd

the standard deviation of the normal distribution (on the log scale).

a

lower domain of the number of cell counts.

b

upper domain of the number of cell counts.

f

final dilution factor.

u

amount put on the plate.

USL

upper specification limit.

n_sim

number of simulations (large simulations provide a more precise estimation).

mu_low

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

mu_high

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

Value

true concentration level when sample collected from a heterogeneous batch.

Details

Let Y be the count of microorganisms and C be the true concentration level (in counts per ml). When diluted sample collected from heterogeneous (non-homogeneous) batch, Y can be modelled by Poisson lognormal distribution with parameter \(\mu, \sigma\). Let X be the count of microorganisms on a plate, and it can be modelled by truncated Poisson lognormal distribution with parameters \(\mu_d,\sigma, a, b\). Also, \(\lambda_d\) can be written in terms of \(\mu\),f and u. It is given by $$\mu_d = \mu + \log(f)+ \log(u)$$ And the true concentration level is given by $$C = \frac{X}{f*u}$$ where \(f\) is final dilution factor and \(u\) is amount of diluted sample on plate. Based on the literatures, we used \(\sigma = 0.2\) in these dilution process; see Gonzales-Barron et al. (2013, p. 370) and Schothorst et al. (2009).

References

  • Gonzales-Barron, U.A., Pilão Cadavez, V.A., Butler, F., 2013. Statistical approaches for the design of sampling plans for microbiological monitoring of foods, in: Mathematical and Statistical Methods in Food Science and Technology. Wiley, Chichester, UK, pp.363–384.

  • Schothorst, M. van, Zwietering, M.H., Ross, T., Buchanan, R.L., Cole, M.B., 2009. Relating microbiological criteria to food safety objectives and performance objectives. Food Control 20, 967–979.

Examples

mu_low <- 0
mu_high <- 10
sd <- 0.2
a <- 0
b <- 300
f <- c(0.01,0.1)
u <- c(0.1,0.1)
USL <- 1000
n_sim <- 50000
true_concentration_curves_heterogeneous(mu_low, mu_high, sd, a, b, f, u, USL, n_sim)