BioMed Research International

Volume 2014 (2014), Article ID 503920, 8 pages

http://dx.doi.org/10.1155/2014/503920

Research Article

Chaotic Analysis of the Electroretinographic Signal for Diagnosis

Surya S. Nair and K. Paul Joseph

Electrical Engineering Department, National Institute of Technology, Calicut, Kerala 673601, India

Received 24 February 2014; Accepted 23 May 2014; Published 15 June 2014

Academic Editor: Kevin Gregory-Evans

Copyright © 2014 Surya S. Nair and K. Paul Joseph. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Electroretinogram (ERG) is a time-varying potential which arises from different layers of retina. To be specific, all the physiological signals may contain some useful information which is not visible to our naked eye. However this subtle information is difficult to monitor directly. Therefore the ERG signal features which are extracted and analyzed using computers are highly useful for diagnosis. This work discusses the chaotic aspect of the ERG signal for the controls, congenital stationary night blindness (CSNB), and cone-rod dystrophy (CRD) classes. In this work, nonlinear parameters like Hurst exponent (HE), the largest Lyapunov exponent (LLE), Higuchi’s fractal dimension (HFD), and approximate entropy (ApEn) are analyzed for the three different classes. It is found that the measures like HE dimension and ApEn are higher for controls as compared to the other two classes. But LLE shows no distinguishable variation for the three cases. We have also analyzed the recurrence plots and phase-space plots which shows a drastic variation among the three groups. The results obtained show that the ERG signal is highly complex for the control groups and less complex for the abnormal classes with value less than 0.05.

1. Introduction

All physiological signals exhibit complex behavior which reflects the nonlinear dynamic properties of a biological system. Considering this, the use of nonlinear tools to exhibit the chaotic behavior may be a better approach to explore the nature of the electroretinographic signal. The randomness of the ERG signal does not allow any form of time-series prediction. The study of nonlinear dynamics can contribute to the understanding of the ERG signal and the underlying retinal processes [1, 2].

1.1. Electroretinogram (ERG)

Electroretinogram is the time-varying potential which arises from different retinal layers and is elicited by a brief flash of light stimuli. Contact lens type electrode which carries a silver chloride wire is used to record ERG clinically. The electrode is placed on the cornea and is in the shape of a cup filled with saline. The reference electrode is placed either on the earlobe, temple, or forehead. The amplitude of the ERG waveform is in the range of tenths of millivolt which depends upon the stimulating and physiological conditions [3].

Figure 1 shows the cross-sectional view of human retina and the origin of signals from different layers of the retina. An ERG signal comprises of early receptor potential (ERP) and , -wave ( and ), oscillatory potentials, -wave (cone -wave and rod -wave), and -wave. The initial changes in the photo pigment molecules of the photoreceptor cells (cone cells and rod cells) due to the brief flash of light stimulus will give rise to a positive deflection followed by deflection which forms the ERP. It is then followed by a late receptor potential (LRP) after a 2 ms delay which forms the main portion of the -wave and is a corneo-negative wave. The -wave comprises of two negative dips and which shows the contribution of cone cells and rod cells, respectively. It lasts for about 30 ms. These cone cells and rod cells can be separated by applying appropriate stimuli. A dim-blue light with the dark background extracts a rod ERG and a bright red light with the light adapted background will give a cone ERG. The Muller cells in the inner retina which contribute to the -wave are a corneo-positive wave. Muller responses can be obtained either from the cone cells or from the rod cells separately. The oscillatory potentials which occur in the rising edge of the -wave are small amplitude wavelets that reflect the activity of amacrine cells of the inner retinal layers. The -wave which is generated by the retinal pigment epithelium (RPE) as a result of interaction with the rod cells is a slower positive wave [3].

Figure 1: The cells of retina and the standard ERG waveform (courtesy: Jaakko Malmimo and Robert Plonsey, “Bioelectro-magnetism-Principles and Applications of Bioelectric and Biomagnetic fields”).

1.2.  Recording of the ERG Signal

International Standards for Clinical Electrophysiology of Vision (ISCEV) sets the recording of the ERG signal in five different steps.

1.2.1. Rod Response

