Throughout, PXβ and PYβ are probability measures on (S,B(S)) and (T,B(T)), hence Ο-finite, the constant sequence of the whole space witnessing it, so that PXββPYβ and Tonelli's and Fubini's theorems are available.
Step 1: the joint law. The function (X,Y):Ξ©βSΓT is F/(B(S)βB(T))-measurable (The Product Ο-Algebra: Rectangles, Generators, Sections, Projections, Pairings and the Swap (vii)), so {(X,Y)βC}=(X,Y)β1(C) is an event for every CβB(S)βB(T), and Q=Pβ(X,Y)β1 is a probability measure on (SΓT,B(S)βB(T)) (Image Measures and Measures with a Density: Measures, Integrals and Change of Variables (i)). For AβB(S) and BβB(T) we have (X,Y)β1(AΓB)=Xβ1(A)β©Yβ1(B), so by The Law of a Random Element (i) and Product Measure (i)
Q(AΓB)=P(Xβ1(A)β©Yβ1(B)),(PXββPYβ)(AΓB)=P(Xβ1(A))P(Yβ1(B)).(1)
Step 2: independence on rectangles. As {XβA,Β YβB}=Xβ1(A)β©Yβ1(B), {XβA}=Xβ1(A) and {YβB}=Yβ1(B) (Probability Spaces: Notation (iii)), Independent Random Elements: Preimages, Functions and Blocks (ii) shows that X and Y are independent if and only if
P(Xβ1(A)β©Yβ1(B))=P(Xβ1(A))P(Yβ1(B))(AβB(S),Β BβB(T)).(2)
Proof of (i). By (1), condition (2) says that Q and PXββPYβ agree on the measurable rectangles. If it holds, the two measures agree on Rect(B(S),B(T)), a Ο-system containing SΓT (The Product Ο-Algebra: Rectangles, Generators, Sections, Projections, Pairings and the Swap (i)) that generates B(S)βB(T) (Measurable Rectangles and the Product Ο-Algebra (ii)); both give the whole space the finite mass 1, so Q=PXββPYβ by Uniqueness of Measures Agreeing on a Generating Ο-System (i), which is the displayed identity. Conversely, the displayed identity for C=AΓB is (2) by (1).
Proof of (ii). As X and Y are independent, Q=PXββPYβ by (i); and g(X,Y)=gβ(X,Y). Change of variables (Image Measures and Measures with a Density: Measures, Integrals and Change of Variables (ii)) gives gβ(X,Y)βM+(Ξ©,F), and with Tonelli's theorem (Tonelli's and Fubini's Theorems (i))
Eg(X,Y)=β«Ξ©βgβ(X,Y)dP=β«SΓTβgd(PXββPYβ)=β«Sβ(β«TβgxβdPYβ)dPXβ(x).
Proof of (iii). Again Q=PXββPYβ. By change of variables (Image Measures and Measures with a Density: Measures, Integrals and Change of Variables (iii)), g(X,Y)=gβ(X,Y) is a random variable, it lies in L1(P) if and only if gβL1(PXββPYβ), and then Eg(X,Y)=β«SΓTβgd(PXββPYβ). Fubini's theorem (Tonelli's and Fubini's Theorems (ii)) gives the assertions on N and I and the equality of this integral with β«SβIdPXβ.