Physiologically Guided Modeling for EEG Multichannel Signals

TL;DR

This paper introduces a new method for analyzing EEG signals that incorporates physiological principles to better understand brain activity.

Contribution

The novel contribution is a physiologically guided model using a parametric formulation to capture EEG signal structure with fewer sources.

Findings

The proposed model explains interchannel EEG correlations using few latent sources.

The framework enables statistical testing of fixed dipole orientation assumptions.

Results show high accuracy in capturing temporal and spatial EEG signal structures.

Abstract

Electroencephalographic (EEG) recordings often exhibit strong interchannel correlations due to scalp potentials reflecting electric fields generated by localized neural sources commonly modeled as current dipoles. Despite this physiological basis, many widely used approaches (e.g., independent component analysis) are largely data-driven, may require many components, and offer limited interpretability and limited support for assessing physiological assumptions. We propose a physiologically guided multichannel model built on a parametric Frequency-Modulated Möbius (FMM) formulation, where each electrode signal is expressed as a linear combination of a small number of latent dipole-related sources with parametrically described trajectories. The resulting framework captures temporal and spatial structure with high accuracy and enables likelihood-based testing of the fixed dipole-orientation assumption through a reduced-rank formulation. Using publicly available EEG data, we illustrate that the proposed approach explains interchannel dependence with few latent sources and provides statistical evidence regarding the plausibility of fixed orientation.

Figures

Full text

1. Introduction

EEGis a non-invasive technique that measures voltage fluctuations generated by synchronized activity of cortical pyramidal neurons using multiple electrodes placed on the scalp [1]. EEG is widely used in clinical and research settings, supporting applications such as epilepsy diagnosis [2], sleep staging [3], anesthesia monitoring [4], and the assessment of neurodegenerative and psychiatric disorders [5,6,7]. Beyond these classical clinical applications, EEG has become increasingly relevant in computational psychiatry and affective computing. In major depressive disorder (MDD), recent studies leverage advanced feature engineering, time-frequency representations, and machine/deep learning methods to improve diagnostic classification and to assess the practical feasibility and limitations of EEG-based biomarkers [8,9,10,11,12]. In parallel, EEG-based emotion recognition has grown rapidly due to its potential in human–computer interaction, with modern approaches combining signal transformations and deep neural architectures, alongside recent reviews synthesizing methodological trends and evaluation practices [13,14,15,16,17].

Physiologically, EEG activity arises when groups of pyramidal neurons become synchronously active. Their combined currents can be approximated at a large scale by an equivalent current dipole, a vector that summarizes the overall strength and direction of the neural activity within a small cortical region [18]. The voltage measured at each electrode is the result of the superposition of the fields produced by many dipoles active across the cortex. Consequently, each EEG channel reflects a mixture of these underlying contributions, forming a multivariate time series whose oscillatory patterns encode the dynamics of the neural generators [19,20,21,22,23]. Accurate modeling is essential for identifying rhythms, tracking temporal changes, and deriving clinically useful biomarkers; this calls for models that link cortical sources to the observed signals [7,21,24,25].

Most existing analytical approaches, such as independent component analysis (ICA) or principal component analysis (PCA), decompose the multichannel signal into a set of data-driven components, e.g., [26,27,28]. Although these methods have proved effective for artifact removal and exploratory analysis, they offer limited physiological interpretability and no explicit link to the dipolar generators that underlie EEG activity; moreover, recent benchmarking studies have shown that reasonable preprocessing choices can influence ERP estimates and downstream conclusions [29,30,31]. In ICA, the most popular approach, the observed channels are modeled as linear mixtures of statistically independent sources, a simple statistical representation that, however, ignores the spatial and biophysical structure. Moreover, in EEG, there is no fixed or predetermined relationship between the location of a scalp electrode and the underlying cortical source or sources that generate the recorded signals. The field of source analysis encompasses a range of methods developed to estimate how many sources are active and where they are located within the brain. However, this area remains surrounded by significant uncertainty and controversy, as these methods rely on a number of debatable assumptions. One particularly problematic assumption concerns the orientation of the dipoles representing the neural sources [20,22,32]. Despite its central role in electrophysiological modeling, there remains considerable debate over whether the orientation of cortical dipoles should be assumed fixed or allowed to vary across time and experimental conditions. To date, this assumption has rarely been examined through formal statistical testing, leaving a long-standing hypothesis in EEG analysis largely unverified [5,18,19,33].

