For a cell i, j
CountWhenIndependent(i,j) = Total(i) * Total(j) / dBase
PrcInd (i,j) = CountWhenIndependent(i,j) / dBase
PrcObs (i,j) = Observe(i,j) / dBase
Sigma = (PrcObs (i,j) - PrcInd (i,j)) /sqrt(PrcInd (i,j) * (1. - PrcInd (i,j)) / dBase)
N1 = Total(j)
N2 = dBase - N1
If N1 and N2 >
P1 = Observe(i,j) / N1
P2 = (Total(i) - Observe(i,j) ) / N2
Denominator = Squared root ( P1 * (1.- P1) / N1 + P2 * (1.-P2) / N2)
Sigma = (P1-P2) / Denominator
N1 = Total(i)
N2 = dBase - N1
If N1 and N2 >
P1 = Observe(i,j) / N1
P2 = (Total(j) - Observe(i,j) ) / N2
Denominator = Squared root ( P1 * (1.- P1) / N1 + P2 * (1.-P2) / N2)
Sigma = (P1-P2) / Denominator