We study the response to perturbation of non-Poisson dichotomous fluctuations that generate super-diffusion. We adopt the Liouville perspective and with it a quantum-like approach based on splitting the density distribution into a symmetric and an anti-symmetric component. To accomodate the equilibrium condition behind the stationary correlation function, we study the time evolution of the anti-symmetric component, while keeping the symmetric component at equilibrium. For any realistic form of the perturbed distribution density we expect a breakdown of the Onsager principle, namely, of the property that the subsequent regression of the perturbation to equilibrium is identical to the corresponding equilibrium correlation function. We find the directions to follow for the calculation of higher-order correlation functions, an unsettled problem, which has been addressed in the past by means of approximations yielding quite different physical effects.