The framework developed in this work addresses this need by providing a physiologically motivated way to represent dipolar neural generators and their temporal evolution, offering a clearer and more interpretable connection between cortical activity and multichannel EEG recordings. We propose a physiologically guided multivariate FMM model as a reduced-rank variant within the FMM framework [34,35]. The model represents scalp potentials as the superposition of a small number of latent FMM sources. Each source corresponds to a cortical dipole whose moment vector has a parametric temporal evolution and projects onto channel-specific directions, producing the observed voltages.

This additive structure reflects the physical laws of electric field propagation. Furthermore, the FMM model is fully parametric, providing a statistically efficient and interpretable representation of multichannel EEG. Its physiologically guided structure, which explicitly incorporates interchannel relationships, yields stable parameter estimation and improves robustness to noise.

Beyond signal representation, the proposed framework allows a formal statistical comparison between models assuming fixed and variable dipole orientations. By embedding this hypothesis within a likelihood-based reduced-rank formulation, we provide a unified framework for estimating dipole dynamics and formally comparing variable- and fixed-orientation models through rank constraints, within a single likelihood-based approach.

Empirical analyses on publicly available datasets demonstrate that our model yields a more parsimonious and physiologically grounded representation of oscillatory activity. This capability bridges statistical modeling and physiological validation in EEG analysis, offering a unified framework for both representation and inference.

The remainder of the paper is organized as follows. Section 2 recalls the biophysical formulation of EEG potentials in terms of dipole position and orientation, reviews the main elements of the FMM model and introduces the concept of reduced-rank multivariate modeling. Section 3 presents the proposed restricted FMM model, describing how physical constraints are incorporated, the estimation algorithm and how the model permits formal validation of the fixed-orientation hypothesis. Section 4 reports empirical analyses using publicly available EEG data, and Section 5 concludes with a discussion of implications and future research directions.

2. Background

2.1. Biophysical Formulation

The electric potentials measured on the scalp by an EEG are generated by the summed activity of large populations of cortical pyramidal neurons, which can be approximated by a set of equivalent current dipoles, e.g., [5,18,19].

The equivalent dipole moment $[eqn]$ represents the net current flow generated by a population of cortical pyramidal neurons within a localized region. Its magnitude and orientation jointly reflect the balance between excitatory and inhibitory synaptic activity and the degree of temporal synchronization among neurons. At rest, when neuronal firing is largely uncorrelated and synaptic currents are statistically balanced, the net dipole moment fluctuates near zero. Following stimulation or spontaneous synchronization, a subset of pyramidal neurons becomes transiently active, generating coherent current flow along their dendritic axes. Consequently, $[eqn]$ increases in magnitude, reaches a peak corresponding to maximal population activation, and subsequently declines as the cortical tissue returns to baseline. The temporal trajectory depends on physiological and experimental conditions.

Traditionally, EEG modeling assumes that the orientation of the dipole remains constant and only its magnitude varies in time [36]. Under this classical model, $[eqn]$ evolves linearly along a fixed direction, representing a localized and stable cortical generator whose activity is purely amplitude modulated. This assumption is supported by the near-parallel organization of pyramidal neurons in a cortical column, which results in a macroscopic current flow along a nearly fixed axis. However, real cortical activity is often more complex and spatially distributed than this idealization suggests [37,38]. Recent work indicates that dipole orientation can vary over time in ways that may reflect underlying physiological dynamics, and the combined activity of nearby cortical regions may produce gradual changes in the effective dipole direction [33]. In such situations, the trajectory of $[eqn]$ can be reasonably approximated by an elliptical path within a plane, representing the smooth evolution of the resultant dipole orientation. This generalized formulation does not imply a physical rotation of a single current element, but rather provides a compact geometrical description of distributed or anisotropic cortical activation.

