The level statistics in the transition between delocalized and localized {phases of} many body interacting systems is {considered}. We recall the joint probability distribution for eigenvalues resulting from the statistical mechanics for energy level dynamics as introduced by Pechukas and Yukawa. The resulting single parameter analytic distribution is probed numerically {via Monte Carlo method}. The resulting higher order spacing ratios are compared with data coming from different {quantum many body systems}. It is found that this Pechukas-Yukawa distribution compares favorably with {$beta$--Gaussian ensemble -- a single parameter model of level statistics proposed recently in the context of disordered many-body systems.} {Moreover, the Pechukas-Yukawa distribution is also} only slightly inferior to the two-parameter $beta$-h ansatz shown {earlier} to reproduce {level statistics of} physical systems remarkably well.