7: compute Sj = SN = -log10 ℙ(N ≥ Nobs);