In this contribution we analyze modulation electron cyclotron resonance heating (MECH) experiment and discuss higher harmonic frequency dependence of transport coefficients. We use the bidirectional heat pulse propagation method, in which both inward propagating heat pulse and outward propagating heat pulse are analyzed at a radial range, in order to distinguish frequency dependence of transport coefficients due to hysteresis from that due to other reasons, such as radially dependent transport coefficients, a finite damping term, or boundary effects. The method is applied to MECH experiments performed in various helical and tokamak devices, i.e. Large Helical Device (LHD), TJ-II, Korea Superconducting Tokamak Advanced Research (KSTAR), and Doublet III-D (DIII-D) with different plasma conditions. The frequency dependence of transport coefficients are clearly observed, showing a possibility of existence of transport hysteresis in flux-gradient relation.