s(a,b) = int(kxlog(M(a,b)/C(a,b)))
dove int sta per intero del valore ottenuto moltiplicando per k il logaritmo decimale del rapporto fra M(a,b) e C(a,b).
M(a,b) è la frequenza di sostituzione dell’amminoacido a nell’amminoacido b, osservata nei blocchi di proteine omologhe considerate.
C(a,b) è la frequenza di sostituzione attesa, stimata come prodotto delle frequenze degli amminoacidi a e b nella totalità dei blocchi di proteine omologhe considerate.