Skip to content

foapy.characteristics.ma.identifying_information

identifying_information(intervals, dtype=None)

Calculates identifying informations (amount of information) of the intervals grouped by congeneric sequence.

\[ \left[ H_j \right]_{1 \le j \le m} = \left[ \log_2 { \left(\frac{1}{n_j} * \sum_{i=1}^{n_j} \Delta_{ij} \right) } \right]_{1 \le j \le m} \]

where \( \Delta_{ij} \) represents \(i\)-th interval of \(j\)-th congeneric intervals array, \( n_j \) is the total number of intervals in \(j\)-th congeneric intervals array and \(m\) is number of congeneric intervals arrays.

Parameters:

Name Type Description Default

intervals

array_like

An array of congeneric intervals array

 required

dtype

dtype

The dtype of the output

 None

Returns:

Type Description
array

An array of the identifying information of congeneric intervals.

Examples:

Calculate the identifying information of a sequence.

1
2
3
4
5
6
7
8
9
import foapy
import numpy as np

source = np.array(['a', 'b', 'a', 'c', 'a', 'd'])
order = foapy.ma.order(source)
intervals = foapy.ma.intervals(order, foapy.binding.start, foapy.mode.normal)
result = foapy.characteristics.ma.identifying_information(intervals)
print(result)
# [0.73696559 1.         2.         2.5849625 ]

Calculate the identifying information of congeneric intervals of a sequence.

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
import foapy

X = []
X.append([1, 1, 4, 4])
X.append([3, 1, 3])
X.append([5, 3, 1])

result = foapy.characteristics.ma.identifying_information(X)
print(result)
# [1.32192809 1.22239242 1.5849625 ]
Source code in .tox/docs-deploy/lib/python3.11/site-packages/foapy/characteristics/ma/_identifying_information.py 
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
def identifying_information(intervals, dtype=None):
    """
    Calculates identifying informations (amount of information) of the intervals
    grouped by congeneric sequence.

    $$
    \\left[ H_j \\right]_{1 \\le j \\le m} =
    \\left[
    \\log_2 { \\left(\\frac{1}{n_j} * \\sum_{i=1}^{n_j} \\Delta_{ij} \\right) }
    \\right]_{1 \\le j \\le m}
    $$

    where \\( \\Delta_{ij} \\) represents $i$-th interval of $j$-th
    congeneric intervals array, \\( n_j \\) is the total
    number of intervals in $j$-th congeneric intervals array
    and $m$ is number of congeneric intervals arrays.

    Parameters
    ----------
    intervals : array_like
        An array of congeneric intervals array
    dtype : dtype, optional
        The dtype of the output

    Returns
    -------
    : array
        An array of the identifying information of congeneric intervals.

    Examples
    --------

    Calculate the identifying information of a sequence.

    ``` py linenums="1"
    import foapy
    import numpy as np

    source = np.array(['a', 'b', 'a', 'c', 'a', 'd'])
    order = foapy.ma.order(source)
    intervals = foapy.ma.intervals(order, foapy.binding.start, foapy.mode.normal)
    result = foapy.characteristics.ma.identifying_information(intervals)
    print(result)
    # [0.73696559 1.         2.         2.5849625 ]
    ```

    Calculate the identifying information of congeneric intervals of a sequence.

    ``` py linenums="1"
    import foapy

    X = []
    X.append([1, 1, 4, 4])
    X.append([3, 1, 3])
    X.append([5, 3, 1])

    result = foapy.characteristics.ma.identifying_information(X)
    print(result)
    # [1.32192809 1.22239242 1.5849625 ]
    ```
    """  # noqa: W605

    return np.asanyarray(
        [
            np.log2(np.mean(line, dtype=dtype), dtype=dtype) if len(line) != 0 else 0
            for line in intervals
        ],
        dtype=dtype,
    )