2: compute S = -log10 ℙ(N ≥ Nobs);