# American Institute of Mathematical Sciences

• Previous Article
Role of white-tailed deer in geographic spread of the black-legged tick Ixodes scapularis : Analysis of a spatially nonlocal model
• MBE Home
• This Issue
• Next Article
Optimal design for dynamical modeling of pest populations
August  2018, 15(4): 1011-1032. doi: 10.3934/mbe.2018045

## Probability of Escherichia coli contamination spread in ground beef production

 1 School of Mathematical Sciences, Rochester Institute of Technology, Rochester, NY 14623, USA 2 Department of Mathematics & Statistics, University of Guelph, Guelph ON N1G 2W1, Canada

* Corresponding author: Allan R. Willms

Received  September 07, 2017 Accepted  October 27, 2017 Published  March 2018

Human illness due to contamination of food by pathogenic strains of Escherichia coli is a serious public health concern and can cause significant economic losses in the food industry. Recent outbreaks of such illness sourced from ground beef production motivates the work in this paper. Most ground beef is produced in large facilities where many carcasses are butchered and various pieces of them are ground together in sequential batches. Assuming that the source of contamination is a single carcass and that downstream from the production facility ground beef from a particular batch has been identified as contaminated by E. coli, the probability that previous and subsequent batches are also contaminated is modelled. This model may help the beef industry to identify the likelihood of contamination in other batches and potentially save money by not needing to cook or recall unaffected batches of ground beef.

Citation: Petko M. Kitanov, Allan R. Willms. Probability of Escherichia coli contamination spread in ground beef production. Mathematical Biosciences & Engineering, 2018, 15 (4) : 1011-1032. doi: 10.3934/mbe.2018045
##### References:

show all references

##### References:
Schematic diagram showing spread of meat in ground beef production. Each carcass is spread in a like manner within a region of a raw source. The centres of regions from sequential carcasses are shifted forward in the raw source. The model can account for carcasses being spread across raw sources (dashed line). Material from the raw sources is sequentially input to consecutive batches of ground beef
Spread of carcasses in a raw source. (a): Example base probability density function for pieces from a carcass. The density is symmetric, and piece-wise linear. In this example there are $K = 2$ linear segments both to the left and right of zero, and there are discontinuities at $\pm N_2$. The parameters $K$, $N$, and $H$ must be chosen so that the area under the curve equals one. (b): Distributions of carcasses in a raw source. The centres, $\mu_c$, of the distributions of consecutive carcasses are evenly spaced along the raw source. The individual distributions overlap (more than depicted in the figure). For carcasses in the middle of the source, the probability density function $G_{sc}$ is just a shifted version of the base base probability $F_s$, as illustrated by $G_{s5}$. At the ends, distributions that extend beyond the boundary are reflected back, as indicated by the arrow and dashed line at the bottom left, and the reflected portion is added to the distribution already there yielding the dotted line distribution $G_{s1}$.
Mass input from raw source $s$ to batches. Batch $b$ receives a mass of $m_{sb}$ from source $s$. This mass is located in the interval $B_{sb} = (M_{sb}, M_{sb}+m_{sb}]$, where $M_{sb}$ is the total mass used from this source in batches prior to $b$
Probability (%) that a particular carcass in a hot source is hot given that batch $h$ is the hot batch, Equation (10), for the synthetic data set. The four separate curves in each plot correspond to hot batch $h = 5, 8, 10$, and $12$, left to right, respectively
Probability (%) that carcass number 100 in Source Ⅵ is present in other batches given that it is present in Batch 2, Equation (12), for the synthetic data set
Probability of contamination for each batch given that a fixed batch is contaminated (hot). The five separate curves correspond to Batches 2, 5, 8, 11, and 14 being the hot batch
Probability of contamination for each batch given that a fixed batch is contaminated (hot). The fifteen separate curves correspond to Batches 2, 5, 8, $\dotsc$, 41 and 44 being the hot batch
List of Symbols
 $S$ Total number of raw sources. $B$ Total number of ground beef batches. $h$ The contaminated (hot) batch. $C_s$ Number of carcasses in raw source $s$. $p_s$ Number of pieces supplied by each carcass in raw source $s$. $a_s$ Average size of pieces from each carcass in raw source $s$. $M_s$ Total mass in raw source $s$. $x$ Mass location in raw source. $\mu_c$ Mid point of piece distribution for carcass $c$ (in source $s$). $F_s$ Base probability density function for piece distribution in source $s$. $G_{sc}(x)$ Probability density function for piece distribution for carcass $c$ in source $s$. $Q_{sc}(R)$ Probability that a piece from carcass $c$ is located in region $R$ in source $s$. $K$ Half the number of piece-wise linear segments of $F_s$. $N_i$ Boundaries of piece-wise linear segments of $F_s$, measured in number of carcasses from centre, $\mu_c$, of distribution. $H^\pm_i$ Values of $F_s$ at boundary $N_i$, approaching from the left, $-$, or the right, $+$. $m_{sb}$ Mass from source $s$ input to batch $b$. $M_{sb}$ Mass from source $s$ input to batches $1$ through $b-1$. $B_{sb}$ Interval of mass locations in source $s$ input to batch $b$. $A_{sc}(B_{sb})$ Probability that carcass $c$ is absent from the set $B_{sb}$, that is, carcass $c$ contributes no pieces to batch $b$ through source $s$. $f_s$ Fraction of fat in raw source $s$. $g_s$ Relative susceptibility to contamination factor for source $s$. $V_{s_1 s_2}$ Fraction of carcasses present in both raw sources $s_1$ and $s_2$.
 $S$ Total number of raw sources. $B$ Total number of ground beef batches. $h$ The contaminated (hot) batch. $C_s$ Number of carcasses in raw source $s$. $p_s$ Number of pieces supplied by each carcass in raw source $s$. $a_s$ Average size of pieces from each carcass in raw source $s$. $M_s$ Total mass in raw source $s$. $x$ Mass location in raw source. $\mu_c$ Mid point of piece distribution for carcass $c$ (in source $s$). $F_s$ Base probability density function for piece distribution in source $s$. $G_{sc}(x)$ Probability density function for piece distribution for carcass $c$ in source $s$. $Q_{sc}(R)$ Probability that a piece from carcass $c$ is located in region $R$ in source $s$. $K$ Half the number of piece-wise linear segments of $F_s$. $N_i$ Boundaries of piece-wise linear segments of $F_s$, measured in number of carcasses from centre, $\mu_c$, of distribution. $H^\pm_i$ Values of $F_s$ at boundary $N_i$, approaching from the left, $-$, or the right, $+$. $m_{sb}$ Mass from source $s$ input to batch $b$. $M_{sb}$ Mass from source $s$ input to batches $1$ through $b-1$. $B_{sb}$ Interval of mass locations in source $s$ input to batch $b$. $A_{sc}(B_{sb})$ Probability that carcass $c$ is absent from the set $B_{sb}$, that is, carcass $c$ contributes no pieces to batch $b$ through source $s$. $f_s$ Fraction of fat in raw source $s$. $g_s$ Relative susceptibility to contamination factor for source $s$. $V_{s_1 s_2}$ Fraction of carcasses present in both raw sources $s_1$ and $s_2$.
