+

WO2018235076A1 - Procédé et système permettant de prédire une réponse à un traitement pharmacologique à partir d'un eeg - Google Patents

Procédé et système permettant de prédire une réponse à un traitement pharmacologique à partir d'un eeg Download PDF

Info

Publication number
WO2018235076A1
WO2018235076A1 PCT/IL2018/050675 IL2018050675W WO2018235076A1 WO 2018235076 A1 WO2018235076 A1 WO 2018235076A1 IL 2018050675 W IL2018050675 W IL 2018050675W WO 2018235076 A1 WO2018235076 A1 WO 2018235076A1
Authority
WO
WIPO (PCT)
Prior art keywords
data
feature
machine learning
eeg
covariance
Prior art date
Application number
PCT/IL2018/050675
Other languages
English (en)
Inventor
Shahar ARZY
Original Assignee
Hadasit Medical Research Services And Development Ltd.
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hadasit Medical Research Services And Development Ltd. filed Critical Hadasit Medical Research Services And Development Ltd.
Publication of WO2018235076A1 publication Critical patent/WO2018235076A1/fr

Links

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems

Definitions

  • the present invention in some embodiments thereof, relates to neurophysiology and, more particularly, but not exclusively, to method and system for predicting a response to pharmacological treatment from EEG.
  • EEG excitatory post synaptic potentials
  • U.S. Patent No. 7,725,174 discloses a method of indicating a subject's reaction to different agents administered to induce anesthesia using two sets of EEG signals are received from two sets of electrodes located on two regions of a head of a subject.
  • U.S. Patent No. 8,632,750 discloses a method for recommending therapy for a patient having a behaviorally diagnosed psychiatric condition, using patient univariate measures extracted from recorded EEG.
  • a method of analyzing EEG data collected from a brain of a subject diagnosed with a neuropsychiatric disorder comprises: extracting a feature from the EEG data; receiving information pertaining to a candidate pharmacological treatment for the neuropsychiatric disorder; executing a machine learning procedure wherein the extracted feature and the received information are inputs of the machine learning procedure; and determining a predicted response of the subject to the pharmacological treatment based on an output of the machine learning procedure.
  • the execution of the machine learning and the determination of the predicted response is by a server system and the method comprises receiving the EEG data and the information by the server system over a communication network.
  • the method comprises transmitting the predicted response to a client computer.
  • the method comprises collecting EEG signals from the brain of the subject, and digitizing the EEG signals to provide the EEG data.
  • the feature comprises covariance between data channels.
  • the method comprises separating the data into an absolute component and a phase component prior to the feature extraction, wherein the covariance is calculated separately for at least one of the components.
  • the method comprises decomposing the data into frequency bands, wherein the covariance is calculated separately for at least one of the frequency bands.
  • the feature comprises a coherence level of the data.
  • the machine learning procedure comprises at least one procedure selected from the group consisting of k-nearest neighbors analysis, support vector machine, clustering, linear modeling, decision tree learning, ensemble learning procedure, neural networks, probabilistic model, graphical model, Bayesian network, and association rule learning.
  • the neuropsychiatric disorder is a disorder affecting the central nervous system.
  • the neuropsychiatric disorder comprises at least one disorder selected from the group consisting of dementia, movement disorder, depression, anxiety disorder, bipolar disorder, schizophrenia, obsessive-compulsive disorder, Parkinson's disease, and chronic pain.
  • the neuropsychiatric disorder comprises epilepsy.
  • the candidate pharmacological treatment comprises administration of a predetermined therapeutic amount of an antiepileptic drug.
  • the candidate pharmacological treatment comprises administration of a predetermined therapeutic amount of a combination of at least two antiepileptic drugs.
  • the candidate pharmacological treatment comprises administration of a predetermined therapeutic amount of a drug selected from the group consisting of valproic acid, levetiracetam, lamotrigine, ezogabine, pregabalin, topiramate, rufinamide, gabapentin, clonazepam, oxcarbazepine, phenobarbital, phenytoin, sultiam and felbamate, or salts or prodrugs or derivatives of any purity or sustained release formulation thereof.
  • a drug selected from the group consisting of valproic acid, levetiracetam, lamotrigine, ezogabine, pregabalin, topiramate, rufinamide, gabapentin, clonazepam, oxcarbazepine, phenobarbital, phenytoin, sultiam and felbamate, or salts or prodrugs or derivatives of any purity or sustained release formulation thereof.
  • the candidate pharmacological treatment comprises administration of a predetermined therapeutic amount of carbamazepine or salts or prodrugs or derivatives of any purity or sustained release formulation thereof.
  • the neuropsychiatric disorder is epilepsy
  • the candidate pharmacological treatment comprises administration of a predetermined therapeutic amount of valproic acid
  • the method comprises decomposing the data into frequency bands
  • the extracting the feature comprises calculating covariance between data channels separately for at least two of the frequency bands
  • the feature comprises a mean of the covariance over the at least two of the frequency bands
  • the machine learning procedure comprises k-nearest neighbors analysis.
  • the neuropsychiatric disorder is epilepsy
  • the candidate pharmacological treatment comprises administration of a predetermined therapeutic amount of levetiracetam
  • the method comprises, prior to the feature extraction, decomposing the data into frequency bands, and separating the data into an absolute component and a phase component, wherein the feature comprises covariance between data channels calculated separately for each component and each frequency band, and wherein the machine learning procedure comprises k-nearest neighbors analysis.
  • the neuropsychiatric disorder is epilepsy
  • the candidate pharmacological treatment comprises administration of a predetermined therapeutic amount of levetiracetam
  • the method comprises, prior to the feature extraction, decomposing the data into frequency bands, and extracting a phase component of the data for each frequency band, wherein the feature comprises covariance between data channels calculated separately for at least two of the frequency bands, and wherein the machine learning procedure comprises k-nearest neighbors analysis.
  • the neuropsychiatric disorder is epilepsy
  • the candidate pharmacological treatment comprises administration of a predetermined therapeutic amount of lamotrigine
  • the method comprises, prior to the feature extraction, decomposing the data into frequency bands, and extracting an absolute component of the data for each frequency band, wherein the extracting the feature comprises calculating covariance between data channels for the absolute component separately for at least two of the frequency bands, wherein the feature comprises a mean of the covariance over the at least two of the frequency bands, and wherein the machine learning procedure comprises k- nearest neighbors analysis.
  • the EEG data are collected after at least one session of the candidate pharmacological treatment, and the method comprises determining a response of the subject to the at least one session of the pharmacological treatment and generating output including a recommendation for future sessions of the pharmacological treatment based on an output of the machine learning procedure.
  • a computer software product comprises a computer-readable medium in which program instructions are stored, which instructions, when read by a data processor, cause the data processor to receive EEG data collected from a brain of a subject diagnosed with a neuropsychiatric disorder, and information pertaining to a candidate pharmacological treatment for the neuropsychiatric disorder, and to execute the method as delineated above and optionally as further detailed below.
  • a system for analyzing EEG data comprises a processor arranged to perform code instructions, comprises: code instructions for receiving EEG data collected from a brain of a subject diagnosed with a neuropsychiatric disorder; code instructions for extracting a feature from the EEG data; code instructions for receiving information pertaining to a candidate pharmacological treatment for the neuropsychiatric disorder; code instructions for executing a machine learning procedure wherein the extracted feature and the received information are inputs of the machine learning procedure; and code instructions for determining a predicted response of the subject to the pharmacological treatment based on an output of the machine learning procedure.
  • the system comprises a transceiver arranged to receive and transmit information on a communication network; wherein the processor is arranged to communicate with the transceiver over the communication network; and wherein the code instructions comprise code instructions for transmitting the predicted response to a client computer.
  • the determining the predicted response comprises classifying whether the treatment is effective, partially effective or not effective.
  • the code instructions comprise code instructions for generating an output indicating a recommended change in the treatment.
  • the code instructions comprise code instructions for generating an output pertaining to potential side effects cause by the treatment.
  • the code instructions comprise code instructions for receiving blood test data of the subject, following the treatment, analyzing the blood test data to determine a response of the subject to the treatment, and comparing the response to the predicted response.
  • the system wherein the analysis comprises determining vitamin deficiency.
  • a data processor such as a computing platform for executing a plurality of instructions.
  • the data processor includes a volatile memory for storing instructions and/or data and/or a non-volatile storage, for example, a magnetic hard-disk and/or removable media, for storing instructions and/or data.
  • a network connection is provided as well.
  • a display and/or a user input device such as a keyboard or mouse are optionally provided as well.
  • FIG. 1 is a flowchart diagram of a method of analyzing EEG data collected from a brain of a subject diagnosed with a neuropsychiatric disorder, according to various exemplary embodiments of the present invention
  • FIG. 2 is a schematic illustration of a machine learning system suitable for use according to some embodiments of the present invention.
  • FIG. 3 is a schematic illustration of a computing platform that can be used according to some embodiments of the present invention. DESCRIPTION OF SPECIFIC EMBODIMENTS OF THE INVENTION
  • the present invention in some embodiments thereof, relates to neurophysiology and, more particularly, but not exclusively, to method and system for predicting a response to pharmacological treatment from EEG.
  • the present inventor devised a technique suitable for determining a predicted response of a subject diagnosed with a neuropsychiatric disorder to a pharmacological treatment.
  • neuropsychiatric disorder refers to any disorder affecting the central nervous system (CNS), and may include, but is not limited to, one or more of epilepsy, movement disorder, depression, dementia, anxiety disorder, bipolar disorder, schizophrenia, obsessive- compulsive disorder, Parkinson's disease, and chronic pain.
  • Epilepsy is a chronic neurological condition that is characterized by recurrent seizures .
  • the primary indication of epilepsy is synchronized electrical activity between large numbers of brain neurons.
  • a seizure a group of neurons in the brain demonstrate a highly synchronized firing pattern.
  • Epileptic seizures are characterized by symptoms ranging from a major motor convulsion to a brief period of lack of awareness.
  • Epileptic seizures are typically classified as being either generalized and or partial.
  • a partial seizure commences from a single location within one of the hemispheres and may result from a brain lesion.
  • Generalized seizures involve the whole cortex simultaneously and are assumed to have a genetic basis.
  • AEDs antiepileptic drugs
  • AEDs antiepileptic drugs
  • modulation of voltage-gated ion channels e.g., sodium, GAB A, and calcium
  • inhibition of synaptic excitation e.g., sodium, GAB A, and calcium
  • enhancement of synaptic inhibition e.g., enhancement of synaptic inhibition
  • the technique of the present embodiments utilizes EEG data for predicting a response of a subject to pharmacological treatment and can therefore provide healthcare providers, e.g. , physicians and patients, as well as to healthcare reimbursement providers, e.g., providers of prescription drug insurance, with insightful information useful for clinical decision making.
  • the method utilizes EEG data for determining a response of the subject to previous pharmacological treatments and so as to allow correcting a dosage and/or drug regimen for future treatments.
  • the technique of the present embodiments can optionally and preferably classify whether the treatment is effective, partially effective or not effective. Based on this classification, an output indicating a recommended change in the treatment can be generated.
  • an output pertaining to potential side effects cause by the treatment is generated.
  • a subject is treated by the treatment, and blood test data of the subject is obtained after the treatment.
  • the blood test data can optionally and preferably be analyzed (e.g., to identify vitamin deficiency) so as to determine a response of the subject to the treatment. This response can be compared to the predicted response.
  • FIG. 1 is a flowchart diagram of a method of analyzing EEG data collected from a brain of a subject diagnosed with a neuropsychiatric disorder, according to various exemplary embodiments of the present invention. It is to be understood that, unless otherwise defined, the operations described hereinbelow can be executed either contemporaneously or sequentially in many combinations or orders of execution. Specifically, the ordering of the flowchart diagrams is not to be considered as limiting. For example, two or more operations, appearing in the following description or in the flowchart diagrams in a particular order, can be executed in a different order (e.g., a reverse order) or substantially contemporaneously. Additionally, several operations described below are optional and may not be executed.
  • At least part of the operations can be can be implemented by a data processing system, e.g. , a dedicated circuitry or a general purpose computer, configured for receiving data and executing the operations described below.
  • At least part of the operations can be implemented by a cloud-computing facility at a remote location.
  • the data processing system or cloud -computing facility can serve, at least for part of the operations as an image processing system, wherein the data received by the data processing system or cloud-computing facility include image data.
  • Computer programs implementing the method of the present embodiments can commonly be distributed to users on a distribution medium such as, but not limited to, a floppy disk, a CD- ROM, a flash memory device and a portable hard drive. From the distribution medium, the computer programs can be copied to a hard disk or a similar intermediate storage medium. The computer programs can be run by loading the computer instructions either from their distribution medium or their intermediate storage medium into the execution memory of the computer, configuring the computer to act in accordance with the method of this invention. All these operations are well-known to those skilled in the art of computer systems.
  • the method of the present embodiments can be embodied in many forms. For example, it can be embodied in on a tangible medium such as a computer for performing the method operations. It can be embodied on a computer readable medium, comprising computer readable instructions for carrying out the method operations. In can also be embodied in electronic device having digital computer capabilities arranged to run the computer program on the tangible medium or execute the instruction on a computer readable medium.
  • the method begins at 10 optionally and preferably continues to 11 at which EEG data are received.
  • the data can be received by recording EEG signals directly from the subject by placing a set of EEG electrodes on the scalp of the subject, and digitizing the EEG signals to provide the data.
  • the EEG data can be collected from the brain of the subject before or after the subject is treated by a pharmacological treatment of interest.
  • EEG signals refers to recordings of cerebral electrical activity, or other types of brain electrically activity, recorded from any location on the cranium.
  • the portion of the data that corresponds to signals acquired by each EEG electrode is referred to as a "data channel.”
  • EEG electrode refers to any electrode placed on a person's head and capable of detecting brain electrical activity. EEG electrodes can be placed according to known positioning systems, such as but not limited to, the expanded international 10/20 placement system.
  • EEG can include cerebral electrical activity or other types of brain electrical activity, and it is to be understood that the technique of the present disclosure can be applied to any type of brain electrical activity signal.
  • the EEG data can alternatively be received from an external source, such as a computer readable memory medium on which the EEG data are stored, in which case the EEG signals are recorded and digitized by a separate system.
  • the method can proceed to 12 at which information pertaining to a candidate pharmacological treatment for the neuropsychiatric disorder is received.
  • the information typically includes a candidate drug or a combination of two or more drugs to be administered to the subject, and optionally and preferably also full or partial information regarding a candidate dosage regimen associated with the candidate drug or drugs (e.g. , formulation, route of administration, unit dose, frequency of treatment, duration of treatment, etc.).
  • the candidate drug is optionally and preferably, but not necessarily, an antiepileptic drug or a combination of two or more antiepileptic drugs. These embodiments are particularly useful when the neuropsychiatric disorder includes epilepsy, but are also useful for treating other neuropsychiatric disorders.
  • Other drugs not necessarily drugs considered antiepileptic, are also contemplated in some embodiments of the present invention.
  • the method proceeds to 13 at which subject profile is received.
  • the subject profile typically includes one or more of weight, body mass index, gender, age, ethnicity, race, allergies, clinical history, family history.
  • subject profile also includes a genetic profile, which can encompass the genes in an entire genome of the subject, or it can encompass a specific subset of genes.
  • the genetic profile may include genomic profile, a proteomic profile, an epigenomic profile and/or a transcriptomic profile.
  • Any of operations 11-13 can be executed digitally by transmitting digital data to a data processor that executes the method, or by instructing data processor to access a computer readable medium and retrieve the data digitally from the medium.
  • the method proceeds to 14 at which the EEG data are preprocessed.
  • the preprocessing operation or set of operations are preferably applied digitally after the EEG signals are digitized.
  • Frequency bands suitable for the present embodiments including, without limitation, the EEG ⁇ band, which is typically less than 4 Hz, but not less than 1 Hz, the EEG ⁇ band, which is typically from about 4 Hz to about 7 Hz, the EEG a band, which is typically from about 8 Hz to about 15 Hz, the EEG ⁇ band, which is typically from about 16 Hz to about 31 Hz, and the EEG ⁇ band, which is typically from about 32 Hz to about 80 Hz.
  • the EEG ⁇ band which is typically less than 4 Hz, but not less than 1 Hz
  • the EEG ⁇ band which is typically from about 4 Hz to about 7 Hz
  • the EEG a band which is typically from about 8 Hz to about 15 Hz
  • the EEG ⁇ band which is typically from about 16 Hz to about 31 Hz
  • the EEG ⁇ band which is typically from about 32 Hz to about 80 Hz.
  • frequency bands which may or may not overlap with these bands (e.g., the EEG ⁇ range, typically from about 8 Hz to about 12 Hz, and/or high and very high frequencies, for example, up to 1000 Hz or more) are also contemplated.
  • the EEG ⁇ range typically from about 8 Hz to about 12 Hz, and/or high and very high frequencies, for example, up to 1000 Hz or more
  • Each of the above bands can be subdivided into two or more sub-bands. Frequency bands that do not necessarily match the aforementioned EEG bands are also contemplated.
  • Another type of preprocessing operation that can be employed according to preferred embodiments of the invention includes removal of artifacts from the EEG data.
  • artifacts refers to any electrical potential recorded while obtaining an EEG or other brain electrical activity signal that is not of cerebral origin or is the result of abnormal brain activity. EEG artifacts are oftentimes also referred to as "noise.”
  • Non-physiological artifacts such as, but not limited to, cable or electrode movement causing over-range artifacts, and impulse artifacts (for example due to due to an abrupt change in the impedance of an electrode); non-brain-generated physiological artifacts, such as, but not limited to, horizontal/lateral eye movements, vertical eye movements (e.g., blinks), and electromyographic activity; and brain-generated artifacts, such as, but not limited to, significantly low amplitude signal, and atypical electrical activity pattern (for example due to paroxysmal brain activity).
  • non-physiological artifacts such as, but not limited to, cable or electrode movement causing over-range artifacts, and impulse artifacts (for example due to due to an abrupt change in the impedance of an electrode)
  • non-brain-generated physiological artifacts such as, but not limited to, horizontal/lateral eye movements, vertical eye movements (e.g., blinks), and electromyographic activity
  • brain-generated artifacts such as
  • Lateral eye movements can be identified, as waveforms of 1 Hz or less that have opposite polarity at EEG channels F7 and F8, and can be removed using a FIR filter with passband of about 0.5-3 Hz.
  • Vertical eye movement can be identified by locating large excursions on the Fpl and Fp2 leads, and can be removed a low-pass filter in the range of about 0.5-5 Hz.
  • Cable or electrode movement artifacts can be identified when the signal amplitude is greater than a threshold which can be, for example, between about 50 ⁇ and about 250 ⁇ .
  • Impulse artifacts can be identified as follows. A frontal EEG channel is first high-pass filtered with cutoff frequency to remove an alpha band from the signal in that channel.
  • the cut-off frequency is 15 Hz.
  • high-frequency variations of signal amplitude in successive segments with overlap are examined.
  • a segment is identified as containing an impulse artifact when the difference between the maximal and minimal levels value exceeds a threshold.
  • the threshold can be between 25 ⁇ and 250 ⁇ .
  • Electromyographic activity can be identified when a relative energy between two sub-bands of an EEG ⁇ band exceeds a threshold.
  • Data segments containing significantly low amplitude signal can be identified by searching for signal epochs with mean-square energy below a threshold, e.g., between 1 ⁇ 2 and 25 ⁇ 2.
  • Atypical electrical activity pattern can be identified using a combination of wavelet analysis and fractal dimension computation.
  • Another type of preprocessing operation that can be employed according to preferred embodiments of the invention includes dimensionality reduction. This can be achieved using any dimensionality reduction technique known in the art, including, without limitation, Principal Component Analysis (PCA), Linear Discriminant Analysis (LDA), Independent Component Analysis (ICA), kernel PCA, and feature selection.
  • PCA Principal Component Analysis
  • LDA Linear Discriminant Analysis
  • ICA Independent Component Analysis
  • kernel PCA feature selection
  • the method optionally and preferably extracts one or more features from the EEG data.
  • feature of EEG refers to any quantitative measure that can be extracted from EEG data.
  • Representative examples include, without limitation, source localization, cross- correlation, covariance, coherence, coupling, cross-covariance, cross-power, symmetry, correlation, autoregressive coefficients, entropy, spectral entropy, energy, energy derivative, filtered amplitude squared, fractal dimension, mean frequency, nonlinear decorrelation lag, nonlinear energy operator, number of zero crossings, Pisarenko harmonic decomposition, power distribution in frequency bands, principal components, principal Lyapunov exponent, cepstrum, spike (occurrence, amplitude, curvature), wavelet subband energy, wavelet compression coefficients, epileptiform discharge complexity (a measure of number of peaks, amplitude, frequency content and morphology of spike waveforms), amount of background disruption (amount of deviation from baseline time and frequency characteristics of electrical signals), regional coherence (coherence of activity in a focal brain region compared to that of
  • Source localization is a feature that identifies the sources in the brain for each activity as measured by the EEG electrodes.
  • the source potential variable x is an ordered set corresponding to a number of discrete sources of electromagnetic radiation, and likewise the measure potential variable y is an ordered set corresponding to a number of EEG electrodes placed at discrete locations.
  • the transformation operator A is a matrix that specifies, generally, the physics of propagating electromagnetic radiation from points inside the body to points outside the body.
  • the operator A can be represented as a matrix and can be determined by modeling or measurement.
  • inverse problem In source localization the inverse of the forward problem is solved ("inverse problem").
  • the solution to the inverse problem works backward from the sensed potentials to infer the set of source potentials that produced them.
  • additional criteria are typically provided in the form of one or more assumptions regarding the expected behavior of the source variable x.
  • a typical assumption is that the source variable x can exhibit spatial smoothness. For example, in what is known as the "minimum norm" approach, the source variable is constrained to minimum variance.
  • LORETA Low Resolution Electromagnetic Tomographic Analysis
  • Utilizing additional criteria allows selection of a solution to the inverse problem in which the criteria are best met, and rejection of alternative solutions in which the criteria are less well met, or are not met.
  • the results of the source localization are supplemented to, or replace the EEG data.
  • the identified sources of brain activity can be used as channels of the EEG data.
  • one or more features of EEG can be extracted for one or more particular frequency bands and/or a particular frequency sub-bands, and more preferably separately for each frequency band and sub- band.
  • covariance between data channels is used as a feature of the EEG data.
  • the covariance is optionally and preferably calculated with respect to one or more components of the data.
  • the data of each channel are separated into an absolute component and a phase component, and the covariance is calculated separately for one of these components or both these components.
  • the covariance of one or both of the absolute and phase components is optionally and preferably calculated separately for at least one of the frequency bands and sub-bands.
  • Use of covariance between data channels calculated separately for each component and each frequency band is particularly useful when the neuropsychiatric disorder is epilepsy, and the pharmacological drug comprises levetiracetam.
  • Use of covariance between channels calculated for the absolute component of the data for each frequency band is particularly useful when the neuropsychiatric disorder is epilepsy, and the pharmacological drug comprises lamotrigine.
  • a covariance quantity, or a covariance quantity of one or both of the absolute and phase components is/are calculated separately for more than one frequency band and/or sub-band, and the mean of these calculated covariance quantities over these frequency bands and/or sub-bands.
  • the mean of the calculated covariance quantities is referred to a mean covariance.
  • the mean covariance is used as a feature of the EEG data.
  • the mean of the calculated covariance quantities is useful, for example, when the neuropsychiatric disorder is epilepsy, and the candidate drug comprises valproic acid.
  • Another feature that can be used according to some embodiments of the present invention is a coherence level of the data. This feature can be calculated separately for one or more of the frequency bands and sub-bands, or for the full range of frequencies.
  • the coherence level typically characterizes the level by which the same frequency bands or sub-band of two data channels preserve their relative phase.
  • the method proceeds to 16 at which a machine learning procedure is executed.
  • machine learning refers to a procedure embodied as a computer program configured to induce patterns, regularities, or rules from previously collected data to develop an appropriate response to future data, or describe the data in some meaningful way.
  • Machine learning system 20 suitable for use according to some embodiments of the present invention is schematically illustrated in FIG. 2.
  • Machine learning system 20 comprises a processing system 22, which can be a part of a central processing unit (CPU) circuit, an input 24 and an output 26.
  • Input 24 can be fed by more than one source.
  • input 24 can receive the extracted EEG features, the information that pertains to the candidate pharmacological treatment, and optionally also the subject profile.
  • Machine learning system 20 is preferably pre-trained in a manner that the received inputs are processes by processing system 22 to provide, at 28, output that is indicative of the predicted response of the subject to the pharmacological treatment.
  • the machine learning procedure comprises, or is, a supervised learning procedure.
  • supervised learning global or local goal functions are used to optimize the structure of the learning system.
  • supervised learning there is a desired response, which is used by the system to guide the learning.
  • the machine learning procedure comprises, or is, an unsupervised learning procedure.
  • unsupervised learning there are typically no goal functions.
  • the learning system is not provided with a set of rules.
  • One form of unsupervised learning according to some embodiments of the present invention is unsupervised clustering in which the data objects are not class labeled, a priori.
  • the machine learning procedure comprises, or is, a semi- supervised learning procedure.
  • a semi- supervised learning the learning system is first provided with set of an annotated data on which a desired response is expected. This set is referred to as supervised data. Then the learning system uses the set of rules extracted from the supervised data to learn another set of data to which set of rules is not provided.
  • the output provided by machine learning system 20 can be a binary output that confirms or disconfirms a predicted existence of a response.
  • output can include one of two values, e.g., a first output value (e.g., " 1") when machine learning system 20 predicts that the subject will be responsive to the candidate pharmacological treatment, and second output value (e.g., "0") when machine learning system 20 predicts that the subject will be non-responsive to the pharmacological treatment.
  • the output can be a non- binary output. In this case, it can include one of a plurality of discreet values or one of a range of values indicative of the predicted level or intensity by which the subject will be responsive.
  • the output can include a set of values, each pertaining to different property of the predicted response, such as, but not limited to, level or intensity of the response, polarization of the response, speed of the response, etc.
  • the set of output values can include one or more values for each of the conditions.
  • machine learning procedures suitable for the present embodiments, include, without limitation, clustering, association rule algorithms, feature evaluation algorithms, subset selection algorithms, support vector machines, classification rules, cost-sensitive classifiers, vote algorithms, stacking algorithms, Bayesian networks, decision trees, neural networks, instance-based algorithms, linear modeling algorithms, k-nearest neighbors (KNN) analysis, ensemble learning algorithms, probabilistic models, graphical models, logistic regression methods (including multinomial logistic regression methods), gradient ascent methods, singular value decomposition methods and principle component analysis.
  • the self-organizing map and adaptive resonance theory are commonly used unsupervised learning algorithms.
  • the adaptive resonance theory model allows the number of clusters to vary with problem size and lets the user control the degree of similarity between members of the same clusters by means of a user-defined constant called the vigilance parameter.
  • Support vector machines are algorithms that are based on statistical learning theory.
  • a support vector machine (SVM) according to some embodiments of the present invention can be used for classification purposes and/or for numeric prediction.
  • a support vector machine for classification is referred to herein as “support vector classifier”
  • support vector machine for numeric prediction is referred to herein as “support vector regression”.
  • An SVM is typically characterized by a kernel function, the selection of which determines whether the resulting SVM provides classification, regression or other functions.
  • the SVM maps input vectors into high dimensional feature space, in which a decision hyper- surface (also known as a separator) can be constructed to provide classification, regression or other decision functions.
  • a decision hyper- surface also known as a separator
  • the surface is a hyper-plane (also known as linear separator), but more complex separators are also contemplated and can be applied using kernel functions.
  • the data points that define the hyper- surface are referred to as support vectors.
  • the support vector classifier selects a separator where the distance of the separator from the closest data points is as large as possible, thereby separating feature vector points associated with objects in a given class from feature vector points associated with objects outside the class.
  • a high-dimensional tube with a radius of acceptable error is constructed which minimizes the error of the data set while also maximizing the flatness of the associated curve or function.
  • the tube is an envelope around the fit curve, defined by a collection of data points nearest the curve or surface.
  • An advantage of a support vector machine is that once the support vectors have been identified, the remaining observations can be removed from the calculations, thus greatly reducing the computational complexity of the problem.
  • An SVM typically operates in two phases: a training phase and a testing phase.
  • a training phase a set of support vectors is generated for use in executing the decision rule.
  • decisions are made using the decision rule.
  • a support vector algorithm is a method for training an SVM. By execution of the algorithm, a training set of parameters is generated, including the support vectors that characterize the SVM.
  • a representative example of a support vector algorithm suitable for the present embodiments includes, without limitation, sequential minimal optimization.
  • KNN analysis the affinity or closeness of objects is determined. The affinity is also known as distance in a feature space between objects.
  • the KNN analysis is a technique to find distance-based outliers based on the distance of an object from its kth-nearest neighbors in the feature space. Specifically, each object is ranked on the basis of its distance to its kth-nearest neighbors. The farthest away object is declared the outlier. In some cases the farthest objects are declared outliers. That is, an object is an outlier with respect to parameters, such as, a k number of neighbors and a specified distance, if no more than k objects are at the specified distance or less from the object.
  • the KNN analysis is a classification technique that uses supervised learning.
  • An item is presented and compared to a training set with two or more classes.
  • the item is assigned to the class that is most common amongst its k-nearest neighbors. That is, compute the distance to all the items in the training set to find the k nearest, and extract the majority class from the k and assign to item.
  • KNN is useful, for example, when the neuropsychiatric disorder is epilepsy, the candidate drug comprises valproic acid, and the extracted feature is a mean of covariance values over two or more frequency bands or sub-bands.
  • KNN is also useful when the neuropsychiatric disorder is epilepsy, the candidate drug comprises levetiracetam, and the EEG feature comprises covariance between channels calculated separately for each component and each frequency band.
  • KNN is also useful when the neuropsychiatric disorder is epilepsy, the candidate drug comprises lamotrigine and the EEG feature comprises covariance between channels calculated for the absolute component of the data for each frequency band.
  • Association rule algorithm is a technique for extracting meaningful association patterns among features.
  • association in the context of machine learning, refers to any interrelation among features, not just ones that predict a particular class or numeric value. Association includes, but it is not limited to, finding association rules, finding patterns, performing feature evaluation, performing feature subset selection, developing predictive models, and understanding interactions between features.
  • association rules refers to elements that co-occur frequently within the datasets. It includes, but is not limited to association patterns, discriminative patterns, frequent patterns, closed patterns, and colossal patterns.
  • a usual primary step of association rule algorithm is to find a set of items or features that are most frequent among all the observations. Once the list is obtained, rules can be extracted from them.
  • the aforementioned self-organizing map is an unsupervised learning technique often used for visualization and analysis of high-dimensional data. Typical applications are focused on the visualization of the central dependencies within the data on the map.
  • the map generated by the algorithm can be used to speed up the identification of association rules by other algorithms.
  • the algorithm typically includes a grid of processing units, referred to as "neurons". Each neuron is associated with a feature vector referred to as observation.
  • the map attempts to represent all the available observations with optimal accuracy using a restricted set of models. At the same time the models become ordered on the grid so that similar models are close to each other and dissimilar models far from each other. This procedure enables the identification as well as the visualization of dependencies or associations between the features in the data.
  • Feature evaluation algorithms are directed to the ranking of features or to the ranking followed by the selection of features based on their impact.
  • Information gain is one of the machine learning methods suitable for feature evaluation.
  • the definition of information gain requires the definition of entropy, which is a measure of impurity in a collection of training instances.
  • the reduction in entropy of the target feature that occurs by knowing the values of a certain feature is called information gain.
  • Information gain may be used as a parameter to determine the effectiveness of a feature in explaining the response to the treatment.
  • Symmetrical uncertainty is an algorithm that can be used by a feature selection algorithm, according to some embodiments of the present invention. Symmetrical uncertainty compensates for information gain's bias towards features with more values by normalizing features to a [0,1] range.
  • Subset selection algorithms rely on a combination of an evaluation algorithm and a search algorithm. Similarly to feature evaluation algorithms, subset selection algorithms rank subsets of features. Unlike feature evaluation algorithms, however, a subset selection algorithm suitable for the present embodiments aims at selecting the subset of features with the highest impact on the subject's response to the treatment, while accounting for the degree of redundancy between the features included in the subset.
  • the benefits from feature subset selection include facilitating data visualization and understanding, reducing measurement and storage requirements, reducing training and utilization times, and eliminating distracting features to improve classification.
  • Two basic approaches to subset selection algorithms are the process of adding features to a working subset (forward selection) and deleting from the current subset of features (backward elimination).
  • forward selection is done differently than the statistical procedure with the same name.
  • the feature to be added to the current subset in machine learning is found by evaluating the performance of the current subset augmented by one new feature using cross-validation.
  • subsets are built up by adding each remaining feature in turn to the current subset while evaluating the expected performance of each new subset using cross-validation.
  • the feature that leads to the best performance when added to the current subset is retained and the process continues.
  • Backward elimination is implemented in a similar fashion. With backward elimination, the search ends when further reduction in the feature set does not improve the predictive ability of the subset.
  • the present embodiments contemplate search algorithms that search forward, backward or in both directions.
  • Representative examples of search algorithms suitable for the present embodiments include, without limitation, exhaustive search, greedy hill-climbing, random perturbations of subsets, wrapper algorithms, probabilistic race search, schemata search, rank race search, and Bayesian classifier.
  • a decision tree is a decision support algorithm that forms a logical pathway of steps involved in considering the input to make a decision.
  • decision tree refers to any type of tree-based learning algorithms, including, but not limited to, model trees, classification trees, and regression trees.
  • a decision tree can be used to classify the datasets or their relation hierarchically.
  • the decision tree has tree structure that includes branch nodes and leaf nodes.
  • Each branch node specifies an attribute (splitting attribute) and a test (splitting test) to be carried out on the value of the splitting attribute, and branches out to other nodes for all possible outcomes of the splitting test.
  • the branch node that is the root of the decision tree is called the root node.
  • Each leaf node can represent a classification (e.g., whether a particular portion of the group dataset matches a particular portion of the subject-specific dataset) or a value.
  • the leaf nodes can also contain additional information about the represented classification such as a confidence score that measures a confidence in the represented classification (i.e., the likelihood of the classification being accurate).
  • the confidence score can be a continuous value ranging from 0 to 1, which a score of 0 indicating a very low confidence (e.g., the indication value of the represented classification is very low) and a score of 1 indicating a very high confidence (e.g., the represented classification is almost certainly accurate).
  • Decision tree can be applied to the extracted EEG feature(s) and/or the possible clinical decisions that may be taken by the clinician. This allows the technique of the present embodiments to aid the clinician in determining the desired AED and/or its dosage.
  • Regression techniques which may be used in accordance with the present invention include, but are not limited to linear Regression, Multiple Regression, logistic regression, probit regression, ordinal logistic regression ordinal Probit-Regression, Poisson Regression, negative binomial Regression, multinomial logistic Regression (MLR) and truncated regression.
  • a logistic regression or logit regression is a type of regression analysis used for predicting the outcome of a categorical dependent variable (a dependent variable that can take on a limited number of values, whose magnitudes are not meaningful but whose ordering of magnitudes may or may not be meaningful) based on one or more predictor variables. Logistic regression may also predict the probability of occurrence for each data point. Logistic regressions also include a multinomial variant. The multinomial logistic regression model is a regression model which generalizes logistic regression by allowing more than two discrete outcomes.
  • a Bayesian network is a model that represents variables and conditional interdependencies between variables.
  • variables are represented as nodes, and nodes may be connected to one another by one or more links.
  • a link indicates a relationship between two nodes.
  • Nodes typically have corresponding conditional probability tables that are used to determine the probability of a state of a node given the state of other nodes to which the node is connected.
  • a Bayes optimal classifier algorithm is employed to apply the maximum a posteriori hypothesis to a new record in order to predict the probability of its classification, as well as to calculate the probabilities from each of the other hypotheses obtained from a training set and to use these probabilities as weighting factors for future predictions of the subject's response to treatment.
  • An algorithm suitable for a search for the best Bayesian network includes, without limitation, global score metric -based algorithm.
  • Markov blanket can be employed. The Markov blanket isolates a node from being affected by any node outside its boundary, which is composed of the node's parents, its children, and the parents of its children.
  • Instance-based techniques generate a new model for each instance, instead of basing predictions on trees or networks generated (once) from a training set.
  • instance in the context of machine learning, refers to an example from a dataset. Instance-based techniques typically store the entire dataset in memory and build a model from a set of records similar to those being tested. This similarity can be evaluated, for example, through nearest-neighbor or locally weighted methods, e.g. , using Euclidian distances. Once a set of records is selected, the final model may be built using several different techniques, such as the naive Bayes.
  • Neural networks are a class of algorithms based on a concept of inter-connected "neurons.”
  • neurons contain data values, each of which affects the value of a connected neuron according to connections with pre-defined strengths, and whether the sum of connections to each particular neuron meets a pre-defined threshold.
  • connection strengths and threshold values a process also referred to as training
  • a neural network can achieve efficient recognition of images and characters.
  • these neurons are grouped into layers in order to make connections between groups more obvious and to each computation of values.
  • Each layer of the network may have differing numbers of neurons, and these may or may not be related to particular qualities of the input data.
  • each of the neurons in a particular layer is connected to and provides input value to those in the next layer. These input values are then summed and this sum compared to a bias, or threshold. If the value exceeds the threshold for a particular neuron, that neuron then holds a positive value which can be used as input to neurons in the next layer of neurons. This computation continues through the various layers of the neural network, until it reaches a final layer. At this point, the output of the neural network routine can be read from the values in the final layer.
  • convolutional neural networks operate by associating an array of values with each neuron, rather than a single value. The transformation of a neuron value for the subsequent layer is generalized from multiplication to convolution.
  • the method proceeds to 17 at which the predicted response of the subject to the pharmacological treatment is determined based on the output of the machine learning procedure. This can be done by associating the output value or values with a human readable descriptor or descriptors.
  • the method can display, print or store in a computer readable medium a report regarding the predicted response.
  • the method optionally and preferably determines a response of the subject to these previous sessions based on the output of the machine learning procedure.
  • the method can generate an output including a recommendation to correcting future sessions of the treatment, based on the output.
  • the output can include a recommendation to increase the dosage of the drug or to change a dosage ratio among two or more drugs, or the like.
  • FIG. 3 illustrates a computing platform that can be used according to some embodiments of the present invention for executing method 10.
  • client computer 30 having a hardware processor 32, which typically comprises an input/output (I/O) circuit 34, a hardware central processing unit (CPU) 36 (e.g., a hardware microprocessor), and a hardware memory 38 which typically includes both volatile memory and non-volatile memory.
  • CPU 36 is in communication with I/O circuit 34 and memory 38.
  • Client computer 30 preferably comprises a graphical user interface (GUI) 42 in communication with processor 32.
  • I/O circuit 34 preferably communicates information in appropriately structured form to and from GUI 42.
  • An EEG system 66 configured for collecting EEG signals using an arrangement of EEG electrodes (not shown), optionally and preferably communicates with client computer 30, particularly with hardware processor.
  • a server computer 50 which can similarly include a hardware processor 52, an I/O circuit 54, a hardware CPU 56, a hardware memory 58.
  • I/O circuits 34 and 54 of client 30 and server 50 computers preferable operate as transceivers that communicate information with each other via a wired or wireless communication.
  • client 30 and server 50 computers can communicate via a network 40, such as a local area network (LAN), a wide area network (WAN) or the Internet.
  • Server computer 50 can be in some embodiments be a part of a cloud computing resource of a cloud computing facility in communication with client computer 30 over the network 40.
  • GUI 42 and processor 32 can be integrated together within the same housing or they can be separate units communicating with each other.
  • GUI 42 can optionally and preferably be part of a system including a dedicated CPU and I/O circuits (not shown) to allow GUI 42 to communicate with processor 32.
  • Processor 32 issues to GUI 42 graphical and textual output generated by CPU 36.
  • Processor 32 also receives from GUI 42 signals pertaining to control commands generated by GUI 42 in response to user input.
  • GUI 42 can be of any type known in the art, such as, but not limited to, a keyboard and a display, a touch screen, and the like.
  • GUI 42 is a GUI of a mobile device such as a smartphone, a tablet, a smartwatch and the like.
  • GUI 42 is a GUI of a mobile device
  • processor 32 the CPU circuit of the mobile device can serve as processor 32 and can execute the code instructions described herein.
  • Client 30 and server 50 computers can further comprise one or more computer-readable storage media 44, 64, respectively.
  • Media 44 and 64 are preferably non-transitory storage media storing computer code instructions as further detailed herein, and processors 32 and 52 execute these code instructions.
  • the code instructions can be run by loading the respective code instructions into the respective execution memories 38 and 58 of the respective processors 32 and 52.
  • Storage media 64 preferably also store a library of reference data as further detailed hereinabove.
  • EEG system 66 optionally and preferably collects EEG signals from the brain of the subject, and transmits the signals to processor 32.
  • Processor 32 can digitize the EEG signals to provide the EEG data as further detailed hereinabove.
  • the EEG data can be stored on storage 44, in which case processor 32 reads the EEG data from storage 44 instead of receiving it from EEG system 66.
  • Processor 32 of client computer 30 receives, via GUI 42, information pertaining to a candidate pharmacological treatment for the neuropsychiatric disorder, and optionally also a subject profile, as further detailed hereinabove.
  • Processor 32 can transmit the EEG data and the input from GUI 42 over the network 40 to server computer 50.
  • Computer 50 can extract a feature from the EEG data, and executes a machine learning procedure wherein the extracted feature, the received information and optionally the subject profile are inputs of the machine learning procedure, as further detailed hereinabove.
  • Computer 50 can also determining the predicted response of the subject to the pharmacological treatment based on the output of the machine learning procedure, and transmit the predicted response to client computer 32.
  • one or more of the operations, or all the operations, can be executed by computer 32.
  • compositions, method or structure may include additional ingredients, steps and/or parts, but only if the additional ingredients, steps and/or parts do not materially alter the basic and novel characteristics of the claimed composition, method or structure.
  • a compound or “at least one compound” may include a plurality of compounds, including mixtures thereof.
  • range format is merely for convenience and brevity and should not be construed as an inflexible limitation on the scope of the invention. Accordingly, the description of a range should be considered to have specifically disclosed all the possible subranges as well as individual numerical values within that range. For example, description of a range such as from 1 to 6 should be considered to have specifically disclosed subranges such as from 1 to 3, from 1 to 4, from 1 to 5, from 2 to 4, from 2 to 6, from 3 to 6 etc., as well as individual numbers within that range, for example, 1, 2, 3, 4, 5, and 6. This applies regardless of the breadth of the range.
  • Epilepsy is a chronic disease, affecting approximately 1% of the population. Most patients start suffering from the disease in childhood and adolescence or when they are over 65. The majority of epileptic seizures are controlled pharmaceutically by antiepileptic drugs (AEDs), which reduce neuronal burst firing and synchronization in localized neuronal ensembles and inhibit the spread of abnormal firing to distant brain sites via three main molecular mechanisms: (i) modulation of voltage-gated ion channels (mostly calcium and sodium); (ii) inhibition of synaptic excitation; and (iii) enhancement of synaptic inhibition. There are several different types of AEDs, four of which are considered primary, as they cover together roughly 80% of the cases.
  • EEG recordings of patients from diverse age groups have been collected, and different machine learning procedures were employed to classify the recordings into two groups: recordings of patients successfully medicated by one AED, and recordings of patients that were not medicated at the time of the EEG recording.
  • the task of predicting the effectiveness of an AED is formulated into a supervised binary classification problem.
  • a machine learner was given a training set with example from two classes, along with labels that indicates from which class each example was sampled.
  • the task of the learner was to find a prediction rule that given an unlabeled example predicts its label.
  • the performance of the learner was measured on a new unseen data set, referred to as the test set.
  • the two classes were epileptic patients with no medicine, referred to herein as class DO, and epileptic patients that are stable with a treatment of a single AED, where "stable" was defined as either having no seizure over the past year, referred to herein as subclass D4, or stable clinically with no symptoms in the EEG but for less than a year, referred to herein as subclass D3.
  • LamictalTM For each AED a different learner was trained and its performance was evaluated separately from the other two. Because of the difficulty in collecting a large data set of patients, a leave-one-out method was employed where at each iteration a learner is trained on the entire dataset minus one patient that might have one or more records, and then tested on that patient.
  • the average accuracy of the learner on the patients was calculated to provide an estimation of the learner's performance.
  • a preprocessing was applied the EEG data.
  • the preprocessing consisted of high pass filter, which passed only frequencies higher than 0.1 Hz; a low pass filter which passed only frequencies lower than 70 Hz; a notch filter at 50 Hz; and an electrode interpolation using spherical spline interpolating method. Good segments, at different lengths, were selected from the data by a human expert.
  • covariance and coherence Two different features of the EEG data were tested: covariance and coherence.
  • covariance feature the covariance of the phase and the power in different frequency bands, in each pair of electrodes.
  • a Fourier Transform of the signal was computed at each electrode by FFT, and the data were separated into the following frequency bands 1-4 Hz, 4-8 Hz, 8-13 Hz, 13-30 Hz and 30-50 Hz, referred to in this example as delta, theta, alpha, beta and gamma, respectively.
  • the power covariance of the FFT result the covariance applied to the absolute component
  • phase covariance the covariance applied to the angle component
  • the cross spectrum G uv (fn) estimate based on the complex valued coefficients F u k (f n ), F u k ifn) for channel pair ( ⁇ , ⁇ ) is:
  • the coherence was calculated between all pairs of data channels.
  • the KNN learner stored, at the training stage, all the training data into the memory. When a new example was presented, the learner found the K closest examples from the training set according to a similarity measure, and predicts the most common label among those examples. Both LI and L2 norms were used for distance measure, with similar results. The SVM learner found a maximal margin hyperplane that divided the training data into two classes, where the margin was defines as the minimal distance of an example from the hyperplane.
  • Each patient in the data set had one or more records, which might have different labels. From each record, a number of clean segments were extracted by the human expert. The learner was trained on all of the segments from the current training set. When a new record was presented for classification, the learner predicted a label for each of the record's segments. In order to predict the label of the record, voting over the labels of the different segments was employed. Several voting methods were tested, all providing similar results. In this example, a simple majority vote is presented.
  • Weighted accuracy was used as performance measure.
  • the weighted accuracy was defined as (TP/No. of Positive + TN/No. of Negative)/2, where TP and TN are true positive and true negative, respectively.
  • phase and power frequency covariance were used and the performances of SVM and KNN were compared on the different data set.
  • the K parameter for KNN was chosen empirically.
  • the obtained predicted accuracy values are summarized in Table 2.
  • Table 4 demonstrates one or two frequency bands with a better prediction accuracy than the others for DepaleptTM and LamictalTM, but a small difference between the prediction accuracy for each frequency band for KeppraTM.
  • PCA was employed to reduce the data dimension of the frequency representation (both power and phase) from 3610 to 100.
  • the obtained predicted accuracy values for band frequency covariance are summarized in Table 5.
  • Table 5 demonstrates that PCA improves the prediction accuracy of the frequency band covariance for DepaleptTM and LamictalTM, which is in line with the results presented in Table 4, above.
  • PCA was also applied for the mean of covariance values as calculated per frequency band. In this case, the dimension was reduced from 722 to 100.
  • the predicted accuracy values for the mean of covariance values are summarized in Table 6.
  • Table 6 demonstrates that the prediction accuracy of the mean of covariance values is not affected by the dimension of the data.
  • PDF conditioned probability density function
  • the PDFs were estimated using kernel density estimation which allows estimating the PDFs without making assumptions about the distribution of the data.
  • the following two Kullback Leibler divergences were calculated for Q + and Q " , D_KL(Q " II Q + ) and D_KL(Q + H Q " ).
  • the average of these divergences was computed and used as a score for selecting the features. This allowed selecting the feature or features that best separated between the two classes, since the Kullback Leibler divergence between the conditional probability densities can serve as a measure of separateness.
  • the obtained predicted accuracy values of the forward greedy feature selection for the band frequency covariance are summarized in Table 7.
  • S denotes the number of features (with highest score) that were selected.
  • the present Example demonstrates the ability of the technique of the present embodiments to classify EEG data into two classes: stable patients medicated by a single medicine, and unstable patients, with no medicine.
  • KNN analysis was found to be a better learner for this task than SVM.
  • Some other future direction might include features with larger interpretability - for instance, it might be useful to look at the signal at the brain sources rather than the electrodes. This can be done by finding the sources of the activity recorded in the EEG using inverse solution methods.

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Pathology (AREA)
  • Epidemiology (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

L'invention concerne un procédé d'analyse de données d'EEG collectées dans un cerveau d'un sujet chez qui un trouble neuropsychiatrique a été diagnostiqué. Le procédé consiste à extraire une caractéristique des données d'EEG et à recevoir des informations concernant un traitement pharmacologique candidat pour le trouble neuropsychiatrique. Le procédé consiste également à exécuter une procédure d'apprentissage automatique, la caractéristique extraite et les informations reçues étant des entrées de la procédure d'apprentissage automatique, et à déterminer une réponse prédite du sujet au traitement pharmacologique sur la base d'une sortie de la procédure d'apprentissage automatique.
PCT/IL2018/050675 2017-06-21 2018-06-18 Procédé et système permettant de prédire une réponse à un traitement pharmacologique à partir d'un eeg WO2018235076A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US201762522780P 2017-06-21 2017-06-21
US62/522,780 2017-06-21

Publications (1)

Publication Number Publication Date
WO2018235076A1 true WO2018235076A1 (fr) 2018-12-27

Family

ID=63077929

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/IL2018/050675 WO2018235076A1 (fr) 2017-06-21 2018-06-18 Procédé et système permettant de prédire une réponse à un traitement pharmacologique à partir d'un eeg

Country Status (1)

Country Link
WO (1) WO2018235076A1 (fr)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109754886A (zh) * 2019-01-07 2019-05-14 广州达美智能科技有限公司 治疗方案智能生成系统、方法及可读存储介质、电子设备
CN110765978A (zh) * 2019-11-04 2020-02-07 西安邮电大学 一种基于分形维数的通道选择方法
CN111000555A (zh) * 2019-11-29 2020-04-14 中山大学 一种癫痫脑电信号的训练数据生成方法、自动识别模型建模方法和自动识别方法
CN111938594A (zh) * 2020-08-13 2020-11-17 山东大学 基于层次图模型的癫痫脑电信号异常监测系统及设备
CN112349371A (zh) * 2020-11-18 2021-02-09 南通市第一人民医院 一种化疗病人药物记录评估方法和装置
US20210401356A1 (en) * 2018-04-26 2021-12-30 The Penn State Research Foundation Biological marker and methods
US20240047040A1 (en) * 2021-02-10 2024-02-08 Hitachi, Ltd. Information processing system and information processing method
CN118380126A (zh) * 2024-06-25 2024-07-23 武汉飞宇益克科技有限公司 一种医用物品的流转管理方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7725174B2 (en) 2006-06-28 2010-05-25 The University Of Utah Research Foundation Distinguishing different drug effects from the electroencephalogram
US8632750B2 (en) 1997-09-06 2014-01-21 Cns Response, Inc. Methods for recommending neurophysiological disorder therapy

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8632750B2 (en) 1997-09-06 2014-01-21 Cns Response, Inc. Methods for recommending neurophysiological disorder therapy
US7725174B2 (en) 2006-06-28 2010-05-25 The University Of Utah Research Foundation Distinguishing different drug effects from the electroencephalogram

Non-Patent Citations (11)

* Cited by examiner, † Cited by third party
Title
ARZY ET AL., EUROPEAN JOURNAL OF NEUROLOGY, vol. 17, no. 10, 2010, pages 1308 - 1312
KHODAYARI-ROSTAMABAD A ET AL: "A pilot study to determine whether machine learning methodologies using pre-treatment electroencephalography can predict the symptomatic response to clozapine therapy", CLINICAL NEUROPHYSIOLOGY, ELSEVIER SCIENCE, IE, vol. 121, no. 12, 1 December 2010 (2010-12-01), pages 1998 - 2006, XP027448967, ISSN: 1388-2457, [retrieved on 20100617], DOI: 10.1016/J.CLINPH.2010.05.009 *
KHODAYARI-ROSTAMABAD A ET AL: "Using pre-treatment EEG data to predict response to SSRI treatment for MDD", 2010 ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY : (EMBC 2010) ; BUENOS AIRES, ARGENTINA, 31 AUGUST - 4 SEPTEMBER 2010, IEEE, PISCATAWAY, NJ, USA, 31 August 2010 (2010-08-31), pages 6103 - 6106, XP032109523, ISBN: 978-1-4244-4123-5, DOI: 10.1109/IEMBS.2010.5627823 *
KHODAYARI-ROSTAMABAD AHMAD ET AL: "Latent Variable Dimensionality Reduction Using a Kullback-Leibler Criterion and Its Application to Predict Antidepressant Treatment Response", 2013 INTERNATIONAL WORKSHOP ON PATTERN RECOGNITION IN NEUROIMAGING, IEEE, 22 June 2013 (2013-06-22), pages 148 - 151, XP032486013, DOI: 10.1109/PRNI.2013.46 *
KHODAYARI-ROSTAMABAD ET AL., CLINICAL NEUROPHYSIOLOGY, vol. 121, no. 12, 2010, pages 1998 - 2006
KHODAYARI-ROSTAMABAD ET AL.: "A machine learning approach using EEG data to predict response to SSRI treatment for major depressive disorder", CLINICAL NEUROPHYSIOLOGY, vol. 124, 2013, pages 1975 - 1985, XP028703547, DOI: doi:10.1016/j.clinph.2013.04.010
RAVAN MARYAM ET AL: "A machine learning approach using auditory odd-ball responses to investigate the effect of Clozapine therapy", CLINICAL NEUROPHYSIOLOGY, ELSEVIER SCIENCE, IE, vol. 126, no. 4, 27 August 2014 (2014-08-27), pages 721 - 730, XP029203467, ISSN: 1388-2457, DOI: 10.1016/J.CLINPH.2014.07.017 *
SINISA COLIC ET AL: "Prediction of antiepileptic drug treatment outcomes using machine learning", JOURNAL OF NEURAL ENGINEERING, INSTITUTE OF PHYSICS PUBLISHING, BRISTOL, GB, vol. 14, no. 1, 30 November 2016 (2016-11-30), pages 16002, XP020312943, ISSN: 1741-2552, [retrieved on 20161130], DOI: 10.1088/1741-2560/14/1/016002 *
SONJA SIMPRAGA ET AL: "EEG machine learning for accurate detection of cholinergic intervention and Alzheimer's disease", SCIENTIFIC REPORTS, vol. 7, no. 1, 18 July 2017 (2017-07-18), pages 5775, XP055503684, DOI: 10.1038/s41598-017-06165-4 *
SONJA SIMPRAGA: "EEG biomarker integration for better decision making in clinical trials", INTERNATIONAL PHARMACO-EEG SOCIETY, 26 October 2016 (2016-10-26), http://www.ipeg-society.org/userfiles/files/IPEG%202016-.pdf, pages 41 - 42, XP055503693 *
WAJID MUMTAZ ET AL: "A wavelet-based technique to predict treatment outcome for Major Depressive Disorder", PLOS ONE, vol. 12, no. 2, 2 February 2017 (2017-02-02), pages e0171409, XP055503720, DOI: 10.1371/journal.pone.0171409 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20210401356A1 (en) * 2018-04-26 2021-12-30 The Penn State Research Foundation Biological marker and methods
US11744508B2 (en) * 2018-04-26 2023-09-05 The Penn State Research Foundation Biological marker and methods
CN109754886A (zh) * 2019-01-07 2019-05-14 广州达美智能科技有限公司 治疗方案智能生成系统、方法及可读存储介质、电子设备
CN110765978A (zh) * 2019-11-04 2020-02-07 西安邮电大学 一种基于分形维数的通道选择方法
CN110765978B (zh) * 2019-11-04 2022-08-16 西安邮电大学 一种基于分形维数的通道选择方法
CN111000555A (zh) * 2019-11-29 2020-04-14 中山大学 一种癫痫脑电信号的训练数据生成方法、自动识别模型建模方法和自动识别方法
CN111000555B (zh) * 2019-11-29 2022-09-30 中山大学 一种癫痫脑电信号的训练数据生成方法、自动识别模型建模方法和自动识别方法
CN111938594A (zh) * 2020-08-13 2020-11-17 山东大学 基于层次图模型的癫痫脑电信号异常监测系统及设备
CN111938594B (zh) * 2020-08-13 2021-07-02 山东大学 基于层次图模型的癫痫脑电信号异常监测系统及设备
CN112349371A (zh) * 2020-11-18 2021-02-09 南通市第一人民医院 一种化疗病人药物记录评估方法和装置
US20240047040A1 (en) * 2021-02-10 2024-02-08 Hitachi, Ltd. Information processing system and information processing method
CN118380126A (zh) * 2024-06-25 2024-07-23 武汉飞宇益克科技有限公司 一种医用物品的流转管理方法及系统

Similar Documents

Publication Publication Date Title
WO2018235076A1 (fr) Procédé et système permettant de prédire une réponse à un traitement pharmacologique à partir d'un eeg
Siuly et al. Designing a robust feature extraction method based on optimum allocation and principal component analysis for epileptic EEG signal classification
Ahmadlou et al. Fuzzy synchronization likelihood with application to attention-deficit/hyperactivity disorder
Alvi et al. Neurological abnormality detection from electroencephalography data: a review
Soni et al. Graphical representation learning-based approach for automatic classification of electroencephalogram signals in depression
Sameer et al. CNN based framework for detection of epileptic seizures
Agarwal et al. Classification of alcoholic and non-alcoholic EEG signals based on sliding-SSA and independent component analysis
Ibrahim et al. Deep‐learning‐based seizure detection and prediction from electroencephalography signals
Baghdadi et al. A channel-wise attention-based representation learning method for epileptic seizure detection and type classification
Abenna et al. EEG-based BCI: A novel improvement for EEG signals classification based on real-time preprocessing
Karthick et al. Automated detection of subthalamic nucleus in deep brain stimulation surgery for Parkinson’s disease using microelectrode recordings and wavelet packet features
Li et al. A between-subject fNIRS-BCI study on detecting self-regulated intention during walking
Moussa et al. Explainable computer-aided detection of obstructive sleep apnea and depression
Grubov et al. Two-stage approach with combination of outlier detection method and deep learning enhances automatic epileptic seizure detection
Shen et al. GA-SVM modeling of multiclass seizure detector in epilepsy analysis system using cloud computing
Fatma et al. Survey on Epileptic Seizure Detection on Varied Machine Learning Algorithms
Umair et al. Emotion fusion-sense (Emo Fu-sense)–a novel multimodal emotion classification technique
Hosseini Brain-computer interface for analyzing epileptic big data
Huang et al. Analysis and Classification of Sleep Stages based on common frequency pattern from a single-channel EEG signal
Badiei et al. A novel approach for sleep arousal disorder detection based on the interaction of physiological signals and metaheuristic learning
Chakrabarti et al. A review on epileptic seizure detection and prediction using soft computing techniques
Baghersalimi et al. M2SKD: Multi-to-Single Knowledge Distillation of Real-Time Epileptic Seizure Detection for Low-Power Wearable Systems
Sidaoui et al. Epilepsy Seizure Prediction from EEG Signal Using Machine Learning Techniques.
Ghasimi et al. Assessment of cognitive workload using simultaneous EEG and fNIRS: A comparison of feature combinations
Borhade et al. ResneXt-Lenet: A hybrid deep learning for epileptic seizure prediction

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 18749161

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 18749161

Country of ref document: EP

Kind code of ref document: A1

点击 这是indexloc提供的php浏览器服务，不要输入任何密码和下载