Bi-spectrum estimation