Model parameters for the raw sources in the synthetic data set
 Source $g_s$ $f_s$ $p_s$ $a_s$ (kg) $N_{1s}$ $C_s$ $M_s$ (kg) Ⅰ(frozen lean) 0.2 0.05 25 0.5 15 160 2000 Ⅱ (frozen lean) 0.2 0.09 25 0.5 15 160 2000 Ⅲ (frozen lean) 0.2 0.07 25 0.5 15 160 2000 Ⅳ(fresh lean) 0.8 0.10 20 0.25 20 500 2500 Ⅴ (fresh lean) 0.8 0.08 20 0.25 20 600 3000 Ⅵ (fresh fat) 1.0 0.40 40 0.2 30 250 2000 Ⅶ (fresh fat) 1.0 0.45 40 0.2 30 250 2000
 Source $g_s$ $f_s$ $p_s$ $a_s$ (kg) $N_{1s}$ $C_s$ $M_s$ (kg) Ⅰ(frozen lean) 0.2 0.05 25 0.5 15 160 2000 Ⅱ (frozen lean) 0.2 0.09 25 0.5 15 160 2000 Ⅲ (frozen lean) 0.2 0.07 25 0.5 15 160 2000 Ⅳ(fresh lean) 0.8 0.10 20 0.25 20 500 2500 Ⅴ (fresh lean) 0.8 0.08 20 0.25 20 600 3000 Ⅵ (fresh fat) 1.0 0.40 40 0.2 30 250 2000 Ⅶ (fresh fat) 1.0 0.45 40 0.2 30 250 2000
Source input mass, $m_{sb}$, (kg) and total fat percentage for the synthetic data set
 Source frozen lean fresh lean fresh fat Batch Ⅰ Ⅱ Ⅲ Ⅳ Ⅴ Ⅵ Ⅶ fat % 1 312 136 552 25 2 384 52 564 25 3 114 404 260 222 25 4 262 239 231 268 25 5 201 205 89 293 212 25 6 320 180 292 100 108 15 7 407 105 284 204 15 8 390 456 154 15 9 300 205 325 170 15 10 209 211 543 37 10 11 293 132 536 39 10 12 318 94 540 48 10 13 479 454 67 10 14 701 226 73 10
 Source frozen lean fresh lean fresh fat Batch Ⅰ Ⅱ Ⅲ Ⅳ Ⅴ Ⅵ Ⅶ fat % 1 312 136 552 25 2 384 52 564 25 3 114 404 260 222 25 4 262 239 231 268 25 5 201 205 89 293 212 25 6 320 180 292 100 108 15 7 407 105 284 204 15 8 390 456 154 15 9 300 205 325 170 15 10 209 211 543 37 10 11 293 132 536 39 10 12 318 94 540 48 10 13 479 454 67 10 14 701 226 73 10