Distinguishing between linear and elliptical trajectories of $[eqn]$ thus provides a methodological means to evaluate the adequacy of the source model and to characterize differences in the spatial organization of cortical activity across conditions.

2.2. Geometric Representation

Let $[eqn]$ denote the vector from the source to the channel, and let $[eqn]$ be its length. The scalar potential $[eqn]$ generated by an electric dipole in a homogeneous, infinite, and isotropic medium of conductivity $[eqn]$ is given by

$$[eqn]$$

The dipole vector $[eqn]$ does not represent a position but a directed quantity (the dipole moment).

We assume that each dipole is described by an elliptical trajectory lying in an arbitrary plane. An elliptical trajectory provides a natural geometric description of a dipole that is not confined to a single direction and is therefore adopted as a flexible geometric model of dipolar dynamics. Under this assumption, the dipole evolves along an ellipse contained in a plane $[eqn]$ . Consider an orthonormal basis $[eqn]$ adapted to that plane, where $[eqn]$ . The trajectory of the dipole can be written as

$$[eqn]$$

Importantly, the phase $[eqn]$ is an internal parameter of the dipole trajectory and should not be confused with the angle $[eqn]$ which denotes the instantaneous angle between $[eqn]$ and the source–channel direction $[eqn]$ ; see Figure 1a.

A purely linear trajectory arises as a degenerate case of this parametrization: when one of the semi-axes vanishes and, without loss of generality, we assume $[eqn]$ , the ellipse collapses into a line segment; see Figure 1b. The same phase function still provides a valid representation, although the resulting curve no longer explores a two-dimensional elliptical geometry.

From now on, the superscript d denotes quantities associated with dipole d, and the subscript g denotes quantities associated with channel g. For each dipole–channel pair, let $[eqn]$ denote the corresponding source–channel direction.

The scalar product of the dipole moment with this direction can be written as

$$[eqn]$$

Therefore, in the ideal homogeneous infinite-medium case, the contribution of dipole d to the signal measured at channel g, denoted by $[eqn]$ , is given by $[eqn]$ up to a time-independent constant $[eqn]$ .

Using the identity $[eqn]$ , $[eqn]$ can be rewritten as

$$[eqn]$$

where

$$[eqn]$$

and the amplitude $[eqn]$ and phase offset $[eqn]$ are given by

$$[eqn]$$

Finally, since the recorded voltage at channel g is the superposition of the contributions from D dipoles, the observation model takes the form

$$[eqn]$$

where $[eqn]$ and $[eqn]$ represents additive Gaussian noise.

2.3. FMM Representation

The dipole is assumed to traverse its trajectory without reversals. Since the motion is fully determined by the phase function $[eqn]$ , this assumption requires $[eqn]$ to be monotonic. From a modeling perspective, explicitly parametrizing $[eqn]$ is advantageous, as it yields a small and interpretable set of parameters and promotes parsimony. Such parsimony is essential because the number of active dipoles is rarely known a priori; without structural assumptions of this kind, the EEG inverse problem becomes severely under-determined and non-unique [22,39].

The simplest monotonic choice is a linear phase, $[eqn]$ , corresponding to uniform angular motion. A slightly richer alternative is provided by Möbius-type reparameterizations, which preserve monotonicity while allowing for non-uniform traversal of the ellipse. These functions represent the simplest nonlinear extensions compatible with circular geometry and are defined as

$$[eqn]$$

