Accumulated source imaging of brain activity with both low and high-frequency neuromagnetic signals

Jing Xiang, Qian Luo, Rupesh Kotecha, Abraham Korman, Fawen Zhang, Huan Luo, Hisako Fujiwara, Nat Hemasilpin, Douglas F Rose, Jing Xiang, Qian Luo, Rupesh Kotecha, Abraham Korman, Fawen Zhang, Huan Luo, Hisako Fujiwara, Nat Hemasilpin, Douglas F Rose

Abstract

Recent studies have revealed the importance of high-frequency brain signals (>70 Hz). One challenge of high-frequency signal analysis is that the size of time-frequency representation of high-frequency brain signals could be larger than 1 terabytes (TB), which is beyond the upper limits of a typical computer workstation's memory (<196 GB). The aim of the present study is to develop a new method to provide greater sensitivity in detecting high-frequency magnetoencephalography (MEG) signals in a single automated and versatile interface, rather than the more traditional, time-intensive visual inspection methods, which may take up to several days. To address the aim, we developed a new method, accumulated source imaging, defined as the volumetric summation of source activity over a period of time. This method analyzes signals in both low- (1~70 Hz) and high-frequency (70~200 Hz) ranges at source levels. To extract meaningful information from MEG signals at sensor space, the signals were decomposed to channel-cross-channel matrix (CxC) representing the spatiotemporal patterns of every possible sensor-pair. A new algorithm was developed and tested by calculating the optimal CxC and source location-orientation weights for volumetric source imaging, thereby minimizing multi-source interference and reducing computational cost. The new method was implemented in C/C++ and tested with MEG data recorded from clinical epilepsy patients. The results of experimental data demonstrated that accumulated source imaging could effectively summarize and visualize MEG recordings within 12.7 h by using approximately 10 GB of computer memory. In contrast to the conventional method of visually identifying multi-frequency epileptic activities that traditionally took 2-3 days and used 1-2 TB storage, the new approach can quantify epileptic abnormalities in both low- and high-frequency ranges at source levels, using much less time and computer memory.

Keywords: brain; high-frequency oscillations; magnetic source imaging; magnetoencephalography; multi-frequency.

Figures

