Skip to content

Synamic Phase-Locking

dPL(N, Tmax, phases)

Dynamic Phase-Locking.

This fuction returns the dynamic Phase-Locking for all parcels/voxels of the input.

Parameters:

Name Type Description Default
N int

Number of parcels/voxels of the input array.

 required
Tmax int

BOLD signal samples count.

 required
phases ndarray

Phases signal array for all parcels/voxels in the format [N, Tmax].

 required

Returns:

Name Type Description
tuple

syncConnAux : Synchronicity matrix for all parcels/voxels in the format [N, N], leidaArrayAux : Leading eigenvector of synchornicity matrix [Tmax, N].

References

.. [1] Cabral, J. et al. (2017). ‘Cognitive performance in healthy older adults relates to spontaneous switching between states of functional connectivity during rest’, Scientific Reports. Nature Publishing Group, 7(1), p. 5135. doi: 10.1038/s41598-017-05425-7.

.. [2] Lord et al,. (2019). Dynamical exploration of the repertoire of brain networks at rest is modulated by psilocybin. NeuroImage, 199(April), 127–142. https://doi.org/10.1016/j.neuroimage.2019.05.060

Source code in dynfc/dPL.py 
 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
def dPL(N, Tmax, phases):
    """Dynamic Phase-Locking.

    This fuction returns the dynamic Phase-Locking for all parcels/voxels
    of the input.

    Args:
        N (int): Number of parcels/voxels of the input array.
        Tmax (int): BOLD signal samples count.
        phases (ndarray): Phases signal array for all parcels/voxels in the format [N, Tmax].

    Returns:
        tuple: 
            syncConnAux : Synchronicity matrix for all parcels/voxels in the format [N, N], 
            leidaArrayAux : Leading eigenvector of synchornicity matrix [Tmax, N].

    References
    ----------

    .. [1] 
    Cabral, J. et al. (2017). ‘Cognitive performance in healthy 
    older adults relates to spontaneous switching between states 
    of functional connectivity during rest’, Scientific Reports. 
    Nature Publishing Group, 7(1), p. 5135. 
    doi: 10.1038/s41598-017-05425-7.

    .. [2] 
    Lord et al,. (2019). Dynamical exploration of the 
    repertoire of brain networks at rest is 
    modulated by psilocybin. NeuroImage, 199(April), 127–142. 
    https://doi.org/10.1016/j.neuroimage.2019.05.060

    """

    T = arange(10, Tmax - 10)

    syncConnAux = zeros([N, N, Tmax-20])
    syncAux = zeros([len(T), 1])

    for t in range(0, len(T)):

        for j in range(0, N):

            for k in range(0, j + 1):

                syncConnAux[j, k, t] = phDiff(phases[j, T[t]],
                                                       phases[k, T[t]])
                syncConnAux[k, j, t] = syncConnAux[j, k, t]

    leidaArrayAux = get_LEiDA(syncConnAux)

    return syncConnAux, leidaArrayAux