Пакет Python, поддерживающий вычисление взвешенной ковариации

Существует ли статистический пакет Python, который поддерживает вычисление взвешенной ковариации (, т. е. каждое наблюдение имеет вес )? К сожалению, numpy.cov не поддерживает веса.

Желательно работать в среде numpy/scipy (, т. е. иметь возможность использовать массивы numpy для ускорения вычислений ).

Большое спасибо!