- Research
- Open Access
- Published:

# A comparative study of a theoretical neural net model with MEG data from epileptic patients and normal individuals

*Theoretical Biology and Medical Modelling***volume 2**, Article number: 37 (2005)

## Abstract

### Objective

The aim of this study was to compare a theoretical neural net model with MEG data from epileptic patients and normal individuals.

### Methods

Our experimental study population included 10 epilepsy sufferers and 10 healthy subjects. The recordings were obtained with a one-channel biomagnetometer SQUID in a magnetically shielded room.

### Results

Using the method of x^{2}-fitting it was found that the MEG amplitudes in epileptic patients and normal subjects had Poisson and Gauss distributions respectively. The Poisson connectivity derived from the theoretical neural model represents the state of epilepsy, whereas the Gauss connectivity represents normal behavior. The MEG data obtained from epileptic areas had higher amplitudes than the MEG from normal regions and were comparable with the theoretical magnetic fields from Poisson and Gauss distributions. Furthermore, the magnetic field derived from the theoretical model had amplitudes in the same order as the recorded MEG from the 20 participants.

### Conclusion

The approximation of the theoretical neural net model with real MEG data provides information about the structure of the brain function in epileptic and normal states encouraging further studies to be conducted.

## Introduction

Epilepsy is a disorder involving recurrent unprovoked seizures: episodes of abnormally synchronized and high-frequency firing of neurons in the brain that result in abnormal behaviors or experiences. This is a fairly common disorder, affecting close to 1% of the population. The lifetime risk of having a seizure is even higher, with estimates ranging from 10 to 15% of the population. Epilepsy can be caused by genetic, structural, metabolic or other abnormalities. Epileptic disorders can be generalized, partial (focal) or undetermined. A primary generalized seizure starts as a disturbance in both hemispheres synchronously, without evidence of a localized onset. Partial forms of epilepsy start in a focal area of the brain and may remain localized without alteration of consciousness.

MEG is a noninvasive imaging technique, applicable to the human brain with temporal resolution approximately ~1 ms [1]. Several authors during the last decade have demonstrated the importance of MEG in the investigation of normal and pathological brain conditions [2, 3]. The major advantage of MEG over electroencephalography (EEG) is that MEG has higher localization accuracy. This is because the different structures of the head (brain, liquor cerebrospinalis, skull and scalp) influence the magnetic fields less than the volume current flow that causes the EEG. Also, MEG is reference free, so that the localization of sources with a given precision is easier for MEG than it is for EEG [4].

The goal of this study is to compare the theoretical model that follows Poisson or Gauss distributed connectivity [5–12] with experimental MEG data from epileptic patients and healthy volunteers.

## Methods

### Description of the model

Neural nets are assumed to be constructed of discrete sets of randomly interconnected neurons of similar structure and function. The neural connections are set up by means of chemical markers carried by the individual cells. Thus, the neural population of the net is treated as a set of subpopulations of neurons, each of them characterized by a specific chemical marker. We attribute the appropriate Poisson or Gauss distribution law to each subsystem to describe connectivity.

The elementary unit, the neuron, is bistable. It can be either in the resting or in an active (firing) state. The transition from the resting to the firing state occurs when the sum of postsynaptic potentials (PSPs) arriving at the cell exceeds the firing threshold θ of the neuron. PSPs may be excitatory (EPSPs) or inhibitory (IPSPs), shifting the membrane potential closer to or further away from θ, respectively. Each neuron may carry an electrical potential of a few millivolts, which it passes on to the neurons to which it is connected.

In this model, a net with N markers is assumed to be constructed of A formal neurons. A fraction h (0<h<l) of these are inhibitory with all the axon branches generating IPSPs, while the rest are excitatory with all their axon branches generating EPSPs. Each neuron receives, on average, μ^{+} EPSPs and μ^{-} IPSPs. The size of the PSP produced by an excitatory (inhibitory) unit is K^{+} (K^{-}). The neurons are also characterized by the absolute refractory period and the synaptic delay τ. If a neuron fires at time t, it produces the appropriate PSP after a fixed time interval τ, the synaptic delay. PSPs arriving at a neuron are summed instantly, and if this sum is greater or equal to θ, then the neuron will fire immediately, otherwise it will be idle. PSPs (if below θ) will persist with or without decrement for a period called the summation time, which is assumed to be less than the synaptic delay. Firing is momentary and causes the neuron to be insensitive to further stimulation for a time interval called the (absolute) refractory period [5–12].

