<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD with MathML3 v1.1d2 20140930//EN" "JATS-journalpublishing1-mathml3.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="1.1d2" xml:lang="en">
  <front>
    <journal-meta>
      <journal-id journal-id-type="nlm-ta">BISH</journal-id>
      <journal-id journal-id-type="publisher-id">ICCK</journal-id>
      <journal-title-group>
        <journal-title>Biomedical Informatics and Smart Healthcare</journal-title>
      </journal-title-group>
      <issn pub-type="ppub" publication-format="print">pending</issn>
      <issn pub-type="epub" publication-format="electronic">pending</issn>
      <publisher>
        <publisher-name>Institute of Central Computation and Knowledge Inc</publisher-name>
        <publisher-loc>522 W RIVERSIDE AVE STE N, SPOKANE, WA, 99201, UNITED STATES</publisher-loc>
      </publisher>
    </journal-meta>
    <article-meta>
      <article-id pub-id-type="doi">10.62762/BISH.2025.346171</article-id>
      <article-categories>
        <subj-group subj-group-type="heading">
          <subject>Research Article</subject>
        </subj-group>
      </article-categories>
      <title-group>
        <article-title>Exploring the Temporal Dynamics and Causal Interactions Between the Amygdala and vmPFC: A Multidimensional Approach Using Time-Lagged Correlation, Granger Causality, and Entropy Analysis</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <name>
            <surname>Su</surname>
            <given-names>Wynn</given-names>
          </name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <contrib-id contrib-id-type="orcid">https://orcid.org/0000-0002-7457-0060</contrib-id>
          <name>
            <surname>Mishra</surname>
            <given-names>Satyam</given-names>
          </name>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <contrib-id contrib-id-type="orcid">https://orcid.org/0000-0002-9612-9073</contrib-id>
          <name>
            <surname>Bijalwan</surname>
            <given-names>Anchit</given-names>
          </name>
          <xref ref-type="aff" rid="aff3">3</xref>
        </contrib>
        <aff id="aff1"><label>1</label>Department of Science, Concordia International School, Hanoi, Vietnam</aff>
        <aff id="aff2"><label>2</label>Faculty of Engineering, Vietnam National University, Hanoi, Vietnam</aff>
        <aff id="aff3"><label>3</label>Computing, Games, and Data Science Academic Research Cluster (CGDS ARC), British University Vietnam, Hanoi, Vietnam</aff>
      </contrib-group>
      <author-notes>
        <corresp id="cor1">Corresponding Author: Wynn Su. Email: <email>wynnsu9487@gmail.com</email></corresp>
      </author-notes>
      <pub-date date-type="pub" pub-type="epub" publication-format="online">
        <day>24</day>
        <month>6</month>
        <year>2025</year>
      </pub-date>
      <volume>1</volume>
      <issue>1</issue>
      <fpage>27</fpage>
      <lpage>34</lpage>
      <history>
        <date date-type="received">
          <day>21</day>
          <month>4</month>
          <year>2025</year>
        </date>
        <date date-type="accepted">
          <day>07</day>
          <month>5</month>
          <year>2025</year>
        </date>
      </history>
      <permissions>
        <copyright-statement>© 2025 by the Authors. Published by Institute of Central Computation and Knowledge. This is an open access article under the CC BY license (https://creativecommons.org/licenses/by/4.0/).</copyright-statement>
        <license xlink:href="https://creativecommons.org/licenses/by/4.0/">
        <license-p>This work is licensed under a <ext-link ext-link-type="uri" xlink:type="simple" xlink:href="https://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution 4.0 International License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</license-p>
        </license>
      </permissions>
      <self-uri xlink:href="https://www.icck.org/article/abs/bish.2025.346171">This article is available from https://www.icck.org/article/abs/bish.2025.346171</self-uri>
      <abstract>
        <p>Understanding the dynamic interaction between the amygdala and the ventromedial prefrontal cortex (vmPFC) is essential for unraveling the neural mechanisms underlying emotion regulation. This study introduces a multidimensional analytical framework that integrates functional connectivity, time-lagged correlation, Granger causality analysis (GCA), and entropy-based complexity metrics to explore the amygdala–vmPFC relationship during emotionally aversive tasks using intracranial EEG (iEEG) data. Our findings reveal a significant negative functional correlation (<inline-formula><mml:math alttext="r=-0.009" display="inline"><mml:mrow><mml:mi>r</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mo>−</mml:mo><mml:mn>0.009</mml:mn></mml:mrow></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math alttext="p&lt;0.05" display="inline"><mml:mrow><mml:mi>p</mml:mi><mml:mo>&lt;</mml:mo><mml:mn>0.05</mml:mn></mml:mrow></mml:math></inline-formula>) between the amygdala and vmPFC, indicating an inhibitory relationship. Time-lagged correlation analysis further uncovers a temporal delay ( 30 time points), suggesting that amygdala activity precedes vmPFC modulation. Entropy analysis shows that the vmPFC exhibits higher signal complexity (H = 4.26) than the amygdala (H = 3.98), consistent with its role in higher-order emotional regulation. However, GCA results yielded no statistically significant directional influence, highlighting the non-linear and multifaceted nature of these interactions. This study is among the first to combine entropy, time-lag, and causality metrics for analyzing vmPFC–amygdala dynamics in human iEEG. Our integrated framework offers deeper insights into brain region interplay and lays the groundwork for future research in affective neuroscience and neuroeconomics.</p>
      </abstract>
      <kwd-group kwd-group-type="author" xml:lang="en">
        <kwd>amygdala-vmPFC connectivity</kwd>
        <kwd>time-lagged correlation</kwd>
        <kwd>granger causality</kwd>
        <kwd>EEG signal analysis</kwd>
        <kwd>entropy metrics</kwd>
        <kwd>neural dynamics</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="S1">
      <label>1.</label>
      <title>Introduction</title>
      <sec id="S1.SS1">
        <label>1.1</label>
        <title>Background</title>
        <p id="S1.SS1.p1">Research has identified a negative correlation between amygdala and ventromedial prefrontal cortex (vmPFC) functional connectivity, which plays a crucial role in emotion regulation. This relationship is influenced by factors such as inflammation, development, and neurofeedback training. Mehta et al. [<xref rid="ref001" ref-type="bibr">1</xref>] found that increased inflammation predicted decreased amygdala-vmPFC connectivity in depressed patients, correlating with anxiety symptoms. Motzkin et al. [<xref rid="ref002" ref-type="bibr">2</xref>] demonstrated that vmPFC lesions led to heightened amygdala activity, supporting its regulatory role. Gee et al. [<xref rid="ref003" ref-type="bibr">3</xref>] observed a developmental shift from positive to negative amygdala-prefrontal connectivity during the transition to adolescence, potentially explaining improvements in emotional regulation. Paret et al. [<xref rid="ref004" ref-type="bibr">4</xref>] showed that voluntary amygdala down-regulation through fMRI neurofeedback increased amygdala-vmPFC connectivity, with a bottom-up information flow from amygdala to vmPFC. These findings highlight the complex interplay between amygdala and vmPFC in emotional processing and regulation, offering insights into potential therapeutic approaches for mood and anxiety disorders.</p>
      </sec>
      <sec id="S1.SS2">
        <label>1.2</label>
        <title>Literature Review</title>
        <p id="S1.SS2.p1">Understanding the neural dynamics among brain regions has provided neuroscientists with valuable insights into constructing a robust model for how the brain regulates emotion and makes economic decisions [<xref rid="ref005" ref-type="bibr">5</xref>, <xref rid="ref006" ref-type="bibr">6</xref>]. Such study is especially crucial in the fields of affective neuroscience and neuroeconomics, notably the investigation of the relationship between vmPFC and amygdala in rational adults [<xref rid="ref007" ref-type="bibr">7</xref>, <xref rid="ref008" ref-type="bibr">8</xref>, <xref rid="ref002" ref-type="bibr">2</xref>, <xref rid="ref009" ref-type="bibr">9</xref>, <xref rid="ref010" ref-type="bibr">10</xref>]. These literatures have shown a negative functional connectivity between these two regions, indicating a negative temporal correlation across the two distant regions. However, correlation in neuroscience alone does not necessarily prove causation [<xref rid="ref011" ref-type="bibr">11</xref>]. Most of these studies primarily rely on using static correlation, which only records activity at a specific moment in time, to identify the functional connectivity between regions.</p>
        <p id="S1.SS2.p2">Other studies, however, have attempted to explore causality between vmPFC and amygdala using Granger causality analysis (GCA). For example, Seth et al. [<xref rid="ref013" ref-type="bibr">13</xref>] used GCA to interpret Functional Magnetic Resonance Imaging (fMRI) data, revealing its potential for inferring causal interaction in neural networks. Still, these analysis methods rely on low-temporal, high-spatial resolution data from fMRI, which might undermine the rapid, dynamic causal interplay between regions (in this case, vmPFC and amygdala) involved in emotion regulation and economic decision-making [<xref rid="ref012" ref-type="bibr">12</xref>].</p>
        <p id="S1.SS2.p3">As Dosenbach et al. [<xref rid="ref014" ref-type="bibr">14</xref>] suggested, a multidimensional approach is essential to robustly capture the complexity of neural interactions. Yet, most studies still tend to utilize one or two methods for analysis. For instance, Gold et al. [<xref rid="ref015" ref-type="bibr">15</xref>] concentrated on connectivity but failed to examine signal complexity or causal relationships between vmPFC and amygdala, limiting their understanding of the vmPFC-amygdala interaction. Hence, in this paper, a multidimensional approach is utilized to study amygdala-vmPFC interaction from EEG data.</p>
        <p id="S1.SS2.p4">First, we move beyond static connectivity by applying time-lagged correlation to explore how these regions influence each other over time, a technique not commonly used in previous research. Second, we employ GCA to investigate predictive relationships, something often missing in EEG-based studies of these regions. Third, the entropy metric is implemented into our analysis to offer insights into the complexity of signals in both regions. This provides a unique and dynamic perspective on how the amygdala and vmPFC might handle different levels of information processing during tasks involving emotion regulation or decision-making. Lastly, by combining these techniques, we offer a holistic analysis that moves beyond what single-method studies have achieved, filling critical gaps in literature.</p>
        <p id="S1.SS2.p5">In other words, the following questions are addressed within this paper:</p>
        <p>
          <list list-type="order" id="S1.I1">
            <list-item id="S1.I1.i1">
              <p id="S1.I1.i1.p1">To what extent are the amygdala and vmPFC functionally connected?</p>
            </list-item>
            <list-item id="S1.I1.i2">
              <p id="S1.I1.i2.p1">Does the interaction between amygdala and vmPFC abide with a time-lagged pattern?</p>
            </list-item>
            <list-item id="S1.I1.i3">
              <p id="S1.I1.i3.p1">Which brain region displays more complex activities?</p>
            </list-item>
            <list-item id="S1.I1.i4">
              <p id="S1.I1.i4.p1">Can neural activity of one region predict that of another region?</p>
            </list-item>
          </list>
        </p>
        <p id="S1.SS2.p7">In essence, our solution not only maps the connectivity between the amygdala and vmPFC but also explores the temporal dynamics, predictability, and complexity of their interactions in a way that previous studies have largely overlooked.</p>
      </sec>
    </sec>
    <sec id="S2">
      <label>2.</label>
      <title>Methodology</title>
      <sec id="S2.SS1">
        <label>2.1</label>
        <title>Dataset Collection</title>
        <p id="S2.SS1.p1">The dataset used in this study comprises intracranial EEG (iEEG) and neuronal data recorded from the amygdalae of nine human subjects during exposure to emotionally aversive and neutral visual stimuli [<xref rid="ref016" ref-type="bibr">16</xref>]. Data collection was performed in a controlled clinical setting, and the dataset was obtained from OpenNeuro in the BIDS standard format.</p>
      </sec>
      <sec id="S2.SS2">
        <label>2.2</label>
        <title>Dataset Analysis</title>
        <p id="S2.SS2.p1">The electrophysiological data include:</p>
        <p>
          <list list-type="bullet" id="S2.I1">
            <list-item id="S2.I1.i1">
              <p id="S2.I1.i1.p1"><bold>Intracranial EEG (iEEG)</bold>: High temporal resolution signals from electrodes implanted in the amygdala. Spectral analysis validated the iEEG signals.</p>
            </list-item>
            <list-item id="S2.I1.i2">
              <p id="S2.I1.i2.p1"><bold>Neuronal Spike Times and Waveforms</bold>: Neuronal firing data, including spike sorting quality metrics.</p>
            </list-item>
            <list-item id="S2.I1.i3">
              <p id="S2.I1.i3.p1"><bold>Metadata</bold>: Task-related information, subject details, session information, and electrode specifications.</p>
            </list-item>
          </list>
        </p>
        <p id="S2.SS2.p3">This dataset is publicly available through OpenNeuro under the title "Dataset of neurons and intracranial EEG from human amygdala during aversive dynamic visual stimulation" [<xref rid="ref016" ref-type="bibr">16</xref>]. The dataset underwent rigorous technical validation, including:</p>
        <p>
          <list list-type="bullet" id="S2.I2">
            <list-item id="S2.I2.i1">
              <p id="S2.I2.i1.p1">Spike sorting quality assessment for neuronal firing data.</p>
            </list-item>
            <list-item id="S2.I2.i2">
              <p id="S2.I2.i2.p1">Spectral analysis of iEEG signals to ensure data integrity.</p>
            </list-item>
          </list>
        </p>
        <p id="S2.SS2.p4">This comprehensive dataset provides a unique opportunity to study the amygdala's role in processing emotional stimuli at both macroscopic and microscopic levels. It was instrumental in analyzing the functional connectivity and causal interactions between the amygdala and vmPFC in the present study.</p>
      </sec>
      <sec id="S2.SS3">
        <label>2.3</label>
        <title>Functional Connectivity and Correlation Analysis</title>
        <p id="S2.SS3.p1">Functional connectivity demonstrates the temporal relationship between spatially distant brain regions. Functional connectivity, which is usually determined through EEG and fMRI [<xref rid="ref017" ref-type="bibr">17</xref>]. In our study, we used the Pearson correlation coefficient to examine functional connectivity between the neural activities of amygdala and vmPFC, as it indicates how synchronized these two regions are. However, Rousselet and Pernet [<xref rid="ref020" ref-type="bibr">20</xref>] have pointed out that the Pearson correlation might be an unrobust metric due to its high sensitivity to outliers. This indicates that relying solely on Pearson correlation can yield unreliable conclusions for understanding the relation between brain and behavior. Despite this limitation, we included the Pearson correlation as an initial measure, utilizing it as a preliminary reference to support further validation of functional connectivity between the amygdala and vmPFC. This approach allows us to capture basic connectivity patterns, potentially helping future research to refine more robust analytical methods.</p>
      </sec>
      <sec id="S2.SS4">
        <label>2.4</label>
        <title>Time-Lagged Correlation</title>
        <p id="S2.SS4.p1">Time-lagged correlation measures how one variable temporally influences another, allowing us to examine the stability of the Pearson correlation coefficient in representing functional connectivity. In this study, time-lagged correlation is used to test the robustness of the Pearson correlation by assessing whether the observed causality between the amygdala and vmPFC is significant over time. While this approach can hint at causation across functionally connected regions, it falls short of providing conclusive evidence of temporal causation. According to Rykhlevskaia et al. [<xref rid="ref019" ref-type="bibr">19</xref>], time-lagged correlation does not provide adequate information to prove the existence of such temporal causation. To address this drawback, we used Granger Causality Analysis (GCA) to corroborate the directional influences between the amygdala and vmPFC.</p>
      </sec>
      <sec id="S2.SS5">
        <label>2.5</label>
        <title>Granger Causality Analysis (GCA)</title>
        <p id="S2.SS5.p1">Granger Causality Analysis (GCA) is a statistical method used to evaluate the predictive relationship between two time-series variables by assessing whether the past activity of one variable can predict the future activity of another. In this study, GCA was selected as the primary method for exploring the directional influence between the amygdala and the ventromedial prefrontal cortex (vmPFC). By leveraging the high temporal resolution of intracranial EEG (iEEG) signals, GCA provides insights into how these regions dynamically influence each other during emotionally aversive tasks [<xref rid="ref013" ref-type="bibr">13</xref>].</p>
        <boxed-text id="S2.SS5.SSS0.Px1">
          <caption>
            <p>Why GCA Over Other Statistical Tests?</p>
          </caption>
          <p id="S2.SS5.SSS0.Px1.p1">While several statistical methods exist for analyzing relationships between time-series data, GCA offers unique advantages that align with our research objectives:</p>
          <p>
            <list list-type="bullet" id="S2.I3">
              <list-item id="S2.I3.i1">
                <p id="S2.I3.i1.p1"><bold>SSR-Based F-Test</bold>: Evaluates whether adding the past activity of one region significantly improves the prediction of another's future activity. While useful, it primarily identifies reduction in residual error but does not explicitly model temporal influence.</p>
              </list-item>
              <list-item id="S2.I3.i2">
                <p id="S2.I3.i2.p1"><bold>SSR-Based Chi-Square Test</bold>: Measures variance explained by including lagged predictors but lacks the time-lag-specific modeling offered by GCA.</p>
              </list-item>
              <list-item id="S2.I3.i3">
                <p id="S2.I3.i3.p1"><bold>Likelihood Ratio Test</bold>: Compares the likelihood of the data under two models (with and without lagged predictors) but does not address the directionality of influence.</p>
              </list-item>
              <list-item id="S2.I3.i4">
                <p id="S2.I3.i4.p1"><bold>Parameter F-Test</bold>: Assesses individual regression coefficients for significance but is less suited for identifying overall causal patterns in a time-lagged context.</p>
              </list-item>
            </list>
          </p>
          <p id="S2.SS5.SSS0.Px1.p3">Although these methods provide robust tools for testing relationships, they are either limited in their ability to model temporal directionality or focus narrowly on specific aspects of causality. GCA, in contrast, integrates temporal lags directly into its framework and evaluates the directional predictive power between time-series variables, making it particularly relevant for examining the dynamic interplay between the amygdala and vmPFC.</p>
          <p id="S2.SS5.SSS0.Px1.p4">Using GCA, we tested whether including the past activity of one region significantly improved predictions of the other's future activity. This was achieved through statistical measures like the F-test and p-values, where a significant p-value (e.g., below 0.05) would indicate a meaningful Granger causal relationship.</p>
        </boxed-text>
      </sec>
      <sec id="S2.SS6">
        <label>2.6</label>
        <title>Entropy Metrics</title>
        <p id="S2.SS6.p1">To quantify the complexity and randomness of neural activity in the amygdala and vmPFC, we employed Shannon Entropy, a widely used metric for assessing signal variability. Based on Carhart-Harris's entropy brain hypothesis, entropy provides a mathematical framework to evaluate the degree of randomness within a system, which in this context reflects the unpredictability of neural signals in each brain region.</p>
        <boxed-text id="S2.SS6.SSS0.Px1">
          <caption>
            <p>Parameters and Calculation:</p>
          </caption>
          <p>
            <list list-type="bullet" id="S2.I4">
              <list-item id="S2.I4.i1">
                <p id="S2.I4.i1.p1"><bold>Signal Input</bold>: The entropy calculation was performed on pre-processed iEEG signals recorded from the amygdala and vmPFC. Pre-processing steps included noise and artifact removal.</p>
              </list-item>
              <list-item id="S2.I4.i2">
                <p id="S2.I4.i2.p1"><bold>Time Windows</bold>: Entropy was computed using a sliding window approach (e.g., 500 ms) to capture temporal variations.</p>
              </list-item>
              <list-item id="S2.I4.i3">
                <p id="S2.I4.i3.p1"><bold>Probability Distribution</bold>: Signal amplitudes within each time window were discretized into bins. The probability <inline-formula><mml:math alttext="p_{i}" display="inline"><mml:msub><mml:mi>p</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:math></inline-formula> for each bin was calculated.</p>
              </list-item>
              <list-item id="S2.I4.i4">
                <p id="S2.I4.i4.p1"><bold>Shannon Entropy Formula</bold>:</p>
                <p>
                  <disp-formula id="S2.Ex1">
                    <mml:math alttext="H=-\sum_{i=1}^{k}p_{i}\log_{2}(p_{i})" display="block">
                      <mml:mrow>
                        <mml:mi>H</mml:mi>
                        <mml:mo>=</mml:mo>
                        <mml:mrow>
                          <mml:mo>−</mml:mo>
                          <mml:mrow>
                            <mml:munderover>
                              <mml:mo movablelimits="false">∑</mml:mo>
                              <mml:mrow>
                                <mml:mi>i</mml:mi>
                                <mml:mo>=</mml:mo>
                                <mml:mn>1</mml:mn>
                              </mml:mrow>
                              <mml:mi>k</mml:mi>
                            </mml:munderover>
                            <mml:mrow>
                              <mml:msub>
                                <mml:mi>p</mml:mi>
                                <mml:mi>i</mml:mi>
                              </mml:msub>
                              <mml:mo lspace="0.167em">⁢</mml:mo>
                              <mml:mrow>
                                <mml:msub>
                                  <mml:mi>log</mml:mi>
                                  <mml:mn>2</mml:mn>
                                </mml:msub>
                                <mml:mo>⁡</mml:mo>
                                <mml:mrow>
                                  <mml:mo stretchy="false">(</mml:mo>
                                  <mml:msub>
                                    <mml:mi>p</mml:mi>
                                    <mml:mi>i</mml:mi>
                                  </mml:msub>
                                  <mml:mo stretchy="false">)</mml:mo>
                                </mml:mrow>
                              </mml:mrow>
                            </mml:mrow>
                          </mml:mrow>
                        </mml:mrow>
                      </mml:mrow>
                    </mml:math>
                  </disp-formula>
                </p>
                <p>where <inline-formula><mml:math alttext="H" display="inline"><mml:mi>H</mml:mi></mml:math></inline-formula> is the entropy, <inline-formula><mml:math alttext="p_{i}" display="inline"><mml:msub><mml:mi>p</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:math></inline-formula> is the probability of the <inline-formula><mml:math alttext="i" display="inline"><mml:mi>i</mml:mi></mml:math></inline-formula>-th amplitude bin, and <inline-formula><mml:math alttext="k" display="inline"><mml:mi>k</mml:mi></mml:math></inline-formula> is the number of bins.</p>
              </list-item>
            </list>
          </p>
        </boxed-text>
        <boxed-text id="S2.SS6.SSS0.Px2">
          <caption>
            <p>Outputs:</p>
          </caption>
          <p>
            <list list-type="bullet" id="S2.I5">
              <list-item id="S2.I5.i1">
                <p id="S2.I5.i1.p1"><bold>Regional Entropy</bold>: Separate entropy values were computed for the amygdala and vmPFC.</p>
              </list-item>
              <list-item id="S2.I5.i2">
                <p id="S2.I5.i2.p1"><bold>Temporal Trends</bold>: Changes in entropy over time were analyzed to assess variability.</p>
              </list-item>
            </list>
          </p>
          <p id="S2.SS6.SSS0.Px2.p2">Entropy metrics provided critical insights into the functional dynamics between the amygdala and vmPFC:</p>
          <p>
            <list list-type="bullet" id="S2.I6">
              <list-item id="S2.I6.i1">
                <p id="S2.I6.i1.p1"><bold>Complexity Analysis</bold>: Comparing entropy across regions highlighted differences in signal unpredictability.</p>
              </list-item>
              <list-item id="S2.I6.i2">
                <p id="S2.I6.i2.p1"><bold>Dynamic Interactions</bold>: Aligned temporal entropy patterns helped infer excitatory or inhibitory interactions.</p>
              </list-item>
              <list-item id="S2.I6.i3">
                <p id="S2.I6.i3.p1"><bold>Causal Insights</bold>: Synchrony in entropy fluctuations hinted at potential directional influence.</p>
              </list-item>
            </list>
          </p>
          <p id="S2.SS6.SSS0.Px2.p3">Incorporating entropy into our methodology allowed us to go beyond conventional connectivity metrics, offering a nuanced view of how the amygdala and vmPFC manage information under emotional stimuli. Combined with time-lagged correlation and GCA, this creates a multidimensional framework for analyzing brain region interactions.</p>
        </boxed-text>
      </sec>
    </sec>
    <sec id="S3">
      <label>3.</label>
      <title>Results</title>
      <sec id="S3.SS1">
        <label>3.1</label>
        <title>Functional Connectivity and Correlational Analysis</title>
        <p id="S3.SS1.p1">The Pearson correlation analysis revealed a significant negative correlation of <inline-formula><mml:math alttext="r=-0.009" display="inline"><mml:mrow><mml:mi>r</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mo>−</mml:mo><mml:mn>0.009</mml:mn></mml:mrow></mml:mrow></mml:math></inline-formula> (<inline-formula><mml:math alttext="p" display="inline"><mml:mi>p</mml:mi></mml:math></inline-formula>-value &lt; 0.05) between the amygdala and the vmPFC, as shown in Figure <xref ref-type="fig" rid="F1">1</xref>. This suggests an inverse relationship: when one region becomes more active, the other tends to become less active. Such an inhibitory interaction supports the vmPFC's role in regulating the amygdala's emotional reactivity.</p>
        <p id="S3.SS1.p2">This negative correlation confirms functional connectivity between the two regions, which concludes our first research objective successfully, i.e., <italic>Are these regions functionally connected?</italic></p>
        <p>
          <fig id="F1">
            <label>Figure 1.</label>
            <caption>
              <p>Correlation Calculation between Amygdala and vmPFC. The negative correlation suggests that as the activity of one region increases, the other decreases.</p>
            </caption>
            <graphic xlink:href="correlation_analysis.png"/>
          </fig>
        </p>
      </sec>
      <sec id="S3.SS2">
        <label>3.2</label>
        <title>Time-Lagged Correlation</title>
        <p id="S3.SS2.p1">Time-lagged correlation analysis revealed the strongest negative correlation at a lag of approximately <inline-formula><mml:math alttext="-30" display="inline"><mml:mrow><mml:mo>−</mml:mo><mml:mn>30</mml:mn></mml:mrow></mml:math></inline-formula> time points, indicating that the amygdala's activity leads the vmPFC by this delay, as seen in Figure <xref ref-type="fig" rid="F2">2</xref>. This temporal pattern suggests a sequential process, with the amygdala initiating neural responses and the vmPFC following to modulate emotional processing.</p>
        <p id="S3.SS2.p2">This time lag indicates a sequential interaction with the amygdala leading, which concludes our second research objective successfully, i.e., <italic>Do their interactions follow a time-lagged pattern?</italic></p>
        <p>
          <fig id="F2">
            <label>Figure 2.</label>
            <caption>
              <p>Time-lagged Correlation: Plot showing the correlation values at different time lags between the amygdala and vmPFC. The strongest negative correlation occurs at <inline-formula><mml:math alttext="\sim" display="inline"><mml:mo>∼</mml:mo></mml:math></inline-formula>30 time points, indicating a delay in the vmPFC's response to the amygdala.</p>
            </caption>
            <graphic xlink:href="timelagged_correlation.png"/>
          </fig>
        </p>
      </sec>
      <sec id="S3.SS3">
        <label>3.3</label>
        <title>Entropy Metrics</title>
        <p id="S3.SS3.p1">Entropy analysis revealed that the vmPFC had higher entropy (<inline-formula><mml:math alttext="H=4.26" display="inline"><mml:mrow><mml:mi>H</mml:mi><mml:mo>=</mml:mo><mml:mn>4.26</mml:mn></mml:mrow></mml:math></inline-formula>) compared to the amygdala (<inline-formula><mml:math alttext="H=3.98" display="inline"><mml:mrow><mml:mi>H</mml:mi><mml:mo>=</mml:mo><mml:mn>3.98</mml:mn></mml:mrow></mml:math></inline-formula>), indicating more complex and unpredictable activity in the vmPFC, as shown in Figure <xref ref-type="fig" rid="F3">3</xref>. This aligns with its role in integrating higher-order cognitive functions, such as decision-making and emotional regulation.</p>
        <p id="S3.SS3.p2">These outputs conclude that the vmPFC's higher entropy reflects its more complex neural activity, which concludes our third research objective successfully, i.e., <italic>Is one region's activity more complex than the other?</italic></p>
        <p>
          <fig id="F3">
            <label>Figure 3.</label>
            <caption>
              <p>Entropy Values of Amygdala and vmPFC. The vmPFC exhibits higher entropy, reflecting its greater complexity in processing and regulating neural activity.</p>
            </caption>
            <graphic xlink:href="entropy_results.png"/>
          </fig>
        </p>
      </sec>
      <sec id="S3.SS4">
        <label>3.4</label>
        <title>Granger Causality</title>
        <p id="S3.SS4.p1">Granger causality testing did not reveal statistically significant predictive relationships between the amygdala and vmPFC (all <inline-formula><mml:math alttext="p" display="inline"><mml:mi>p</mml:mi></mml:math></inline-formula>-values &gt; 0.05), as shown in Figure <xref ref-type="fig" rid="F4">4</xref>. This suggests that the regions' interactions may not follow a simple linear causal pattern but could involve more complex dynamics influenced by external factors or other brain regions.</p>
        <p id="S3.SS4.p2">These outputs conclude that the weak evidence for predictive influence between the amygdala and vmPFC suggests that their interactions may involve more complex, non-linear dynamics. This addresses our fourth research objective, i.e., <italic>Can we predict one region's activity based on the other?</italic>, by indicating that a straightforward causal relationship is unlikely in this context.</p>
        <p>
          <fig id="F4">
            <label>Figure 4.</label>
            <caption>
              <p>Granger Causality Results: Statistical results of Granger causality tests at different lags. The lack of significant <inline-formula><mml:math alttext="p" display="inline"><mml:mi>p</mml:mi></mml:math></inline-formula>-values indicates weak evidence for direct causal relationships between the amygdala and vmPFC.</p>
            </caption>
            <graphic xlink:href="granger_results.png"/>
          </fig>
        </p>
      </sec>
    </sec>
    <sec id="S4">
      <label>4.</label>
      <title>Discussion</title>
      <p id="S4.p1">This study aimed to clarify the dynamic relationship between the amygdala and the ventromedial prefrontal cortex (vmPFC) in regulating emotions. By employing a multidimensional analysis that integrated functional connectivity, time-lagged correlation, Granger causality, and entropy metrics, we provided a nuanced understanding of how these regions interact during emotionally aversive tasks. Below, we discuss our findings in the context of prior literature and highlight their implications for neuroscience and neuroeconomics.</p>
      <sec id="S4.SS1">
        <label>4.1</label>
        <title>Functional Connectivity and Non-Linear Dynamics</title>
        <p id="S4.SS1.p1">Consistent with previous studies [<xref rid="ref007" ref-type="bibr">7</xref>, <xref rid="ref021" ref-type="bibr">21</xref>, <xref rid="ref023" ref-type="bibr">23</xref>], our results confirmed a significant functional connection between the amygdala and vmPFC, evidenced by a statistically significant negative correlation. This aligns with the established role of the vmPFC in inhibiting amygdala-driven emotional reactivity. However, Pearson correlation, which only captures linear relationships, may oversimplify the complexity of brain interactions, as neural communication often involves non-linear dynamics [<xref rid="ref018" ref-type="bibr">18</xref>, <xref rid="ref024" ref-type="bibr">24</xref>]. This limitation underscores the need for advanced methods that account for such complexity, as functional connectivity alone cannot fully explain the causal dynamics between brain regions.</p>
      </sec>
      <sec id="S4.SS2">
        <label>4.2</label>
        <title>Time-Lagged Correlation and Causality</title>
        <p id="S4.SS2.p1">The time-lagged correlation analysis revealed a temporal delay of approximately <inline-formula><mml:math alttext="-30" display="inline"><mml:mrow><mml:mo>−</mml:mo><mml:mn>30</mml:mn></mml:mrow></mml:math></inline-formula> to <inline-formula><mml:math alttext="-40" display="inline"><mml:mrow><mml:mo>−</mml:mo><mml:mn>40</mml:mn></mml:mrow></mml:math></inline-formula> time points, with the amygdala leading and the vmPFC following. This finding supports the hypothesis that the amygdala initiates emotional responses, which are subsequently modulated by the vmPFC. However, Granger causality analysis provided weak evidence for a direct causal relationship between the regions, contradicting some previous findings [<xref rid="ref002" ref-type="bibr">2</xref>]. This discrepancy may stem from GCA's inherent limitations, including its reliance on linear assumptions and sensitivity to noise in high-variance neural systems [<xref rid="ref022" ref-type="bibr">22</xref>]. These results highlight the possibility that the amygdala-vmPFC interaction involves more intricate, non-linear dynamics than can be captured by GCA alone.</p>
      </sec>
      <sec id="S4.SS3">
        <label>4.3</label>
        <title>A Novel Approach to Causality Analysis</title>
        <p id="S4.SS3.p1">Entropy metrics introduced a novel dimension to our understanding of the amygdala-vmPFC relationship. The vmPFC exhibited higher entropy values compared to the amygdala, reflecting its greater complexity and integration of diverse neural inputs. This finding complements prior studies showing that the vmPFC's activity correlates with emotional regulation [<xref rid="ref025" ref-type="bibr">25</xref>]. By combining entropy metrics with time-lagged correlation and functional connectivity, our analysis suggests a regulatory mechanism where the vmPFC modulates amygdala activity over time. This integrated approach offers a robust framework for exploring causality in neural systems and addresses the limitations of linear causality methods like GCA.</p>
      </sec>
      <sec id="S4.SS4">
        <label>4.4</label>
        <title>Limitations and Future Directions</title>
        <p id="S4.SS4.p1">While our study achieved significant and novel insights into the dynamics of the amygdala-vmPFC relationship, there remain opportunities to expand upon this work. The integration of functional connectivity, time-lagged correlation, Granger causality, and entropy metrics provides a powerful and innovative framework that advances the understanding of brain region interactions. However, the complexities of neural systems and the limitations of current methodologies invite further exploration to build upon our novel findings.</p>
      </sec>
    </sec>
    <sec id="S5">
      <label>5.</label>
      <title>Conclusion</title>
      <p id="S5.p1">This study investigated the dynamic neural relationship between the amygdala and the ventromedial prefrontal cortex (vmPFC) during emotionally aversive tasks using a multidimensional analysis framework. We confirmed a statistically significant negative functional correlation (<inline-formula><mml:math alttext="r=-0.009" display="inline"><mml:mrow><mml:mi>r</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mo>−</mml:mo><mml:mn>0.009</mml:mn></mml:mrow></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math alttext="p&lt;0.05" display="inline"><mml:mrow><mml:mi>p</mml:mi><mml:mo>&lt;</mml:mo><mml:mn>0.05</mml:mn></mml:mrow></mml:math></inline-formula>), indicating inhibitory connectivity. Time-lagged correlation analysis revealed that amygdala activity precedes vmPFC response by approximately 30 time points, suggesting sequential emotional regulation. Entropy analysis showed that the vmPFC exhibited greater signal complexity (H = 4.26) than the amygdala (H = 3.98), supporting its role in higher-order cognitive processing. Granger causality analysis, however, did not detect a significant predictive relationship, implying the presence of more complex or nonlinear dynamics.</p>
      <p id="S5.p2">
        <bold>Key Contributions:</bold>
      </p>
      <p>
        <list list-type="bullet" id="S5.I1">
          <list-item id="S5.I1.i1">
            <p id="S5.I1.i1.p1">Proposed a novel, multidimensional methodology combining functional connectivity, time-lagged correlation, Granger causality, and entropy metrics on human iEEG data.</p>
          </list-item>
          <list-item id="S5.I1.i2">
            <p id="S5.I1.i2.p1">Demonstrated that vmPFC–amygdala interactions are temporally ordered and functionally connected, but not linearly causal.</p>
          </list-item>
          <list-item id="S5.I1.i3">
            <p id="S5.I1.i3.p1">Introduced entropy-based complexity analysis as a valuable metric to complement traditional connectivity and causality tools.</p>
          </list-item>
          <list-item id="S5.I1.i4">
            <p id="S5.I1.i4.p1">Provided a reproducible framework for studying emotion-related neural circuits with implications for neuroscience, affective computing, and neuroeconomics.</p>
          </list-item>
        </list>
      </p>
      <p id="S5.p3">These findings offer a richer understanding of the amygdala-vmPFC axis and emphasize the need to move beyond linear models in neuroscience. Future work may explore nonlinear causal inference techniques, integrate multimodal data, or apply this framework to clinical populations to assess psychiatric or decision-making disorders.</p>
    </sec>
  </body>
  <back>
    <ack>
      <title>Acknowledgments</title>
      <p id="ack.p1">This work was supported without any funding.</p>
    </ack>
    <sec id="sec0100" sec-type="COI-statement">
      <title>Conflict of interest</title>
      <p>The authors declare no conflicts of interest.</p>
    </sec>
    <ref-list>
      <title>References</title>
      <ref id="ref001">
        <label>[1]</label>
        <mixed-citation> Mehta, N. D., Haroon, E., Xu, X., Woolwine, B. J., Li, Z., &amp; Felger, J. C. (2018). Inflammation negatively correlates with amygdala-ventromedial prefrontal functional connectivity in association with anxiety in patients with depression: Preliminary results. <italic>Brain, behavior, and immunity, 73</italic>, 725-730. [<uri>https://doi.org/10.1016/j.bbi.2018.07.026</uri>] </mixed-citation>
      </ref>
      <ref id="ref002">
        <label>[2]</label>
        <mixed-citation> Motzkin, J. C., Philippi, C. L., Wolf, R. C., Baskaya, M. K., &amp; Koenigs, M. (2015). Ventromedial prefrontal cortex is critical for the regulation of amygdala activity in humans. <italic>Biological Psychiatry, 77</italic>(3), 276–284. [<uri>https://doi.org/10.1016/j.biopsych.2014.02.014</uri>] </mixed-citation>
      </ref>
      <ref id="ref003">
        <label>[3]</label>
        <mixed-citation> Gee, D. G., Humphreys, K. L., Flannery, J., Goff, B., Telzer, E. H., Shapiro, M., … &amp; Tottenham, N. (2013). A developmental shift from positive to negative connectivity in human amygdala–prefrontal circuitry. <italic>Journal of Neuroscience, 33</italic>(10), 4584–4593. [<uri>https://doi.org/10.1523/JNEUROSCI.3446-12.2013</uri>] </mixed-citation>
      </ref>
      <ref id="ref004">
        <label>[4]</label>
        <mixed-citation> Paret, C., Ruf, M., Gerchen, M. F., Kluetsch, R., Demirakca, T., Jungkunz, M., … &amp; Schmahl, C. (2016). fMRI neurofeedback of amygdala response to aversive stimuli enhances prefrontal–limbic brain connectivity. <italic>NeuroImage, 125</italic>, 182–188. [<uri>https://doi.org/10.1016/j.neuroimage.2015.10.027</uri>] </mixed-citation>
      </ref>
      <ref id="ref005">
        <label>[5]</label>
        <mixed-citation> Vuilleumier, P. (2005). How brains beware: Neural mechanisms of emotional attention. <italic>Trends in Cognitive Sciences, 9</italic>(12), 585–594. [<uri>https://doi.org/10.1016/j.tics.2005.10.011</uri>] </mixed-citation>
      </ref>
      <ref id="ref006">
        <label>[6]</label>
        <mixed-citation> Wager, T. D., Davidson, M. L., Hughes, B. L., Lindquist, M. A., &amp; Ochsner, K. N. (2008). Prefrontal–subcortical pathways mediating successful emotion regulation. <italic>Neuron, 59</italic>(6), 1037–1050. [<uri>https://doi.org/10.1016/j.neuron.2008.09.006</uri>] </mixed-citation>
      </ref>
      <ref id="ref007">
        <label>[7]</label>
        <mixed-citation> Etkin, A., Egner, T., &amp; Kalisch, R. (2011). Emotional processing in anterior cingulate and medial prefrontal cortex. <italic>Trends in Cognitive Sciences, 15</italic>(2), 85–93. [<uri>https://doi.org/10.1016/j.tics.2010.11.004</uri>] </mixed-citation>
      </ref>
      <ref id="ref008">
        <label>[8]</label>
        <mixed-citation> Frontiers in Human Neuroscience. (2017). Prefrontal–amygdala connectivity and state anxiety during fear extinction recall in adolescents. <italic>Frontiers in Human Neuroscience, 11</italic>, 587. [<uri>https://doi.org/10.3389/fnhum.2017.00587</uri>] </mixed-citation>
      </ref>
      <ref id="ref009">
        <label>[9]</label>
        <mixed-citation> Phillips, M. L., Drevets, W. C., Rauch, S. L., &amp; Lane, R. (2003). Neurobiology of emotion perception I: The neural basis of normal emotion perception. <italic>Biological Psychiatry, 54</italic>(5), 504–514. [<uri>https://doi.org/10.1016/S0006-3223(03)00168-9</uri>] </mixed-citation>
      </ref>
      <ref id="ref010">
        <label>[10]</label>
        <mixed-citation> Silvers, J. A., Insel, C., Powers, A., Franz, P., Helion, C., Martin, R. E., … &amp; Ochsner, K. N. (2017). vlPFC–vmPFC–amygdala interactions underlie age-related differences in cognitive regulation of emotion. <italic>Cerebral cortex, 27</italic>(7), 3502-3514. [<uri>https://doi.org/10.1093/cercor/bhw073</uri>] </mixed-citation>
      </ref>
      <ref id="ref011">
        <label>[11]</label>
        <mixed-citation> Harnad, S. (2000). <italic>Correlation vs. Causality: How/Why the Mind/Body Problem Is Hard.[Invited Commentary of Humphrey, N." How to Solve the Mind-Body Problem"] Journal of Consciousness Studies 7 (4): 54-61</italic>. </mixed-citation>
      </ref>
      <ref id="ref012">
        <label>[12]</label>
        <mixed-citation> Bressler, S. L., &amp; Seth, A. K. (2011). Wiener–Granger causality: a well established methodology. <italic>Neuroimage, 58</italic>(2), 323-329. [<uri>https://doi.org/10.1016/j.neuroimage.2010.02.059</uri>] </mixed-citation>
      </ref>
      <ref id="ref013">
        <label>[13]</label>
        <mixed-citation> Seth, A. K., Barrett, A. B., &amp; Barnett, L. (2015). Granger causality analysis in neuroscience and neuroimaging. <italic>Journal of Neuroscience, 35</italic>(8), 3293–3297. [<uri>https://doi.org/10.1523/JNEUROSCI.4399-14.2015</uri>] </mixed-citation>
      </ref>
      <ref id="ref014">
        <label>[14]</label>
        <mixed-citation> Dosenbach, N. U. F., Fair, D. A., Miezin, F. M., Cohen, A. L., Wenger, K. K., Dosenbach, R. A. T., … &amp; Petersen, S. E. (2007). Distinct brain networks for adaptive and stable task control in humans. <italic>Proceedings of the National Academy of Sciences, 104</italic>(26), 11073–11078. [<uri>https://doi.org/10.1073/pnas.0704320104</uri>] </mixed-citation>
      </ref>
      <ref id="ref015">
        <label>[15]</label>
        <mixed-citation> Gold, A. L., Shechner, T., Farber, M. J., Spiro, C. N., Leibenluft, E., Pine, D. S., &amp; Britton, J. C. (2016). Amygdala–cortical connectivity: associations with anxiety, development, and threat. <italic>Depression and anxiety, 33</italic>(10), 917-926. [<uri>https://doi.org/10.1002/da.22470</uri>] </mixed-citation>
      </ref>
      <ref id="ref016">
        <label>[16]</label>
        <mixed-citation> Fedele, T., Boran, E., Chirkov, V., Hilfiker, P., Grunwald, T., &amp; Stieglitz, L. (2020). Dataset of neurons and intracranial EEG from human amygdala during aversive dynamic visual stimulation. <italic>OpenNeuro doi, 10</italic>. [<uri>https://doi.org/10.12751/g-node.270z59</uri>] </mixed-citation>
      </ref>
      <ref id="ref017">
        <label>[17]</label>
        <mixed-citation> Feinsod, M., Kreinin, B., Chistyakov, A., &amp; Klein, E. (1998). Preliminary evidence for a beneficial effect of low-frequency, repetitive transcranial magnetic stimulation in patients with major depression and schizophrenia. <italic>Depression and Anxiety, 7</italic>(2), 65–68. [<uri>https://doi.org/10.1002/(SICI)1520-6394(1998)7:2&lt;65::AID-DA2&gt;3.0.CO;2-4</uri>] </mixed-citation>
      </ref>
      <ref id="ref018">
        <label>[18]</label>
        <mixed-citation> Faghiri, A., Stephen, J. M., Wang, Y. P., Wilson, T. W., &amp; Calhoun, V. D. (2019). Brain development includes linear and multiple nonlinear trajectories: a cross-sectional resting-state functional magnetic resonance imaging study. <italic>Brain connectivity, 9</italic>(10), 777-788. [<uri>https://doi.org/10.1089/brain.2018.0641</uri>] </mixed-citation>
      </ref>
      <ref id="ref019">
        <label>[19]</label>
        <mixed-citation> Rykhlevskaia, E., Fabiani, M., &amp; Gratton, G. (2006). Lagged covariance structure models for studying functional connectivity in the brain. <italic>Neuroimage, 30</italic>(4), 1203-1218. [<uri>https://doi.org/10.1016/j.neuroimage.2005.11.019</uri>] </mixed-citation>
      </ref>
      <ref id="ref020">
        <label>[20]</label>
        <mixed-citation> Rousselet, G. A., &amp; Pernet, C. R. (2012). Improving standards in brain–behavior correlation analyses. <italic>Frontiers in Human Neuroscience, 6</italic>, 119. [<uri>https://doi.org/10.3389/fnhum.2012.00119</uri>] </mixed-citation>
      </ref>
      <ref id="ref021">
        <label>[21]</label>
        <mixed-citation> Ghashghaei, H. T., &amp; Barbas, H. (2002). Pathways for emotion: Interactions of prefrontal and anterior temporal pathways in the amygdala of the rhesus monkey. <italic>Neuroscience, 115</italic>(4), 1261–1279. [<uri>https://doi.org/10.1016/S0306-4522(02)00446-3</uri>] </mixed-citation>
      </ref>
      <ref id="ref022">
        <label>[22]</label>
        <mixed-citation> Stokes, P. A., &amp; Purdon, P. L. (2017). A study of problems encountered in Granger causality analysis from a neuroscience perspective. <italic>Proceedings of the national academy of sciences, 114</italic>(34), E7063-E7072. [<uri>https://doi.org/10.1073/pnas.1704663114</uri>] </mixed-citation>
      </ref>
      <ref id="ref023">
        <label>[23]</label>
        <mixed-citation> McDonald, A. J., Mascagni, F., &amp; Guo, L. (1996). Projections of the medial and lateral prefrontal cortices to the amygdala: A Phaseolus vulgaris leucoagglutinin study in the rat. <italic>Neuroscience, 71</italic>(1), 55–75. [<uri>https://doi.org/10.1016/0306-4522(95)00417-3</uri>] </mixed-citation>
      </ref>
      <ref id="ref024">
        <label>[24]</label>
        <mixed-citation> Rusakov, D. A. (2023). A misadventure of the correlation coefficient. <italic>Trends in Neurosciences, 46</italic>(2), 94–96. [<uri>https://doi.org/10.1016/j.tins.2022.09.009</uri>] </mixed-citation>
      </ref>
      <ref id="ref025">
        <label>[25]</label>
        <mixed-citation> Urry, H. L., van Reekum, C. M., Johnstone, T., Kalin, N. H., Thurow, M. E., Schaefer, H. S., … &amp; Davidson, R. J. (2006). Amygdala and ventromedial prefrontal cortex are inversely coupled during regulation of negative affect and predict the diurnal pattern of cortisol secretion among older adults. <italic>Journal of Neuroscience, 26</italic>(16), 4415–4425. [<uri>https://doi.org/10.1523/JNEUROSCI.3215-05.2006</uri>] </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>