This phase function is used to describe the temporal trajectory of the dipole vector. The parameter $[eqn]$ sets the temporal shift and thus the latency of the response, $[eqn]$ controls how fast the trajectory is traversed around its maximum, and $[eqn]$ determines the relative timing of the rising and falling phases of the dipole magnitude.

The projection of an elliptical trajectory therefore yields a unimodal but potentially asymmetric waveform, as a direct geometric consequence, whose skewness depends on the geometry of the ellipse and the projection direction. Waveforms based on this parametrization are known as FMM modulations. Previous studies [34,35] have shown that FMM models provide a compact and flexible parametric representation for oscillatory signals.

In multichannel settings such as EEG, where data on multiple electrodes are recorded, the observed signals arise as linear mixtures of the same dipole sources. Each electrode records the projection of the same elliptical trajectory onto its own direction. Accordingly, all channels are assumed to share the same parameters $[eqn]$ governing the instantaneous phase $[eqn]$ of each dipole source, while electrode-specific differences are captured through $[eqn]$ and $[eqn]$ parameters. Multichannel FMM representations have been successfully applied to multivariate time series, including ECG and EEG, enabling accurate modeling of frequency variation and waveform asymmetry [40,41].

2.4. Reduced-Rank Regression

In the classical multivariate regression model

$$[eqn]$$

where $[eqn]$ is an $[eqn]$ matrix of response variables, $[eqn]$ is an $[eqn]$ matrix of predictors, $[eqn]$ is a $[eqn]$ matrix of regression coefficients, and $[eqn]$ denotes the $[eqn]$ matrix of errors, a reduced-rank constraint is imposed by requiring that

$$[eqn]$$

This assumption implies that the responses depend on a limited number of latent linear combinations of the predictors. Equivalently, the coefficient matrix can be factorized as

$$[eqn]$$

where $[eqn]$ is a $[eqn]$ matrix and $[eqn]$ is a $[eqn]$ matrix. The model can thus be written as

$$[eqn]$$

so that the p response variables are explained through r latent predictor combinations $[eqn]$ .

Estimation of $[eqn]$ and $[eqn]$ is based on minimizing the sum of squared residuals

$$[eqn]$$

The least-squares solution is obtained from the singular value decomposition (SVD) of the sample cross-covariance between predictors and responses, which provides the best rank-r approximation to the multivariate regression surface in the least-squares sense [42,43]. The singular vectors associated with the largest r canonical correlations define the subspaces that maximize the explained cross-covariance between $[eqn]$ and $[eqn]$ , establishing a close connection between reduced-rank regression and canonical correlation analysis.

3. Proposed Model

Let $[eqn]$ denote the sampling points, with $[eqn]$ and $[eqn]$ collect the nonlinear parameters of the D sources. Using the reparametrization $[eqn]$ , for each $[eqn]$ , models (5) and (6) can be written as

$$[eqn]$$

where $[eqn]$ is the vector of voltages across channels, $[eqn]$ is the channel-wise baseline, $[eqn]$ is the spatial mixing matrix with

$$[eqn]$$

and $[eqn]$ collects the temporal regressors

$$[eqn]$$

where $[eqn]$ ; and $[eqn]$ is variance-covariance matrix of the residuals.

3.1. Reduced-Rank Model Representing the Fixed-Orientation Hypothesis

The restriction of fixed dipole orientations translates into linear constraints on the structure of the coefficient matrix $[eqn]$ . Going back to the geometric representation given in (3) and (4), a linear trajectory corresponds to $[eqn]$ which implies:

$$[eqn]$$

Therefore, when the dipole direction does not change over time, the phase shift $[eqn]$ collapses to $[eqn]$ . The effect of adding or subtracting $[eqn]$ to the argument of a cosine or a sine is given by the identities

$$[eqn]$$