Figure 1
Figure 1
Workflow for computing accumulated spectrogram (left) and the basic principle of computing accumulated spectrogram (right). Since the analysis of high-frequency MEG signals requires high-sampling rate MEG data, MEG data are digitized in a high-frequency range. To improve the performance and optimize the use of computer memory for analyzing both low- and high-frequency MEG signals, the new method can re-sample MEG data dynamically according to the analysis frequency ranges. If the data points of the recorded data are smaller than the minimum data point of wavelet transform in frequency range, the “Data Padding” function can pad some data points so as to meet the requirements of wavelet transform. The “Thresholding” indicates that a spectral value can be rejected or accepted by the accumulated spectrogram according to a threshold value. MEG data recorded are waveforms, which are divided to segments (e.g., “Waveform 1,” “Waveform N”) to minimize the use of memory for wavelet transform (“Wavelet transform”). In the new method, wavelet transform transfers each segment of waveform data to a spectrum (e.g., “Spectrum 1,” “Spectrum N”). Of note, “N” indicates the total number of segments or spectra, which can be theoretically infinitely large. The “+” indicates the process of accumulation, which add all spectra together to produce an accumulated spectrum (“Accumulated Spectrum”). The left view of the sensor distribution of our MEG system is shown on the top right.
Figure 2
Figure 2
Workflow for computing accumulated source images. The workflow includes two main computing pipelines. One computing pipeline processes MEG data with filter or wavelet transforms so as to generate multi-frequency datasets. MEG signals in multi-frequency datasets are in a set of frequency ranges. Another computing pipeline works on several tasks, which included the creation of a three-dimensional source grid (3D grid), performing forward solution by calculating lead fields, ranking the norm for each source for all sensors, and performing SVD. The node-beam lead field is completed by selecting a group of sensors which have a larger lead field norm (or weights). Of note, each location in accumulated source imaging can have multiple parameters (e.g., “Frequency Index,” “Source Strength”). Some processes are optional (red lines) and additional parameters can also be added to the workflow.
Figure 3
Figure 3
Accumulated global spectrograms in three conditions. “Magnetic Noise” was computed with MEG data recorded without subjects. “Control Subject” was computed with MEG data recorded from a healthy child. “Epilepsy Subject” was computed with MEG data recorded from a child with epilepsy between seizures (interictal). The sampling rate of all MEG recordings was 6000 Hz. An accumulated global spectrogram represents the “spatial summation” of the entire MEG sensor array accumulated spectrograms. The three spectrograms show that the epilepsy subject has elevated spectral power as compared to the control subject.
Figure 4
Figure 4
Accumulated spectrograms show the well-known alpha activity in a healthy subject and an epilepsy subject. Noticeably, healthy subject (“Healthy Subject”) has a clear activity around 8–12 Hz (alpha activity). However, the epilepsy subject (“Epilepsy Subject”) has incrased activity in 2–4 Hz (low-frequency activity). The color bar shows the color coding of spectral power.
Figure 5
Figure 5
Global accumulated spectrograms from 10 epilepsy subjects and 1 healthy subject show the main frequency components of neuromagnetic signals in 70–200 Hz in epilepsy patients. Noticeably, the activity patterns vary across patients. The color bar shows the color coding of spectral power for all the global accumulated spectrograms.
Figure 6
Figure 6
Accumulated spectral contour maps from 10 epilepsy subjects and 1 healthy subject show the spatial distributions. Noticeably, the spectral distribution varies across patients. All the contour maps have the same orientation defined by the arrows: the “L” indicates the left side of the head and the “R” indicates the right side of the head. The “F” indicates that the upper part of the contour map represents the frontal region of the head; the “B” indicates that the lower part of the contour map represents the posterior region of the head. Each small circle represents one MEG sensor. The color bar shows the color coding of spectral power for all the contour maps.
Figure 7
Figure 7
An illustration of the basic principle of accumulated source imaging. The top waveforms show MEG data at sensor levels. The bottom images show individual structural magnetic resonance image and the region of interest (ROI, blue lines) for source scanning. MEG sensor data are firstly divided into small segments (e.g., “Sensor Data Segment 1,” “Sensor Data Segment 2,” “Sensor Data Segment N”). Volumetric sources are then produced by scanning the entire ROI with each segment of sensor data. The red, yellow and white small cubes indicate the sources (or voxels) identified. For illustration purposes, a very low resolution (12 millimeter) spatial resolution was used. An accumulated source image is generated by spatially adding all volumetric sources together. Of note, only sources reach certain thresholds (in this case, 75%) are added to accumulated source images, which differentiate this accumulating process from averaging. The color bar indicates the color coding of the source strength.
Figure 8
Figure 8
Accumulated source imaging shows low frequency brain activity in 8–12 Hz (alpha) in an epilepsy subject (“Epilepsy Subject”) and a healthy subject (“Control Subject”). Alpha activity is localized to the occipital cortex in the healthy subject. However, alpha activity is overshadowed by epileptic activity in the epilepsy subject. The epileptic activity is localized to the left and right parietal cortices in the epilepsy subject, which is concordant with clinical findings.
Figure 9
Figure 9
A digital photo of intracranial recording (“ECoG”) and accumulated source imaging (“ASI”) show the concordance of the two methods. The two images are placed in the similar orientation. “Frontal” indicate the frontal cortex; “Temporal” indicates temporal lobe. The left green arrow points to the epileptic area invasively defined by intracranial recording; the right blue arrow points to the epileptic area noninvasively localized with high-frequency neuromagnetic signals. Noticeably, the areas are matched in gyrus level. The color bar shows the color coding of accumulated source imaging. The value of the source voxel is normalized T value (no unit).

