We consider a class of nonequilibrium systems of interacting agents with pairwise interactions and quenched disorder in the dynamics featuring, in the thermodynamic limit, phase transitions. We identify mathematical conditions on the microscopic interaction structure, namely the separability of the interaction kernel, that lead to a dimension reduction of the system in terms of a finite number of reaction coordinates (RCs). Such RCs prove to be proper nonequilibrium thermodynamic variables as they carry information on correlation, memory and resilience properties of the system. Phase transitions can be identified and quantitatively characterised as singularities of the complex valued susceptibility functions associated to the RCs. We provide analytical and numerical evidence of how the singularities affect the physical properties of finite size systems.