Probability (%) that sources are hot, given hot batch $h$, for the synthetic data set. ${\rm Prob}(s_1\text{ is hot }| \;h)$ is computed from Equation (9) and the data from Tables 2 and 3. Blank entries indicate zero probability due to no mass input
 Source hot frozen lean fresh lean fresh fat Batch Ⅰ Ⅱ Ⅲ Ⅳ Ⅴ Ⅵ Ⅶ 1 1.3 4.6 94.0 2 1.6 1.8 96.6 3 0.5 13.6 43.8 42.1 4 1.1 8.2 39.4 51.4 5 0.9 1.6 3.2 52.0 42.3 6 2.7 2.7 19.7 33.8 41.1 7 3.4 1.6 18.9 76.2 8 6.2 32.3 61.4 9 4.5 13.8 17.5 64.2 10 5.2 23.4 48.2 23.1 11 7.8 15.6 50.7 25.9 12 9.1 2.1 54.7 34.2 13 10.2 44.1 45.7 14 17.2 25.3 57.5
 Source hot frozen lean fresh lean fresh fat Batch Ⅰ Ⅱ Ⅲ Ⅳ Ⅴ Ⅵ Ⅶ 1 1.3 4.6 94.0 2 1.6 1.8 96.6 3 0.5 13.6 43.8 42.1 4 1.1 8.2 39.4 51.4 5 0.9 1.6 3.2 52.0 42.3 6 2.7 2.7 19.7 33.8 41.1 7 3.4 1.6 18.9 76.2 8 6.2 32.3 61.4 9 4.5 13.8 17.5 64.2 10 5.2 23.4 48.2 23.1 11 7.8 15.6 50.7 25.9 12 9.1 2.1 54.7 34.2 13 10.2 44.1 45.7 14 17.2 25.3 57.5
Probability of contamination for each batch in percent. Each column corresponds to a different hot batch
 hot batch Batch 1 2 3 4 5 6 7 8 9 10 11 12 13 14 1 100 46 3 0 0 0 0 1 0 1 0 0 0 0 2 60 100 33 12 0 1 0 1 0 1 1 0 0 0 3 4 46 100 63 32 1 1 1 1 1 1 1 1 0 4 1 20 75 100 70 40 16 1 1 1 1 1 1 0 5 1 1 32 63 100 78 41 16 1 1 1 1 0 0 6 1 1 1 29 61 100 63 28 10 1 0 0 0 0 7 1 1 1 10 23 44 100 58 31 7 5 4 1 0 8 1 1 1 1 10 20 61 100 56 14 13 15 15 11 9 1 1 1 1 1 8 35 58 100 48 24 26 29 29 10 1 1 1 1 1 1 16 29 62 100 53 28 31 32 11 1 1 1 1 1 1 11 25 43 47 100 50 35 36 12 1 1 1 1 1 1 7 21 39 20 41 100 56 42 13 1 1 1 1 1 1 2 17 34 18 22 47 100 67 14 0 0 0 0 0 0 0 10 27 14 18 27 57 100
 hot batch Batch 1 2 3 4 5 6 7 8 9 10 11 12 13 14 1 100 46 3 0 0 0 0 1 0 1 0 0 0 0 2 60 100 33 12 0 1 0 1 0 1 1 0 0 0 3 4 46 100 63 32 1 1 1 1 1 1 1 1 0 4 1 20 75 100 70 40 16 1 1 1 1 1 1 0 5 1 1 32 63 100 78 41 16 1 1 1 1 0 0 6 1 1 1 29 61 100 63 28 10 1 0 0 0 0 7 1 1 1 10 23 44 100 58 31 7 5 4 1 0 8 1 1 1 1 10 20 61 100 56 14 13 15 15 11 9 1 1 1 1 1 8 35 58 100 48 24 26 29 29 10 1 1 1 1 1 1 16 29 62 100 53 28 31 32 11 1 1 1 1 1 1 11 25 43 47 100 50 35 36 12 1 1 1 1 1 1 7 21 39 20 41 100 56 42 13 1 1 1 1 1 1 2 17 34 18 22 47 100 67 14 0 0 0 0 0 0 0 10 27 14 18 27 57 100
