cocoatree.statistics.position.compute_entropy

cocoatree.statistics.position.compute_entropy(aa_freq)

Computes Shannon’s entropy for each position in the alignment

\[H(a) = -\sum_i f_{ia} \log f_{ia}\]

where H(a) is the relative entropy of amino acid a,

fia is the frequency of amino acid a at position i

Parameters

aa_freqnp.ndarray,

amino acid frequencies per position

Returns

s: array of shape (N_pos)