The patient is dark adapted for at least 20 minutes. The standard stimulus of dim-white flash (2 seconds between the flash) of 2.5 log units or blue stimulus is given. It is the first signal measured during the ERG recording [4].

1.2.2. Maximal Response

The combined response of both cone cells and rod cells gives maximal response and is produced by the white standard flash (10 seconds between the flashes) [4].

1.2.3. Oscillatory Potentials

Oscillatory potentials or OP are obtained either from light adapted eye (1.5 seconds between flashes) or from dark adapted eye (15 seconds between flashes) using white standard flash. The frequency of interest is set by the band pass filter and the lowest frequency range is 75–100 Hz and 300 Hz and above at the higher end [4].

1.2.4. Cone Response

Before recording the cone response, the patient should be light adapted for at least 10 minutes. The flash used is white or bright red in colour and 0.5 seconds between the flashes is the interval given [4].

1.2.5. 30 Hz Flicker Response

Here the flicker type stimulus is used and flashes are given at the rate of 30 stimuli per second. Rods follow the flickering light up to 12–17 Hz and cones will follow up to 60–70 Hz [4].

1.3. Literature Study

The first remarkable work done in the field of ERG signal analysis by Bornschein et al. was the study of electroretinography in normal, colour-blind, and night-blind subjects under various states of adaptation with varying stimuli [5]. A series of work followed by Barraco et al. shows the three-frequency range of occurrence between 20 and 200 Hz [1] and also analyzed the time-frequency characteristics of the -wave in congenital stationary night blindness (CSNB) patients [2, 6]. Study on the basis of principal component analysis and wavelet analysis was used to visualise the time domain features and wavelet features (Rogala and Brykalski) [7]. Another work reported in the area of ophthalmology is multifocal ERG analysis using wavelet transform by Miguel et al. for the diagnosis of glaucoma [8, 9]. Nair and Paul Joseph have analysed the ERG signal using wavelets and entropy analysis [10]. Study of Crevier and Meister showed that period-doubling occurs in the nonlinear dynamical system and observation of chaotic behavior in the nervous system is the period-doubling route to chaos in flicker vision of the ERG [11]. Molaie et al. showed that the parameters like flash frequency and contrast have a greater impact on the recorded ERG signals which cause bifurcations resulting in a period-doubling and the work defines neural network to be a powerful tool for modeling highly chaotic behavior in the nervous system [12].

To the best of our knowledge, no work is done on the nonlinear aspect of the ERG signal. Here, in this work, we are analysing the ERG signal of controls, congenital stationary night blindness (CSNB), and cone-rod dystrophy (CRD) with the nonlinear chaotic perspective of analysis.

1.4. Organisation of Work

The organisation of work includes the following subsections. In “Data Acquisition” Section we present the details of the acquisition of ERG recordings and various pathological diseases affecting the eye. “Analysis Methods” provides a brief overview of various characteristic measures like Hurst exponent, Lyapunov exponent, approximate entropy, fractal dimension, and recurrence plots for various cases. Following that the results are included. Final conclusion and discussion of the study are also reported.

2. Data Acquisition

Acquisition of the recordings of the ERG signal is performed using TOMEY EP 1000 version 3.0.4 from 15 control subjects, 20 subjects with congenital stationary night blindness (CSNB) type I, 15 subjects with CSNB type II, 35 subjects with cone-rod dystrophy also called as retinitis pigmentosa, among which 15 subjects are of typical RP (retinitis pigmentosa), 10 subjects of early onset RP, and 10 subjects of late stage RP. From the five steps of the recording of ERG signals, we are analyzing only the maximum response and 30 Hz flicker response from the above subjects. Patient data were collected from Little Flower Hospital and Research Centre, India, with the proper consent from the clinicians. Brief descriptions of the pathologies analyzed are given below.

2.1. Controls

The normal amplitude range of rod response is 140–250 μV with the implicit time of 80–90 ms, maximum response ranges from 250 to 500 μV amplitude, and implicit time is 45 ms. Cone response amplitude is 100–180 μV with the implicit time of 32 ms and flicker response amplitude is 50 μV approximately.

2.2. Congenital Stationary Night Blindness (CSNB)

