US6832110B2 - Method for analysis of ongoing and evoked neuro-electrical activity - Google Patents
Method for analysis of ongoing and evoked neuro-electrical activity Download PDFInfo
- Publication number
- US6832110B2 US6832110B2 US10/235,163 US23516302A US6832110B2 US 6832110 B2 US6832110 B2 US 6832110B2 US 23516302 A US23516302 A US 23516302A US 6832110 B2 US6832110 B2 US 6832110B2
- Authority
- US
- United States
- Prior art keywords
- deflections
- stimuli
- neuro
- time
- subject
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related, expires
Links
- 238000000034 method Methods 0.000 title claims abstract description 68
- 230000000763 evoking effect Effects 0.000 title claims abstract description 64
- 238000004458 analytical method Methods 0.000 title claims abstract description 33
- 230000000694 effects Effects 0.000 title claims description 82
- 238000009826 distribution Methods 0.000 claims abstract description 72
- 230000004044 response Effects 0.000 claims abstract description 51
- 230000000008 neuroelectric effect Effects 0.000 claims abstract description 41
- 230000004069 differentiation Effects 0.000 claims abstract description 12
- 230000001953 sensory effect Effects 0.000 claims abstract description 11
- 210000004556 brain Anatomy 0.000 claims abstract description 10
- 238000005259 measurement Methods 0.000 claims abstract description 10
- 230000000007 visual effect Effects 0.000 claims description 11
- 208000027418 Wounds and injury Diseases 0.000 claims description 4
- 230000006378 damage Effects 0.000 claims description 4
- 201000010099 disease Diseases 0.000 claims description 4
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims description 4
- 208000014674 injury Diseases 0.000 claims description 4
- 239000003814 drug Substances 0.000 claims description 2
- 230000008569 process Effects 0.000 claims description 2
- 238000011084 recovery Methods 0.000 claims description 2
- 230000000246 remedial effect Effects 0.000 claims description 2
- 238000011002 quantification Methods 0.000 claims 1
- 230000001054 cortical effect Effects 0.000 description 10
- 238000012935 Averaging Methods 0.000 description 5
- 230000001537 neural effect Effects 0.000 description 4
- 230000003321 amplification Effects 0.000 description 3
- 238000013459 approach Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000007689 inspection Methods 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 230000003340 mental effect Effects 0.000 description 3
- 238000003199 nucleic acid amplification method Methods 0.000 description 3
- 210000004761 scalp Anatomy 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000004424 eye movement Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 238000012109 statistical procedure Methods 0.000 description 2
- 230000000638 stimulation Effects 0.000 description 2
- 206010002091 Anaesthesia Diseases 0.000 description 1
- 230000036626 alertness Effects 0.000 description 1
- 230000037005 anaesthesia Effects 0.000 description 1
- 230000003925 brain function Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 230000036648 cognitive speed Effects 0.000 description 1
- 230000036992 cognitive tasks Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 210000001061 forehead Anatomy 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000001339 gustatory effect Effects 0.000 description 1
- 210000001595 mastoid Anatomy 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 210000005195 occipital scalp Anatomy 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000001936 parietal effect Effects 0.000 description 1
- 210000005199 parietal scalp Anatomy 0.000 description 1
- 210000000977 primary visual cortex Anatomy 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000008521 reorganization Effects 0.000 description 1
- 230000000284 resting effect Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000010972 statistical evaluation Methods 0.000 description 1
- 230000002459 sustained effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/369—Electroencephalography [EEG]
- A61B5/377—Electroencephalography [EEG] using evoked responses
- A61B5/378—Visual stimuli
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/369—Electroencephalography [EEG]
- A61B5/377—Electroencephalography [EEG] using evoked responses
- A61B5/38—Acoustic or auditory stimuli
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/398—Electrooculography [EOG], e.g. detecting nystagmus; Electroretinography [ERG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
- A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
Definitions
- the present invention relates to the field of measurement and analysis of brain function and, more particularly, to a method of analysis of EEG and evoked potential (EP) generated data which allows differentiation between ongoing cortical neuro-electric activity and evoked cortical neuro-electric activity of a subject.
- EEG EEG and evoked potential
- the frequency analysis technique has also been suggested to study the evoked potentials in single trials.
- the evoked potentials can be considered to result from a reorganization of the phases of the ongoing EEG or from selective amplification of EEG in specific frequency bands.
- these spectral (frequency domain) methods do not relate to the local time properties of the responses (e.g. latencies of the particular EP components), and the methods do not allow study of the local time variations of the EP.
- the present invention describes a novel method for the analysis of both the neural electrical activity initiated in response to a series of sensory stimuli and of ongoing neural electrical activity.
- the positive and negative deflections in the EEG activity are detected and the time distributions and the amplitude as a function of time distributions of these deflections are obtained.
- This provides a statistical description of the appearance of positive and negative deflections or fluctuations in the ongoing EEG before, during and after a series of sensory stimuli.
- An algorithm based upon the statistical description is used to quantify the differentiation between the neural electrical activity evoked by the stimuli and the background EEG, both with respect to the number of these deflections and with respect to their amplitude. This differentiation provides insights into the possible relations between evoked and background EEG activity, leading to a better understanding of how evoked potentials are generated.
- a method of analyzing EEG and EP data in order to differentiate between ongoing neuro-electric activity of the brain of a subject and evoked neuro-electric activity of the brain of the subject, the analysis being based upon the frequency and the amplitude of deflections from the base line of the EEG and EP data comprising a. recording by EEG measurement the ongoing neuro-electric activity; b. administering a series of sensory stimuli to the subject; c. recording by EEG and EP measurement neuro-electric activity evoked in response to the stimuli; d. marking the time of occurrence, polarity and amplitude of each of the deflections recorded; e.
- FIG. 1 is a flow chart illustrating the method provided in accordance with the present invention.
- FIGS. 2 to 10 are graphic representations of data relevant to a described study of the implementation of the method of FIG. 1 .
- the present invention is of a method of EEG and EP analysis which preserves the moment-to-moment variation in scalp recorded activity, helps differentiate between evoked and ongoing electrical activity and provides a constructive basis for the study of the single evoked responses. Moreover, the variability of the responses can be evaluated by this method.
- FIG. 1 shows a flow chart illustrating a method of analysis of EEG and EP data according to the teachings of the present invention, which is referred to hereinbelow as method 10 .
- Step 1 consists of recording ongoing neuro-electric activity within the frequency range of interest of a subject by common equipment at the subject's scalp using a set of electroencephalograph electrodes and amplifier and analog/digital(A/D) converters to provide a set of digital data representing the subject's neuro-electric activity prior to the provision of stimuli.
- the duration of the measurement may vary according to the quantity of data and degree of precision of differentiation desired.
- Step 2 consists of administering a series of sensory stimuli to said subject.
- the stimuli may consist of any sensory stimulus, including but not limited to visual, auditory, tactile, gustatory or olfactory and/or may include performing an attention demanding task or engaging the memory process or other cognitive tasks. It is appreciated that the nature of the stimulus is not a material element of the invention and the method of the present invention may be employed with respect to any stimulus that evokes neuro-electrical activity.
- Step 3 consists of recording the neuro-electric activity evoked in response to the stimuli. The measurement is taken as described with respect to step 1 above.
- Steps 4 through 8 are best illustrated by example and are therefore described more fully hereinafter with respect to the study set forth hereinafter. Such steps include the following:
- Step 4 consists of marking the time of occurrence, polarity and amplitude of each of the deflections recorded by the EEG device.
- the deflections may be both positive and negative deflections from the baseline of the EEG data.
- Step 5 consists of generating a histogram of time periods of the deflections for a defined duration prior to, during and subsequent to each stimulus.
- the histogram will reflect the deflections within time periods ranging from 3 to 10 ms and preferably 6 ms. Said time periods may vary in duration as desired by the Examiner, depending on the frequency range of EEG to be analyzed.
- the histogram may span any time period desired by the examiner, for example a duration of 1300 ms, 300 ms prior to the stimulus and 1000 ms after the stimulus.
- Step 6 consists of generating a deflection time distribution based upon data from the histogram.
- Step 7 consists of generating an amplitude profile of the deflections.
- Step 8 consists of using the deflection time distribution and the amplitude profile along with statistical procedures to evaluate the ongoing and the evoked neuro-electric activity by assigning a value to each time period.
- the statistical procedures can be modified to take into account all frequency bands of recorded brain electrical activity. As described in the following study, synchronization of the alpha-beta frequency deflections with onset of the stimulus seems to be limited to latencies up to 200 ms., while later components (e.g. P300) might be generated by EEG waves lower in frequency content. It is understood that the algorithm may preferably be modifiable to allow analysis of evoked neuro-electric activity of different frequencies and amplitudes, including shorter and longer latency evoked responses and lower and higher frequency deflections.
- the analysis may be used for examining the timing of the neuro-electric activity, the variability of the neuro-electric activity and the correlation of the neuro-electric activity and therefore may also be helpful in identifying and/or analyzing the stimuli evoking the neuro-electric activity.
- the algorithm is preferably actuatable by a computer and, like all modern EEG measurement algorithms, is capable of being programmed.
- Method 10 may be used for all customary purposes carried out or facilitated by EEG and EP examination, including but not limited to measuring the effect of a stimulus on a subject compared to the effect of the stimulus on other subjects, measuring the effect of a stimulus on a subject over a period of time, measuring the effect of medicine or anesthesia on a response of a subject to a stimulus, measuring the effectiveness of a remedial program on a response of a subject to a stimulus, measuring the effect of fatigue on a response of a subject to a stimulus, measuring the effect of injury or disease on a response of a subject to a stimulus and measuring the effect of recovery from injury or disease on a response of a subject to a stimulus.
- Method 10 may also be useful for evaluating a subject's level of alertness, mental efforts and brain utilization, sustained focused attention, neurocognitive strategy, cognitive speed, transient focused attention, response to changes in mental workload and quickness to adapt.
- Auditory and visual stimuli were presented in alternating order at an overall rate of 1/sec.
- the duration of a stimulation/recording session was 600 sec (300 auditory and 300 visual stimuli: Aud-Vis-Aud-Vis . . . ).
- Six sessions on two experimental days were conducted with each subject. To insure subject attentiveness, the subject was instructed to perform mental counting of each perceived auditory stimulus.
- the auditory stimulus was a binaural click of 30 dB above the subject's subjective threshold (i.e. SL).
- the visual stimulus consisted of an alternating checkerboard with fill field of 11° (25.7 cm.) and the angle for check viewing of 1.35° (a rather low basic spatial frequency of visual stimuli). Subjects were sitting 140 cm. from a 27 ⁇ 21 cm. CTX computer screen, which provided the visual stimuli. In 3 of these subjects a session was conducted during which stimuli were not presented. This report describes the analysis of the responses to the visual stimuli only.
- Scalp activity was recorded with an Electrocap (Neuroscan, Inc., U.S.A.) using sixteen monopolar electrode sites: Fz, F 3 , F 4 , Cz, C 3 , C 4 , T 3 , T 4 , T 5 , T 6 , Pz, P 3 , P 4 , O 1 , O 2 , and Oz of the international 10-20 electrode placement system with the reference electrode at the mastoid, and the ground electrode was located on the forehead.
- Vertical and horizontal electrooculograms (EOGs) were recorded in order to control for eye movement artifacts.
- the stimulus presentations and data acquisition were provided by a STIM/SCAN EEG/EP workstation (Neuroscan, Inc., U.S.A.).
- smoothing of traces was conducted across every 5 data points. Eye movement artifact reduction was performed for all data using a NeuroScan procedure on the basis of the vertical electro-oculogram (VEOG).
- Computer algorithms were used to identify the positive and negative deflections and the time of their peaks (FIG. 2) during the 300 ms pre-stimulus and the 1,000 ms post-stimulus period. The base to peak amplitude of each such deflection was also measured (FIG. 3 ).
- the time coordinates of the deflections with respect to the onset of the stimulus trigger were determined separately for positive and negative deflections for each trial.
- each single trial could be described by a sequence of events (deflections), each of which had a time coordinate (with respect to the stimulus), polarity and amplitude.
- the number of positive and negative deflections in bins of 6 ms each calculated for a series of up to 300 trials was defined as “Deflection Time Distribution” for the given subject, given number of trials and given session (FIG. 2 ).
- the Deflection Time Distributions obtained over the six different recording sessions from the same subject could be combined (averaged or summated, wherever appropriate) into a composite Deflection Time Distributions. Deflection time distributions of the times of the positive and negative (separately) deflections were obtained for each subject.
- the amplitudes (base to peak) of the same deflections within each time bin as a function of post stimulus time could be used for the calculation of the Amplitude Profile for up to 300 trials over six sessions for each subject (FIG. 3 ). In general all 300 trials from a single session were combined into an Amplitude Profile.
- These deflection time distributions and amplitude profiles were the basis for further analysis which included, for example, statistical evaluation of the proportion of positive and negative deflections and their amplitudes in specific time periods with respect to the stimulus and analysis of the sequence of positive and then negative deflections and vice versa, in relation to the stimulus.
- the EEG recorded in three of the same subjects in the sessions without sensory stimuli was analyzed with the new technique in order to specify the time and amplitude distributions of true background activity. This would help in evaluating the time period which is truly “response ” in nature.
- the 300 ms of prestimulus time and 300 ms of late poststimulus time (interval 700-1000 ms) in the deflection time and amplitude distributions were compared. Similar distributions of uniform low level of number and amplitude of deflections were obtained in these two time intervals. In addition, similar uniform distributions were obtained in three subjects when the recordings were also conducted in the absence of any sensory stimuli. This helps define the 300 ms prestimulus period and the 300 ms (700-1000 ms) period post stimulus as the stationary (background) activity.
- a parameter R was defined as the normalized difference between the numbers of positive and negative deflections obtained in the particular time bin over the set of trials:
- #p number of positive deflections obtained in particular time bin,—obtained from Positive Deflection Time Distribution
- #n number of negative deflections obtained in particular time bin,—obtained from Negative Deflection Time Distribution
- the R value would be distributed approximately as a standard normal variance N(0,1) if the probability to obtain positive deflections in the given time bin over N trials is equal to the probability of obtaining negative deflections (i.e. background activity).
- FIG. 4 shows the R-values obtained during the resting EEG (absence of stimulus) and in the prestimulus intervals in subject SB, along with the distribution N(0,1) superimposed. This figure illustrates that the experimental data during background activity is in agreement with a standard normal variance distribution.
- Absolute values of R calculated from the recordings in subject SB are presented in FIG. 5 .
- the mean time interval between successive deflections and the reciprocal of this value is one of the principal characteristics of the background activity and can be estimated in the following way: For example, there were 2-6 positive deflections in each of the 300 ms prestimulus time intervals recorded in subject SB, with a mean across all 300 trials of 4.24 ⁇ 1.07 deflections. A mean number of 4.24 deflections in 300 ms gives a mean rate in subject SB of 14.13 ⁇ 3.56 Hz (range 6.67-20 Hz), and a mean time between successive deflections (period) of 70.8 ⁇ 24.1 ms.
- time interval histograms show that more than 95% of the positive EEG deflections detected by our algorithm appear with time intervals shorter than 110 ms between successive deflections, i.e. with a rate greater than 9 Hz, mostly in the Alpha and Beta range.
- the Fourier spectrum of the raw EEG segments from the same data had a much broader frequency content, including lower frequencies.
- the relatively high frequency content of the deflections (without low frequencies) used in this analysis is obviously the result of the low pass filtering at 30 Hz and the subject algorithm which detects all peaks and then takes the time interval between successive peaks from which to calculate rate. It therefore concentrates on the highest frequencies present in the raw EEG and does not take into account the broader (lower frequency) content.
- the R-values for the description of the background activity can also be used to define the time range of the onset, termination and duration of the visual evoked activity. These are presented in Table 1 below.
- the onset of the response at the Oz electrode was significantly earlier probably because this is the region of the primary visual cortex. Duration of the overall responses at these electrodes did not differ significantly.
- Table 1 presents data for the whole set of EP-components: N 1 -P 1 -N 2 .
- each of the components can be statistically defined as different peaks using the criteria R-value.
- positive and negative components can be presented separately, depending on the sign of the nominator used to calculate the R-value (#p ⁇ #n), (formula 1).
- FIG. 7 shows both the R-value definition of evoked activity based on the deflection time distributions and the conventionally averaged VEP from the same data.
- Inspection of the R-value definition in FIG. 7 reveals the principal components of VEP (N 1 , P 1 , N 2 ) and in addition the possible presence of evoked activity which is not apparent in the conventionally averaged VEP (positive components at 60 ms and 180 ms).
- the time segment (T 1 , T 2 ) has to include the time bin (T peak ) which corresponds to the latency P 1 VEP and contains a peak value of (positive) deflections;
- the border bins T 1 and T 2 may have R-values equal to zero.
- That part of the positive deflection time distribution related to the time bins [T 1 , T 2 ] was defined as the partial deflection time distribution corresponding to the VEP component P 1 .
- the partial deflection time distributions corresponding to VEP components N 1 , N 2 can be described in a similar way.
- the partial distributions corresponding to the VEP N 1 , P 1 , N 2 components were determined for the Oz electrode in all 14 subjects.
- An example of the partial distributions corresponding to the P 1 component obtained in the recordings in subject SB (one session, 300 trials; using the procedure presented in FIG. 8A is presented in FIG. 8 B.
- time interval [T 1 , T 2 ] where positive deflections related to the P 1 appear has a duration of 30 ms (5 time bins of 6 ms). Within this time range over 300 trials there were 204 positive deflections not uniformly distributed over these 5 bins, and the maximal number of deflections (54 out of 300 trials) was obtained in the 110-116 ms bin (i.e. T peak , the mode of the distribution, was 113 ms).
- the partial deflection time distributions corresponding to the EP components can be used as a measure of the variability of the evoked response across single trials, this is illustrated in FIG. 8 B.
- the mean latency and standard deviation (s) of these deflections was: 118.3 ⁇ 7.3 ms.
- the parameter Standard Deviation [ms] may serve as a measure of the “time jitter” of the given EP-component, i.e. trial to trial variability of latency.
- An additional measure of the variability of the latency of responses may be provided by the possible asymmetry of the partial time distribution of a wave, i.e. if the slopes on either side of the peak differ. In that case the modal latency and the mean latency will differ. For example, in FIG.
- Synchronization (time-locking) of deflections is reflected in the relative number of positive and negative deflections obtained in the time bins, and the R-value (formula 1) has been applied in order to separate between background and evoked activity.
- a high R-value is produced in a particular time bin when this bin contains more deflections with a specific polarity than some criterion. This means that a deflection of that polarity (positive or negative) was generated in the given time bin in many trials, while the opposite polarity deflections appeared only rarely.
- a significantly high R-value in the particular bin, or for the set of bins i.e.
- N 1 and P 1 components of the VEP were produced by the time locking of EEG deflections (synchronization in response to the successive stimuli) probably those in the high alpha—and beta range (7-30 Hz). These deflections can be considered as the single stimulus trial responses which contribute to the final averaged response.
- This suggestion for generation of the VEP components may not be appropriate for component N 2 because it seems that its duration (time interval of the corresponding partial deflection distribution) is greater than 30 ms and there may be more than one deflection in this time range.
- the data obtained over periods of background activity can be combined.
- the amplitudes of the deflections in the 300 ms prestimulus periods in 300 trials in one session for subject AG were calculated.
- FIG. 9 shows the distribution of their amplitudes.
- the two distributions are very similar as expected, since, by definition, during background activity the EEG is symmetrical both with respect to the rate of deflections (see R-value, described above) and their amplitude. These amplitude distributions could be approximated by a normal (Gaussian) variance symmetrically around the mean, with identical SDs for both positive and negative deflections (see FIG. 9 ).
- Such mean value of the amplitude of the “background” deflections presents a level around that the mean amplitudes calculated for the deflections which contributed to the given (successive) time bins (i.e. values of the Amplitude Profile), fluctuate (FIG. 3 for example).
- these fluctuations are much smaller than deviation of the amplitude of the single deflections, since the values of Amplitude Profiles are the product of averaging.
- the fluctuations of the amplitude of the positive EEG deflections calculated without averaging in time bins across trials were evaluated in subject SB over 6 sessions as 6.48 ⁇ 6.91 uV, and at the same time the values of the Amplitude Profiles (calculated over 6 sessions) fluctuated as 6.48 ⁇ 0.55 uV.
- the variation of the values of the Amplitude Profile for the time period of background activity in a given subject may be used to generate a statistical criteria for the determination of significant amplitude deviations which could serve as a sign of evoked activity (similar to the criteria used to differentiate between background and evoked activity based on the deflection time distributions.
- the critical value for positive deflections is 8.40 uV (6.48+3.5 ⁇ 0.55), and for negative deflections ⁇ 8.38 uV.
- the correlation coefficient ranged for N 1 -component from 0.879 (Oz) to 0.983 (P 4 ), and for P 1 -component from 0.690 (Oz) to 0.914 (P 4 ).
- the slope of regression line was maximal on the occipital electrodes (Oz,O 1 ,O 2 ): these occipital values ranged from 1.38 to 1.23 for N 1 , and from 1.77 to 1.49 for P 1 -component.
- the method of the present invention offers a number of advantages over prior methods of analysis of EEG generated data, as follows:
- the method of the present invention allows for detailed analysis of background and evoked EEG data within different frequency ranges, providing information about the stimulus evoking the response. Such analysis may be useful in the evaluation of cortical electrical activity in clinical cases.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Psychiatry (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Psychology (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Acoustics & Sound (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
A method of analyzing EEG and EP data in order to differentiate between ongoing neuro-electric activity of the brain of a subject and evoked neuro-electric activity of the brain of the subject, the analysis being based upon the frequency and the amplitude of deflections from the base line of the EEG and EP data, the method including a. recording by EEG measurement the ongoing neuro-electric activity; b. administering a series of sensory stimuli to the subject; c. recording by EEG and EP measurement neuro-electric activity evoked in response to the stimuli; d. marking the time of occurrence, polarity and amplitude of each of the deflections recorded; e. generating a histogram of time periods of the deflections for a defined duration prior to, during and subsequent to each of the single stimuli, the histogram reflecting the deflections within the time periods, generating a deflection time distribution based upon data from the histogram, generating an amplitude profile of the deflections, and f using the deflection time distribution and the amplitude profile to statistically analyze and evaluate the ongoing and evoked neuro-electric activity; such that the differentiation between the ongoing and the evoked neuro-electric activity may be quantified.
Description
The present application claims benefit of U.S. Provisional Application No. 60/316,974 filed Sep. 5, 2001.
The present invention relates to the field of measurement and analysis of brain function and, more particularly, to a method of analysis of EEG and evoked potential (EP) generated data which allows differentiation between ongoing cortical neuro-electric activity and evoked cortical neuro-electric activity of a subject.
The conventional analysis of EPs is based on averaged data due to the small signal (response) to noise (background EEG) ratio. Thus, data concerning individual trials and the variation of the evoked responses from trial to trial are lost and the nature of the individual response is unclear. For example, a typical goal in EP analysis in clinical cases is to elucidate the reason for the attenuation of particular EP components: is it due to variations in latency of the components across trials or to lower amplitudes of that component in several individual trials? Therefore, many studies have searched for a technique which could possibly provide single trial analysis of the evoked responses
An alternative approach is based on the assumption that the waveform of the conventionally averaged EP is similar to that of the single response so that the averaged response can be used as a criteria (template) for the selection of the appropriate trials for single trial analysis. Thus, only the components of EP which are obtained after averaging can then be used in this approach. However, it is known that some single trial components can be distorted, highly decreased in amplitude or even lost as a result of the averaging, leading to a loss of possibly important information needed for such a study of single trial responses.
The frequency analysis technique has also been suggested to study the evoked potentials in single trials. Thus the evoked potentials can be considered to result from a reorganization of the phases of the ongoing EEG or from selective amplification of EEG in specific frequency bands. However, these spectral (frequency domain) methods do not relate to the local time properties of the responses (e.g. latencies of the particular EP components), and the methods do not allow study of the local time variations of the EP.
Recent studies have shown that background activity and single trial responses are integrated with each other. Consequently it would be advantageous to use the same method for processing and analyzing both aspects (responses and background) of the recorded activity. However, the problem of differentiation of single trial evoked responses from the background EEG activity has yet to be solved, and reliable methods for the analysis of EP and background variability are needed.
There is thus a widely recognized need for a method of EEG and EP data analysis that is devoid of the above described limitations and that facilitates the differentiation of the ongoing cortical neuro-electric activity and the evoked cortical neuro-electric activity of the brain of a subject.
Accordingly, it is an object of the present invention to provide a method of EEG and EP data analysis that facilitates the differentiation of the ongoing cortical neuro-electric activity and the evoked cortical neuro-electric activity of the brain of a subject.
It is a further object of the present invention to provide a method that may identify responses to single stimuli evoking neuro-electrical activity.
It is a further object of the present invention to provide a method that may be actuated by an algorithm.
It is a further object of the present invention to provide a method that employs an algorithm that is computer programmable.
It is a further object of the present invention to provide a method that employs an algorithm that may be modified to provide for analysis of differing EEG and EP data.
The present invention describes a novel method for the analysis of both the neural electrical activity initiated in response to a series of sensory stimuli and of ongoing neural electrical activity. The positive and negative deflections in the EEG activity are detected and the time distributions and the amplitude as a function of time distributions of these deflections are obtained. This provides a statistical description of the appearance of positive and negative deflections or fluctuations in the ongoing EEG before, during and after a series of sensory stimuli. An algorithm based upon the statistical description is used to quantify the differentiation between the neural electrical activity evoked by the stimuli and the background EEG, both with respect to the number of these deflections and with respect to their amplitude. This differentiation provides insights into the possible relations between evoked and background EEG activity, leading to a better understanding of how evoked potentials are generated.
Accordingly, there is provided a method of analyzing EEG and EP data in order to differentiate between ongoing neuro-electric activity of the brain of a subject and evoked neuro-electric activity of the brain of the subject, the analysis being based upon the frequency and the amplitude of deflections from the base line of the EEG and EP data, the method comprising a. recording by EEG measurement the ongoing neuro-electric activity; b. administering a series of sensory stimuli to the subject; c. recording by EEG and EP measurement neuro-electric activity evoked in response to the stimuli; d. marking the time of occurrence, polarity and amplitude of each of the deflections recorded; e. generating a histogram of time periods of the deflections for a defined duration prior to, during and subsequent to each of the single stimuli, the histogram reflecting the deflections within the time periods, generating a deflection time distribution based upon data from the histogram, generating an amplitude profile of the deflections; and f using the deflection time distribution and the amplitude profile to statistically analyze and evaluate the ongoing and evoked neuro-electric activity; such that the differentiation between the ongoing and the evoked neuro-electric activity may be quantified.
With specific reference now to the drawings in detail, it is stressed that the particulars shown are by way of example and for the purposes of illustrative discussion of the preferred embodiment of the present invention only, and are presented in the cause of providing what is believed to be the most useful and readily understood description of the principles and conceptual aspects of the invention. In this regard, no attempt is made to show details of the invention in more detail than is necessary for a fundamental understanding of the invention, the description taken with the drawings making apparent to those skilled in the art how the several forms of the invention may be embodied in practice.
In the drawings:
FIG. 1 is a flow chart illustrating the method provided in accordance with the present invention; and
FIGS. 2 to 10 are graphic representations of data relevant to a described study of the implementation of the method of FIG. 1.
The present invention is of a method of EEG and EP analysis which preserves the moment-to-moment variation in scalp recorded activity, helps differentiate between evoked and ongoing electrical activity and provides a constructive basis for the study of the single evoked responses. Moreover, the variability of the responses can be evaluated by this method.
Reference is now made to FIG. 1 which shows a flow chart illustrating a method of analysis of EEG and EP data according to the teachings of the present invention, which is referred to hereinbelow as method 10.
Step 7 consists of generating an amplitude profile of the deflections.
The statistical procedures can be modified to take into account all frequency bands of recorded brain electrical activity. As described in the following study, synchronization of the alpha-beta frequency deflections with onset of the stimulus seems to be limited to latencies up to 200 ms., while later components (e.g. P300) might be generated by EEG waves lower in frequency content. It is understood that the algorithm may preferably be modifiable to allow analysis of evoked neuro-electric activity of different frequencies and amplitudes, including shorter and longer latency evoked responses and lower and higher frequency deflections.
The analysis may be used for examining the timing of the neuro-electric activity, the variability of the neuro-electric activity and the correlation of the neuro-electric activity and therefore may also be helpful in identifying and/or analyzing the stimuli evoking the neuro-electric activity. The algorithm is preferably actuatable by a computer and, like all modern EEG measurement algorithms, is capable of being programmed.
Additional objects, advantages, and novel features of the present invention will become apparent to one ordinarily skilled in the art upon examination of the following study, which is not intended to be limiting. Additionally, each of the various aspects of the present invention as delineated herein and/or as claimed in the claims section below finds experimental support in the following study.
Subjects
Fourteen male subjects (age range 18-30 years) were invited to participate in the study as paid participants. All subjects had normal vision (in some cases corrected with glasses) and normal hearing.
Stimuli
Auditory and visual stimuli were presented in alternating order at an overall rate of 1/sec. The duration of a stimulation/recording session was 600 sec (300 auditory and 300 visual stimuli: Aud-Vis-Aud-Vis . . . ). Six sessions on two experimental days (three sessions per day) were conducted with each subject. To insure subject attentiveness, the subject was instructed to perform mental counting of each perceived auditory stimulus. The auditory stimulus was a binaural click of 30 dB above the subject's subjective threshold (i.e. SL). The visual stimulus consisted of an alternating checkerboard with fill field of 11° (25.7 cm.) and the angle for check viewing of 1.35° (a rather low basic spatial frequency of visual stimuli). Subjects were sitting 140 cm. from a 27×21 cm. CTX computer screen, which provided the visual stimuli. In 3 of these subjects a session was conducted during which stimuli were not presented. This report describes the analysis of the responses to the visual stimuli only.
Apparatus and Procedure
Scalp activity was recorded with an Electrocap (Neuroscan, Inc., U.S.A.) using sixteen monopolar electrode sites: Fz, F3, F4, Cz, C3, C4, T3, T4, T5, T6, Pz, P3, P4, O1, O2, and Oz of the international 10-20 electrode placement system with the reference electrode at the mastoid, and the ground electrode was located on the forehead. Vertical and horizontal electrooculograms (EOGs) were recorded in order to control for eye movement artifacts. The stimulus presentations and data acquisition were provided by a STIM/SCAN EEG/EP workstation (Neuroscan, Inc., U.S.A.). EEG and EOG channels were continuously recorded through Synamps Amplifiers (bandpass=0.1 to 100 Hz) and digitized at a rate of 1000 Hz. Data were stored and analyzed off line.
Data Processing
The 1,300 ms EEG segment recorded, beginning at 300 ms before the visual stimulus and continuing 1000 ms after stimulus presentation was defined as a single trial. The frequency content of the recorded activity was limited by the bandpass filter used (1-30 Hz). In addition, smoothing of traces was conducted across every 5 data points. Eye movement artifact reduction was performed for all data using a NeuroScan procedure on the basis of the vertical electro-oculogram (VEOG). Computer algorithms were used to identify the positive and negative deflections and the time of their peaks (FIG. 2) during the 300 ms pre-stimulus and the 1,000 ms post-stimulus period. The base to peak amplitude of each such deflection was also measured (FIG. 3).
The time coordinates of the deflections with respect to the onset of the stimulus trigger were determined separately for positive and negative deflections for each trial. Thus each single trial could be described by a sequence of events (deflections), each of which had a time coordinate (with respect to the stimulus), polarity and amplitude. The number of positive and negative deflections in bins of 6 ms each calculated for a series of up to 300 trials was defined as “Deflection Time Distribution” for the given subject, given number of trials and given session (FIG. 2). The Deflection Time Distributions obtained over the six different recording sessions from the same subject could be combined (averaged or summated, wherever appropriate) into a composite Deflection Time Distributions. Deflection time distributions of the times of the positive and negative (separately) deflections were obtained for each subject.
In addition, the amplitudes (base to peak) of the same deflections within each time bin as a function of post stimulus time could be used for the calculation of the Amplitude Profile for up to 300 trials over six sessions for each subject (FIG. 3). In general all 300 trials from a single session were combined into an Amplitude Profile. These deflection time distributions and amplitude profiles were the basis for further analysis which included, for example, statistical evaluation of the proportion of positive and negative deflections and their amplitudes in specific time periods with respect to the stimulus and analysis of the sequence of positive and then negative deflections and vice versa, in relation to the stimulus.
These same data were also subjected to conventional averaging and the deflection time distributions and the Amplitude Profiles obtained with this new technique were compared to the conventionally averaged EPs obtained from the same subjects, from the same recording sites over the same time period and for the same components. The latency of conventionally averaged VEP P1 components was identified as the most positive point between 75 and 135 ms, N1 as the most negative point between 55 and 95 ms, and N2 as the most negative point between 130 and 206 ms.
The EEG recorded in three of the same subjects in the sessions without sensory stimuli was analyzed with the new technique in order to specify the time and amplitude distributions of true background activity. This would help in evaluating the time period which is truly “response ” in nature.
Results
1. Differentiation Between Background and Evoked Activity on the Basis of the Deflection Time Distributions
Inspection of the time and amplitude distributions of the deflections as a function of time in a typical subject over all 6 recording sessions (FIG. 2) reveals that during the 300 ms prestimulus period and from about 400 ms after the stimulus (time 0), there are no major peaks. This indicates that during these periods, the number of deflections is at a uniform low level and that their amplitude is relatively constant, reflecting background activity. On the other hand, inspection of time and amplitude distribution of the deflections during the period from the stimulus up to about 400 ms after the stimulus (FIGS. 2, 3) reveals several time periods with clear increases in the number of positive and negative deflections and peaks in the amplitude profiles. This obviously represents the stimulus evoked activity reflected in these distributions.
In order to quantitatively determine the time periods during which only background (not time-locked) activity was present, the 300 ms of prestimulus time and 300 ms of late poststimulus time (interval 700-1000 ms) in the deflection time and amplitude distributions were compared. Similar distributions of uniform low level of number and amplitude of deflections were obtained in these two time intervals. In addition, similar uniform distributions were obtained in three subjects when the recordings were also conducted in the absence of any sensory stimuli. This helps define the 300 ms prestimulus period and the 300 ms (700-1000 ms) period post stimulus as the stationary (background) activity.
In order to define a statistical criteria for the presence of an evoked response, one must first define the variance of the background activity for the time period described above. An evoked response will then be indicated by a significant increase in this variance above some critical level.
One may assume that the difference between the number of positive and negative deflections in the time bins representing the background activity across trials would be small and differ only by chance from zero, i.e. the probabilities of detecting positive or negative deflections in the particular time bins (across trials) would be equal. In order to obtain a statistical measure of this variance between the number of positive and negative deflections which would help differentiate between background and evoked activity, a parameter R was defined as the normalized difference between the numbers of positive and negative deflections obtained in the particular time bin over the set of trials:
where
#p—number of positive deflections obtained in particular time bin,—obtained from Positive Deflection Time Distribution;
#n—number of negative deflections obtained in particular time bin,—obtained from Negative Deflection Time Distribution;
SQRT—square root.
The R value would be distributed approximately as a standard normal variance N(0,1) if the probability to obtain positive deflections in the given time bin over N trials is equal to the probability of obtaining negative deflections (i.e. background activity).
FIG. 4 shows the R-values obtained during the resting EEG (absence of stimulus) and in the prestimulus intervals in subject SB, along with the distribution N(0,1) superimposed. This figure illustrates that the experimental data during background activity is in agreement with a standard normal variance distribution.
Entering the values of R obtained during background activity into the table of normal (Gauss) distribution N(0,1) provides the criterion for assessing the presence of an R value significantly greater than that during background activity (which would then be evoked activity). However, one has to keep in mind that the critical level Rcr(p<0.05) is dependent on the duration of the time interval of the analysis, since higher values of R are more probable during longer samples of background activity. Thus a stronger criteria for the detection of a significant deviation of R must be applied during longer time intervals. For example, a duration of 1300 ms (216 time bins of 6 ms) requires a critical level Rcr (p<0.05)=3.5. However, for 100 ms poststimulus time (only 16 time bins of 6 ms), the critical level Rcr (p<0.05) would be 2.73.
Absolute values of R calculated from the recordings in subject SB are presented in FIG. 5. One can see a significant increase in the R-value in the specific time range (e.g. from 50 to about 300 ms) where the conventional evoked potential is usually obtained. The calculated critical level Rcr=3.5 is indicated in FIG. 5 as the line parallel to the time axis and cutting the profile of the R-values. The time interval during which the R-values project above the critical level of 3.5 and then below it, can be used to define the period of appearance of evoked activity.
2. Detailed Analysis of Background Activity on the Basis of the Deflection Time Distributions
Time Intervals Between Successive Deflections and Rate of Deflections
The mean time interval between successive deflections and the reciprocal of this value (termed here rate, or frequency of deflections) is one of the principal characteristics of the background activity and can be estimated in the following way: For example, there were 2-6 positive deflections in each of the 300 ms prestimulus time intervals recorded in subject SB, with a mean across all 300 trials of 4.24±1.07 deflections. A mean number of 4.24 deflections in 300 ms gives a mean rate in subject SB of 14.13±3.56 Hz (range 6.67-20 Hz), and a mean time between successive deflections (period) of 70.8±24.1 ms.
A similar operation was performed on the 300 ms prestimulus periods from all 6 sessions recorded in all 14 subjects. The mean rate across subjects was 15.25±0.78 Hz. Thus the rate of deflections detected in the 300 ms prestimulus period (background activity) in the 14 subjects lies within the beta range of EEG activity (for the Oz electrode).
While the mean rate of deflections calculated from the deflection time distribution algorithm is in the beta range, the frequency content of the raw EEG was obviously broader, also containing oscillations lower in frequency. The histograms of the time intervals between the successive positive deflections calculated by our algorithm in the 300 ms prestimulus EEG in 2 subjects was constructed (data over 300 EEG segments, FIG. 6).
The time interval histograms (FIG. 6) show that more than 95% of the positive EEG deflections detected by our algorithm appear with time intervals shorter than 110 ms between successive deflections, i.e. with a rate greater than 9 Hz, mostly in the Alpha and Beta range. The Fourier spectrum of the raw EEG segments from the same data had a much broader frequency content, including lower frequencies.
The relatively high frequency content of the deflections (without low frequencies) used in this analysis is obviously the result of the low pass filtering at 30 Hz and the subject algorithm which detects all peaks and then takes the time interval between successive peaks from which to calculate rate. It therefore concentrates on the highest frequencies present in the raw EEG and does not take into account the broader (lower frequency) content.
a. Characterization of the Overall Evoked Response
The R-values for the description of the background activity (see formula 1) can also be used to define the time range of the onset, termination and duration of the visual evoked activity. These are presented in Table 1 below. The onset of the response at the Oz electrode was significantly earlier probably because this is the region of the primary visual cortex. Duration of the overall responses at these electrodes did not differ significantly.
Table 1 presents data for the whole set of EP-components: N1-P1-N2. However, one can consider each of the components separately since they can be statistically defined as different peaks using the criteria R-value. Moreover, positive and negative components can be presented separately, depending on the sign of the nominator used to calculate the R-value (#p−#n), (formula 1).
TABLE 1 |
Mean and standard deviations values of the onset, termination, |
and duration of the overall activity evoked by the visual stimulation |
in 14 subjects. Time duration of analysis was 1300 ms with a critical |
Rcr value = 3.5. Data from four electrodes. |
Oz | Pz | Cz | Fz | ||
Onset: | 60.0 ± 7.9 | 77.5 ± 22.0 | 92.0 ± 16.6 | 87.5 ± 14.9 |
[ms] | ||||
Termination: | 214.5 ± 23.4 | 248.0 ± 25.1 | 227.3 ± 20.5 | 234.5 ± 13.7 |
[ms] | ||||
Duration: | 154.5 ± 26.7 | 170.5 ± 38.0 | 135.3 ± 27.0 | 147.0 ± 25.5 |
[ms] | ||||
FIG. 7 shows both the R-value definition of evoked activity based on the deflection time distributions and the conventionally averaged VEP from the same data. One can see good agreement between both curves: the principal VEP-components are present in the time regions where the R-value projects above the critical level Rcr=3.5. Inspection of the R-value definition in FIG. 7 reveals the principal components of VEP (N1, P1, N2) and in addition the possible presence of evoked activity which is not apparent in the conventionally averaged VEP (positive components at 60 ms and 180 ms).
b. Breakdown of the Particular EP-components from the Peaks of Deflection Time Distributions
The correspondence between the known components of conventional EP, the deflection time distributions and R-values permits the separation of the overall deflection time distribution of evoked activity into its local distribution for individual EP components. The separation of the positive deflections associated with VEP P1 component from the deflection time distribution obtained in subject SB is illustrated in the FIG. 8A. The goal of this procedure is the determination of the left (T1) and right (T2) borders of the time segment containing time bins with a high percentage of the positive deflections which can be related to VEP P1. The following conditions have to be fulfilled by this procedure:
1) the time segment (T1, T2) has to include the time bin (Tpeak) which corresponds to the latency P1 VEP and contains a peak value of (positive) deflections;
2) the R-value calculated for this time bin (Tpeak) has to exceed the critical level Rcr;
3) the R-values of the other time bins neighboring the Tpeak bin and within this segment (T1, T2) have to have the same sign as the R-value of the Tpeak bin (i.e. to be greater a zero for the distribution of positive deflections);
4) the border bins T1 and T2 may have R-values equal to zero.
That part of the positive deflection time distribution related to the time bins [T1, T2] was defined as the partial deflection time distribution corresponding to the VEP component P1. The partial deflection time distributions corresponding to VEP components N1, N2 can be described in a similar way. The partial distributions corresponding to the VEP N1, P1, N2 components were determined for the Oz electrode in all 14 subjects. An example of the partial distributions corresponding to the P1 component obtained in the recordings in subject SB (one session, 300 trials; using the procedure presented in FIG. 8A is presented in FIG. 8B. One can see that the time interval [T1, T2] where positive deflections related to the P1 appear, has a duration of 30 ms (5 time bins of 6 ms). Within this time range over 300 trials there were 204 positive deflections not uniformly distributed over these 5 bins, and the maximal number of deflections (54 out of 300 trials) was obtained in the 110-116 ms bin (i.e. Tpeak, the mode of the distribution, was 113 ms).
c. Analysis of the EP-components Variability on the Basis of the Partial Deflection Time Distributions
The partial deflection time distributions corresponding to the EP components can be used as a measure of the variability of the evoked response across single trials, this is illustrated in FIG. 8B. The mean latency and standard deviation (s) of these deflections was: 118.3±7.3 ms. The parameter Standard Deviation [ms] may serve as a measure of the “time jitter” of the given EP-component, i.e. trial to trial variability of latency. An additional measure of the variability of the latency of responses may be provided by the possible asymmetry of the partial time distribution of a wave, i.e. if the slopes on either side of the peak differ. In that case the modal latency and the mean latency will differ. For example, in FIG. 8B (subject SB) the time coordinate of the peak number of deflections (mode) Tpeak is 113 ms and mean value of the latency of the deflections Tmean is 118.3 ms. This gives an asymmetry of the partial distribution (as): as=Tpeak−Tmean=113−118.3=−5.3 ms.
The parameters of the partial time deflection distributions, calculated for VEP-components N1, P1, N2 over 14 subjects by methods described above, are presented hereinafter in Table 2. The calculation of the T1 and T2 for VEP components N1 and P1 (electrode Oz) in all 14 subjects revealed that as a rule, the duration of time interval (T1, T2) for these components is less then 30 ms. Since the time interval between successive deflections of a specific polarity is significantly longer then 30 ms (for example the mean time interval evaluated in background activity over all 14 subjects was 65.6 ms), it is not likely that there can be more then one deflection contribution to the time range of the partial deflection time distribution of N1 or P1 in response to a single stimulus (single trial). In other words, there seems to be a one-to-one correspondence between the deflections which make up the partial time distribution and the trials related to this distribution. The number of trials whose deflections contribute to this distribution is less than the total number of trials used in the analysis (e.g. in FIG. 8B there were only 204 deflections in the P1 time distribution across 300 trials). This means that in some stimulus trials a positive deflection was absent in the P1 time range. The mean percentage of trials containing the “appropriate deflection” can be evaluated from Table 2. For example, a deflection in the P1 component time range was present in our subjects in approximately 224.7 trials out of 300, i.e. in 75% of the trials. These trials containing “appropriate” deflections can be considered as an approach to single trial analysis.
TABLE 2 |
Parameters of VEP-components in Partial deflection time distributions: |
N1 | P1 | N2 | ||
Peak Value of the | 18.9 ± 4.4 | 22.1 ± 4.9 | 15.3 ± 3.3 |
partial distribution: [%] | |||
# deflection/300 | |||
Latency of the peak | 88.8 ± 8.7 | 121.0 ± 4.4 | 157.0 ± 11.0 |
number of deflections: | |||
Tpeak [ms] | |||
T1 [ms]: | 80.4 ± 8.0 | 108.3 ± 5.3 | 140.9 ± 9.3 |
Left border of the | |||
partial distribution | |||
T2 [ms]: | 101.4 ± 6.7 | 132.7 ± 4.6 | 180.3 ± 17.6 |
Right border of the | |||
partial distribution | |||
Width (duration) of the | 21.0 ± 3.0 | 24.4 ± 4.8 | 39.4 ± 13.6 |
partial distribution | |||
T2 − T1 [ms] | |||
Mean Latency of | 90.2 ± 7.1 | 120.2 ± 4.6 | 156.5 ± 12.8 |
deflections contributed | |||
to the partial distribution | |||
Tmean [ms] | |||
Standard Deviation of | 6.9 ± 0.8 | 7.5 ± 1.3 | 9.4 ± 2.6 |
the latency of | |||
deflections contributed | |||
to the partial distribution: | |||
s [ms] | |||
Asymmetry of the | −1.4 ± 2.8 | 0.8 ± 3.2 | 0.5 ± 14.7 |
partial distribution [ms] | |||
Tmode − Tmean | |||
Mean # Deflections | 186.3 ± 30.2 | 224.7 ± 35.4 | 236.5 ± 70.2 |
contributed to the | |||
partial distribution | |||
(from 300 trials) | |||
Mean ± SD in 14 subjects. Electrode Oz. |
d. Application of the Deflection Time Distributions (R-criteria) for Response Determination
The clear correspondence between the time distributions of deflections detected by the algorithm used here and the components of the conventional averaged EP with respect to N1 and P1 leads to the suggestion that these principal components of the VEP are generated by the EEG deflections themselves. During background EEG, the deflections are dispersed randomly in time, at a rate of about 15 Hz, as shown. Upon repeated stimulus presentations (trials) there is a tendency for the deflections to appear within narrow time periods after the stimulus. That is, when examining these periods following repeated stimulus trials, there is a “synchronization” of the deflections over many, but not all, stimulus trials.
Synchronization (time-locking) of deflections is reflected in the relative number of positive and negative deflections obtained in the time bins, and the R-value (formula 1) has been applied in order to separate between background and evoked activity. A high R-value is produced in a particular time bin when this bin contains more deflections with a specific polarity than some criterion. This means that a deflection of that polarity (positive or negative) was generated in the given time bin in many trials, while the opposite polarity deflections appeared only rarely. Thus a significantly high R-value in the particular bin, or for the set of bins (i.e. for the some time interval),—is an indication of synchronization (time-locking) of the deflections of that polarity in response to the stimulus. This measure of synchronization (high magnitude of R) is exclusively dependent on the presence of the deflections and not on their amplitude. Therefore the correspondence between VEP parameters and R-value provide evidence for the important role of synchronization of deflections in producing the conventionally averaged EP.
Thus it is proposed that the N1 and P1 components of the VEP were produced by the time locking of EEG deflections (synchronization in response to the successive stimuli) probably those in the high alpha—and beta range (7-30 Hz). These deflections can be considered as the single stimulus trial responses which contribute to the final averaged response. This suggestion for generation of the VEP components may not be appropriate for component N2 because it seems that its duration (time interval of the corresponding partial deflection distribution) is greater than 30 ms and there may be more than one deflection in this time range.
3. Amplitude of Deflections and the Amplitude Profile in Background Activity
Due to the stationarity of background activity, the data obtained over periods of background activity can be combined. In order to illustrate the amplitudes and. variation of the amplitude of the EEG deflections measured during background activity, the amplitudes of the deflections in the 300 ms prestimulus periods in 300 trials in one session for subject AG were calculated. The mean (±SD) amplitude of the positive deflections was 6.06±7.28 uV (n=1483 deflections) and for the negative deflections was −6.05±7.24 uV (n=1482). FIG. 9 shows the distribution of their amplitudes. The two distributions (positive and negative deflections) are very similar as expected, since, by definition, during background activity the EEG is symmetrical both with respect to the rate of deflections (see R-value, described above) and their amplitude. These amplitude distributions could be approximated by a normal (Gaussian) variance symmetrically around the mean, with identical SDs for both positive and negative deflections (see FIG. 9). The distributions of the amplitude of deflections during background activity differed in different subjects: the means ranged from 2 to 8 uV over the 14 subjects, with a mean amplitude over all subjects of 4.73±1.80 uV for the positive deflections, and −5.03±1.70 uV for negative deflections (n=14).
Such mean value of the amplitude of the “background” deflections (positive or negative) presents a level around that the mean amplitudes calculated for the deflections which contributed to the given (successive) time bins (i.e. values of the Amplitude Profile), fluctuate (FIG. 3 for example). One can note that these fluctuations are much smaller than deviation of the amplitude of the single deflections, since the values of Amplitude Profiles are the product of averaging. For example, the fluctuations of the amplitude of the positive EEG deflections calculated without averaging in time bins across trials were evaluated in subject SB over 6 sessions as 6.48±6.91 uV, and at the same time the values of the Amplitude Profiles (calculated over 6 sessions) fluctuated as 6.48±0.55 uV.
4. Differentiation Between Background and Evoked Activity on the Basis of the Peak Amplitude of Deflections
The variation of the values of the Amplitude Profile for the time period of background activity in a given subject may be used to generate a statistical criteria for the determination of significant amplitude deviations which could serve as a sign of evoked activity (similar to the criteria used to differentiate between background and evoked activity based on the deflection time distributions.
For example, in subject SB the mean and standard deviation of the amplitude of the deflections during background activity over 6 sessions was 6.48±0.55 uV for positive deflections and −6.52±0.53 uV for negative. By how much should amplitudes have to differ from this mean in order to be considered to have surpassed it and therefore defined as evoked activity; i.e. what is the criterion magnitude of this critical deviation?
As noted above a significance of p<0.05 may be achieved in our case if one uses a criteria of ±3.5 SDs from the mean value obtained during background activity. Thus, with respect to subject SB, the critical value for positive deflections is 8.40 uV (6.48+3.5×0.55), and for negative deflections −8.38 uV. (In addition the criteria for determination of amplitudes significantly lower than background can also be calculated: for positive deflections=4.56 uV, for negative=−4.66 uV). These criteria for the positive and negative deflections are indicated in FIG. 10 as the shaded area. The Amplitude Profiles reach and project above these critical levels during the time regions corresponding to the evoked components, and this seems to be clearer for the positive waves.
This leads to the suggestion that there are two different mechanisms of response generation acting in different subjects in different proportions: first, one responsible for the occurrence of the wave (the peak of neuronal activity, EEG-deflection) at a certain moment after the onset of the stimulus (latency), and the second one, responsible for its magnitude. Thus, appearance of a peak (component) in the conventional averaged EP in response to repeated stimuli is probably determined by two factors: (1) the synchronization of the deflections (the time locking of deflections) and (2) an increase in the amplitude of these deflections, (i.e. the amplitude of the deflections during the specific period is greater than that during background activity). The second one can be called “amplification” of deflections and it is really due to the enhanced synchronization of synaptic-dendritic potentials initiated by stimuli.
5. Correlation Between the Mean Amplitude of Deflections Corresponding to the EP Components and the EEG Deflections in Background Activity
There was a strong positive correlation between the peak value of the Amplitude Profiles corresponding to the N1 and P1 components and the mean amplitude of deflections in the background activity at the parietal and occipital scalp electrodes (Pz,P3,P4, Oz,O1,O2,T5,T6) where response (conventional VEP) was maximal, and at the central electrodes (Cz,C3,C4). This, it seems, that subjects with larger amplitudes of deflections during background activity also had larger amplitudes of their N1 and P1 VEPs. The parameters of this dependence (correlation coefficient, slope of regression line and intercept) differed slightly between N1 and P1 components and were different on electrode. The correlation coefficient ranged for N1-component from 0.879 (Oz) to 0.983 (P4), and for P1-component from 0.690 (Oz) to 0.914 (P4). The slope of regression line was maximal on the occipital electrodes (Oz,O1,O2): these occipital values ranged from 1.38 to 1.23 for N1, and from 1.77 to 1.49 for P1-component.
The results of this part of the study on the interrelations between evoked and background activity leads to the suggestion that different components of the EP might be differently related to and originating from the background activity.
The method of the present invention offers a number of advantages over prior methods of analysis of EEG generated data, as follows:
1. It provides a form of single trial analysis with a greater degree of distinction between ongoing cortical electrical activity and evoked cortical electrical activity.
2. It allows individual analysis of both the timing and the amplitude of deflections, the two mechanisms involved in generation of the conventionally averaged EEG data for evoked responses.
3. It provides insight into possible mechanisms involved in the generation of evoked cortical electrical activity based on the time-locking (EEG-synchronization) and amplitude amplification of EEG deflections.
4. It allows detailed evaluation of the temporal properties of EP activity using methods similar to the study of trains of spikes (for example, analogous to post-stimulus time histograms—deflection time distributions and time interval histograms which were used in the present study).
5. It allows study of correlation between the time structure of the occurrence of deflections and their amplitude.
In summary, the method of the present invention allows for detailed analysis of background and evoked EEG data within different frequency ranges, providing information about the stimulus evoking the response. Such analysis may be useful in the evaluation of cortical electrical activity in clinical cases.
It is appreciated that certain features of the invention, which are, for clarity, described in the context of separate embodiments, may also be provided in combination in a single embodiment. Conversely, various features of the invention which are, for brevity, described in the context of a single embodiment, may also be provided separately or in any suitable subcombination.
It will be appreciated by persons skilled in the art that the present invention is not limited to what has been particularly shown and described hereinabove. Rather, the scope of the present invention is defined by the appended claims and includes both combinations and subcombinations of the various features described hereinabove as well as variations and modifications thereof which would occur to persons skilled in the art upon reading the foregoing description. Accordingly, it is intended to embrace all such alternatives, modifications and variations that fall within the spirit and broad scope of the appended claims.
Claims (19)
1. A method of analyzing EEG data in order to differentiate between ongoing neuro-electric activity of the brain of a subject and evoked neuro-electric activity of the brain of said subject, said analysis being based upon the time, polarity and amplitude distributions of deflections from the base line of said EEG and EP data, said method comprising
a. recording by EEG measurement said ongoing neuro-electric activity;
b. administering a series of sensory stimuli to said subject;
c. recording by EEG and EP measurement neuro-electric activity evoked in response to said stimuli;
d. marking the time of occurrence, polarity and amplitude of each of said deflections recorded;
e. generating a histogram of time periods of said deflections for a defined duration prior to, during and subsequent to each of said stimuli, said histogram reflecting the deflections within said time periods, generating a deflection time distribution based upon data from said histogram, generating an amplitude profile of said deflections; and
f. using said deflection time distribution and said amplitude profile to statistically analyze and evaluate the ongoing and evoked neuro-electric activity,
such that the differentiation between said ongoing neuro-electric activity and said evoked neuro-electric activity may be quantified.
2. The method of claim 1 , wherein said method is modifiable to allow quantification of evoked neuro-electric activity of different frequencies and amplitudes.
3. The method of claim 1 , further comprising identifying at least one of said stimuli evoking said neuro-electric activity.
4. The method of claim 1 , further comprising analyzing at least one of said stimuli evoking said neuro-electric activity, and examining said evoked neuro-electric activity, said examination including one or more of the group consisting of the timing of said neuro-electric activity; the variability of said neuro-electric activity; and the correlation of said neuro-electric activity.
5. The method of claim 1 , wherein said stimuli are sensory stimuli.
6. The method of claim 5 , wherein said sensory stimuli are visual stimuli.
7. The method of claim 1 , wherein said stimuli include performing an attention demanding task.
8. The method of claim 1 , wherein said stimuli include engaging the memory process.
9. The method of claim 1 , wherein said time periods range between 3 ms and 10 ms.
10. The method of claim 1 , wherein said time period is 6 ms.
11. The method of claim 1 , wherein said deflections are both positive and negative deflections.
12. The method of claim 1 , wherein said method is actuatable and operatable by a computer.
13. The method of claim 1 , including using the method to measure the effect of said stimuli on said subject compared to the effect of said stimuli on other subjects.
14. The method of claim 1 , including using the method to measure the effect of said stimuli on said subject over a period of time.
15. The method of claim 1 , including using the method to measure the effect of medicine on a response of said subject to said stimuli.
16. The method of claim 1 , including using the method to measure the effectiveness of a remedial program on a response of said subject to said stimuli.
17. The method of claim 1 , including using the method to measure the effect of fatigue on a response of said subject to said stimuli.
18. The method of claim 1 , including using the method to measure the effect of injury or disease on a response of said subject to said stimuli.
19. The method of claim 1 , including using the method to measure the effect of recovery from injury or disease on a response of said subject to said stimuli.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US10/235,163 US6832110B2 (en) | 2001-09-05 | 2002-09-05 | Method for analysis of ongoing and evoked neuro-electrical activity |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US31697401P | 2001-09-05 | 2001-09-05 | |
US10/235,163 US6832110B2 (en) | 2001-09-05 | 2002-09-05 | Method for analysis of ongoing and evoked neuro-electrical activity |
Publications (2)
Publication Number | Publication Date |
---|---|
US20030073921A1 US20030073921A1 (en) | 2003-04-17 |
US6832110B2 true US6832110B2 (en) | 2004-12-14 |
Family
ID=26928627
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US10/235,163 Expired - Fee Related US6832110B2 (en) | 2001-09-05 | 2002-09-05 | Method for analysis of ongoing and evoked neuro-electrical activity |
Country Status (1)
Country | Link |
---|---|
US (1) | US6832110B2 (en) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050075832A1 (en) * | 2003-09-22 | 2005-04-07 | Ikeguchi Edward F. | System and method for continuous data analysis of an ongoing clinical trial |
US20130216985A1 (en) * | 2010-06-28 | 2013-08-22 | The Regents Of The University Of California | Method of Suppressing of Irrelevant Stimuli |
CN104871160A (en) * | 2012-09-28 | 2015-08-26 | 加利福尼亚大学董事会 | Systems and methods for sensory and cognitive profiling |
US9814426B2 (en) | 2012-06-14 | 2017-11-14 | Medibotics Llc | Mobile wearable electromagnetic brain activity monitor |
US9940844B2 (en) | 2010-11-11 | 2018-04-10 | The Regents Of The University Of California | Enhancing cognition in the presence of distraction and/or interruption |
US10130277B2 (en) | 2014-01-28 | 2018-11-20 | Medibotics Llc | Willpower glasses (TM)—a wearable food consumption monitor |
US10182736B2 (en) | 2012-10-12 | 2019-01-22 | The Regents Of The University Of California | Configuration and spatial placement of frontal electrode sensors to detect physiological signals |
US10258291B2 (en) | 2012-11-10 | 2019-04-16 | The Regents Of The University Of California | Systems and methods for evaluation of neuropathologies |
Families Citing this family (55)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10331800A1 (en) * | 2003-07-14 | 2005-02-03 | Daimlerchrysler Ag | Load probing for drivers |
US8064993B2 (en) * | 2005-11-14 | 2011-11-22 | General Electric Company | Measurement of EEG reactivity |
US8145297B2 (en) | 2005-11-14 | 2012-03-27 | General Electric Company | Measurement of EEG reactivity |
DK2139390T3 (en) * | 2006-12-22 | 2017-11-27 | Neuro-Insight Pty Ltd | Procedure for Evaluating the Effectiveness of Commercial Communication |
AU2008212131A1 (en) * | 2007-02-09 | 2008-08-14 | Agency For Science, Technology And Research | A system and method for processing brain signals in a BCI system |
US7844324B2 (en) * | 2007-02-14 | 2010-11-30 | The General Electric Company | Measurement of EEG reactivity |
US8484081B2 (en) | 2007-03-29 | 2013-07-09 | The Nielsen Company (Us), Llc | Analysis of marketing and entertainment effectiveness using central nervous system, autonomic nervous system, and effector data |
US8386312B2 (en) | 2007-05-01 | 2013-02-26 | The Nielsen Company (Us), Llc | Neuro-informatics repository system |
WO2008137581A1 (en) * | 2007-05-01 | 2008-11-13 | Neurofocus, Inc. | Neuro-feedback based stimulus compression device |
US20090024449A1 (en) * | 2007-05-16 | 2009-01-22 | Neurofocus Inc. | Habituation analyzer device utilizing central nervous system, autonomic nervous system and effector system measurements |
US8392253B2 (en) | 2007-05-16 | 2013-03-05 | The Nielsen Company (Us), Llc | Neuro-physiology and neuro-behavioral based stimulus targeting system |
US8494905B2 (en) | 2007-06-06 | 2013-07-23 | The Nielsen Company (Us), Llc | Audience response analysis using simultaneous electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) |
EP2152155A4 (en) * | 2007-06-06 | 2013-03-06 | Neurofocus Inc | Multi-market program and commercial response monitoring system using neuro-response measurements |
US20090030287A1 (en) * | 2007-06-06 | 2009-01-29 | Neurofocus Inc. | Incented response assessment at a point of transaction |
WO2009018374A1 (en) | 2007-07-30 | 2009-02-05 | Neurofocus, Inc. | Neuro-response stimulus and stimulus attribute resonance estimator |
US20090036755A1 (en) * | 2007-07-30 | 2009-02-05 | Neurofocus, Inc. | Entity and relationship assessment and extraction using neuro-response measurements |
US8386313B2 (en) | 2007-08-28 | 2013-02-26 | The Nielsen Company (Us), Llc | Stimulus placement system using subject neuro-response measurements |
EP2180825A4 (en) * | 2007-08-28 | 2013-12-04 | Neurofocus Inc | Consumer experience assessment system |
US8635105B2 (en) | 2007-08-28 | 2014-01-21 | The Nielsen Company (Us), Llc | Consumer experience portrayal effectiveness assessment system |
US8392255B2 (en) | 2007-08-29 | 2013-03-05 | The Nielsen Company (Us), Llc | Content based selection and meta tagging of advertisement breaks |
US9058473B2 (en) * | 2007-08-29 | 2015-06-16 | International Business Machines Corporation | User authentication via evoked potential in electroencephalographic signals |
US8494610B2 (en) | 2007-09-20 | 2013-07-23 | The Nielsen Company (Us), Llc | Analysis of marketing and entertainment effectiveness using magnetoencephalography |
US20090083129A1 (en) | 2007-09-20 | 2009-03-26 | Neurofocus, Inc. | Personalized content delivery using neuro-response priming data |
WO2009089532A1 (en) * | 2008-01-11 | 2009-07-16 | Oregon Health & Science University | Rapid serial presentation communication systems and methods |
EP2306890B1 (en) * | 2008-07-09 | 2015-12-02 | McKinley, Laurence M. | Optic function monitoring process and apparatus |
US8270814B2 (en) | 2009-01-21 | 2012-09-18 | The Nielsen Company (Us), Llc | Methods and apparatus for providing video with embedded media |
US9357240B2 (en) | 2009-01-21 | 2016-05-31 | The Nielsen Company (Us), Llc | Methods and apparatus for providing alternate media for video decoders |
US8464288B2 (en) | 2009-01-21 | 2013-06-11 | The Nielsen Company (Us), Llc | Methods and apparatus for providing personalized media in video |
US20100250325A1 (en) | 2009-03-24 | 2010-09-30 | Neurofocus, Inc. | Neurological profiles for market matching and stimulus presentation |
US20110046473A1 (en) * | 2009-08-20 | 2011-02-24 | Neurofocus, Inc. | Eeg triggered fmri signal acquisition |
US20110046502A1 (en) * | 2009-08-20 | 2011-02-24 | Neurofocus, Inc. | Distributed neuro-response data collection and analysis |
US8655437B2 (en) * | 2009-08-21 | 2014-02-18 | The Nielsen Company (Us), Llc | Analysis of the mirror neuron system for evaluation of stimulus |
US10987015B2 (en) | 2009-08-24 | 2021-04-27 | Nielsen Consumer Llc | Dry electrodes for electroencephalography |
US9560984B2 (en) | 2009-10-29 | 2017-02-07 | The Nielsen Company (Us), Llc | Analysis of controlled and automatic attention for introduction of stimulus material |
US8209224B2 (en) | 2009-10-29 | 2012-06-26 | The Nielsen Company (Us), Llc | Intracluster content management using neuro-response priming data |
US20110106750A1 (en) | 2009-10-29 | 2011-05-05 | Neurofocus, Inc. | Generating ratings predictions using neuro-response data |
US8335715B2 (en) * | 2009-11-19 | 2012-12-18 | The Nielsen Company (Us), Llc. | Advertisement exchange using neuro-response data |
US8335716B2 (en) | 2009-11-19 | 2012-12-18 | The Nielsen Company (Us), Llc. | Multimedia advertisement exchange |
US8684742B2 (en) | 2010-04-19 | 2014-04-01 | Innerscope Research, Inc. | Short imagery task (SIT) research method |
US8655428B2 (en) | 2010-05-12 | 2014-02-18 | The Nielsen Company (Us), Llc | Neuro-response data synchronization |
US8392251B2 (en) | 2010-08-09 | 2013-03-05 | The Nielsen Company (Us), Llc | Location aware presentation of stimulus material |
US8392250B2 (en) | 2010-08-09 | 2013-03-05 | The Nielsen Company (Us), Llc | Neuro-response evaluated stimulus in virtual reality environments |
US8396744B2 (en) | 2010-08-25 | 2013-03-12 | The Nielsen Company (Us), Llc | Effective virtual reality environments for presentation of marketing materials |
CN103098493A (en) * | 2010-11-12 | 2013-05-08 | 松下电器产业株式会社 | Sound pressure evaluation system, and method and program therefor |
KR101347606B1 (en) * | 2011-11-25 | 2014-01-07 | 한국과학기술연구원 | Device and the method thereo for diagnosis of autism based on robot using eeg |
US9569986B2 (en) | 2012-02-27 | 2017-02-14 | The Nielsen Company (Us), Llc | System and method for gathering and analyzing biometric user feedback for use in social media and advertising applications |
US9292858B2 (en) | 2012-02-27 | 2016-03-22 | The Nielsen Company (Us), Llc | Data collection system for aggregating biologically based measures in asynchronous geographically distributed public environments |
US9451303B2 (en) | 2012-02-27 | 2016-09-20 | The Nielsen Company (Us), Llc | Method and system for gathering and computing an audience's neurologically-based reactions in a distributed framework involving remote storage and computing |
US9060671B2 (en) | 2012-08-17 | 2015-06-23 | The Nielsen Company (Us), Llc | Systems and methods to gather and analyze electroencephalographic data |
US9320450B2 (en) | 2013-03-14 | 2016-04-26 | The Nielsen Company (Us), Llc | Methods and apparatus to gather and analyze electroencephalographic data |
EP4241830A3 (en) * | 2013-11-07 | 2023-10-04 | SafeOp Surgical, Inc. | Systems and methods for detecting nerve function |
US9622702B2 (en) | 2014-04-03 | 2017-04-18 | The Nielsen Company (Us), Llc | Methods and apparatus to gather and analyze electroencephalographic data |
US9936250B2 (en) | 2015-05-19 | 2018-04-03 | The Nielsen Company (Us), Llc | Methods and apparatus to adjust content presented to an individual |
CN108348187A (en) * | 2015-06-26 | 2018-07-31 | 歌迪·沙哈夫 | Method and system for determining the state of mind |
US20210338140A1 (en) * | 2019-11-12 | 2021-11-04 | San Diego State University (SDSU) Foundation, dba San Diego State University Research Foundation | Devices and methods for reducing anxiety and treating anxiety disorders |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3901215A (en) * | 1971-08-20 | 1975-08-26 | Erwin Roy John | Method of testing the senses and cognition of subjects |
US4493327A (en) * | 1982-07-20 | 1985-01-15 | Neurometrics, Inc. | Automatic evoked potential detection |
US4672978A (en) * | 1984-10-17 | 1987-06-16 | Joseph Danto | Brain stimulation of balance response |
US5392788A (en) * | 1993-02-03 | 1995-02-28 | Hudspeth; William J. | Method and device for interpreting concepts and conceptual thought from brainwave data and for assisting for diagnosis of brainwave disfunction |
US6067467A (en) * | 1994-02-07 | 2000-05-23 | New York University | EEG operative and post-operative patient monitoring method |
-
2002
- 2002-09-05 US US10/235,163 patent/US6832110B2/en not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3901215A (en) * | 1971-08-20 | 1975-08-26 | Erwin Roy John | Method of testing the senses and cognition of subjects |
US4493327A (en) * | 1982-07-20 | 1985-01-15 | Neurometrics, Inc. | Automatic evoked potential detection |
US4672978A (en) * | 1984-10-17 | 1987-06-16 | Joseph Danto | Brain stimulation of balance response |
US5392788A (en) * | 1993-02-03 | 1995-02-28 | Hudspeth; William J. | Method and device for interpreting concepts and conceptual thought from brainwave data and for assisting for diagnosis of brainwave disfunction |
US6067467A (en) * | 1994-02-07 | 2000-05-23 | New York University | EEG operative and post-operative patient monitoring method |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050075832A1 (en) * | 2003-09-22 | 2005-04-07 | Ikeguchi Edward F. | System and method for continuous data analysis of an ongoing clinical trial |
US20080046469A1 (en) * | 2003-09-22 | 2008-02-21 | Ikeguchi Edward F | System and method for continuous data analysis of an ongoing clinical trial |
US7752057B2 (en) | 2003-09-22 | 2010-07-06 | Medidata Solutions, Inc. | System and method for continuous data analysis of an ongoing clinical trial |
US20110161101A1 (en) * | 2003-09-22 | 2011-06-30 | Medidata Solutions, Inc. | System and method for continuous data analysis of an ongoing clinical trial |
US20130216985A1 (en) * | 2010-06-28 | 2013-08-22 | The Regents Of The University Of California | Method of Suppressing of Irrelevant Stimuli |
US11139066B2 (en) | 2010-06-28 | 2021-10-05 | The Regents Of The University Of California | Method of suppressing of irrelevant stimuli |
US10672292B2 (en) * | 2010-06-28 | 2020-06-02 | The Regents Of The University Of California | Method of suppressing of irrelevant stimuli |
US9940844B2 (en) | 2010-11-11 | 2018-04-10 | The Regents Of The University Of California | Enhancing cognition in the presence of distraction and/or interruption |
US11049408B2 (en) | 2010-11-11 | 2021-06-29 | The Regents Of The University Of California | Enhancing cognition in the presence of distraction and/or interruption |
US11830379B2 (en) | 2010-11-11 | 2023-11-28 | The Regents Of The University Of California | Enhancing cognition in the presence of distraction and/or interruption |
US9814426B2 (en) | 2012-06-14 | 2017-11-14 | Medibotics Llc | Mobile wearable electromagnetic brain activity monitor |
US9886493B2 (en) | 2012-09-28 | 2018-02-06 | The Regents Of The University Of California | Systems and methods for sensory and cognitive profiling |
CN104871160B (en) * | 2012-09-28 | 2018-12-07 | 加利福尼亚大学董事会 | System and method for feeling and recognizing anatomy |
US10891313B2 (en) | 2012-09-28 | 2021-01-12 | The Regents Of The University Of California | Systems and methods for sensory and cognitive profiling |
CN104871160A (en) * | 2012-09-28 | 2015-08-26 | 加利福尼亚大学董事会 | Systems and methods for sensory and cognitive profiling |
US10182736B2 (en) | 2012-10-12 | 2019-01-22 | The Regents Of The University Of California | Configuration and spatial placement of frontal electrode sensors to detect physiological signals |
US10258291B2 (en) | 2012-11-10 | 2019-04-16 | The Regents Of The University Of California | Systems and methods for evaluation of neuropathologies |
US10130277B2 (en) | 2014-01-28 | 2018-11-20 | Medibotics Llc | Willpower glasses (TM)—a wearable food consumption monitor |
Also Published As
Publication number | Publication date |
---|---|
US20030073921A1 (en) | 2003-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6832110B2 (en) | Method for analysis of ongoing and evoked neuro-electrical activity | |
Twomey et al. | The classic P300 encodes a build‐to‐threshold decision variable | |
Wilson et al. | Spike detection: a review and comparison of algorithms | |
Ayoubian et al. | Automatic seizure detection in SEEG using high frequency activities in wavelet domain | |
Anderer et al. | An E-health solution for automatic sleep classification according to Rechtschaffen and Kales: validation study of the Somnolyzer 24× 7 utilizing the Siesta database | |
Newton et al. | Cognitive event-related potentials in multiple sclerosis | |
van de Velde et al. | Detection of muscle artefact in the normal human awake EEG | |
Eklund et al. | Auditory awareness negativity is an electrophysiological correlate of awareness in an auditory threshold task | |
Retter et al. | All-or-none face categorization in the human brain | |
Correa-Jaraba et al. | Involuntary capture and voluntary reorienting of attention decline in middle-aged and old participants | |
Shahaf et al. | Comprehensive analysis suggests simple processes underlying EEG/ERP–demonstration with the go/no-go paradigm in ADHD | |
Loehr et al. | It's not just my fault: Neural correlates of feedback processing in solo and joint action | |
Hooper | Comparison of the distributions of classical and adaptively aligned EEG power spectra | |
Eklund et al. | Is auditory awareness negativity confounded by performance? | |
Rahne et al. | Sorted averaging—application to auditory event-related responses | |
Beretta-Piccoli et al. | Motor unit synchronization and firing rate correlate with the fractal dimension of the surface EMG: A validation study | |
Huang et al. | Selection of valid and reliable EEG features for predicting auditory and visual alertness levels | |
Ross et al. | Frequency-specific threshold determination with the CERAgram method: basic principle and retrospective evaluation of data | |
CN112587136B (en) | Taste sensory evaluation method and system | |
Lange et al. | Challenging perceptual tasks require more attention: The influence of task difficulty on the N1 effect of temporal orienting | |
Xu et al. | Approximate entropy analysis of event-related potentials in patients with early vascular dementia | |
Boashash et al. | Time-frequency methodology for newborn electroencephalographic seizure detection | |
Suter et al. | Spatial-frequency-tuned attenuation and enhancement of the steady-state VEP by grating adaptation | |
Huupponen et al. | Fuzzy detection of EEG alpha without amplitude thresholding | |
AU8588082A (en) | Brain electrical activity mapping (beam) |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FEPP | Fee payment procedure |
Free format text: PAYOR NUMBER ASSIGNED (ORIGINAL EVENT CODE: ASPN); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
FPAY | Fee payment |
Year of fee payment: 4 |
|
REMI | Maintenance fee reminder mailed | ||
LAPS | Lapse for failure to pay maintenance fees | ||
STCH | Information on status: patent discontinuation |
Free format text: PATENT EXPIRED DUE TO NONPAYMENT OF MAINTENANCE FEES UNDER 37 CFR 1.362 |
|
FP | Lapsed due to failure to pay maintenance fee |
Effective date: 20121214 |