Number of profiles (on the diagonal) and profile matches between batches
 Batch 1 2 3 4 5 6 7 8 9 10 11 12 13 15 16 20 21 22 23 27 28 29 30 33 40 41 42 43 44 45 1 87 7 15 12 4 7 2 2 1 1 1 0 1 1 1 0 0 1 1 0 0 3 0 0 0 0 1 0 0 0 2 56 10 2 3 2 1 1 0 0 0 0 1 1 0 0 1 1 0 0 0 1 0 0 0 0 1 1 0 0 3 64 6 0 2 4 2 0 0 0 0 0 0 1 0 0 0 0 0 0 1 1 0 0 0 1 1 0 0 4 61 4 8 3 4 0 2 2 1 0 0 0 2 0 0 0 0 0 2 2 0 0 0 0 0 0 0 5 57 14 1 4 2 1 4 3 0 0 1 1 1 0 0 0 0 1 0 0 0 0 0 0 0 0 6 62 5 6 3 3 1 2 0 0 0 0 0 0 0 0 0 6 2 0 0 1 0 1 0 0 7 66 12 6 4 3 3 0 0 0 0 4 1 0 0 0 0 0 0 0 1 0 1 0 0 8 50 3 2 3 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 2 0 0 9 56 9 6 5 1 1 0 0 1 0 0 0 0 0 0 0 0 0 1 1 1 0 10 50 4 4 0 0 0 1 1 0 0 0 0 0 0 0 0 1 0 0 0 1 11 58 9 2 0 2 1 1 0 0 0 0 0 1 0 0 0 0 0 0 0 12 61 3 1 3 0 0 0 0 0 1 0 2 0 0 0 2 2 0 0 13 67 3 1 0 0 1 0 0 2 2 2 0 1 0 0 0 0 0 15 63 11 0 3 2 2 1 0 4 1 0 2 0 0 0 0 2 16 57 0 6 3 2 0 1 3 3 0 0 0 0 0 0 0 20 49 7 7 0 1 1 1 0 0 0 0 0 0 0 2 21 87 11 4 0 2 3 0 5 0 0 1 0 0 0 22 45 3 1 2 1 0 0 0 0 1 0 0 0 23 56 2 1 1 1 0 4 0 0 0 0 0 27 52 3 4 2 2 4 0 0 0 0 1 28 51 3 3 0 0 0 0 0 0 1 29 63 10 5 1 1 0 0 0 0 30 71 0 0 1 0 0 0 0 33 62 0 1 1 0 0 1 40 53 5 4 4 8 3 41 50 7 3 6 8 42 78 27 13 5 43 81 15 11 44 65 14 45 62
 Batch 1 2 3 4 5 6 7 8 9 10 11 12 13 15 16 20 21 22 23 27 28 29 30 33 40 41 42 43 44 45 1 87 7 15 12 4 7 2 2 1 1 1 0 1 1 1 0 0 1 1 0 0 3 0 0 0 0 1 0 0 0 2 56 10 2 3 2 1 1 0 0 0 0 1 1 0 0 1 1 0 0 0 1 0 0 0 0 1 1 0 0 3 64 6 0 2 4 2 0 0 0 0 0 0 1 0 0 0 0 0 0 1 1 0 0 0 1 1 0 0 4 61 4 8 3 4 0 2 2 1 0 0 0 2 0 0 0 0 0 2 2 0 0 0 0 0 0 0 5 57 14 1 4 2 1 4 3 0 0 1 1 1 0 0 0 0 1 0 0 0 0 0 0 0 0 6 62 5 6 3 3 1 2 0 0 0 0 0 0 0 0 0 6 2 0 0 1 0 1 0 0 7 66 12 6 4 3 3 0 0 0 0 4 1 0 0 0 0 0 0 0 1 0 1 0 0 8 50 3 2 3 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 2 0 0 9 56 9 6 5 1 1 0 0 1 0 0 0 0 0 0 0 0 0 1 1 1 0 10 50 4 4 0 0 0 1 1 0 0 0 0 0 0 0 0 1 0 0 0 1 11 58 9 2 0 2 1 1 0 0 0 0 0 1 0 0 0 0 0 0 0 12 61 3 1 3 0 0 0 0 0 1 0 2 0 0 0 2 2 0 0 13 67 3 1 0 0 1 0 0 2 2 2 0 1 0 0 0 0 0 15 63 11 0 3 2 2 1 0 4 1 0 2 0 0 0 0 2 16 57 0 6 3 2 0 1 3 3 0 0 0 0 0 0 0 20 49 7 7 0 1 1 1 0 0 0 0 0 0 0 2 21 87 11 4 0 2 3 0 5 0 0 1 0 0 0 22 45 3 1 2 1 0 0 0 0 1 0 0 0 23 56 2 1 1 1 0 4 0 0 0 0 0 27 52 3 4 2 2 4 0 0 0 0 1 28 51 3 3 0 0 0 0 0 0 1 29 63 10 5 1 1 0 0 0 0 30 71 0 0 1 0 0 0 0 33 62 0 1 1 0 0 1 40 53 5 4 4 8 3 41 50 7 3 6 8 42 78 27 13 5 43 81 15 11 44 65 14 45 62
