On this topic, apparently the main (?) approach to non-equilibrium stat mech is in replacing n-particle density matrices (distribution functions) with 'partial density matrices' (partial distributions) which involve only single particle and two-particle interactions, setting up
chains of eom between them, and trying to get 'statistical' concepts from these things