Вычисление функции взаимной корреляции?

В R я использую ccf или acf для вычисления пары функция взаимной корреляции, чтобы я мог узнать, какой сдвиг дает мне максимальное значение. Судя по всему, R дает мне нормализованную последовательность значений. Есть ли что-то подобное в scipy Python, или я должен сделать это с помощью модуля fft ? В настоящее время я делаю это следующим образом:

xcorr = lambda x,y : irfft(rfft(x)*rfft(y[::-1]))
x = numpy.array([0,0,1,1])
y = numpy.array([1,1,0,0])
print xcorr(x,y)