It is an X-linked retinal disorder with abnormal nocturnal vision. It has two forms depending upon the severity: complete form CSNB type I and an incomplete form CSNB type II. The main difference between the two is in the complete form; there are no measurable rod cells, whereas in the incomplete form some response is obtained due to the rod cells. In the complete form, cone activity is also affected [13, 14].

2.3. Cone-Rod Dystrophy (CRD) or Retinitis Pigmentosa (RP)

CRD also called as retinitis pigmentosa is an inherited retinal dystrophy with retinal pigment deposits visible on fundus examinations. The -wave and -wave amplitude are reduced. It is again classified into three types, namely, typical RP, early onset RP, and late stage RP. In the first type, the symptom is night blindness. In the second type, macular involvement occurs early and there is an involvement of rod cells which supports the diagnosis. In the third type, there is a decrease in the visual acuity and also macular involvement. The appearance of night blindness or loss of central vision supports the diagnosis [15, 16].

3. Analysis Methods

In this work various characteristic measures like Hurst exponent, Largest Lyapunov exponent, Higuchi’s fractal dimension, approximate entropy, and recurrence and phase-space plots are analyzed. A brief description of each of the parameters is given below.

3.1. Hurst Exponent

Hurst exponent is used to evaluate the long range dependence of data and its degree in a time series. Hurst exponent is the measure of smoothness of a fractal time series. It can also be defined as

where is the duration of the ERG sample of data and is the value of rescaled range. If , the time series acts as a random walk. If , the time series covers less distance than a random walk. If , the time series covers larger distance than a random walk.

To estimate the Hurst exponent, the dependence of the rescaled range on the time span is first estimated. A time series of length is divided into shorter time series of length , , . The average rescaled range is computed for each “” [17]. Step-by-step explanation of rescaled range calculation is as follows.

(i)Calculate the mean:

(ii)Create a mean-adjusted time series:

(iii)Compute the cumulative deviate series say :

(iv)Compute the range :

(v)Compute , standard deviation,

(vi)Calculate the rescaled range and average over all time series of “.”

3.2. The Largest Lyapunov Exponent (LLE)

Lyapunov exponent is used to distinguish between the periodic and chaotic signals. In phase space, the trajectories of chaotic dynamics follow typical patterns. It is the rate at which the neighboring trajectories separate from each other. A zero exponent indicates that the orbits maintain their relative positions. A negative exponent shows that the orbits approach a common point and the positive exponent shows that they are on chaotic attractor.

For any two points in a space say and , each of the points generates their own orbit in the space using a set of equations, where is the separation between the two orbits. This separation is the function of initial value . Then the Lyapunov exponent is given by

The largest Lyapunov exponent is computed by the least square fit to average line and is defined as

where is distance between th phase-space point and its nearest neighbors at th time and is the average overall phase-space points [18–20].

3.3. Approximate Entropy (ApEn)

Approximate entropy which is applied to the relatively short and noisy data is the logarithmic likelihood that the sample points which are close to each other will be same for the next comparison with a longer pattern. Smaller ApEn value shows that the signal is deterministic and higher ApEn value shows that the signal is random [21, 22]:

where is correlation integral, is pattern length, and is effective filter [23, 24].

3.4. Fractal Dimension (FD)

In traditional geometry, the Euclidean dimension of an object is referred to as the number of directions each differential of the object occupies in a space. The FD is used to provide a measure of how much space is occupied by an object between the Euclidean dimensions. In this work, we are using Higuchi’s algorithm for the analysis [25].

Let , be the time series to be analyzed. Let be the new time series:

where and is initial time value and is the discrete time interval.

For each , (length) is computed by

where is the total length of data , is normalization factor, and is .

This procedure is repeated for different values of ranging from 1 to , obtaining the average length. Fractal dimension is the slope of the least square linear best fit of the graph versus .

3.5. Recurrence Plots (RP)

Recurrence plot is a visualization technique which is used to detect hidden dynamical patterns and correlations in the data. In general, the RP reveals all those times at which the phase-space trajectory visits roughly the same area in the phase space. Suppose we have the time series representing the trajectory in the phase space with . RP is based on the following equation:

, .   Heaviside function. = norm. = predefined threshold.