References

    1. Agirre-Arrizubieta Z., Huiskamp G. J., Ferrier C. H., Van Huffelen A. C., Leijten F. S. (2009). Interictal magnetoencephalography and the irritative zone in the electrocorticogram. Brain 132, 3060–3071 10.1093/brain/awp137
    1. Andrade-Valenca L., Mari F., Jacobs J., Zijlmans M., Olivier A., Gotman J., et al. (2012). Interictal high frequency oscillations (HFOs) in patients with focal epilepsy and normal MRI. Clin. Neurophysiol. 123, 100–105 10.1016/j.clinph.2011.06.004
    1. Ball T., Kern M., Mutschler I., Aertsen A., Schulze-Bonhage A. (2009). Signal quality of simultaneously recorded invasive and non-invasive EEG. Neuroimage 46, 708–716 10.1016/j.neuroimage.2009.02.028
    1. Benar C. G., Chauviere L., Bartolomei F., Wendling F. (2010). Pitfalls of high-pass filtering for detecting epileptic oscillations: a technical note on “false” ripples. Clin. Neurophysiol. 121, 301–310 10.1016/j.clinph.2009.10.019
    1. Benbadis S. R., Lafrance W. C., Jr., Papandonatos G. D., Korabathina K., Lin K., Kraemer H. C., et al. (2009). Interrater reliability of EEG-video monitoring. Neurology 73, 843–846 10.1212/WNL.0b013e3181b78425
    1. Blanco J. A., Stead M., Krieger A., Stacey W., Maus D., Marsh E., et al. (2011). Data mining neocortical high-frequency oscillations in epilepsy and controls. Brain 134, 2948–2959 10.1093/brain/awr212
    1. Brinkmann B. H., Bower M. R., Stengel K. A., Worrell G. A., Stead M. (2009). Large-scale electrophysiology: acquisition, compression, encryption, and storage of big data. J. Neurosci. Methods 180, 185–192 10.1016/j.jneumeth.2009.03.022
    1. Chen Y., Xiang J., Kirtman E. G., Wang Y., Kotecha R., Liu Y. (2010). Neuromagnetic biomarkers of visuocortical development in healthy children. Clin. Neurophysiol. 121, 1555–1562 10.1016/j.clinph.2010.03.029
    1. Crepon B., Navarro V., Hasboun D., Clemenceau S., Martinerie J., Baulac M., et al. (2010). Mapping interictal oscillations greater than 200 Hz recorded with intracranial macroelectrodes in human epilepsy. Brain 133, 33–45 10.1093/brain/awp277
    1. Dalal S. S., Guggisberg A. G., Edwards E., Sekihara K., Findlay A. M., Canolty R. T., et al. (2008). Five-dimensional neuroimaging: localization of the time-frequency dynamics of cortical activity. Neuroimage 40, 1686–1700 10.1016/j.neuroimage.2008.01.023
    1. De Gooijer-Van De Groep K. L., Leijten F. S., Ferrier C. H., Huiskamp G. J. (2013). Inverse modeling in magnetic source imaging: comparison of MUSIC, SAM(g2), and sLORETA to interictal intracranial EEG. Hum. Brain Mapp. 34, 2032–2044 10.1002/hbm.22049
    1. De Munck J. C., Bijma F. (2009). Three-way matrix analysis, the MUSIC algorithm and the coupled dipole model. J. Neurosci. Methods 183, 63–71 10.1016/j.jneumeth.2009.06.040
    1. Dumpelmann M., Jacobs J., Kerber K., Schulze-Bonhage A. (2012). Automatic 80–250 Hz “ripple” high frequency oscillation detection in invasive subdural grid and strip recordings in epilepsy by a radial basis function neural network. Clin. Neurophysiol. 123, 1721–1731 10.1016/j.clinph.2012.02.072
    1. Engel J., Jr., Bragin A., Staba R., Mody I. (2009). High-frequency oscillations: what is normal and what is not? Epilepsia 50, 598–604 10.1111/j.1528-1167.2008.01917.x
    1. Fatima Z., Quraan M. A., Kovacevic N., McIntosh A. R. (2013). ICA-based artifact correction improves spatial localization of adaptive spatial filters in MEG. Neuroimage 78, 284–294 10.1016/j.neuroimage.2013.04.033
    1. Ghuman A. S., Mcdaniel J. R., Martin A. (2011). A wavelet-based method for measuring the oscillatory dynamics of resting-state functional connectivity in MEG. Neuroimage 56, 69–77 10.1016/j.neuroimage.2011.01.046
    1. Gotman J. (2010). High frequency oscillations: the new EEG frontier? Epilepsia 51(Suppl. 1), 63–65 10.1111/j.1528-1167.2009.02449.x
    1. Gramfort A., Luessi M., Larson E., Engemann D.A., Strohmeier D., Brodbeck C., et al. (2014). MNE software for processing MEG and EEG data. Neuroimage 86, 446–460 10.1016/j.neuroimage.2013.10.027
    1. Gross J., Baillet S., Barnes G. R., Henson R. N., Hillebrand A., Jensen O., et al. (2013). Good practice for conducting and reporting MEG research. Neuroimage 65, 349–363 10.1016/j.neuroimage.2012.10.001
    1. Guggisberg A. G., Dalal S. S., Findlay A. M., Nagarajan S. S. (2007). High-frequency oscillations in distributed neural networks reveal the dynamics of human decision making. Front. Hum. Neurosci 1:14 10.3389/neuro.09.014.2007
    1. Gummadavelli A., Wang Y., Guo X., Pardos M., Chu H., Liu Y., et al. (2013). Spatiotemporal and frequency signatures of word recognition in the developing brain: a magnetoencephalographic study. Brain Res. 1498, 20–32 10.1016/j.brainres.2013.01.001
    1. Haegelen C., Perucca P., Chatillon C. E., Andrade-Valenca L., Zelmann R., Jacobs J., et al. (2013). High-frequency oscillations, extent of surgical resection, and surgical outcome in drug-resistant focal epilepsy. Epilepsia 54, 848–857 10.1111/epi.12075
    1. Hamalainen M. S., Sarvas J. (1987). Feasibility of the homogeneous head model in the interpretation of neuromagnetic fields. Phys. Med. Biol. 32, 91–97 10.1088/0031-9155/32/1/014
    1. Huang M. X., Mosher J. C., Leahy R. M. (1999). A sensor-weighted overlapping-sphere head model and exhaustive head model comparison for MEG. Phys. Med. Biol. 44, 423–440 10.1088/0031-9155/44/2/010
    1. Huang M. X., Shih J. J., Lee R. R., Harrington D. L., Thoma R. J., Weisend M. P., et al. (2004). Commonalities and differences among vectorized beamformers in electromagnetic source imaging. Brain Topogr. 16, 139–158 10.1023/B:BRAT.0000019183.92439.51
    1. Jacobs J., Staba R., Asano E., Otsubo H., Wu J. Y., Zijlmans M., et al. (2012). High-frequency oscillations (HFOs) in clinical epilepsy. Prog. Neurobiol. 98, 302–315 10.1016/j.pneurobio.2012.03.001
    1. Jacobs J., Zelmann R., Jirsch J., Chander R., Dubeau C. E., Gotman J. (2009). High frequency oscillations (80–500 Hz) in the preictal period in patients with focal seizures. Epilepsia 50, 1780–1792 10.1111/j.1528-1167.2009.02067.x
    1. Jerbi K., Freyermuth S., Dalal S., Kahane P., Bertrand O., Berthoz A., Lachaux J. P. (2009). Saccade related gamma-band activity in intracerebral EEG: dissociating neural from ocular muscle activity. Brain Topogr. 22, 18–23 10.1007/s10548-009-0078-5
    1. Jirsch J. D., Urrestarazu E., Levan P., Olivier A., Dubeau F., Gotman J. (2006). High-frequency oscillations during human focal seizures. Brain 129, 1593–1608 10.1093/brain/awl085
    1. Kirsch H. E., Robinson S. E., Mantle M., Nagarajan S. (2006). Automated localization of magnetoencephalographic interictal spikes by adaptive spatial filtering. Clin. Neurophysiol. 117, 2264–2271 10.1016/j.clinph.2006.06.708
    1. Kotecha R., Xiang J., Wang Y., Huo X., Hemasilpin N., Fujiwara H., et al. (2009). Time, frequency and volumetric differences of high-frequency neuromagnetic oscillation between left and right somatosensory cortices. Int. J. Psychophysiol. 72, 102–110 10.1016/j.ijpsycho.2008.10.009
    1. Kovach C. K., Tsuchiya N., Kawasaki H., Oya H., Howard M. A., 3rd., Adolphs R. (2011). Manifestation of ocular-muscle EMG contamination in human intracranial recordings. Neuroimage 54, 213–233 10.1016/j.neuroimage.2010.08.002
    1. Le Van Quyen M., Staba R., Bragin A., Dickson C., Valderrama M., Fried I., et al. (2010). Large-scale microelectrode recordings of high-frequency gamma oscillations in human cortex during sleep. J. Neurosci. 30, 7770–7782 10.1523/JNEUROSCI.5049-09.2010
    1. Levesque M., Bortel A., Gotman J., Avoli M. (2011). High-frequency (80–500 Hz) oscillations and epileptogenesis in temporal lobe epilepsy. Neurobiol. Dis. 42, 231–241 10.1016/j.nbd.2011.01.007
    1. Matsumoto A., Brinkmann B. H., Matthew Stead S., Matsumoto J., Kucewicz M., Marsh W. R., et al. (2013). Pathological and physiological high-frequency oscillations in focal human epilepsy. J. Neurophysiol. 110, 1958–1964 10.1152/jn.00341.2013
    1. Montazeri N., Shamsollahi M. B., Hajipour S. (2009). MEG based classification of wrist movement. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2009, 986–989 10.1109/IEMBS.2009.5334472
    1. Mosher J. C., Leahy R. M. (1998). Recursive MUSIC: a framework for EEG and MEG source localization. IEEE Trans. Biomed. Eng. 45, 1342–1354 10.1109/10.725331
    1. Mosher J. C., Leahy R. M., Lewis P. S. (1999). EEG and MEG: forward solutions for inverse methods. IEEE Trans. Biomed. Eng. 46, 245–259 10.1109/10.748978
    1. Ou W., Hamalainen M. S., Golland P. (2009). A distributed spatio-temporal EEG/MEG inverse solver. Neuroimage 44, 932–946 10.1016/j.neuroimage.2008.05.063
    1. Pail M., Halamek J., Daniel P., Kuba R., Tyrlikova I., Chrastina J., et al. (2013). Intracerebrally recorded high frequency oscillations: simple visual assessment versus automated detection. Clin. Neurophysiol. 124, 1935–1942 10.1016/j.clinph.2013.03.032
    1. Papadelis C., Poghosyan V., Fenwick P. B., Ioannides A. A. (2009). MEG's ability to localise accurately weak transient neural sources. Clin. Neurophysiol. 120, 1958–1970 10.1016/j.clinph.2009.08.018
    1. Pulvermuller F., Birbaumer N., Lutzenberger W., Mohr B. (1997). High-frequency brain activity: its possible role in attention, perception and language processing. Prog. Neurobiol. 52, 427–445 10.1016/S0301-0082(97)00023-3
    1. Rau R., Raschka C., Koch H. J. (2002). Uniform decrease of alpha-global field power induced by intermittent photic stimulation of healthy subjects. Braz. J. Med. Biol. Res. 35, 605–611 10.1590/S0100-879X2002000500014
    1. Restuccia D., Del Piero I., Martucci L., Zanini S. (2011). High-frequency oscillations after median-nerve stimulation do not undergo habituation: a new insight on their functional meaning? Clin. Neurophysiol. 122, 148–152 10.1016/j.clinph.2010.06.008
    1. Robinson S. E. (2004). Localization of event-related activity by SAM(erf). Neurol. Clin. Neurophysiol. 2004, 109
    1. Robinson S. E., Nagarajan S. S., Mantle M., Gibbons V., Kirsch H. (2004). Localization of interictal spikes using SAM(g2) and dipole fit. Neurol. Clin. Neurophysiol. 2004, 74
    1. Sarvas J. (1987). Basic mathematical and electromagnetic concepts of the biomagnetic inverse problem. Phys. Med. Biol. 32, 11–22 10.1088/0031-9155/32/1/004
    1. Srejic L. R., Valiante T. A., Aarts M. M., Hutchison W. D. (2013). High-frequency cortical activity associated with postischemic epileptiform discharges in an in vivo rat focal stroke model. J. Neurosurg. 118, 1098–1106 10.3171/2013.1.JNS121059
    1. Staba R. J., Bragin A. (2011). High-frequency oscillations and other electrophysiological biomarkers of epilepsy: underlying mechanisms. Biomark. Med. 5, 545–556 10.2217/bmm.11.72
    1. Stacey W. C., Kellis S., Greger B., Butson C. R., Patel P. R., Assaf T., et al. (2013). Potential for unreliable interpretation of EEG recorded with microelectrodes. Epilepsia 54, 1391–1401 10.1111/epi.12202
    1. Tadel F., Baillet S., Mosher J. C., Pantazis D., Leahy R. M. (2011). Brainstorm: a user-friendly application for MEG/EEG analysis. Comput. Intell. Neurosci. 2011:879716 10.1155/2011/879716
    1. Tort A. B., Scheffer-Teixeira R., Souza B. C., Draguhn A., Brankack J. (2013). Theta-associated high-frequency oscillations (110–160 Hz) in the hippocampus and neocortex. Prog. Neurobiol. 100, 1–14 10.1016/j.pneurobio.2012.09.002
    1. Uhlhaas P. J., Pipa G., Neuenschwander S., Wibral M., Singer W. (2011). A new look at gamma? High- (>60 Hz) gamma-band activity in cortical networks: function, mechanisms and impairment. Prog. Biophys. Mol. Biol. 105, 14–28 10.1016/j.pbiomolbio.2010.10.004
    1. Uhlhaas P. J., Singer W. (2013). High-frequency oscillations and the neurobiology of schizophrenia. Dialogues Clin. Neurosci. 15, 301–313
    1. Van Essen D. C., Ugurbil K., Auerbach E., Barch D., Behrens T. E., Bucholz R., et al. (2012). The Human Connectome Project: a data acquisition perspective. Neuroimage 62, 2222–2231 10.1016/j.neuroimage.2012.02.018
    1. Vrba J., Robinson S. E. (2001). Signal processing in magnetoencephalography. Methods 25, 249–271 10.1006/meth.2001.1238
    1. Weiss S. A., Banks G. P., McKhann G. M., Jr., Goodman R. R., Emerson R. G., Trevelyan A. J., et al. (2013). Ictal high frequency oscillations distinguish two types of seizure territories in humans. Brain 136, 3796–3808 10.1093/brain/awt276
    1. Worrell G. A., Jerbi K., Kobayashi K., Lina J. M., Zelmann R., Le Van Quyen M. (2012). Recording and analysis techniques for high-frequency oscillations. Prog. Neurobiol. 98, 265–278 10.1016/j.pneurobio.2012.02.006
    1. Xiang J., Degrauw X., Korostenskaja M., Korman A. M., O'Brien H. L., Kabbouche M. A., et al. (2013). Altered cortical activation in adolescents with acute migraine: a magnetoencephalography study. J. Pain 14, 1553–1563 10.1016/j.jpain.2013.04.009
    1. Xiang J., Holowka S., Qiao H., Sun B., Xiao Z., Jiang Y., et al. (2004). Automatic localization of epileptic zones using magnetoencephalography. Neurol. Clin. Neurophysiol. 2004, 98
    1. Xiang J., Liu Y., Wang Y., Kirtman E. G., Kotecha R., Chen Y., et al. (2009a). Frequency and spatial characteristics of high-frequency neuromagnetic signals in childhood epilepsy. Epileptic Disord. 11, 113–125 10.1684/epd.2009.0253
    1. Xiang J., Liu Y., Wang Y., Kotecha R., Kirtman E. G., Chen Y., et al. (2009b). Neuromagnetic correlates of developmental changes in endogenous high-frequency brain oscillations in children: a wavelet-based beamformer study. Brain Res. 1274, 28–39 10.1016/j.brainres.2009.03.068
    1. Xiang J., Wang Y., Chen Y., Liu Y., Kotecha R., Huo X., et al. (2010). Noninvasive localization of epileptogenic zones with ictal high-frequency neuromagnetic signals. J. Neurosurg. Pediatr. 5, 113–122 10.3171/2009.8.PEDS09345
    1. Xiang J., Xiao Z. (2009). Spatiotemporal and frequency signatures of noun and verb processing: a wavelet-based beamformer study. J. Clin. Exp. Neuropsychol. 31, 648–657 10.1080/13803390802448651
    1. Zafeiriou D. I., Vargiami E. (2012). Noninvasive ultra high-frequency (1kHz) oscillations' recording: high-fidelity over somatosensory cortex. Clin. Neurophysiol. 123, 2323–2324 10.1016/j.clinph.2012.05.010
    1. Zijlmans M., Huiskamp G. M., Cremer O. L., Ferrier C. H., Van Huffelen A. C., Leijten F. S. (2012a). Epileptic high-frequency oscillations in intraoperative electrocorticography: the effect of propofol. Epilepsia 53, 1799–1809 10.1111/j.1528-1167.2012.03650.x
    1. Zijlmans M., Jiruska P., Zelmann R., Leijten F. S., Jefferys J. G., Gotman J. (2012b). High-frequency oscillations as a new biomarker in epilepsy. Ann. Neurol. 71, 169–178 10.1002/ana.22548

Source: PubMed

3
Subscribe