Fit carcass distribution parameters for the combined raw sources. The overlap fractions $V_{24}$ and $V_{34}$ are equal to $V_{23}$
 $p$ $a$ (kg) $N_1$ $N_2$ $H^\pm_0$ $H^\pm_1$ $V_{23}$ $V_{56}$ $V_{78}$ 8 0.045 27 6383 $4.351\times 10^{-2}$ $1.177\times 10^{-5}$ 0.08 0 0.20
 $p$ $a$ (kg) $N_1$ $N_2$ $H^\pm_0$ $H^\pm_1$ $V_{23}$ $V_{56}$ $V_{78}$ 8 0.045 27 6383 $4.351\times 10^{-2}$ $1.177\times 10^{-5}$ 0.08 0 0.20
Estimated source input mass, $m_{sb}$, (kg) for real data set
 Source frozen lean fresh lean fresh fat Batch Ⅰ Ⅱ Ⅲ Ⅳ Ⅴ Ⅵ Ⅶ Ⅷ 1-10 355 355 264 255 11-13 355 355 264 255 14-17 355 355 132 132 255 18-20 355 355 264 255 21-22 355 355 264 255 23-40 709 264 255 41-45 709 264 255
 Source frozen lean fresh lean fresh fat Batch Ⅰ Ⅱ Ⅲ Ⅳ Ⅴ Ⅵ Ⅶ Ⅷ 1-10 355 355 264 255 11-13 355 355 264 255 14-17 355 355 132 132 255 18-20 355 355 264 255 21-22 355 355 264 255 23-40 709 264 255 41-45 709 264 255
Estimated and derived model parameters for the raw sources
 Source $g_s$ $f_s$ $C_s$ $M_s$ (kg) Ⅰ (frozen lean) 0.2 0.05 1,950 3,545 Ⅱ (frozen lean) 0.2 0.05 4,290 7,800 Ⅲ (frozen lean) 0.2 0.05 9,360 17,018 Ⅳ (frozen lean) 0.2 0.05 1,950 3,545 Ⅴ(fresh lean) 0.8 0.10 2,175 3,955 Ⅵ (fresh lean) 0.8 0.10 4,350 7,909 Ⅶ (fresh fat) 1.0 0.24 2,800 5,091 Ⅷ (fresh fat) 1.0 0.24 3,500 6,364
 Source $g_s$ $f_s$ $C_s$ $M_s$ (kg) Ⅰ (frozen lean) 0.2 0.05 1,950 3,545 Ⅱ (frozen lean) 0.2 0.05 4,290 7,800 Ⅲ (frozen lean) 0.2 0.05 9,360 17,018 Ⅳ (frozen lean) 0.2 0.05 1,950 3,545 Ⅴ(fresh lean) 0.8 0.10 2,175 3,955 Ⅵ (fresh lean) 0.8 0.10 4,350 7,909 Ⅶ (fresh fat) 1.0 0.24 2,800 5,091 Ⅷ (fresh fat) 1.0 0.24 3,500 6,364
