Displacement distributions and Non-gaussian parameter

GROMACS version: 2018.1
GROMACS modification: No

Is there a way on gromacs to compute the non-gaussian parameter of displacement distribution for a group of molecules? Or a way to compute the k-th moment of displacement distribution to latter compute this parameter?

If not, do you know a external tool to do it using the trayectory.pdb file?