Throughout we use Linearity and Monotonicity of the Lebesgue Integral (linearity and monotonicity of the integral) and Monotone Convergence Theorem, together with the notation and [0,∞] conventions of the statement.
Step 0 (dyadic approximation). For a natural number L≥1 let GL={j2−L:j=0,1,…,L2L}, a finite set of nonnegative reals containing 0, and for t∈[0,∞] let φL(t) be the largest element of GL not exceeding t (this exists, GL being finite and containing 0≤t). Then:
(i) φL(t)≤φL+1(t)≤t for all t, since GL⊆GL+1: indeed j2−L=(2j)2−L−1, and j≤L2L implies 2j≤(L+1)2L+1.
(ii) φL(t)→t as L→∞ for every t∈[0,∞]: for real t≥0 and every natural number L≥t the elements of GL have spacing 2−L throughout [0,L]⊇[0,t], so t−φL(t)≤2−L; and φL(∞)=L→∞.
(iii) If u:Z→[0,∞] is measurable on a measurable space (Z,Z), then φL∘u takes finitely many values, and for c=j2−L with j<L2L one has {φL∘u=c}={u≥c}∖{u≥c+2−L}, while {φL∘u=L}={u≥L}; each set {u≥c}=⋂n≥1{u>c−1/n} is measurable. A real-valued function with finitely many values whose level sets are measurable is measurable (every preimage of a set is the finite union of the relevant level sets), so φL∘u is a nonnegative simple function. By (i) and (ii), φL∘u↑u pointwise.
Step 1 (claim 1). μT(∅)=μ(∅)=0. If B1,B2,⋯∈S are pairwise disjoint, the preimages T−1(Bm) are pairwise disjoint members of F with T−1(⋃mBm)=⋃mT−1(Bm), so countable additivity of μ gives countable additivity of μT. Finally μT(S)=μ(T−1(S))=μ(X), and a probability measure is exactly a measure of total mass 1 (Measure, Measure Space, and Probability Measure).
Step 2 (claim 2). For B∈S one has 1B∘T=1T−1(B) pointwise, so by the integral of simple functions,
∫S1BdμT=μT(B)=μ(T−1(B))=∫X1B∘Tdμ.
For a nonnegative simple s=∑m=1ncm1Bm the identity follows by linearity, since s∘T=∑mcm1T−1(Bm) is again simple. For measurable g:S→[0,∞]: g∘T is measurable, since for real a one has {g∘T>a}=T−1({g>a})∈F; Step 0 gives nonnegative simple φL∘g↑g with (φL∘g)∘T=φL∘(g∘T)↑g∘T; and two applications of Monotone Convergence Theorem give
∫SgdμT=Lsup∫SφL∘gdμT=Lsup∫X(φL∘g)∘Tdμ=∫Xg∘Tdμ.
For measurable g:S→R, apply this to the positive and negative parts g±, noting (g∘T)±=g±∘T pointwise. The two resulting identities show ∫g±dμT and ∫g±∘Tdμ are finite together, which is the stated equivalence of integrability (Integrable Function and the Lebesgue Integral), and in that case subtracting them gives the display in R.
Step 3 (claim 3). νh(∅)=0, the integrand being 0. For pairwise disjoint A1,A2,⋯∈F, pointwise 1⋃mAmh=supn∑m=1n1Amh with nondecreasing partial sums, so Monotone Convergence Theorem and linearity give
νh(m⋃Am)=nsupm=1∑n∫X1Amhdμ=m∑νh(Am),
so νh is a measure. Here and below, sums and products of real-valued measurable functions are measurable by Sequentially Continuous Functions of Measurable Euclidean Maps are Measurable, the maps (u,v)↦u+v and (u,v)↦uv being continuous on R2. Next, for nonnegative simple s=∑mcm1Bm (Bm∈F), linearity gives
∫Xsdνh=m∑cmνh(Bm)=m∑cm∫X1Bmhdμ=∫Xshdμ.
For measurable f:X→[0,∞], Step 0 gives simple φL∘f↑f; then (φL∘f)h↑fh pointwise: where f is finite this is continuity of multiplication, where f=∞ and h>0 one has (φL∘f)h=Lh↑∞=fh (convention ∞⋅a=∞ for a>0), and where f=∞ and h=0 both sides vanish (convention ∞⋅0=0). Hence fh=supL(φL∘f)h pointwise, so for real a the set {fh>a} equals X for a<0 and ⋃L{(φL∘f)h>a} for a≥0 (the sequence being nondecreasing), and fh is measurable, each (φL∘f)h being a product of real-valued measurable functions. Monotone Convergence Theorem on both sides, with the simple case just proved, gives
∫Xfdνh=Lsup∫XφL∘fdνh=Lsup∫X(φL∘f)hdμ=∫Xfhdμ.
Finally, for measurable f:X→R: pointwise (fh)±=f±h since h≥0, so ∫f±dνh=∫f±hdμ; both sides are finite together, giving the integrability equivalence, and subtraction gives the final display. ■