British Journal of Anaesthesia 1991; 66: 81-87
SENSORY INFORMATION PROCESSING DURING GENERAL ANAESTHESIA: EFFECT OF ISOFLURANE ON AUDITORY EVOKED NEURONAL OSCILLATIONS C. MADLER, I. KELLER, D. SCHWENDER AND E. POPPEL
There is evidence from neuropsychological and psychophysical measurements that sensory information is processed in discrete time segments. The segmentation process may be described as neuronal oscillation at a frequency of 30-40 Hz. Stimulus-induced neuronal oscillations of this frequency are found in the middle latency range of the auditory evoked potential (AEP). We have studied the effect of different end-tidal concentrations of isoflurane on auditory evoked 30-40 Hz neuronal oscillations. We studied 13 patients undergoing intra-abdominal urological and gynaecological procedures. AEP were recorded in the awake state and during endexpiratory steady state isoflurane concentrations of 0.3, 0.6 and 1.2vol%. These incremental doses of isoflurane caused a stepwise decrease in frequency of oscillations. The decrease in oscillation frequency and sometimes the disappearance of oscillatory components may be interpreted as suppression of sensory information processing. The measurement of auditory evoked neuronal oscillations in the AEP appears to be a promising tool to monitor both sensory information processing capacity and depth of anaesthesia. KEY WORDS Anaesthesia: general. Anaesthetics, Brain: auditory evoked potentials.
volatile:
isoflurane.
The aim of general anaesthesia is to suppress consciousness during surgery. However, it is not clear what mechanisms are responsible for conscious awareness, although several theories have been proposed [1—4]. Sensory information is processed in spatially separate areas within the
brain. This concept of segregation of function is particularly well established for vision [5,6]. Spatial segregation of function has an apparent parallel in the time domain, because it is suggested by several studies that information processing is not continuous in time and sensory input is processed in discrete successive steps. If the temporal order threshold is determined for visual, auditory and tactile stimuli, similar values are obtained for each modality (20-30 ms [7, 8]). It is interesting to note that these three modalities show almost identical order thresholds, although fusion thresholds for stimuli are very different [9]. Further examples of discontinuous processing of sensory information come from histograms of choice reaction time and latencies of saccadic and pursuit eye movements. Such histograms are often multimodal with intermodal distances of 30-40 ms [10-13]. The multimodalities may be interpreted as an expression of preferred time zones set up by neuronal relaxation oscillations [12], within which reactions are located—that is, temporally segregated. In accordance with these psychophysical and neuropsychological findings, extra- and intracellular recordings in neurones demonstrate the existence of a 40-Hz oscillatory activity in the brain [14—16], which can be observed also in the spontaneous electroencephalogram (EEG) [17-19] and in sensory evoked potentials of different modalities [20-23]. The dominance of the 40-Hz oscillation is represented prominently in the midlatency response of auditory evoked potentials [20, 24]. We C. MADLEH,
M.D.; D. SCHWENDER,
M.D.;
Institute
of
Anaesthesiology, University of Munich, Marchioninistr. 15, 8000 Munich, Germany. I. KELLER, PH.D; E. P6PPEL, PH.D.;
Institute of Medical Psychology, University of Munich, Goethestr. 31, 8000 Munich, Germany. Accepted for Publication: May 3, 1990.
Downloaded from http://bja.oxfordjournals.org/ at Karolinska Institutet on July 22, 2015
SUMMARY
BRITISH JOURNAL OF ANAESTHESIA
82
have found suppression of the 40-Hz oscillatory component after induction of anaesthesia with etomidate [24]. These results indicate a potential link between the oscillatory components in the auditory evoked potential as an expression of information processing and the functional state of the brain. The aim of this study was to evaluate qualitative and dose-dependent effects of the volatile anaesthetic isoflurane on neuronal oscillations monitored in the midlatency range of auditory evoked potentials (AEP).
Patients and anaesthesia
Following institutional Ethics Committee approval and informed consent, we studied 13 patients (aged 28-65 yr, ASA class II) undergoing intra-abdominal urological or gynaecological procedures. None of the patients had either otolaryngological or neurological disorders or was receiving any medication. Patients were premedicated orally with flunitrazepam 2 mg, Analysis of auditory evoked potentials 45—60 min before the beginning of the study. To For each situation, we analysed the AEP of eliminate nociceptive stimuli during the operative individual patients and an interindividual grand procedure, continuous lumbar extradural an- average (GA), calculated from the AEP of all 13 algesia (EA) was administered using 0.5% patients. Off-line data analysis was as follows: bupivacaine. The level and effectiveness of EA latencies of the peaks V, Na, and Pa and relative was tested by pinprick. Anaesthesia was induced amplitudes from peak Na to peak Pa were with thiopentone 5 mg kg"1 and pancuronium measured. Latencies of all peaks were compared O.lmgkg" 1 was given to produce muscular re- with the Wilcoxon test (two-tailed, minimum laxation. After tracheal intubation, the patient's level of significance 5%, including Bonferroni's lungs were ventilated with 50% oxygen and correction). For frequency analysis of the poisoflurane in air to maintain normocapnia. The tential, in addition to the classical Fast Fourier end-expiratory concentration of isoflurane was Transformation (FFT), a rank-correlation promonitored continuously with a multigas analyser cedure (RCF) was used. In contrast to other (Siemens). AEP were measured after 15 min frequency analysing methods, the ranksteady state concentrations of isoflurane 1.2, 0.6 correlation method provides also a high frequency and 0.3 vol % after a minimum time interval of resolution for short epoch signals such as the 45 min after anaesthetic induction. Systemic AEP. The method is based on a non-parametric arterial pressure was monitored invasively, and procedure for the detection of oscillations in time mean arterial pressure was maintained greater series [25, 26] and has been used successfully in than 80 mm Hg. Core temperature was measured another study [27] to determine main oscillation with an oesophageal probe and maintained greater frequencies perioperatively by recorded AEP. than 34.5 °C using heated blankets and infusion of The RCF was calculated for a time window of warmed solutions. approximately 80 ms, always beginning with peak Na. Recording of auditory evoked potentials
The electrodiagnostic system Pathfinder I (Nicolet Instr.) was used for acoustic stimulation, registration and analysis of evoked potentials. Rarefaction clicks of 0.1 ms and 70 dB above the
RESULTS
Auditory evoked potentials were recorded successfully in each patient. Figure 1 shows a rep-
Downloaded from http://bja.oxfordjournals.org/ at Karolinska Institutet on July 22, 2015
PATIENTS AND METHODS
normal hearing level were presented binaurally with a stimulation frequency of 9.3 Hz using acoustically shielded earphones (TDH 39). For recording, silver-silver chloride electrodes were positioned at Cz and Ax with Fpz as ground (according to the international 10-20 system). The impedance of all electrodes was kept less than 0.5 kfi. An epoch of 100 ms (bin width 0.2 ms) was bandpass filtered (1-1500 Hz) with an analog Butterworth-filter (roll-off 6 dB octave"1) and averaged across 1000 stimulus presentations. The recording procedure was controlled visually on a monitor, and an automatic artefact detector rejected signals greater than 96% of full scale. To guarantee reliability of the signal and correct transmission and transduction of the acoustic stimuli, evoked potentials without a brainstem response (peak V) were rejected also. AEP were measured in the awake patient and then at steady state end-expiratory concentrations. One AEP represents the averaged activity during approximately 1 min.
SENSORY INFORMATION PROCESSING DURING ANAESTHESIA
83
latency, brainstem-generated component served as a reference for a correct and successful transduction and transmission of the auditory stimulation. The brainstem potential was followed by the so-called No/Po complex in a latency range between 10 and 15 ms after stimulation. The No/Po potentials were contaminated occasionally with sonomotor potentials in the awake patients. In the MLAEP, beginning 15-20 ms after stimulus application, a sequence of negative and positive peaks depicts a stimulus-dependent, oscillatory process. The peaks Na and Pa were identified easily in all cases. The following components, termed Nb, PI, Nl and P2 were detectable only in the awake state or during low concentrations of isofiurane. 40 60 80 100 AEP and corresponding power spectra (Fast Time (ms) FIG. 1. Auditory evoked potential (AEP) of an awake patient. Fourier Transformation) of the midlatency BAEP = Brainstem-generated auditory evoked potential; oscillations for a characteristic single patient in MLAEP = midlatency auditory evoked potential. Arrows inthe awake state and during the diree different dicate peaks. isofiurane concentrations are shown in figure 2. Midlatency components beginning with peak Na and Pa were obviously reduced in amplitude and resentative AEP of a premedicated, awake patient prolonged in latency with increasing doses of immediately before induction of anaesthesia. It isofiurane, whereas the amplitude and configuration consists of two distinct components: a brainstem- of brainstem potential V did not change with the generated AEP (BAEP) and a midlatency AEP level of anaesthesia. The FFT assigns a dominant (MLAEP). For the BAEP, a prominent peak with frequency of 40 Hz to the midlatency oscillation a latency of 5.8—6.2 ms (corresponding to the of the evoked potential. Under general anaesthesia brainstem potential V of Jewett's nomenclature) using 0.3% isoflurane, the observable increase, in was found in each patient. This robust short- peak latencies resulted in a decrease in the leading AEP
Corresponding power spectra
Awake Isofiurane 0.3 vol% 0.6 vol% 1.2 vol%
r- +
1.25 uV 1
20
i
i
40
I
i
60
Time (ms)
i
i
80
I
I
I
I
I
I
I
10120 30 40 50 60 70 80 Frequency (Hz)
FIG. 2. Auditory evoked potentials (AEP) and corresponding power spectra from a single patient. The graph shows data from the awake state before induction of anaesthesia and during different concentrations of isoflurane. Note the slowing down and disappearance of the oscillatory component from the midlatency range and stability of the brainstem evoked potentials under anaesthesia, the latter indicating intact sensory transduction mechanisms.
Downloaded from http://bja.oxfordjournals.org/ at Karolinska Institutet on July 22, 2015
AEP(n= 1)
BRITISH JOURNAL OF ANAESTHESIA
84 AEP Grand average (n = 13)
Corresponding power spectra
Awake Isoflurane 0.3 vol% 0.6 vol% 1.2 vol% 0.03 20
i l l l 80 40 60 Time (ms)
1 I I I I 1 I 20 30 40 50 60 70 80 Frequency (Hz)
FIG. 3. Grand averages of auditory evoked potentials (AEP) and their power spectra from 13 patients. Data are shown from the awake state before induction of anaesthesia and for different concentrations of isoflurane.
oscillation frequency to 30 Hz. This oscillation frequency of the MLAEP decreased again with 0.6% isoflurane to an oscillation of 10 Hz and remained constant with 1.2% isoflurane. Figure 3 shows the interindividual grand averages of AEP recorded from all 13 patients. The results of the grand averages confirm the suppression of amplitude and reduction of oscillation frequency seen in the single patient. Oscillation frequency again decreased in a stepwise manner. The first suppression in os60 -
|
^
i-
i Q 1-2 vol%
1
1
30-
•til
2 ™ 20 •
10 0 -
i
_E
3.0
0.6 vol%
40 -
E o
*
i—
| Awake
H | 0.3 vol%
50 •Z
cillation frequency was associated with isoflurane concentrations of 0.3%. Isoflurane concentrations of 0.6vol% caused a further reduction of midlatency oscillations, whereas 1.2% isoflurane resulted in their complete abolition. Mean values and SEM for the latencies of the peaks V, Na and Pa, and Na/Pa amplitudes are presented infigure4. The latencies of the peaks V, Na and Pa had a low interindividual variability. A graded and statistically significant (P < 0.05) increase with isoflurane concentration was observed
i
mi
I 491 Na Peak latencies
m
1 m m 1^ /
Pa
- 2.0 CD
•o 3
"a <
Na/Pa Amplitude
FIG. 4. Means and SEM ( • ) of peak latencies V, Na, Pa (left side) and Na/Pa amplitudes (right side). *Significant differences between experimental conditions tested with the Wilcoxon test for matched samples: *P^0.05; * * P < 0 . 0 1 .
Downloaded from http://bja.oxfordjournals.org/ at Karolinska Institutet on July 22, 2015
0.62 uV
SENSORY INFORMATION PROCESSING DURING ANAESTHESIA 60
n
50 40
|
30
2 20 A 10 0 0.3
0.6
1.2
Isoflurane concn (vol%)
FIG. 5. Means, SEM of oscillation frequency measured with the rank-correlation function for each experimental situation. Significant difference between conditions: *P < 0.05 (Wilcoxon test for matched samples).
for both Na and Pa latencies. A significant (P < 0.05) increase in the short latency brainstem potentials (peak V) was detectable only between the awake state and the 0.3 % isoflurane. A further graded change in latency for this component under increasing concentrations of isoflurane was not seen. The Na/Pa amplitudes were significantly (P < 0.05) reduced. Figure 5 shows die means and SEM of the dominant oscillation frequency for all 13 patients ••according to the concentration of isoflurane. Mean oscillation frequency (P < 0.05) decreased significantly from about 45 Hz in the awake state to 29 Hz with 0.3% isoflurane. A further increase of isoflurane concentration to 0.6% resulted in a small and insignificant shift to 28 Hz. With 1.2% isoflurane, oscillation frequency of the MLAEP decreased significantly (P < 0.05), to 23 Hz. DISCUSSION
In all patients, an oscillation frequency of approximately 40 Hz was observed in the awake state. This finding confirms the existence of a stimulus-dependent and interindividual constant oscillatory brain mechanism for sensory information processing. It is not clear if these oscillations are caused by a single generator or by a suprasegmental chain of generators. Studies have suggested, however, that oscillations of the MLAEP represent the activation of several generators situated at different brain levels [1]. Isoflurane is a non-specific acting general
anaesthetic agent and induces a dose-dependent suppression of auditory evoked neuronal oscillations. The suppression of neuronal oscillations was reflected by prolongation of latencies and reduction of amplitudes of the MLAEP components, mainly the peaks Na and Pa. The observed influence of volatile anaesthetic agents on latencies and amplitudes of the MLAEP has been demonstrated previously [28-30]. Changes in these parameters have been proposed, therefore, as an objective criterion for depth of anaesthesia. Our results suggest that direct measurement of neuronal oscillations by frequency analysis seems preferable, for the following reasons. Latencies and amplitudes of MLAEP components are connected only indirectly to these oscillations. Peaks Na and Pa, for instance, are prolonged in latency with each increment of isoflurane concentration. This does not necessarily represent disturbed sensory information processing. Prolongation of midlatency peaks might also be caused by slowed transmission of the acoustic stimulus from the receptor to the primary auditory cortex. Although Na and Pa latencies gradually increased from the awake state to 1.2% isoflurane, the difference in Na-Pa latency for 0.3 and 0.6% was the same, indicating a change in phase but not in frequency of the neuronal oscillation. This observation was confirmed by frequency analysis, which showed that neuronal oscillations were reduced significantly for 0.3 and 1.2 % isoflurane only. The fact that oscillation frequency did not change with increasing concentrations from 0.3 to 0.6% isoflurane indicates that these doses did not suppress significantly those structures involved in the oscillatory process. As oscillation frequency of the MLAEP depends mainly on the Na/Pa complex, it could be expected that the underlying generators of these midlatency peaks are unaffected with concentrations of isoflurane less than or equal to 0.6%. It has been proposed recently [31,32] that the Na/Pa complex is generated in the primary auditory cortex of the temporal lobe. Therefore, despite prolonged peak latencies, it can be assumed that concentrations less than or equal to 0.6% do not suppress sensory information processing of acoustic stimuli at the level of the primary auditory cortex. Amplitudes are reduced gradually by isoflurane also, but have a greater interindividual variability compared with latency and oscillation measures.
Downloaded from http://bja.oxfordjournals.org/ at Karolinska Institutet on July 22, 2015
Patient awake
85
86
General anaesthetic agents seem to be a useful tool to test the neuropsychological hypotheses of consciousness. With the presented neuropsychological model, suppression of the neuronal oscillations should allow the stimulus to enter the brainstem, but any higher cognitive function, such as perception or even retention in memory, should be blocked. The latter has yet to be tested. In this context, it would be of interest to investigate the behaviour of the neuronal oscillator under anaesthetic agents with psychomimetic effects such as ketamine or during anaesthetic regimens associated with intraoperative awareness such as high-dose fentanyl anaesthesia [38]. In the present study, AEP with 0.6% isoflurane still exhibited midlatency oscillations of about 30 Hz, despite prolonged latencies. This corresponds well with clinical observations of inadequate anaesthesia under these concentrations. Therefore, to guarantee suppression of consciousness and sufficient "depth of anaesthesia", isoflurane should be used only in concentrations exceeding 0.6%. REFERENCES 1. Deshmukh VDA. A chain of suprasegmenul neuroscillatory circuits: A human brain theory. Clinical EUctroencephalography 1988; 19: 7-13. 2. Webb A. Consciousness and the cerebral cortex. British Journal of Anaesthesia 1983; 55: 209-217. 3. P6ppel E. Taxonomy of the subjective: an evolutionary perspective. In: Brown JW, ed. Neuropsychology of Visual Perception. Hillsdale: Lawrence Erlbaum Associates, 1989: 219-232. 4. Edclman G. 77K Remembered Present. A Biological Theory of Consciousness. New York: Basic Books 1990. 5. Zeki S. Functional specialisation in the visual cortex of the rhesus monkey. Nature (London) 1978; 274: 423-428. 6. De Yoe EA, van Essen DC. Concurrent processing streams in monkey visual cortex. Trends in Neuroscience 1988;5: 219-226. 7. Hirsh IJ, Sherrick CE jr. Perceived order in different sense modalities. Journal of Experimental Psychology 1961; 62: 423-^32. 8. Ilmberger J. Auditory excitability cycles in choice reaction time and order threshold. Naturwissenschaften 1986; 73: 743-744. 9. POppel E. Time perception. In: Held R, Leibowitz R, Teuber HW, eds. Handbook of Sensory Physiology, Vol. VIII. Berjin, Heidelberg, New York: Springer, 1978. 10. Poppcl E. Oszillatorijche Komponenten in Reaktionszeiten. Naturwissenschaften 1968; 5: 449-450. 11. Harter R, White CT. Periodicity within reaction time distributions and electromyograms. Quarterly Journal of Experimental Psychology 1968; 20: 157-166. 12. P6ppel E. Excitability cycles in central intermittency. Psychologische Forschung 1970; 34: 1-9.
Downloaded from http://bja.oxfordjournals.org/ at Karolinska Institutet on July 22, 2015
This makes it difficult to evaluate adequate anaesthesia with this parameter. The only way to use amplitudes as an indicator of anaesthetic depth would be to compare pre- and intraoperative amplitudes individually. Nevertheless, amplitudes are very sensitive to biological and technical influences. Interindividual variability is caused by subjective hearing level, mass of the generator and distance between generator and electrode. Moreover, intraindividual changes might occur especially during operative procedures such as extracorporal circulation in cardiac surgery. In these situations, alterations of skin resistance caused by sweating or reduced skin perfusion are quite common. A suppression of amplitudes, therefore, seems only a vague predictor of the functional state of the brain. In addition to these practical implications, there are theoretical reasons for preferring measurement of neuronal oscillations as an indicator of the state of the brain under anaesthesia. The proposed neuropsychological model underlying sensory information processing is based on the existence of a time segmentation process. This process synchronizes sensory information from different modalities and by doing so allows the organism to identify sensory events in a time ordered structure [3,33]. The clinical observation that patients often report no lapse of time during anaesthesia is relevant in this context [24]. It can be speculated that the suppression of the neuronal oscillations results in a disturbance of time perception or even in a subjective "loss" of time. This impression is different from the evaluation of time after natural sleep, which may be quite accurate [34]. In agreement with this observation, auditory evoked 40-Hz responses could be recorded during most phases of night sleep [35]. According to the underlying model of sensory information processing and temporal organization of perception [3, 36, 37], not only time perception but also consciousness in general should be interrupted when the oscillating activity is suppressed. One way of testing this hypothesis is to gain insight into the physiological mechanism of the proposed neuronal oscillations using pharmacological intervention. The anatomical localization of the neuronal oscillation generator has not yet been identified, but it is proposed that one or more generators localized in the midbrain and cortical structures produce the 40-Hz response in the EEG and sensory evoked potentials of different modalities [1].
BRITISH JOURNAL OF ANAESTHESIA
SENSORY INFORMATION PROCESSING DURING ANAESTHESIA
27. Keller I, Madler C, Schwender S, Pdppel E. Oscillatory components in perioperative AEP-recordings: A nonparametric procedure for frequency measurement. Clinical Electroencephalography 1990; (In Press). 28. Heneghan CPH, Thornton C, Navaratnarajah M, Jones JG. The effect of isofluranc on the auditory evoked response in man. British Journal of Anaesthesia 1987; 59: 277-282. 29. Jones JG. Use of evoked responses in the EEG to measure depth of anaesthesia. In: Rosen M, Lunn JN, eds. Consciousness Awareness and Pain in General Anaesthesia. London: Butterworth & Co., 1987; 98-139. 30. Thornton C, Heneghan CPH, James MFM, Jones JG. The effect of halothane and enflurane with controlled ventilation on the auditory evoked potential. British Journal of Anaesthesia 1984; 56: 323-325. 31. Deiber MP, Ibanez V, Fischer C, Perrin F, Maugiere F. Sequential mapping favours the hypothesis of distinct generators for Na and Pa middle latency auditory evoked potentials. Electroencephalography and Clinical Neurophysiology 1988; 71: 187-197. 32. Scherg M, Hari R, Hamalainen M. Frequency-specific sources of the auditory N19-P30 detected by a multiple source analysis of evoked magnetic fields and potentials. 1 In: Proceedings of the 7th International Conference on Biomagnetism, New York 1989. 33. P6ppel E, Schill K, von Steinbuchel N. Sensory integration within temporally neutral system states: a hypothesis. Naturwissenschaften 1990; 77: 89-91. 34. dauser G. Die Kopfuhr. Das automatischc Erwachen. Stuttgart: Enke, 1984. 35. Erwin R, Buchwald JS. Midlatency auditory evoked responses: differential effects of sleep in the human. Electroencephalography and Clinical Neurophysiology 1986; 65: 383-392. 36. Poppel E. Time perception. In: Adelman G, ed. Encyclopedia of Neuroscience. Boston: BirkhSuser, 1987; 1215-1216. 37. Poppel E. Mindworks. Time and Conscious Experience. Boston: Harcourt Brace Jovanovich, 1988. 38. Mummaneni N, Rao TLK, Montoya A. Awareness and recall with high-dose fentanyl-oxygen anesthesia. Anesthesia and Analgesia 1980; 59: 948-949.
Downloaded from http://bja.oxfordjournals.org/ at Karolinska Institutet on July 22, 2015
13. Pfippel E, Logothetis N. Neuronal oscillations in the human brain. Naturwissenschaften 1986; 73: 267-268. 14. Llinas RR. Neuronal oscillations in mammalian brain. In: Cohenand M, Sturmwasser F, eds. Comparative Neurology. Chichester: John Wiley & Sons, 1985; 279-290. 15. Ecknorn A, Bauer R, Jordan W, Brosch M, Kruse W, Munk M, Reitboeck HJ. Coherent oscillations: A mechanism of feature unking in the visual cortex? Biological Cybernetics 1988; 60: 121-130. 16. Gray CM, K6nig P, Engel AK, Singer W. Oscillatory responses in cat visual cortex exhibit inter-columnar synchronisation which reflects global stimulus properties. Nature (London) 1989; 338: 334-337. 17. Basar E, Rosen B, Basar-Eroglu C, Greitschus F. The associations between 40 Hz-EEG and the middle latency response of the auditory evoked potential. International Journal of Neuroscience 1987; 33: 103-117. 18. Loring DW, Sheer DE. Laterality of 40 Hz EEG and EMG during cognitive performance. Psychophysiology 1984; 21: 34-38. 19. Sheer DE. Biofeedback training of 40 Hz EEG activity and behaviour. In: Burch N, Altschuler H, eds. Behaviour and Brain Electrical Activity. New York: Plenum Press, 1975. 20. GalambosR, MakeigS,TalmachoffPJ. A 40-Hz auditory potential recorded from the human scalp. Proceedings of the National Academy of Science U.S.A. 1981; 78: 2643-2647. 21. MIkela JP, Hari R. Evidence for the cortical origin of the 40-Hz auditory evoked response in man. Electroencephalography and Clinical Neurophysiology 1987; 66: 539-546. 22. Regan D. A high frequency mechanism which underlies VEP. Electroencephalography and Clinical Neurophysiology 1968; 25: 231-237. 23. Spehlmann R. Evoked Potential Primer. Boston: Butterworth Publishers, 1985. 24. Madler C, Pdppel E. Auditory evoked potentials indicate the loss of neuronal oscillations during general anaesthesia. Naturwissenschaften 1987; 74: 42^13. 25. P6ppel E. Frequency measurement in time series data. Life Sciences and Space Research 1970; 8: 234-238. 26. Miescke KJ, PSppel E. A nonparametric procedure to detect periods in time series. Stochastic Processes and Their Application 1982; 13: 319-325.
87