Thus, a cosine shifted by $[eqn]$ becomes a sine, and a sine shifted by $[eqn]$ becomes a cosine. Therefore, in the linear case the relative combination of the cosine and sine components is identical for every channel: only a single multiplicative factor changes from channel to channel.

To make this constraint explicit, let $[eqn]$ denote the two-column block of the mixing matrix associated with dipole d, i.e., $[eqn]$ . Under the general (elliptical) formulation, $[eqn]$ is complete and typically has rank two. Under the fixed-orientation hypothesis, the cosine and sine spatial patterns must be proportional across channels, which is equivalently expressed as the rank-one restriction

$$[eqn]$$

That is, there exist $[eqn]$ and $[eqn]$ such that $[eqn]$ . This reduced-rank structure provides an explicit statistical encoding of the fixed-orientation assumption at the channel level. As a consequence, the general representation with $[eqn]$ channel-specific coefficients reduces to $[eqn]$ parameters: a single common temporal waveform for dipole d, together with one scalar multiplier per channel.

Hence, the classical fixed-orientation model is obtained as a reduced-rank special case of the general elliptical formulation. Testing these rank conditions provides a principled way to assess whether each cortical source is directionally stable (fixed) or dynamically evolving (elliptical).

3.2. Comparison with PCA, ICA, and Related Decomposition Methods

Classical decomposition methods such as PCA or ICA also represent multichannel signals as linear mixtures of latent components:

$$[eqn]$$

where $[eqn]$ denotes a vector of latent time courses. PCA enforces orthogonality among the trajectories $[eqn]$ , whereas ICA imposes statistical independence. Neither method specifies a particular temporal structure for $[eqn]$ , nor assigns biophysical meaning to the columns of $[eqn]$ .

The proposed FMM formulation (11) preserves the same linear mixing architecture but introduces two key extensions.

First, the latent trajectories are replaced by analytic FMM waveforms whose parameters determine latency, asymmetry, and frequency modulation. Thus, the complete vectors $[eqn]$ of PCA/ICA are replaced by smooth, parametric, and physiologically interpretable temporal dynamics.

Second, the coefficient matrix $[eqn]$ acquires biophysical structure: each dipole contributes two columns corresponding to its cosine and sine projection patterns across channels. Linear relations between these columns encode hypotheses about dipole orientation. A fixed-orientation dipole leads to a rank-one column pair, whereas an elliptical dipole requires a full rank-two pair.

As a result, the FMM model is not a purely algebraic factorization but a reduced-rank nonlinear least-squares framework that integrates structured temporal dynamics with physiological spatial constraints. This yields a representation that is parsimonious, interpretable, and well-suited for EEG analysis and complementary to classical PCA/ICA.

While reduced-rank regression is not a recent technique, it continues to be relevant for analyzing multichannel neural data, including EEG, where it enables low-dimensional representations of inter-signal dependencies [44].

3.3. Parameter Estimation Algorithm

The parameters $[eqn]$ are estimated by maximum likelihood under the Gaussian model (11). Although the likelihood is nonlinear in $[eqn]$ , it is quadratic in $[eqn]$ for fixed $[eqn]$ . Therefore, for any given $[eqn]$ , the optimal coefficient matrix $[eqn]$ is available in closed form via generalized least squares, possibly under rank constraints. This naturally motivates a profile-likelihood strategy in which $[eqn]$ is updated analytically within the optimization over $[eqn]$ .

Alternating Updates

At each iteration, the following steps are performed:

• 1.Update of the nonlinear parameters $[eqn]$ . Given $[eqn]$ , we maximize the profile log-likelihood

$$[eqn]$$

where $[eqn]$ denotes the closed-form GLS (or reduced-rank GLS) estimate of $[eqn]$ for the current value of $[eqn]$ . In practice, this step follows multivariate procedures described in the FMM literature, e.g., [34,35], exploiting the smoothness of the Möbius phase functions and using a profile approach with respect to $[eqn]$ .