Recurrence plots are mainly of four types, homogenous, drift, periodic, and disrupted. Homogenous RP are of stationary type and relaxation time is short with respect to the time spanned by the RP. Oscillating systems have RP (periodic) with diagonal lines and checkerboard structures. Slowly varying parameters with brightened RP at upper left and lower right corners are drift RP. Extreme events or sudden changes in the dynamics produce disrupted recurrence plots [26].

3.6. Surrogate Data Analysis

Surrogate data analysis is used to check the nonlinearity in the original data. Surrogate data is generated by phase randomizing the original dataset. The surrogate data has the same mean, variance, autocorrelation function, and similar spectral properties as of original data but phase relations are different. 15 surrogate series were generated from each original data series. Statistical significance is measured by comparing the experimental data with the surrogate data. If both results differ more than 50%, then the null hypothesis is rejected and it shows that the original data is nonlinear [27].

4. Results

Nonlinear parameters like Hurst exponent, approximate entropy, Higuchi’s fractal dimension, and the largest Lyapunov exponent are calculated for 30 Hz flicker signal and oscillatory potential of ERG signal. Tables 1 and 2 show the unique parameters in different cases along with the values.

Table 1: Nonlinear parameter values in different cases for 30 Hz flicker ERG signal.

Parameters
Controls
CSNB I
CSNB II
Typical RP
Early onset RP
Late stage RP
value

HE
0.5468 ± 0.042
0.3467 ± 0.04
0.2513 ± 0.0214
0.7768 ± 0.012
0.6987 ± 0.0314
0.9144 ± 0.0654
0.0051

ApEn
2.9527 ± 0.0301
2.6553 ± 0.0425
2.4934 ± 0.028
1.9464 ± 0.01
1.2970 ± 0.006
1.0904 ± 0.043
0.0124

HFD
1.8483 ± 0.0401
1.6526 ± 0.0912
1.4366 ± 0.0048
1.3370 ± 0.03
1.2464 ± 0.0062
1.0809 ± 0.012
0.123

LLE
0.4823 ± 0.0031
0.457 ± 0.001
0.41 ± 0.004
0.43 ± 0.0015
0.401 ± 0.007
0.422 ± 0.0017
0.067

Table 2: Nonlinear parameter values in different cases for oscillatory potential of ERG signal.

Parameters
Controls
CSNB I
CSNB II
Typical RP
Early onset RP
Late stage RP
value

HE
0.4551 ± 0.033
0.1120 ± 0.0112
0.1909 ± 0.0241
0.2805 ± 0.0306
0.3548 ± 0.016
0.8168 ± 0.0212
0.002

ApEn
1.8526 ± 0.005
1.4491 ± 0.0041
1.1947 ± 0.0203
0.95608 ± 0.052
0.7522 ± 0.004
0.6405 ± 0.031
0.034

HFD
1.8501 ± 0.0301
1.5448 ± 0.0165
1.4383 ± 0.026
1.0098 ± 0.003
1.3142 ± 0.0046
1.2696 ± 0.011
0.106

LLE
0.7090 ± 0.0015
0.7577 ± 0.0009
0.7134 ± 0.002
0.7043 ± 0.0052
0.741 ± 0.0064
0.732 ± 0.0041
0.079

The Hurst exponent which is in the range of 0-1 is used to measure the long range dependence of a time series. In each of the cases shown in Tables 1 and 2, the HE varies with a value less than 0.05. Approximate entropy which measures the disorder or predictability of the ERG signal shows higher values for control and the range of values reduces with each pathological case. The time lag used is 1 with value of 2. For ApEn value is taken as 15% of the standard deviation. The value is selected on the basis of previous studies which indicate good statistical validity [22]. HFD also showed reduced values for abnormal cases due to reduced rhythmic variation. Dimension value is taken as 9 with the delay of 1 in our study.

Figures 2, 3, and 4 show the Hurst exponent, LLE, and ApEn for the ERG signal. In our analysis, LLE value is a positive exponent which means that the signal is chaotic but the value is not distinguishable for different cases with the value greater than 0.05.

Figure 2: Variation of HE for six cases.

Figure 3: Variation of ApEn for six cases.

Figure 4: Variation of HFD for six cases.

Recurrence plot shows unique pattern in the case of control, CSNB, and CRD groups. In cone-rod dystrophy, also known as retinitis pigmentosa, recurrence plots show more squares in the case of both 30 Hz flicker input and oscillatory potential input (Figures 5 and 6). It indicates that, in CRD groups, a rhythmic variation causes periodicity and more patches of colors indicate variation of the signal.

