Add incidence function

[rogerssam]

Users may want to see how often pairs of treatments co-occur next to each other, or in the same row/column etc.

[wvjgsuhp]

A wrapper around calculate_nb should do it for adjacent pair.