hfscanstat

A python implementation of the hypothesis-free multiple scan statistic with variable window for detecting clusters in point processes (L. Cucala. A hypothesis-free multiple scan statistic with variable window. Biom. J., 50 (2008), pp. 299-310).

Install

pip install hfscanstat

How to use

Using example data from Knox, G. (1959). Secular pattern of congenital oesophageal atresia. British Journal of Preventive Social Medecine 13, 222–226.

knox_1959 = [170, 316, 445, 468, 938, 1034, 1128, 1233, 1248, 1249,
        1252, 1259, 1267, 1305, 1385, 1388, 1390, 1446, 1454,
        1458, 1461, 1491, 1583, 1699, 1702, 1787, 1924, 1974,
        2049, 2051, 2067, 2075, 2108, 2151, 2174, 2191]

clusters = get_clusters(knox_1959)

print(clusters)

        X_i       X_j  i   j         I_H_F            ps
0  0.525977  0.653637  7  21  7.319199e+14 -6.661338e-16