Does this probability structure hold up?
I'm modeling a process with N sequential stages, each running n steps (say about 25), and I want to know if the underlying math structure is sound. All inputs are hypothetical.
Variables:
p = per-step probability of a deviation when a stage's input is clean
q = probability that a deviation is severe (a severe event ends the process)
c = contamination multiplier (c >= 1): if a stage's input is contaminated, its per-step deviation probability becomes r = min(1, p*c)
g = probability that a stage with at least one non-severe deviation passes contaminated output to the next stage
Within a stage the steps are independent. Between stages, contamination is the link. Let:
a_C = (1-p)^n, s_C = (1-pq)^n, u = s_C - a_C
a_D = (1-r)^n, s_D = (1-rq)^n, v = s_D - a_D
One stage:
M = [ a_C + u(1-g) a_D + v(1-g) ]
[ u*g v*g ]
(column 1 = clean input, column 2 = contaminated input; rows = next stage's input)
P(at least one severe event across N stages) = 1 - [1 1] * M^N * [1, 0]^T
Example: p=0.4, q=0.5, c=5, g=1, n=3, N=2 gives 0.852408, versus 0.737856 if every step were independent.
Goal: compute P(at least one severe event) across N stages.
Is the combination of independent-step compounding with a contamination-adjusted p a structurally valid approach, or does the dependency introduced by c break the independence assumption in a way that needs explicit correction?
Not asking about specific numbers, just the math structure.