The mathematical formalism of this study is based on the equations for the expectation values of the activity derived in previous papers [5–12]. A brief mathematical analysis for each case is given below.

#### a) Expectation value of neural activity in noiseless and noisy neural nets with Poisson distributed connectivities

Following the assumptions of previous papers it was shown that the expectation value of the neural activity <α_{n+1}> at t = (n+1) τ, i.e. the average value of α_{n+1} generated by a collection of netlets with identical statistical parameters (μ^{+}, μ^{-}, h, K^{+}, K^{-}, A, θ) and the same α_{n} at t = nτ, is given by:

<α_{n+1}> = (1-α_{n}) P (α_{n}, θ) (1)

where P(α_{n}, θ) is the probability that a neuron receives post synaptic potentials (PSPs) exceeding its threshold at time t = (n+1)τ. Thus:

Here P_{l} and Q_{m} are the probabilities that a neuron will receive l and m EPSPs and IPSPs respectively, and are given by (3):

P_{I} = exp (-α_{n} (1-h) μ^{+}) (α_{n} (1-h) μ^{+})^{l}/l!

Q_{m} = exp (-α_{n} h μ^{-}) (α_{n} hμ^{-})^{m}/m! (3)

In addition, the upper limits in the double sum m_{max} and l_{max} are given by (4):

l_{max} = A α_{n} (1-h) μ^{+}

m_{max} = Aα_{n}hμ^{-} (4)

Taking into account equations (2) and (3), equation (1) takes the form:

Similarly for Poisson nets with noise: if P_{l} and Q_{m} are the probabilities that a given neuron receives I EPSPs and m IPSPs at time t = (n+1)τ, they are given by equation (3). But if T_{δ} (θ) is the probability that the instantaneous threshold value is θ or less than θ, this is given by (6):

Therefore the firing probability per neuron is then given by (7):

where l_{max} and m_{max} are given by equation (4).

Finally, the expectation value of the activity is given by (8):

<α_{n+1}> = (1-α_{n}) P (α_{n}, θ) (8)

#### b) Expectation value of neural activity in neural nets with Poisson distributed connectivities with chemical markers and noise

Similarly, the expectation value of the activity <α_{n+1}> for an isolated neural net with two chemical markers a and b is given by (9):