(a) Control

(b) CSNB

(c) CRD

Figure 5: Recurrence plots for 30 Hz flicker input.

(a) Control

(b) CSNB

(c) CRD

Figure 6: Recurrence plots for oscillatory potential input.

5. Discussion and Conclusion

In this paper we have demonstrated and introduced the nonlinear analysis as a tool to evaluate the electroretinographic signal. All biological signals are random in nature and this randomness does not provide any time domain analysis and prediction. A nonlinear deterministic chaos theory and chaotic indicators are used for the analysis of ERG signal. In our work, quantitative schemes such as Hurst exponent, Lyapunov exponent, Higuchi’s fractal dimension, and approximate entropy are analyzed. Qualitative analysis like recurrence plots and phase-space plots are also computed for controls and different pathological groups. From our work, it can be concluded that all the nonlinear parameter values are higher for the control subjects and smaller values in the case of CSNB and CRD group which clearly shows the reduction of rhythmic variation in the pathological cases. In all the six cases (normal, CSNB I, CSNB II, typical RP, early onset RP, and late stage RP) the parameters like Hurst exponent, ApEn, HFD, and average recurrence show distinguishable numerical information from the ERG signal. From our study the largest Lyapunov exponent (LLE) shows positive value in all cases indicating the confirmation of the chaotic nature of the ERG signal. But LLE is not distinguishable for different groups. Recurrence plots and phase-space plots provide the visual inspection tool to assess the time evolution and the frequency of their recurrences. In our previous work, we have analyzed the ERG signal using wavelet analysis [10]. The method described in this paper is the analysis of the same signal from the nonlinear perspective of view. No work has yet been reported till now on the nonlinear aspect of the ERG signal analysis. Current work provides an excellent method for more advanced studies in the field of ophthalmology. Analysis and methods are individual dependent and the existing literatures prove that nonlinear methods are more effective in the analysis of biomedical signals like EEG (electroencephalogram) and HRV (heart rate variability) analysis.

Conflict of Interests

The authors declare that there is no conflict of interests.

Acknowledgment

The authors would like to thank the entire team of Little Flower Hospital and Research Centre, India, for supplying the data and also for answering the queries as a part of this study.

References

R. Barraco, L. Bellomonte, M. Brai, and D. Perrano Adorno, “Functional analysis of Normal and CSNB a-wave ERG component,” in IFMBE Proceedings, vol. 22, pp. 2668–2671, 2008.

R. Barraco, L. Bellamonte, and M. Brai, “Time-frequency behaviour of the a-wave of the human electroretinogram,” in IFMBE Proceedings, vol. 16, pp. 919–922, 2007.

J. Malmimo and R. Plonsey, Bioelectromagnetism: Principles and Applications of Bioelectric and Biomagnetic Fields, 1995.

M. F. Marmor, G. E. Holder, M. W. Seeliger, and S. Yamamoto, “Standard for clinical electroretinography (2004 update),” Documenta Ophthalmologica, vol. 108, no. 2, pp. 107–114, 2004.

H. Bornschein, G. Goodman, and R. D. Gunkel, “Temporal aspects of the human electroretinogram; a study of the implicit time-amplitude relationship of the b-wave,” A.M.A. Archives of Ophthalmology, vol. 57, no. 3, pp. 386–392, 1957.

R. Barraco, D. Persano Adorno, and M. Brai, “ERG signal analysis using wavelet transform,” Theory in Biosciences, vol. 130, no. 3, pp. 155–163, 2011.

T. Rogala and A. Brykalski, “Wavelet feature space in computer-aided electroretinogram evaluation,” Pattern Analysis and Applications, vol. 8, no. 3, pp. 238–246, 2005.

J. M. Miguel, L. Boquete, S. Ortega, C. A. Cordero, R. Barea, and R. Blanco, “MfERG_LAB: software for processing multifocal electroretinography signals,” Computer Methods and Programs in Biomedicine, vol. 108, no. 1, pp. 377–387, 2012.