• 2.Update of the coefficient matrix $[eqn]$ . Conditional on $[eqn]$ , the model is linear in $[eqn]$ ,

$$[eqn]$$

so we set $[eqn]$ . When orientation constraints are imposed, they are enforced at this stage through the corresponding reduced-rank least-squares update.

• 3.Update of $[eqn]$ . The noise covariance is re-estimated from the current residuals. Under diagonal or structured assumptions this reduces to simple variance estimators; under a full covariance model it corresponds to the empirical covariance of the GLS residuals.

3.3.1. Computational Properties

The alternating scheme exploits the separable structure of the likelihood: the nonlinear parameters $[eqn]$ enter only through the temporal functions, while $[eqn]$ and $[eqn]$ enter linearly. This yields stable iterations and avoids the need for joint optimization in a high-dimensional nonconvex space. In practice, only a small number of outer iterations is required for convergence, as documented in existing multivariate FMM applications [35,40].

3.3.2. Initialization

The noise covariance $[eqn]$ must be initialized. Following standard practice in likelihood- based multivariate models, we use as starting point the empirical covariance matrix of the observed voltages. This choice provides a consistent and data-driven estimate of the overall noise scale while ensuring that the diagonal entries of $[eqn]$ are sufficiently large at initialization. In particular, setting the initial variances to values not smaller than their empirical counterparts prevents numerical degeneracies in the first iterations and guarantees that the profile-likelihood updates of the nonlinear parameters are well conditioned. The initial $[eqn]$ is subsequently refined in step 3 of the algorithm using the residuals of the current FMM fit.

3.3.3. Handling the Linear (Degenerate) Case

The same algorithm accommodates both the general elliptical model and the fixed-orientation (rank-one) case. The difference lies only in step 2: for fixed orientation, the two columns of $[eqn]$ associated with each dipole are constrained to be proportional, whereas the general elliptical case allows a full two-dimensional subspace. The estimation framework thus naturally incorporates hypothesis testing on dipole orientation through rank constraints on $[eqn]$ .

This estimation strategy yields a computationally efficient and statistically robust procedure that exploits the analytic structure of FMM waveforms while respecting the biophysical constraints of dipolar EEG generators.

3.3.4. Model Comparison

Although the fixed- and variable-orientation models are nested, formal likelihood-based tests are not well suited to the present setting. In high-density EEG, both formulations involve a large number of parameters, leading to very large log-likelihood values that are highly sensitive to channel-wise noise. As a result, differences in likelihood may primarily reflect high dimensionality rather than meaningful differences in model adequacy. For similar reasons, information criteria based on the likelihood are also difficult to interpret in this context. We therefore rely primarily on descriptive goodness-of-fit measures, focusing on the coefficient of determination.

For each channel g, the coefficient of determination is defined as

$$[eqn]$$

where $[eqn]$ denotes the model-based fitted value for channel g.

To prevent noise-dominated channels from driving likelihood differences and covariance estimation, we examine the empirical distribution of channel-wise $[eqn]$ values obtained under the complete model and restrict model comparison to channels in its upper tail, defined here by the 95th percentile. This selection is used solely for stability and interpretability and does not constitute a population-level inferential procedure. Throughout the paper, Model A refers to the variable-orientation (complete) formulation in (11), whereas Model B denotes the fixed-orientation formulation under the rank-one restriction (16).

4. Empirical Evidence: Real EEG Data

4.1. Data and Preprocessing

We analyze a publicly available high-density EEG dataset acquired during a visual stimulation task in healthy young adults [45]. The data consist of raw EEG recordings organized according to the EEG-BIDS specification and were collected under ethical approval from the CIUSSS de l’Estrie—CHUS.

In this dataset, visual stimulation produces reproducible and spatially structured EEG responses across trials. Decomposing the EEG into latent components is therefore appropriate for capturing stimulus-related patterns while reducing the dimensionality of the high-density recordings. From a modeling standpoint, we adopt a unified representation in which all spatially coherent contributions to the signal are described within the same parametric framework.