Probability of contamination for Batches 1-24 in percent for hot batches (columns) 1-20. Batches 25-45 (not shown) all had probability of contamination of 2 percent for hot Batches 1-20
 hot batch Batch 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 1 100 15 12 10 8 6 5 3 1 1 0 0 0 0 0 0 0 0 0 0 2 16 100 14 10 8 7 5 3 2 1 0 0 0 0 0 0 0 0 0 0 3 12 14 100 12 8 7 5 4 3 2 1 0 0 0 0 0 0 0 0 0 4 10 10 12 100 10 7 6 5 4 3 2 1 0 0 0 0 0 0 0 0 5 9 8 8 11 100 10 7 6 5 4 3 2 1 0 0 0 0 0 0 0 6 7 7 7 7 10 100 10 7 6 5 4 3 2 1 0 0 0 0 0 0 7 5 5 5 6 7 10 100 10 7 6 5 4 3 2 1 1 0 0 0 0 8 3 3 4 5 6 7 10 100 10 7 5 5 4 3 2 1 1 0 0 0 9 2 2 3 4 5 6 7 10 100 10 6 5 5 3 3 2 1 0 0 0 10 1 1 2 3 4 5 6 7 10 100 9 6 6 4 4 3 2 1 0 0 11 0 0 1 2 3 4 5 5 6 9 100 10 7 6 5 4 3 2 1 1 12 0 0 0 1 2 3 4 5 5 6 10 100 10 6 6 5 4 2 2 1 13 0 0 0 0 1 2 3 4 5 6 7 10 100 9 6 6 5 3 3 2 14 0 0 0 0 0 1 2 3 3 4 6 6 9 100 10 7 6 5 5 5 15 0 0 0 0 0 0 1 2 3 4 5 6 7 10 100 10 7 7 6 6 16 0 0 0 0 0 0 1 1 2 3 4 5 6 7 10 100 11 8 8 7 17 0 0 0 0 0 0 0 1 1 2 3 4 5 6 7 11 100 11 9 9 18 0 0 0 0 0 0 0 0 0 1 2 2 3 5 7 8 11 100 13 11 19 0 0 0 0 0 0 0 0 0 0 1 2 2 5 6 7 9 13 100 14 20 0 0 0 0 0 0 0 0 0 0 0 1 2 5 6 7 8 11 14 100 21 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 3 4 4 5 22 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 3 4 4 23 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 3 3 3 24 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 3
 hot batch Batch 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 1 100 15 12 10 8 6 5 3 1 1 0 0 0 0 0 0 0 0 0 0 2 16 100 14 10 8 7 5 3 2 1 0 0 0 0 0 0 0 0 0 0 3 12 14 100 12 8 7 5 4 3 2 1 0 0 0 0 0 0 0 0 0 4 10 10 12 100 10 7 6 5 4 3 2 1 0 0 0 0 0 0 0 0 5 9 8 8 11 100 10 7 6 5 4 3 2 1 0 0 0 0 0 0 0 6 7 7 7 7 10 100 10 7 6 5 4 3 2 1 0 0 0 0 0 0 7 5 5 5 6 7 10 100 10 7 6 5 4 3 2 1 1 0 0 0 0 8 3 3 4 5 6 7 10 100 10 7 5 5 4 3 2 1 1 0 0 0 9 2 2 3 4 5 6 7 10 100 10 6 5 5 3 3 2 1 0 0 0 10 1 1 2 3 4 5 6 7 10 100 9 6 6 4 4 3 2 1 0 0 11 0 0 1 2 3 4 5 5 6 9 100 10 7 6 5 4 3 2 1 1 12 0 0 0 1 2 3 4 5 5 6 10 100 10 6 6 5 4 2 2 1 13 0 0 0 0 1 2 3 4 5 6 7 10 100 9 6 6 5 3 3 2 14 0 0 0 0 0 1 2 3 3 4 6 6 9 100 10 7 6 5 5 5 15 0 0 0 0 0 0 1 2 3 4 5 6 7 10 100 10 7 7 6 6 16 0 0 0 0 0 0 1 1 2 3 4 5 6 7 10 100 11 8 8 7 17 0 0 0 0 0 0 0 1 1 2 3 4 5 6 7 11 100 11 9 9 18 0 0 0 0 0 0 0 0 0 1 2 2 3 5 7 8 11 100 13 11 19 0 0 0 0 0 0 0 0 0 0 1 2 2 5 6 7 9 13 100 14 20 0 0 0 0 0 0 0 0 0 0 0 1 2 5 6 7 8 11 14 100 21 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 3 4 4 5 22 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 3 4 4 23 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 3 3 3 24 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 3