J. M. Miguel-Jiménez, S. Ortega, L. Boquete, J. M. Rodríguez-Ascariz, and R. Blanco, “Multifocal ERG wavelet packet decomposition applied to glaucoma diagnosis,” BioMedical Engineering Online, vol. 10, article 37, 2011.

S. S. Nair and K. Paul Joseph, “Wavelet based electroretinographic signal analysis for diagnosis,” Biomedical Signal Processing and Control, vol. 9, pp. 37–44, 2014.

D. W. Crevier and M. Meister, “Synchronous period-doubling in flicker vision of salamander and man,” Journal of Neurophysiology, vol. 79, no. 4, pp. 1869–1878, 1998.

M. Molaie, R. Falahian, S. Gharibzadeh, S. Jafari, and J. C. Sprott, “Artificial neural networks: powerful tools for modeling chaotic behavior in the nervous system,” Frontiers in Computational Neuroscience, vol. 8, article 40, 2014.

C. Zeitz, U. Forster, J. Neidhardt et al., “Night blindness-associated mutations in the ligand-binding, cysteine-rich, and intracellular domains of the metabotropic glutamate receptor 6 abolish protein trafficking,” Human Mutation, vol. 28, no. 8, pp. 771–780, 2007.

C. Zeitz, “Molecular genetics and protein function involved in nocturnal vision,” Expert Review of Ophthalmology, vol. 2, no. 3, pp. 467–485, 2007.

C. P. Hamel, “Cone rod dystrophies,” Orphanet Journal of Rare Diseases, vol. 2, no. 1, article 7, 2007.

A. T. Moore, “Cone and cone-rod dystrophies,” Journal of Medical Genetics, vol. 29, no. 5, pp. 289–290, 1992.

S. Dangel, P. F. Meier, H. R. Moser, S. Plibersek, and Y. Shen, “Time series analysis of sleep EEG,” in Computer Assisted Physics, pp. 93–95, 1999.

A. Wolf, J. B. Swift, H. L. Swinney, and J. A. Vastano, “Determining Lyapunov exponents from a time series,” Physica D: Nonlinear Phenomena, vol. 16, no. 3, pp. 285–317, 1985.

A. Das, P. Das, and A. B. Roy, “Applicability of Lyapunov exponent in EEG data analysis,” Complexity International, vol. 9, pp. 1–9, 2002.

M. T. Rosenstein, J. J. Collins, and C. J. De Luca, “A practical method for calculating largest Lyapunov exponents from small data sets,” Physica D: Nonlinear Phenomena, vol. 65, no. 1-2, pp. 117–134, 1993.

U. R. Acharya, K. P. Joseph, N. Kannathal, C. M. Lim, and J. S. Suri, “Heart rate variability: a review,” Medical and Biological Engineering and Computing, vol. 44, no. 12, pp. 1031–1051, 2006.

S. M. Pincus, “Approximate entropy as a measure of system complexity,” Proceedings of the National Academy of Sciences of the United States of America, vol. 88, no. 6, pp. 2297–2301, 1991.

J. S. Richman and J. R. Moorman, “Physiological time-series analysis using approximate and sample entropy,” American Journal of Physiology: Heart and Circulatory Physiology, vol. 278, no. 6, pp. H2039–H2049, 2000.

S. M. Pincus and A. L. Goldberger, “Physiological time-series analysis: what does regularity quantify?” American Journal of Physiology: Heart and Circulatory Physiology, vol. 266, no. 4, pp. H1643–H1656, 1994.

W. Klonowski, E. Olejarczyk, and R. Stepien, “Sleep-EEG analysis using Higuchi’s fractal dimension,” in Proceedings of the International Symposium on Nonlinear Theory and Its Applications (NOLTA '05), pp. 222–225, Bruges, Belgium, 2005.

M. Mohebbi and H. Ghassemian, “Structures of the recurrence plot of heart rate variability signal as a tool for predicting the onset of paroxysmal atrial fibrillation,” Journal of Medical Signals and Sensors, vol. 1, no. 2, pp. 113–121, 2011.

J. Theiler, S. Eubank, A. Longtin, B. Galdrikian, and J. Doyne Farmer, “Testing for nonlinearity in time series: the method of surrogate data,” Physica D: Nonlinear Phenomena, vol. 58, no. 1–4, pp. 77–94, 1992.