In the context of visually evoked EEG responses, the dominant neural activity is assumed to arise from synchronized cortical populations. Over the time window of the evoked response, the resulting effective dipole is therefore expected to exhibit a relatively stable orientation. We accordingly compare a complete multivariate FMM model with a fixed-orientation variant to assess the impact of imposing this constraint on the component decomposition.

4.1.1. Epoch Averaging and Stimulus Selection

For each participant and stimulus condition, EEG epochs were averaged to obtain event-related potentials (ERPs), thereby improving the signal-to-noise ratio and stabilizing the spatial structure of the responses.

Nevertheless, even after averaging, some participant–stimulus ERPs may remain dominated by noise, particularly when the evoked response is weak, temporally misaligned across trials, or based on a limited number of valid epochs. Trial-to-trial latency variability can lead to partial cancellation of stimulus-locked activity, resulting in attenuated or noisy averaged signals [46]. Such responses are not expected to be informative for hypotheses concerning the spatial organization or orientation of neural generators.

Accordingly, the present analysis focuses on participant–stimulus conditions exhibiting a reliable evoked response, as the goal is to enable a stable and meaningful comparison between competing model formulations rather than to characterize all recorded EEG signals. To implement this selection in an objective manner, the complete FMM model was fitted to each averaged ERP using a fixed number of components ( $[eqn]$ ), and the coefficient of determination ( $[eqn]$ ) was computed independently for each channel. Fixing K ensures comparable model complexity across formulations and isolates the effect of the orientation constraint.

Five participants (IDs 3, 9, 13, 19, and 30) were excluded prior to downstream analyses due to pervasive EEG signal-quality failures. Across stimulus conditions, the majority of channels were either near-flat (almost null) or dominated by white-noise–like activity, with only a very small subset, typically fewer than five channels, showing any structured variation. Under these circumstances, it was not possible to ensure that apparent activity reflected genuine stimulus-evoked voltages rather than residual noise or recording artifacts; consequently, these participants were removed from the study.

4.1.2. Qualitative Illustration of Model Fit

To complement the quantitative comparison, we consider a representative participant–stimulus pair retained after the data-selection procedure. The corresponding ERP exhibits clearly structured activity, making it suitable for illustrating both goodness-of-fit and interpretability.

Figure 2 shows a subset of informative channels for waveform-level inspection. At waveform level, both models reconstruct the main temporal features of the ERP with comparable accuracy. Figure 3 extends the comparison to the full montage. Model A may yield components with mixed polarity across distant regions, whereas Model B tends to produce components with more spatially consistent polarity patterns. This visualization highlights that the model preserves the global spatial organization of the response while representing the multichannel ERP with a small number of structured components.

Finally, Figure 4 shows scalp topographic maps derived from the fitted parameters. For each component k, we map the channel-wise amplitude parameter $[eqn]$ from Equation (5). Under the reduced-rank model, the amplitude can take either of two opposite polarities, according to the sign of $[eqn]$ . The resulting signed amplitudes are interpolated over the scalp using template 10–10 electrode coordinates, since subject-specific sensor locations are not available. The resulting model-based topographies provide a compact spatial summary of each latent component, highlight distinct spatial patterns across components, and offer an additional qualitative check that complements the waveform-level fits in Figure 2 and Figure 3.

4.2. Validation of the Hypothesis of Fixed-Orientation

The fixed-orientation hypothesis is formulated at the level of the underlying neural generators. In the context of ERPs, it is implicitly assumed that the dominant generators contributing to the averaged response are consistently engaged across trials, so that their spatial configuration is preserved under epoch averaging. Under this standard assumption, testing the fixed- versus variable-orientation formulations on averaged ERPs provides a meaningful proxy for assessing the orientation stability of the underlying dipolar sources.