where P_{I}, Q_{i}, P'_{l}, Q'_{i'}, are the probabilities that a given neuron will receive l EPSPs, i IPSPs or l'-EPSPs, i'-IPSPs, at time t = (n+1)τ in the subsystems a or b respectively. These probabilities are given by (10):

P_{l} = exp (-α_{n} μ_{a}^{+} (1-h_{a}) m_{a}) (-α_{n} μ_{a}^{+} (1-h_{a}) m_{a})^{l}/l!

Q_{i} = exp (-α_{n} μ_{a}^{-} h_{a} m_{a}) (-α_{n} μ_{a}^{-} h_{a} m_{a})^{i}/i!

P'_{l'} = exp (-α_{n} μ_{b}^{+} (1-h_{b}) (1-m_{a})) (-α_{n} μ_{b}^{+} (1-h_{b}) (1-m_{a}))^{l'}/l'!

Q'_{i'} = exp (-α_{n} μ_{b}^{-} h_{b} (1-m_{a})) (-α_{n} μ_{b}^{-} h_{b} (1-m_{a}))^{i'}/i'! (10)

The upper limits in the sums in equation (9) are given by (11):

l_{max} = A α_{n} μ_{a}^{+} (1-h_{a}) m_{a}

l_{max'} = A α_{n} μ_{b}^{+} (1-h_{b}) (1-m_{a})

i_{max} = A α_{n} μ_{a}^{-} h_{a} m_{a}

i_{max'} = A α_{n} μ_{b}^{-} h_{b} (1-m_{a}) (11)

Finally,

(θ_{a}) and (θ_{b}) are defined as the probabilities that the instantaneous neural thresholds are equal to or less than θ_{a} and θ_{b} in subsystems a and b respectively and are given by (12):

#### b) Expectation value of neural activity in neural nets with Gaussian connectivities in the absence of chemical markers

Let the total PSP of a neuron at t = (n+1)τ be given by:

e_{n+1} = lK^{+} + mK^{-} (13)

where l and m are the numbers of EPSPs and IPSPs respectively. If both l and m are large, their distributions may be approximated by Gaussian distributions about their respective average values

and . The distribution of e_{n+1} is therefore also normal, since the probabilities for l and m are mutually independent, and its variance is the sum of the variances of l and m. Therefore the average PSP will be given by (14):

where K = [μ^{+} (1-h) K^{+} + μ^{-}h K^{-}] (14)

The variance of e_{n+1}, call it , is then given by (15):

= α_{n} [μ^{+} (1-h) (K^{+})^{2} + μ^{-}h (K^{-})^{2}] (15)

The probability that the PSP exceeds a threshold

now becomes:

Equation (16) in conjunction with equation (1) gives values for <α_{n+1}> at t = (n+1)τ.

Let T(θ') be the probability that the instantaneous threshold of a neuron is θ' or less than θ'. This is given by (17):

Here δ is the standard deviation of the Gaussian distribution of the noise. Finally, the probability that a neuron will receive PSPs that will exceed the threshold at time t = (n+1)τ is given by (18):

Since l and m are very large numbers, the double sum can be approximated by

and therefore:

Then the expectation value of <α_{n+1}> of the activity at time t = (n+1)τ will be:

<α_{n+1}> = (1-α_{n}) P(α_{n}, δ_{n+1}, δ) (20)

#### c) Expectation value of neural activity in noisy neural nets with chemical markers and Gaussian distributed connectivities

In a neural netlet of A neurons with two chemical markers a and b, let the fractional numbers corresponding to each chemical marker be m_{a} and m_{b}, and the fractions of inhibitory neurons for each chemical marker be h_{a} and h_{b}, respectively. Also, let α_{n}A be the active neurons in the netlet at t = nτ. Then at t = (n+1)τ the numbers of EPSPs and IPSPs that will appear in the subsystems with a and b markers will be:

l_{a} = A α_{n} μ_{a}^{+} (1-h_{a}) m_{a}

i_{a} = A α_{n} μ_{a}^{-} h_{a} m_{a}

l_{b} = A α_{n} μ_{b}^{+} (1-h_{b}) m_{b}

i_{b} = A α_{n} μ_{b}^{-} h_{b} m_{b} (21)

On the average, the numbers of EPSPs and IPSPs that appear per neuron in subnets with a and b markers will be:

= α_{n} μ_{a}^{+} (1-h_{a}) m_{a}

= α_{n} μ_{a}^{-} h_{a} m_{a}

= α_{n} μ_{b}^{+} (1-h_{b}) m_{b}

= α_{n} μ_{b}^{-} h_{b} m_{b} (22)

As stated in our previous papers [5–12] the total PSP input to a neuron with a and b markers at t = (n+1)τ will be given by (23):

e_{a,n+1} = l_{a}K^{+} + i_{a}K^{-}

e_{b,n+1} = l_{b}K^{+} + i_{b}K^{-} (23)

(Here it is assumed that

)

If the quantities l_{a}, l_{b}, i_{a} and i_{b} are sufficiently large, their distributions may be approximated by Gaussian distributions about their average values, given by (22). Then the average PSPs for the two markers a and b will be given by (24):

and their variances will be given by (25):

Therefore the probability that a neuron with marker a or b will receive a certain number of EPSPs or IPSPs that will shift the membrane potential closer to or further away from the instantaneous threshold will be given by (26):

where:

Thus, the probabilities

and that the instantaneous threshold of a neuron in subsystems a and b is equal to or less than or will be given by (28):

Consequently, as stated in our previous paper [8], the firing probabilities P(α_{n}, δ_{n+1}, δ_{a}) and P'(α_{n}, δ_{n+1}, δ_{b}) that a neuron in subpopulations a and b, respectively, will receive PSPs exceeding threshold at time t = (n+1)τ will be given by (29):

Since the quantities l_{a}, i_{a}, l_{b} and i_{b} are sufficiently large, the double sum in equations (29) will be substituted by the probabilities of the average values of l_{a}, i_{a} and l_{b}, i_{b} for each marker a and b and will be given by (30):

Then according to our previous papers [5–12], the expectation value of activity in this netlet with two markers a and b at time t = (n+1)τ will be given by (31):

The general case for an isolated noisy net with N markers m_{1}, m_{2},..., m_{N}, where m_{i} is the fraction of neurons with the i^{th} marker, is described by an equation analogous to the equation for two markers (31). This general equation for such a netlet at time t = (n+1)τ is:

### Theoretical analysis

#### The electromagnetic fields generated in neural networks with Poisson or Gauss connectivities

Let us consider an isolated neural network with structural parameters A, μ^{+}, μ^{-} and h, and initial activity α_{n} at time t = nτ. The potential generated in this network due to this initial activity will be equal to the summation of all the PSPs [7] and will be given by (33):

V_{n} = α_{n} (A μ^{+} (1-h) - A μ^{-}h) (33)

Similarly, the potential generated by the neural activity α_{n+1} at the next time interval t = (n+1)τ will be given by (34):

V_{n+1} = α_{n+1} (A μ^{+} (1-h) - A μ^{-}h) (34)

By combining equations (33) and (34) and assuming spherical brain symmetry, the potential difference ΔV can be obtained. As is known from classical physics, this generates a magnetic field B_{n} given by (35):

Choosing Δt = 1 ms, the above equation takes the following form:

where μ_{o} and ε_{o} are the magnetic permeability and dielectric constant of the medium respectively.

When the neural network is characterized by two chemical markers a and b, the potentials created at the synapses of the neurons with the a and b markers will be given by (37):

V_{na} = α_{n} (A μ_{a}^{+} (1-h_{a}) m_{a} - A μ_{a}^{-} h_{a} m_{a})

V_{nb} = α_{n} (A μ_{b}^{+} (1-h_{b}) m_{b} - A μ_{b}^{-} h_{b} m_{b}) (37)

On the other hand, the total voltages created at the synapses of the neurons at times t = nτ and t = (n+1)τ will be given by (38):

V_{n} = V_{na} + V_{nb} = α_{n} A [(μ_{a}^{+} (1-h_{a}) m_{a} + μ_{b}^{+} (1-h_{b}) m_{b}) - (μ_{a}^{-} h_{a} m_{a} + μ_{b}^{-} h_{b} m_{b})]

V_{n+1} = α_{n+1} A [(μ_{a}^{+} (1-h_{a}) m_{a} + μ_{b}^{+} (1-h_{b}) m_{b}) - (μ_{a}^{-} h_{a} m_{a} + μ_{b}^{-} h_{b} m_{b})] (38)

Therefore the potential difference between these two time intervals, taking into account equations (38), is given by (39):

ΔV = V_{n+1} - V_{n} = (α_{n+1} - α_{n}) A [(μ_{a}^{+} (1-h_{a}) m_{a} + μ_{b}^{+} (1-h_{b}) m_{b}) - (μ_{a}^{-} h_{a} m_{a} + μ_{b}^{-} h_{b} m_{b})] (39)

Thus, as stated previously, this potential difference will create a magnetic field B_{n}, which is given by (40):

where the neural activity α_{n+1} refers to a Poisson or Gauss distribution of connectivities as given in the previous section.

In the general case, where the neural net has N chemical markers, equation (40) takes the form:

##### Experimental procedure

We compared the theoretical results with the experimental findings obtained using MEG measurements from 10 epileptic patients and 10 healthy volunteers. Informed consent for the methodology and the aim of the study was obtained from all participants prior to the procedure.

Biomagnetic measurements were performed using a second order gradiometer SQUID (Model 601, Biomagnetic Technologies Inc.), which was located in a magnetically shielded room with low magnetic noise. The MEG recordings were performed after positioning the SQUID sensor 3 mm above the scalp of each patient using a reference system. This system is based on the International 10–20 Electrode Placement System [13] and uses any one of the standard EEG recording positions as its origin; in this study we used the P3, P4, T3, T4, F3, and F4 recording positions [14–16]. Around the origin (T3 or T4 for temporal lobes) a rectangular 32-point matrix was used (4 rows × 8 columns, equidistantly spaced in a 4.5 cm × 10.5 cm rectangle) for positioning of the SQUID [14–16]. The MEG was recorded from each temporal lobe at each of the 32 matrix points of the scalp for 32 s and was band-pass filtered with cut-off frequencies of 0.1 and 60 Hz. The MEG recordings were digitized using a 12 bit precision analog-to-digital converter with a sampling frequency of 256 Hz, and were stored in a PC peripheral memory for off-line Fourier statistical analysis. The method, by its nature (i.e. temporal and spatial averaging), eliminates short-term abnormal artifacts in any cortical area, while it retains long-lasting localized activation phenomena. We used the x^{2} – fitting method to analyze the MEG data [17].

This method was based on the following equation (42):

where:

Q_{i}: is the number of elements in the i^{th} interval of the normalized MEG histogram

T_{i}: is the number of elements in the i^{th} interval of the normal distribution with the same mean value and standard deviation as the normalized MEG histogram

k: is the number of intervals

n = k-1: the degrees of freedom of the system

In our case n = 7 and the critical value for distinguishing the Poisson from the Gauss distribution was 14.1 (x_{cr}^{2} = 14.1). If the estimated value of the x^{2} was greater than 14.1, the distribution was Poisson; otherwise it was Gauss.

## Results

Using the x^{2}-fitting method it was found that the MEG recordings from epileptic patients had Poisson distributions whereas those from normal subjects had Gauss distributions. The Poisson connectivity derived from the theoretical model represents the state of epilepsy, whereas the Gauss connectivity represents normal behavior. The magnetic field derived from the theoretical model was approximately in the same order as the recorded MEG in both conditions. Furthermore, the MEG data obtained from epileptic areas had higher amplitudes than those from normal regions and were comparable with the theoretical magnetic fields from Poisson and Gauss distributions.

Figure 1 shows the MEG recorded from an epileptic patient; figure 2 illustrates the MEG recorded from a healthy volunteer.

Figures 3 and 4 show the magnetic fields derived from the theoretical model with Poisson and Gauss distributions respectively.

## Discussion

Over the past three decades, neural nets have been intensively studied from several points of view. An area of considerable importance is that of biological nets, i.e. models of nets designed to imitate the structures and functions of human and other living brains and thus enhance our understanding of learning, memory, understanding etc. Widely used models include the pioneering work of McCulloch and Pitts [18], which treats assemblies of neurons as logical decision elements, the mathematical formalism of Caianiello [19] using the "neuronic equation", and probabilistic neural structures [5, 6] that monitor the net activity, i.e. the fraction of neurons that become active per unit time. All these models have had a measure of success in improving our understanding of functions such as those mentioned above.

The effect of structure on function and dynamic behavior in neural nets has been also a subject of considerable interest in recent years. In the so-called probabilistic nets we have an assembly comprising a large number of neurons, randomly positioned in space, that have only partial connectivity; i.e. each neuron is connected to only a very small fraction of the total number of neurons in the system, randomly chosen. The principal idea is that this connectivity is given by the binomial distribution. In earlier work, probabilistic neural nets were investigated using Poisson or Gauss distributions of interneuronal connectivity; the main conclusion was that when a neuron was connected to a relatively small number of units, a Poisson distribution law was appropriate but if it was connected to a large number of units then a Gaussian law was a fairly good approximation [10–12]. Thus, Poisson neuronal nets may be viewed as approximately Gaussian whenever the number of synaptic connections is relative large.

In this study we measured the MEG of epileptic patients and normal subjects in order to compare the theoretical neural net model [10–12] with real data. Analyzing the MEG data by x^{2}-fitting revealed that the MEG recordings from epileptic areas had Poisson distributions [17]. This finding is consistent with the correspondence between Poisson distributions and low numbers of internal neural connections, and with the synchronization of neural activity during an epileptic seizure [20, 21]. Moreover, the MEG recordings from epileptic areas showed higher amplitudes than those from normal regions, comparable with the results from the theoretical neural model with Poisson and Gauss distributions respectively (Figs. 1, 2, 3, 4).

If a nerve cell is characterized by a given firing threshold which, when exceeded, results in spike discharge, two anatomical situations can be contrasted: one in which only a few synaptic contacts reach the cell in question, and a second in which the cell receives a large number of synaptic inputs. Suppose, in either case, that firing is dependent on the simultaneous excitation of a certain percentage of the total synaptic input (assuming that the ratio of excitatory and inhibitory synapses is the same in both situations so that the inhibitory inputs may be disregarded for the moment). Then it is clear that firing in neurons with a large number of synaptic inputs would require the synchronized activation of a substantial number of synapses; whereas in neurons with few synapses, firing may ensue even from a single excitatory synapse. Thus, a system in which neurons receive small numbers of synaptic connections is likely to exhibit a less "controlled" pattern of activity – and also "spontaneous" discharges [22]. The inverse problem in MEG measurements is the search for unknown sources by analysis of the measured field data. To handle this task one must first study the forward problem, i.e. how the magnetic field and the electrical potential arise from a known source. For practical purposes one also has to choose appropriate models for the source and the biological object as a conductor. Sarvas [23] described basic mathematical and physical concepts relevant to the forward and inverse problems and discussed some new approaches. Especially, he described the forward problem for both homogenous and inhomogenous media. He referred to the Geselowitz's formulae and presented a surface integral equation to handle a piecewise homogenous conductor and a horizontally layered medium. Furthermore, he discussed the non-uniqueness of the solution of the magnetic inverse problem and studied the difficulty caused by the contribution of the electric potential to the magnetic field outside the conductor.

The Poisson distribution corresponds to epileptic areas and the Gauss distribution to normal regions. The approximation of the theoretical neural net model to real MEG data provides a mathematical approach to the structure of brain function and indicates the need for further studies.

## Appendix

The subscript i is a marker label and indicates the properties of a subpopulation in the netlet characterized by the i^{th} marker.

### Structural parameters of the neural net

τ Synaptic delay

A Total number of neurons in the netlet

h_{i} Fraction of inhibitory neurons

The average number of neurons receiving excitatory postsynaptic potentials (EPSPs) from one excitatory neuron

The average number of neurons receiving inhibitory postsynaptic potentials (IPSPs) from one inhibitory neuron

The size of PSP produced by an excitatory neuron of the netlet

The size of PSP produced by an inhibitory neuron of the netlet

m_{i} Fractions of neurons carrying the i^{th} marker in the netlet

θ_{i} Firing thresholds of neurons

### Statistical parameters

δ_{i} Standard deviation of the Gaussian distribution of the neural firing thresholds in the i^{th} subpopulation

### Dynamical parameters

n An integer giving the number of elapsed synaptic delays

α_{n} The activity, i.e. the fractional number of active neurons in the netlet at time t = nτ

## References

- 1.
Mitra PP, Pesaran B: Analysis of dynamic brain imaging data. Biophys J. 1999, 691-708.

- 2.
Timmermann L, Gross J, Dirks M, Volkmann J, Freund HJ, Schnitzler A: The cerebral oscillatory network of parkinsonian resting tremor. Brain. 2003, 126: 199-212. 10.1093/brain/awg022.

- 3.
Volkmann J, Joliot M, Mogilner A, Ioannides AA, Lado F, Fazzini E, Ribary U, Llinas R: Central motor loop oscillations in parkinsonian resting tremor revealed by magnetoencephalography. Neurol. 1996, 46: 1359-1370.

- 4.
Kristeva-Feige R, Rossi S, Feige B, Mergner T, Lucking CH, Rossini PM: The bereitschaftspotential paradigm in investigating voluntary movement organization in humans using magnetoencephalography (MEG). Brain Res Protocol. 1997, 1: 13-22. 10.1016/S1385-299X(97)80327-3.

- 5.
Anninos PA, Beek B, Csermely TJ, Harth EM, Pertile G: Dynamics of neural structures. J Theor Biol. 1970, 26: 121-148.

- 6.
Anninos PA, Kokkinidis M: A neural net model for multiple memory domains. J Theor Biol. 1984, 109: 95-110.

- 7.
Kokkinidis M, Anninos PA: Noisy neural nets exhibiting epileptic features. J Theor Biol. 1985, 113: 559-588.

- 8.
Kotini A, Anninos PA: Dynamics of noisy neural nets with chemical markers and Gauss-distributed connectivity. Connect Sci. 1997, 9: 381-403. 10.1080/095400997116603.

- 9.
Fournou E, Argyrakis P, Kargas B, Anninos PA: Neural nets with chemical markers and firing threshold fluctuations. Trends Biol Cybernet. 1991, 2: 123-131.

- 10.
Fournou E, Argyrakis P, Anninos PA: Neural nets with markers and Gauss-distributed connectivity. Connect Sci. 1993, 5: 77-94.

- 11.
Fournou E, Argyrakis P, Kargas B, Anninos PA: A Gaussian approach to neural nets with multiple memory domains. Connect Sci. 1995, 7: 331-339.

- 12.
Fournou E, Argyrakis P, Kargas B, Anninos PA: Hybrid neural nets with Poisson and Gaussian connectivity. J Stat Phys. 1997, 89: 847-867.

- 13.
Jasper HH: The ten-twenty electrode system of the International Federation. Electroenceph Clin Neurophysiol. 1958, 10: 367-380.

- 14.
Anninos PA, Anogianakis G, Lehnertz K, Pantev C, Hoke M: Biomagnetic measurements using SQUID. Int J Neurosci. 1987, 37: 149-168.

- 15.
Anninos PA, Tsagas N, Jacobson JI, Kotini A: The biological effects of magnetic stimulation in epileptic patients. Panminerva Med. 1999, 41: 207-215.

- 16.
Anninos P, Adamopoulos A, Kotini A, Tsagas N: Nonlinear Analysis of Brain Activity in Magnetic Influenced Parkinson Patients. Brain Topogr. 2000, 13: 135-144. 10.1023/A:1026611219551.

- 17.
Spiegel MR, Schiller J, Srinivasan RA: Schaum's outline of Probability and Statistics. 2000, Mc Graw – Hill Companies Inc, USA, 2

- 18.
McCulloch WS, Pitts W: A logical calculus of the ideas immanent in nervous activity. Bull Math Biol. 1943, 5: 115-133.

- 19.
Caianiello ER: Outline of a theory of thought-processes and thinking machines. J Theor Biol. 1961, 2: 204-235. 10.1016/0022-5193(61)90046-7.

- 20.
Anninos PA, Zenone S, Elul R: Artificial neural nets: dependence of the EEG amplitude's probability distribution on statistical parameters. J Theor Biol. 1983, 103: 339-348. 10.1016/0022-5193(83)90290-4.

- 21.
Leake B, Anninos PA: Effect of connectivity on the activity of neural nets models. J Theor Biol. 1976, 58: 337-363.

- 22.
Anninos PA, Elul R: Effect of structure on function in model nerve nets. Biophys J. 1974, 14: 8-19.

- 23.
Sarvas J: Basic mathematical and electromagnetic concepts of the biomagnetic inverse problem. Phys Med Biol. 1987, 32: 11-22. 10.1088/0031-9155/32/1/004.

## Author information

## Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

## Rights and permissions

## About this article

#### Received

#### Accepted

#### Published

#### DOI

### Keywords

- Poisson distribution
- Gauss distribution
- MEG