Probability of contamination for Batches 1-4 and 16-45 in percent for hot batches (columns) 21-45. Batches 5-15 (not shown) all had probability of contamination of 2 percent for hot batches 21-45
 hot batch Batch 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 3 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 4 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 $\vdots$ 16 3 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 17 3 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 18 4 3 3 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 2 2 2 2 2 19 4 4 3 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 2 2 2 2 2 20 5 4 3 3 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 21 100 14 10 9 7 5 4 3 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 22 14 100 12 9 7 6 4 3 2 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 23 10 12 100 12 8 6 5 4 3 2 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 24 9 9 12 100 11 7 5 4 4 3 2 1 0 0 0 0 0 0 0 0 0 0 0 0 0 25 7 7 8 11 100 10 7 5 4 4 3 2 1 0 0 0 0 0 0 0 0 0 0 0 0 26 6 6 6 7 10 100 10 7 5 4 4 3 2 1 0 0 0 0 0 0 0 0 0 0 0 27 4 4 5 6 7 10 100 10 7 5 4 4 3 2 1 0 0 0 0 0 0 0 0 0 0 28 3 3 4 4 5 7 10 100 10 7 5 4 4 3 2 1 0 0 0 0 0 0 0 0 0 29 1 2 3 4 4 5 7 10 100 10 7 5 4 4 3 2 1 0 0 0 0 0 0 0 0 30 1 1 2 3 4 4 5 7 10 100 10 7 5 4 4 3 2 1 0 0 0 0 0 0 0 31 0 0 1 2 3 4 4 5 7 10 100 10 7 5 4 4 3 2 1 0 0 0 0 0 0 32 0 0 0 1 2 3 4 4 5 7 10 100 10 7 5 4 4 3 2 1 1 0 0 0 0 33 0 0 0 0 1 2 3 4 4 5 7 10 100 10 7 5 4 4 3 2 1 1 0 0 0 34 0 0 0 0 0 1 2 3 4 4 5 7 10 100 10 7 5 4 4 3 2 1 1 0 0 35 0 0 0 0 0 0 1 2 3 4 4 5 7 10 100 10 7 5 4 4 3 2 1 1 0 36 0 0 0 0 0 0 0 1 2 3 4 4 5 7 10 100 10 7 5 4 4 3 2 1 1 37 0 0 0 0 0 0 0 0 1 2 3 4 4 5 7 10 100 10 7 5 5 4 3 2 2 38 0 0 0 0 0 0 0 0 0 1 2 3 4 4 5 7 10 100 10 7 5 5 4 3 3 39 0 0 0 0 0 0 0 0 0 0 1 2 3 4 4 5 7 10 100 11 6 6 5 5 5 40 0 0 0 0 0 0 0 0 0 0 0 1 2 3 4 4 5 7 11 100 9 6 6 6 7 41 0 0 0 0 0 0 0 0 0 0 0 1 1 2 3 4 5 5 6 9 100 12 9 8 8 42 0 0 0 0 0 0 0 0 0 0 0 0 1 1 2 3 4 5 5 6 12 100 13 10 10 43 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 2 3 4 5 6 9 12 100 14 12 44 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 2 3 5 6 8 10 14 100 16 45 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 2 3 5 6 8 9 12 16 100
 hot batch Batch 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 3 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 4 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 $\vdots$ 16 3 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 17 3 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 18 4 3 3 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 2 2 2 2 2 19 4 4 3 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 2 2 2 2 2 20 5 4 3 3 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 21 100 14 10 9 7 5 4 3 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 22 14 100 12 9 7 6 4 3 2 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 23 10 12 100 12 8 6 5 4 3 2 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 24 9 9 12 100 11 7 5 4 4 3 2 1 0 0 0 0 0 0 0 0 0 0 0 0 0 25 7 7 8 11 100 10 7 5 4 4 3 2 1 0 0 0 0 0 0 0 0 0 0 0 0 26 6 6 6 7 10 100 10 7 5 4 4 3 2 1 0 0 0 0 0 0 0 0 0 0 0 27 4 4 5 6 7 10 100 10 7 5 4 4 3 2 1 0 0 0 0 0 0 0 0 0 0 28 3 3 4 4 5 7 10 100 10 7 5 4 4 3 2 1 0 0 0 0 0 0 0 0 0 29 1 2 3 4 4 5 7 10 100 10 7 5 4 4 3 2 1 0 0 0 0 0 0 0 0 30 1 1 2 3 4 4 5 7 10 100 10 7 5 4 4 3 2 1 0 0 0 0 0 0 0 31 0 0 1 2 3 4 4 5 7 10 100 10 7 5 4 4 3 2 1 0 0 0 0 0 0 32 0 0 0 1 2 3 4 4 5 7 10 100 10 7 5 4 4 3 2 1 1 0 0 0 0 33 0 0 0 0 1 2 3 4 4 5 7 10 100 10 7 5 4 4 3 2 1 1 0 0 0 34 0 0 0 0 0 1 2 3 4 4 5 7 10 100 10 7 5 4 4 3 2 1 1 0 0 35 0 0 0 0 0 0 1 2 3 4 4 5 7 10 100 10 7 5 4 4 3 2 1 1 0 36 0 0 0 0 0 0 0 1 2 3 4 4 5 7 10 100 10 7 5 4 4 3 2 1 1 37 0 0 0 0 0 0 0 0 1 2 3 4 4 5 7 10 100 10 7 5 5 4 3 2 2 38 0 0 0 0 0 0 0 0 0 1 2 3 4 4 5 7 10 100 10 7 5 5 4 3 3 39 0 0 0 0 0 0 0 0 0 0 1 2 3 4 4 5 7 10 100 11 6 6 5 5 5 40 0 0 0 0 0 0 0 0 0 0 0 1 2 3 4 4 5 7 11 100 9 6 6 6 7 41 0 0 0 0 0 0 0 0 0 0 0 1 1 2 3 4 5 5 6 9 100 12 9 8 8 42 0 0 0 0 0 0 0 0 0 0 0 0 1 1 2 3 4 5 5 6 12 100 13 10 10 43 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 2 3 4 5 6 9 12 100 14 12 44 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 2 3 5 6 8 10 14 100 16 45 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 2 3 5 6 8 9 12 16 100
 [1] Vincent Calvez, Gaël Raoul, Christian Schmeiser. Confinement by biased velocity jumps: Aggregation of escherichia coli. Kinetic & Related Models, 2015, 8 (4) : 651-666. doi: 10.3934/krm.2015.8.651 [2] Antonio Di Crescenzo, Maria Longobardi, Barbara Martinucci. On a spike train probability model with interacting neural units. Mathematical Biosciences & Engineering, 2014, 11 (2) : 217-231. doi: 10.3934/mbe.2014.11.217 [3] Jie Li, Weinian Zhang. Transition between monostability and bistability of a genetic toggle switch in Escherichia coli. Discrete & Continuous Dynamical Systems - B, 2020, 25 (5) : 1871-1894. doi: 10.3934/dcdsb.2020007 [4] Yinghui Dong, Guojing Wang. Ruin probability for renewal risk model with negative risk sums. Journal of Industrial & Management Optimization, 2006, 2 (2) : 229-236. doi: 10.3934/jimo.2006.2.229 [5] Emilija Bernackaitė, Jonas Šiaulys. The finite-time ruin probability for an inhomogeneous renewal risk model. Journal of Industrial & Management Optimization, 2017, 13 (1) : 207-222. doi: 10.3934/jimo.2016012 [6] Santiago Campos-Barreiro, Jesús López-Fidalgo. KL-optimal experimental design for discriminating between two growth models applied to a beef farm. Mathematical Biosciences & Engineering, 2016, 13 (1) : 67-82. doi: 10.3934/mbe.2016.13.67 [7] H.T. Banks, Jimena L. Davis. Quantifying uncertainty in the estimation of probability distributions. Mathematical Biosciences & Engineering, 2008, 5 (4) : 647-667. doi: 10.3934/mbe.2008.5.647 [8] Jan Lorenz, Stefano Battiston. Systemic risk in a network fragility model analyzed with probability density evolution of persistent random walks. Networks & Heterogeneous Media, 2008, 3 (2) : 185-200. doi: 10.3934/nhm.2008.3.185 [9] Qingwu Gao, Zhongquan Huang, Houcai Shen, Juan Zheng. Asymptotics for random-time ruin probability in a time-dependent renewal risk model with subexponential claims. Journal of Industrial & Management Optimization, 2016, 12 (1) : 31-43. doi: 10.3934/jimo.2016.12.31 [10] Yuebao Wang, Qingwu Gao, Kaiyong Wang, Xijun Liu. Random time ruin probability for the renewal risk model with heavy-tailed claims. Journal of Industrial & Management Optimization, 2009, 5 (4) : 719-736. doi: 10.3934/jimo.2009.5.719 [11] Tianyuan Xu, Shanming Ji, Chunhua Jin, Ming Mei, Jingxue Yin. Early and late stage profiles for a chemotaxis model with density-dependent jump probability. Mathematical Biosciences & Engineering, 2018, 15 (6) : 1345-1385. doi: 10.3934/mbe.2018062 [12] Xiangxiang Huang, Xianping Guo, Jianping Peng. A probability criterion for zero-sum stochastic games. Journal of Dynamics & Games, 2017, 4 (4) : 369-383. doi: 10.3934/jdg.2017020 [13] Eleonora Bardelli, Andrea Carlo Giuseppe Mennucci. Probability measures on infinite-dimensional Stiefel manifolds. Journal of Geometric Mechanics, 2017, 9 (3) : 291-316. doi: 10.3934/jgm.2017012 [14] Michihiro Hirayama. Periodic probability measures are dense in the set of invariant measures. Discrete & Continuous Dynamical Systems - A, 2003, 9 (5) : 1185-1192. doi: 10.3934/dcds.2003.9.1185 [15] Fabio A. C. C. Chalub. An asymptotic expression for the fixation probability of a mutant in star graphs. Journal of Dynamics & Games, 2016, 3 (3) : 217-223. doi: 10.3934/jdg.2016011 [16] Noah H. Rhee, PaweŁ Góra, Majid Bani-Yaghoub. Predicting and estimating probability density functions of chaotic systems. Discrete & Continuous Dynamical Systems - B, 2019, 24 (1) : 297-319. doi: 10.3934/dcdsb.2017144 [17] Subhabrata Samajder, Palash Sarkar. Another look at success probability of linear cryptanalysis. Advances in Mathematics of Communications, 2019, 13 (4) : 645-688. doi: 10.3934/amc.2019040 [18] Benedetto Piccoli, Francesco Rossi. Measure dynamics with Probability Vector Fields and sources. Discrete & Continuous Dynamical Systems - A, 2019, 39 (11) : 6207-6230. doi: 10.3934/dcds.2019270 [19] Rongfei Liu, Dingcheng Wang, Jiangyan Peng. Infinite-time ruin probability of a renewal risk model with exponential Levy process investment and dependent claims and inter-arrival times. Journal of Industrial & Management Optimization, 2017, 13 (2) : 995-1007. doi: 10.3934/jimo.2016058 [20] Veysel Fuat Hatipoğlu. A novel model for the contamination of a system of three artificial lakes. Discrete & Continuous Dynamical Systems - S, 2018, 0 (0) : 0-0. doi: 10.3934/dcdss.2020176

2018 Impact Factor: 1.313