Use a FIR filter, Matlab fir2 maybe, since IIR has nonlinear group delay you will see different time shifts for different frequencies. You can then 'shift' the resulting filtered data by just defining the x-label for plotting purposes or by filtering using filtfilt.