Specifically, we use ERP data to compare two formulations of the proposed model: (A) a variable-orientation (complete) model and (B) fixed-orientation (reduced-rank) model.

Both models are fitted to the same participant–stimulus data using an identical number of components ( $[eqn]$ ), thereby matching model complexity and isolating the effect of the orientation constraint. Table 1 summarizes the 95th percentile of the $[eqn]$ distribution across channels for both models.

Across participants, both model formulations provide a good description of the ERP data on the selected channels, with 95th-percentile $[eqn]$ values typically exceeding 0.80. Differences in goodness of fit between the complete and fixed-orientation models are generally small, as reflected by the near-zero average $[eqn]$ , and vary in sign across individuals. In some participants, the complete model yields slightly higher $[eqn]$ values, whereas in others the fixed-orientation formulation performs comparably or better.

5. Discussion

The empirical results obtained from the visual evoked EEG data indicate that both the complete and the fixed-orientation formulations provide an accurate description of the dominant stimulus-locked activity across participants. The similarity in high-percentile $[eqn]$ values suggests that, for this experimental paradigm, enforcing a fixed dipole orientation provides a parsimonious and interpretable representation of the underlying dynamics. While the complete model captures additional variability in some cases, this added flexibility is not systematically required and may reflect noise in some cases rather than physiologically meaningful signal structure. These findings are consistent with the general observation that increased model complexity does not necessarily translate into improved interpretability or robustness [5,47]. In this context, the proposed framework differs from standard data-driven approaches such as ICA or PCA by explicitly incorporating physiological constraints, allowing orientation-related assumptions to be examined within the model rather than inferred post hoc.

An important aspect of the proposed approach concerns how model complexity is distributed. Instead of increasing the complexity of individual components, expressiveness is controlled through the number of simple, interpretable components. This reflects a common trade-off in EEG analysis, where similar explanatory power may be achieved using either a few complex sources or multiple simpler ones. Favoring the latter supports physiological plausibility while maintaining sufficient flexibility to describe multichannel EEG data.

Uncertainty in EEG recordings represents a key limitation when interpreting the results. EEG signals are affected by noise and artifacts arising from ocular activity, muscle contractions, and environmental factors. Although standard preprocessing steps are applied prior to model estimation, residual artifacts may still influence the fitted parameters. While the proposed framework is not designed as a dedicated artifact-removal method, its emphasis on spatially coherent activity patterns naturally attenuates unstructured noise, which may contribute to the stability of the extracted components.

6. Conclusions

This work introduced a statistically coherent and physiologically motivated framework for modeling multichannel EEG signals using a parametric description of cortical dipolar activity. The proposed FMM-based representation is grounded in biophysical principles and provides interpretable parameters related to latency, temporal asymmetry, and dynamical evolution of stimulus-locked responses.

Using a unified likelihood-based approach, we assessed fixed versus variable dipole-orientation formulations on visual evoked EEG data. Both formulations captured the dominant stimulus-locked activity across participants, and the fixed-orientation specification achieved comparable explanatory performance, supporting its use as a parsimonious and interpretable choice for this paradigm.

More broadly, the framework emphasizes interpretable, physiology-constrained parametrizations as an alternative to purely data-driven dimensionality reduction, and it enables coherent parameter-based spatial summaries of multichannel ERP structure.

Finally, the present study is restricted to a single paradigm and dataset and does not include a systematic quantitative comparison with alternative methods. Future work will extend the analysis to additional experimental settings and incorporate more systematic benchmarking and uncertainty quantification. Another natural continuation, motivated by the fact that the reduced-rank formulation captures the main physiologically plausible features in this setting, is to use it as an explicit ingredient in solving the EEG inverse problem. In particular, the reduced-rank multichannel structure provides a principled bridge to dipole fitting, where the estimated component-wise spatial patterns and temporal parameterizations can inform or constrain the estimation of dipole locations and orientations.