Each result cited is universally quantified over the data in its own statement, and is applied below with the data named at each use.
Fix the metric spaces, the maps q and s, the measure π, the image measures ν1=q#π and ν2=(s∘q)#π=s#ν1, and the kernels κ and λ of the statement. Two facts recalled in the statement are used throughout: for every w∈Y2, λw is a Borel measure on (Y1,d1) with λw(Y1)=1; and for every B∈B(Z), the function κ(⋅,B):Y1→R is measurable with respect to B(Y1) and the Borel σ-algebra of the real line, with values in [0,1] (definition of a probability kernel, applied to κ).
Step 1 (values in [0,1]). Fix w∈Y2 and then B∈B(Z). By claim 6(a) of Borel Measurability and Bounded Integration on a Metric Space, applied to the metric space (Y1,d1) and the Borel measure λw, the constant functions 0 and 1 on Y1 are integrable with respect to λw, with integrals 0⋅λw(Y1)=0 and 1⋅λw(Y1)=1. Since 0≤κ(y,B)≤1 for every y∈Y1, the monotonicity in claim 2 of Linearity and Monotonicity of the Lebesgue Integral, applied on the measure space (Y1,B(Y1),λw), gives 0≤(κ∘λ)(w,B)≤1. Moreover, by claim 6(c) of Borel Measurability and Bounded Integration on a Metric Space (same data, f=κ(⋅,B), M=1), (κ∘λ)(w,B) equals the integral of the nonnegative measurable function κ(⋅,B) in the sense of Lebesgue Integral of a Nonnegative Measurable Function; call this fact (F).
Step 2 (each (κ∘λ)(w,⋅) is a probability measure). Fix w∈Y2 and write ρ(B)=(κ∘λ)(w,B) for B∈B(Z). By Step 1, ρ maps B(Z) into [0,1]⊆[0,∞]. Since each κy is a probability measure on (Z,B(Z)), κ(y,∅)=0 and κ(y,Z)=1 for every y∈Y1, so ρ(∅) and ρ(Z) are the integrals of the constant functions 0 and 1 with respect to λw, that is, ρ(∅)=0 and ρ(Z)=1 by Step 1.
For countable additivity, let (Bj)j∈N be a sequence of pairwise disjoint members of B(Z) with union B. For each y∈Y1, countable additivity of the measure κy (Measure, Measure Space, and Probability Measure) gives κ(y,B)=∑jκ(y,Bj); as this sum is the real number κ(y,B)≤1, the definition of the sum of a sequence in [0,∞] in Measure, Measure Space, and Probability Measure shows that the partial sums are bounded above and κ(y,B) is their least upper bound. For n∈N let gn:Y1→R be the partial sum gn(y)=∑j≤nκ(y,Bj). Each gn is measurable with respect to B(Y1) by claim 2 of Arithmetic, Absolute Values, and Pointwise Limits of Measurable Real-Valued Functions, applied on (Y1,B(Y1)) to the measurable functions κ(⋅,Bj); it is nonnegative, and gn(y)≤gn+1(y) since κ(y,Bn+1)≥0; and supngn(y)=κ(y,B) for every y. By the standing convention on measurability, these real-valued measurable functions are measurable also as [0,∞]-valued functions. Hence Monotone Convergence Theorem, applied on (Y1,B(Y1),λw) to the sequence (gn), gives, for the integrals of nonnegative functions,
∫Y1κ(y,B)λw(dy)=nsup∫Y1gndλw.
By additivity of the nonnegative integral in claim 1 of Linearity and Monotonicity of the Lebesgue Integral (on the same measure space, by induction on n), ∫Y1gndλw=∑j≤n∫Y1κ(y,Bj)λw(dy). By fact (F) of Step 1, applied with B and with each Bj, the nonnegative integrals of κ(⋅,B) and κ(⋅,Bj) are the real numbers ρ(B) and ρ(Bj). So ρ(B) is the least upper bound of the partial sums ∑j≤nρ(Bj), which are real and bounded above by ρ(B); by the definition of the sum in Measure, Measure Space, and Probability Measure, ∑jρ(Bj)=ρ(B). Thus ρ is a measure on (Z,B(Z)) with ρ(Z)=1, a probability measure.
Step 3 (measurability in w). Fix B∈B(Z) and define f:Y2×Y1→R by f(w,y)=κ(y,B), where Y2×Y1 carries the product σ-algebra B(Y2)⊗B(Y1) of Product Sigma-Algebra. Then f=κ(⋅,B)∘p1 with p1(w,y)=y; p1 is measurable with respect to B(Y2)⊗B(Y1) and B(Y1) by claim 5 of Sections of Product-Measurable Sets and Maps Are Measurable, and Insertion Maps into Products Are Measurable, and κ(⋅,B) is measurable with respect to B(Y1) and the Borel σ-algebra of the real line, so f is measurable with respect to B(Y2)⊗B(Y1) by claim 4 of Borel Measurability and Bounded Integration on a Metric Space (with the metric space (Y1,d1), Ω=Y2×Y1, the map p1 in the role of Y, and g=κ(⋅,B)). Also ∣f∣≤1. Apply Integration Against a Probability Kernel: Measurable Sections, the Composite Measure on the Product and the Iterated Integral §bounded with (Y2,B(Y2)) and (Y1,B(Y1)) in place of (Y,Y) and (Z,Z), the probability kernel λ in place of κ, the probability measure ν2 in place of μ, and this f: the function w↦∫Y1f(w,y)λw(dy)=(κ∘λ)(w,B) is measurable with respect to B(Y2). Together with Step 2, the definition of a probability kernel shows that κ∘λ is a probability kernel from (Y2,B(Y2)) to (Z,B(Z)). This is the first assertion.
Step 4 (the composite measure ν2⊗λ is an image of ν1). Let Φ:Y1→Y2×Y1, Φ(y)=(s(y),y). The identity map of Y1 is measurable with respect to B(Y1) and B(Y1), as the preimage of each set is the set itself (Measurable Function and Real-Valued Measurable Function), and s is measurable with respect to B(Y1) and B(Y2). So Pairings into a Product, the Graph of a Measurable Map, and Couplings Concentrated on a Graph §pairing, applied with the measurable space (Y2,B(Y2)) in place of (Y,Y), the metric space (Y1,d1) in place of (Z,dZ), (Ω,O)=(Y1,B(Y1)), F=s and G the identity of Y1, shows that Φ is measurable with respect to B(Y1) and B(Y2)⊗B(Y1). By claim 1 of Image Measures, Measures with Densities, and Change of Variables, applied to (Y1,B(Y1),ν1) and T=Φ, the image measure Φ#ν1, E↦ν1(Φ−1(E)), is a probability measure on (Y2×Y1,B(Y2)⊗B(Y1)).
Let ν2⊗λ be the composite measure of Integration Against a Probability Kernel: Measurable Sections, the Composite Measure on the Product and the Iterated Integral §composite, with the same data as in Step 3. Let A′∈B(Y2) and C∈B(Y1). Then Φ−1(A′×C)={y∈Y1:s(y)∈A′, y∈C}=C∩s−1(A′), so
(Φ#ν1)(A′×C)=ν1(C∩s−1(A′))=∫Y21A′(w)λ(w,C)ν2(dw),
the second equality being the defining identity of a conditional kernel for λ, read with (Y1,d1), (Y2,d2), s, ν1 and s#ν1=ν2 in place of (Z,dZ), (Y,dY), q, π and ν, and with the sets A′ and C in place of A and B. Therefore Integration Against a Probability Kernel: Measurable Sections, the Composite Measure on the Product and the Iterated Integral §rectangles, with the data of Step 3 and the measure Φ#ν1 in the role of π, gives
Φ#ν1=ν2⊗λ.
Step 5 (the conditional-kernel identity). By the definition of a conditional kernel, read with Y2, s∘q (measurable, as recalled in the statement), π and (s∘q)#π=ν2 in place of Y, q, π and ν, and given Step 3, it remains to show that for all A∈B(Y2) and B∈B(Z)
π(B∩(s∘q)−1(A))=∫Y21A(w)(κ∘λ)(w,B)ν2(dw).
Fix A and then B. Define h:Y2×Y1→R by h(w,y)=1A(w)κ(y,B). The indicator 1A is measurable on (Y2,B(Y2)) by claim 1 of Arithmetic, Absolute Values, and Pointwise Limits of Measurable Real-Valued Functions, and the projection p2(w,y)=w is measurable with respect to B(Y2)⊗B(Y1) and B(Y2) by claim 5 of Sections of Product-Measurable Sets and Maps Are Measurable, and Insertion Maps into Products Are Measurable, so 1A∘p2 is measurable by claim 4 of Borel Measurability and Bounded Integration on a Metric Space (with the metric space (Y2,d2)). The function (w,y)↦κ(y,B) is the measurable f of Step 3. Hence h, their pointwise product, is measurable with respect to B(Y2)⊗B(Y1) by claim 3 of Arithmetic, Absolute Values, and Pointwise Limits of Measurable Real-Valued Functions, and ∣h∣≤1.
For each w∈Y2, h(w,⋅)=1A(w)κ(⋅,B) is a real multiple of a function integrable with respect to λw, so by claim 2 of Linearity and Monotonicity of the Lebesgue Integral (with a=1A(w) and b=0) ∫Y1h(w,y)λw(dy)=1A(w)(κ∘λ)(w,B). Thus Integration Against a Probability Kernel: Measurable Sections, the Composite Measure on the Product and the Iterated Integral §bounded, with the data of Step 3 and h in place of f, shows that h is integrable with respect to ν2⊗λ and
∫Y2×Y1hd(ν2⊗λ)=∫Y21A(w)(κ∘λ)(w,B)ν2(dw).
By Step 4 the left-hand side is ∫hd(Φ#ν1). By claim 2 of Image Measures, Measures with Densities, and Change of Variables, applied to (Y1,B(Y1),ν1), T=Φ and the measurable real function g=h, which is integrable with respect to Φ#ν1, the function h∘Φ is integrable with respect to ν1 and
∫Y2×Y1hd(Φ#ν1)=∫Y1h(Φ(y))ν1(dy).
For y∈Y1, h(Φ(y))=1A(s(y))κ(y,B)=1s−1(A)(y)κ(y,B), since s(y)∈A exactly when y∈s−1(A). The set s−1(A) belongs to B(Y1) by measurability of s, so the defining identity of a conditional kernel for κ (given q, with image measure q#π=ν1), applied with the sets s−1(A)∈B(Y1) and B, gives
∫Y11s−1(A)(y)κ(y,B)ν1(dy)=π(B∩q−1(s−1(A)))=π(B∩(s∘q)−1(A)),
the last equality because q−1(s−1(A))=(s∘q)−1(A). Chaining the last three displays with Step 4 proves the required identity, so κ∘λ is a conditional kernel of π given s∘q. □