NON-INVASIVE DETERMINATION OF A PHYSIOLOGICAL STATE OF INTEREST IN A SUBJECT FROM SPECTRAL DATA PROCESSED USING A TRAINED MACHINE LEARNING MODEL

Information

  • Patent Application
  • 20240130690
  • Publication Number
    20240130690
  • Date Filed
    February 11, 2022
    2 years ago
  • Date Published
    April 25, 2024
    10 days ago
Abstract
Methods, systems, and techniques for determining a physiological state of interest of a subject without direct reference to analytes of the subject. Light is directed at a body part of a subject such that the light passes through or is reflected by blood and interstitial fluid of the body part. The light is incident on the body part comprises a range of wavelengths from at least one of the near infrared and visible spectra. A spectrum of the light is measured after the light has one or both of passed through and been reflected by the body part, and the spectrum comprises the range of wavelengths. Determining whether the subject is in the physiological state of interest involves using a trained machine learning model to process the measured spectrum. This machine learning model is trained with reference spectra representative of the physiological state of interest.
Description
TECHNICAL FIELD

The present disclosure relates to methods, systems, and techniques for non-invasive determination of a physiological state of interest in a subject from spectral data processed using a trained machine learning model.


BACKGROUND

The physiological state of a subject may be established by measuring a range of metabolic parameters using a variety of measurement tools. For example, Haslacher H., et. al., (2017, PLoS ONE 12(5): e0177174. doi.org/10.1371/journal. pone.0177174) teaches the measurement of 12 blood constituents to predict physical capability of an elderly subject. However, a sample of blood is required for the analysis and the various blood constituents are measured using a range of techniques, including photometric, enzymatic, enzymatic colourimetry, ELISA, and others.


SUMMARY

According to a first aspect, there is provided a method of non-invasively determining a physiological state of interest in a subject comprising, (a) placing a body part in contact with a receptor; (b) directing a source of electromagnetic radiation (EMR) over a range of wavelengths through the receptor and onto body part so that the EMR reaches the blood and interstitial fluid within the body part; (c) measuring the EMR absorbed by, reflected by, or transmitted through, the blood and interstitial fluid of the body part with a detector to obtain a spectrum over the range of wavelengths; (d) performing a quantitative mathematical analysis of the spectrum using an algorithm to determine an amount of two or more than two analytes within the blood and interstitial fluid of the body part, wherein the two or more than two analytes comprise multiple ghost analytes whose identities are unknown and that are observed to change in response to the state of interest, wherein the state of interest is impairment or intoxication of the subject; (e) comparing the amount of the two or more than two analytes against a reference value of the two or more analytes to derive a biochemical profile; and (f) analyzing the biochemical profile to determine the physiological state of interest in the subject.


In the step of directing, the source of EMR may be provided over a range of wavelengths from about 400 to about 2500 nm.


The physiological state of interest in the subject may be selected from a group of intoxication arising from cannabis, alcohol, a combination of cannabis and alcohol, opiates, fentanyl, amphetamines, phencyclidine, sedatives, anxyolytics, cocaine, caffeine, and nicotine consumption.


The physiological state of interest may be i) cannabis-induced intoxication, and two or more than two analytes are selected from the group of: delta-9-tetrahydrocannabinol (THC), THC glucuronide (THCGlu), 11-nor-9-carboxy-THC (THC-COOH), 11-hydroxy THC (11-OH-THC), THC-COOH/11-OH- THC ratio, 11-nor-9-carboxy-THC glucuronide (THC-COOGlu), cannabidol (CBD), cannbinol (CBN), cannabigerol (CBG), delta-9-tetrahydrocannabivarin (THCV), THCV-carboxylic acid, 11-nor-9-carboxy-delta-tetrahydrocannabivarin (THCV-COOH),total protein, bilirubin, prolactin, triglycerides, creatinine, cortisol, glucose, lactate, Total 4, uric acid, blood urea nitrogen (BUN), blood sugar, calcium, ionized calcium, magnesium, sodium, phosphate, and GABA; ii) alcohol-induced intoxication, and two or more than two analytes are selected from the group of alcohol, aldehyde, lactic acid; or iii) intoxication arising from cannabis, alcohol, opiates, fentanyl, amphetamines, phencyclidine, sedatives, anxyolytics, cocaine, caffeine, and nicotine consumption, then two or more than two analytes may include:, then two or more than two analytes may include: delta-9-tetrahydrocannabinol (THC), THC glucuronide (THCGlu), 11-nor-9-carboxy-THC (THC-COOH), 11-hydroxy THC (11-OH-THC), THC-COOH/11-OH-THC ratio, 11-nor-9-carboxy-THC glucuronide (THC-COOGlu), cannabidol (CBD), cannbinol (CBN), cannabigerol (CBG), delta-9-tetrahydrocannabivarin (THCV), THCV-carboxylic acid, 11-nor-9-carboxy-delta-tetrahydrocannabivarin (THCV-COOH), albumin, apolipoproteins A1 and B (apoA1 and apoB), total protein, bilirubin, prolactin, triglycerides, creatinine, cortisol, glucose, lactate, Total 4, uric acid, blood urea nitrogen (BUN), blood sugar, calcium, ionized calcium, magnesium, sodium, phosphate, GABA, alcohol, aldehyde, and lactic acid.


In the step of analyzing (step f), the physiological state of interest in the subject may be determined by processing a plurality of data sets that are representative of the biochemical profile and that have been obtained from a plurality of subjects, cross validating the plurality of data sets, and training one or more deep neural networks, support vector machines, convolution neural networks, and generalized additive models, to develop a model comprising one or more algorithms used to identify sets of analytes associated with the status of the physiological state of interest, and the model may be used to analyze the biochemical profile of the subject to determine the physiological state of interest in the subject.


The model may be iteratively trained and validated using different data sets, to produce a validated model. The validated model may comprise one or more algorithms used to identify sets of analytes associated with the status of the physiological state of interest, and the model may be used to analyze the biochemical profile of the subject to determine the physiological state of interest in the subject.


According to another aspect, there is provided a method of non-invasively determining a physiological state of interest of a subject comprising, (a) determining one or more than one physiological parameter of the subject; (b) placing a body part in contact with a receptor; (c) directing a source of electromagnetic radiation (EMR) over a range of wavelengths through the receptor and onto body part so that the EMR reaches the blood and interstitial fluid within the body part; (d) measuring the EMR absorbed by, reflected by, or transmitted through, the blood and interstitial fluid of the body part with a detector to obtain a spectrum over the range of wavelengths; (e) performing a quantitative mathematical analysis of the spectrum using an algorithm to determine an amount of two or more than two analytes within the blood and interstitial fluid of the body part, wherein the two or more than two analytes comprise multiple ghost analytes whose identities are unknown and that are observed to change in response to the state of interest, wherein the state of interest is impairment or intoxication of the subject; (f) comparing the amount of the two or more than two analytes against a reference value of the two or more analytes to derive a biochemical profile; and (g) analyzing the biochemical profile and the one or more than one physiological parameter used to determine the physiological state of interest in the subject.


In the step of directing, the source of EMR may be provided over a range of wavelengths from about 400 to about 2500 nm.


The physiological state of interest of the subject may be selected from a group of intoxication arising from cannabis, alcohol, a combination of cannabis and alcohol, opiates, fentanyl, amphetamines, phencyclidine, sedatives, anxyolytics, cocaine, caffeine, and nicotine consumption.


The physiological state of interest may be: i) cannabis induced intoxication, then two or more than two analytes may include: delta-9-tetrahydrocannabinol (THC), THC glucuronide (THCGlu), 11-nor-9-carboxy-THC (THC-COOH), 11-hydroxy THC (11-OH-THC), THC-COOH/11-OH-THC ratio, 11-nor-9-carboxy-THC glucuronide (THC-COOGlu), cannabidol (CBD), cannbinol (CBN), cannabigerol (CBG), delta-9-tetrahydrocannabivarin (THCV), THCV-carboxylic acid, 11-nor-9-carboxy-delta-tetrahydrocannabivarin (THCV-COOH), albumin, apolipoproteins A1 and B (apoA1 and apoB), total protein, bilirubin, prolactin, triglycerides, creatinine, cortisol, glucose, lactate, Total 4, uric acid, blood urea nitrogen (BUN), blood sugar, calcium, ionized calcium, magnesium, sodium, phosphate, and GABA, the physiological parameter may include one or more of: heart rate, pulse rate, body temperature, neuropeptide Y, fatty acid amide hydrolase (FAAH), c reactive protein (cRP), creatine kinase (CK), aspartate amino transferase (AAT), asparate aminotransferase (AST), alanine transaminase (ALT), gamma-glutamyl transpeptidase (GGT), white blood cell count (WBC), red blood cell count (RBC), hemoglobin, hematocrit, neutrophils, lymphocytes, eosinophils, hypoactivity; THC in hair, THC in urine; ii) alcohol-induced intoxication, and two or more than two analytes are selected from the group of alcohol, aldehyde, lactic acid, and the physiological parameter is selected from the group of heart rate, pulse rate, body temperature, neuropeptide Y, aspartate amino transferase (AAT), alanine transaminase (ALT), gamma-glutamyl transpeptidase (GGT); iii) intoxication resulting from a combination of cannabis and alcohol, or opiates, fentanyl, amphetamines, phencyclidine, sedatives, anxyolytics, cocaine, caffeine, and nicotine consumption, then two or more than two analytes may include: then two or more than two analytes may include: delta-9-tetrahydrocannabinol (THC), THC glucuronide (THCGlu), 11-nor-9-carboxy-THC (THC-COOH), 11-hydroxy THC (11-OH-THC), THC-COOH/11-OH-THC ratio, 11-nor-9-carboxy-THC glucuronide (THC-COOGlu), cannabidol (CBD), cannbinol (CBN), cannabigerol (CBG), delta-9-tetrahydrocannabivarin (THCV), THCV-carboxylic acid, 11-nor-9-carboxy-delta-tetrahydrocannabivarin (THCV-COOH), albumin, apolipoproteins A1 and B (apoA1 and apoB), total protein, bilirubin, prolactin, triglycerides, creatinine, cortisol, glucose, lactate, Total 4, uric acid, blood urea nitrogen (BUN), blood sugar, calcium, ionized calcium, magnesium, sodium, phosphate, GABA, alcohol, aldehyde, and lactic acid, the physiological parameter may include one or more of: heart rate, pulse rate, body temperature, neuropeptide Y, fatty acid amide hydrolase (FAAH), c reactive protein (cRP), creatine kinase (CK), aspartate amino transferase (AAT), asparate aminotransferase (AST), alanine transaminase (ALT), gamma-glutamyl transpeptidase (GGT), white blood cell count (WBC), red blood cell count (RBC), hemoglobin, hematocrit, neutrophils, lymphocytes, eosinophils, hypoactivity; THC in hair, THC in urine.


In the step of analyzing (step g), the biochemical profile and the one or more than one physiological parameter used to determine the physiological state of interest in the subject may be determined by processing a plurality of data sets obtained from a plurality of subjects, each data set derived from the biochemical profile and one or more than one physiological parameter, cross validating the plurality of data sets, and training one or more deep neural networks, support vector machines, convolution neural networks, and generalized additive models, to develop a model comprising one or more algorithms used to identify sets of analytes associated with the status of the physiological state of interest, and the model is used to analyze the biochemical profile and the one or more than one physiological parameter of the subject to determine the physiological state of interest in the subject.


The model may be iteratively trained and validated using different data sets, to produce a validated model. The validated model may comprise one or more algorithms used to identify sets of analytes associated with the status of the physiological state of interest, and the validated model may be used to analyze the biochemical profile and the one or more than one physiological parameter of the subject to determine the physiological state of interest in the subject.


According to another aspect, there is provided a device for detecting a physiological state of interest of a subject, comprising: a source of electromagnetic radiation (EMR) that emits a plurality of wavelengths of EMR from about 400 nm to about 2500 nm, the source of EMR being operatively coupled to a power source; a receptor sized to register with, and fit against, a sample, the receptor comprising one or more than one port; one or more than one input radiation guiding element in operable association with the source of EMR, one or more than one output radiation guiding element in operable association with a detector, the one or more than one input radiation guiding element and the one or more than one output radiation guiding element in optical alignment with the one or more than one port located and defining an EMR path within the receptor when the receptor is registered with, and fit against, the sample; the detector for measuring transmitted or reflected EMR received from the sample, the detector operatively coupled to a processing system; the processing system comprising one or more than one algorithm for determining a concentration for two or more than two analytes in the sample, wherein the two or more than two analytes comprise multiple ghost analytes whose identities are unknown and that are observed to change in response to the state of interest, and using the one or more than one algorithm to derive the physiological state of interest of the sample, wherein, the physiological state of interest is: i) cannabis induced intoxication, and two or more than two analytes are selected from the group of: delta-9-tetrahydrocannabinol (THC), THC glucuronide (THCGlu), 11-nor-9-carboxy-THC (THC-COOH), 11-hydroxy THC (11-OH-THC), THC-COOH/11-OH-THC ratio, 11-nor-9-carboxy-THC glucuronide (THC-COOGlu), cannabidol (CBD), cannbinol (CBN), cannabigerol (CBG), delta-9-tetrahydrocannabivarin (THCV), THCV-carboxylic acid, 11-nor-9-carboxy-delta-tetrahydrocannabivarin (THCV-COOH), albumin, apolipoproteins A1 and B (apoA1 and apoB), total protein, bilirubin, prolactin, triglycerides, creatinine, cortisol, glucose, lactate, Total 4, uric acid, blood urea nitrogen (BUN), blood sugar, calcium, ionized calcium, magnesium, sodium, phosphate, and GABA; ii) alcohol induced intoxication, and two or more than two analytes are selected from the group of alcohol, aldehyde, lactic acid; or iii) intoxication arising from cannabis, alcohol, opiates, fentanyl, amphetamines, phencyclidine, sedatives, anxyolytics, cocaine, caffeine, and nicotine consumption, then two or more than two analytes may include:, then two or more than two analytes may include: delta-9-tetrahydrocannabinol (THC),


THC glucuronide (THCGlu), 11-nor-9-carboxy-THC (THC-COOH), 11-hydroxy THC (11-OH-THC), THC-COOH/11-OH-THC ratio, 11-nor-9-carboxy-THC glucuronide (THC-COOGlu), cannabidol (CBD), cannbinol (CBN), cannabigerol (CBG), delta-9-tetrahydrocannabivarin (THCV), THCV-carboxylic acid, 11-nor-9-carboxy-delta-tetrahydrocannabivarin (THCV-COOH), albumin, apolipoproteins A1 and B (apoA1 and apoB), total protein, bilirubin, prolactin, triglycerides, creatinine, cortisol, glucose, lactate, Total 4, uric acid, blood urea nitrogen (BUN), blood sugar, calcium, ionized calcium, magnesium, sodium, phosphate, GABA, alcohol, aldehyde, and lactic acid.


According to another aspect, there is provided a method comprising: directing light at a body part of a subject such that the light passes through or is reflected by blood and interstitial fluid of the body part, wherein the light incident on the body part comprises a range of wavelengths from at least one of the near infrared and visible spectra; measuring a spectrum of the light after the light has one or both of passed through and been reflected by the body part, wherein the spectrum comprises the range of wavelengths; comparing the measured spectrum against a reference spectrum representative of a known physiological state of interest; and determining whether the subject is in the known physiological state of interest from a similarity between the measured spectrum to the reference spectrum and without direct reference to analytes of the subj ect.


The light incident on the body part may comprise a range of wavelengths from both of the near infrared and visible spectra.


The light incident on the body part may comprise a range of wavelengths from the near infrared spectrum.


The light incident on the body part may comprise a range of wavelengths from the visible spectrum.


The spectrum may be measured on the light that has passed through the body part.


The spectrum may be measured on the light that has been reflected by the body part.


The spectrum may be measured on the light that has been through the body part and that has been reflected by the body part.


The measured spectrum may comprise a light reference sample, a dark reference sample, a light sample, and a dark sample, and the comparing may comprise correcting for sensor bias using the light reference sample, the dark reference sample, the light sample, and the dark sample.


The comparing may comprise removing outliers from the measured spectrum and generating a mean centered version of the measured spectrum.


The comparing may comprise applying a transform to the mean centered version of the measured spectrum.


The method may further comprise: applying multiple transforms to the mean centered version of the measured spectrum, wherein the transforms are selected from the group consisting of standard normal variate (SNV), multiplicative scatter correction (MSC), L1 normalization (L1N), L2 normalization (L2N), Savitzky-Golay smoothing (SGS), convolution smoothing (CS), and signal derivative (SD); evaluating performance of each of the multiple transforms to the mean centered version of the measured spectrum; and selecting, from a result of the evaluating, a transformed spectrum, wherein the transformed spectrum is a transformed version of the mean centered version of the measured spectrum.


The method may further comprise selecting at least one range of wavelengths that is a subset of a total wavelength range of the transformed spectrum.


The method may further comprise applying partial least squares to obtain PLS-derived components for the at least one range of wavelengths.


The determining may comprise applying a linear regression model to the fit.


The determining may comprise applying a neural additive model to the fit.


The determining may comprise applying a neural additive model to the at least one range of wavelengths that is the subset of the total wavelength range of the transformed spectrum.


The determining may comprise applying an artificial deep neural network to the PLS-derived components.


The determining may comprise applying an artificial deep neural network to the at least one range of wavelengths that is the subset of the total wavelength range of the transformed spectrum.


The determining may comprise applying a convolutional neural network to the at least one range of wavelengths that is the subset of the total wavelength range of the transformed spectrum.


The method may further comprise applying partial least squares to obtain PLS-derived components for the measured spectrum.


The determining may comprise applying a linear regression model to the fit. The determining may comprise applying a neural additive model to the fit.


The determining may comprise applying a neural additive model to the measured spectrum.


The determining may comprise applying an artificial deep neural network to the PLS-derived components.


The determining may comprise applying an artificial deep neural network to the measured spectrum.


The determining may comprise applying a convolutional neural network to the measured spectrum.


The determining may comprise receiving a sensitivity target and a specificity target, and outputting the physiological state of interest in accordance with the sensitivity and specificity targets.


The physiological state of interest may comprise whether the subject has COVID-19.


The measuring may be performed using a Fourier Transform Near Infrared spectrometer.


The spectrometer may comprise a platform for receiving a sample container, and the measuring may be performed directly on a finger of an individual.


According to another aspect, there is provided a method comprising: directing light at a body part of a subject such that the light passes through or is reflected by blood and interstitial fluid of the body part, wherein the light incident on the body part comprises a range of wavelengths from at least one of the near infrared and visible spectra; measuring a spectrum of the light after the light has one or both of passed through and been reflected by the body part, wherein the spectrum comprises the range of wavelengths; and determining whether the subject is in a physiological state of interest without direct reference to analytes of the subject, wherein the determining comprises using a trained machine learning model to process the measured spectrum and wherein the trained machine learning model is trained with reference spectra representative of the physiological state of interest.


The light incident on the body part may comprise a range of wavelengths from both of the near infrared and visible spectra.


The spectrum may be measured on the light that has passed through the body part.


The spectrum may be measured on the light that has been through the body part and that has been reflected by the body part.


The measured spectrum may comprise a light reference sample, a dark reference sample, a light sample of the subject, and a dark sample of the subject, and the comparing may comprise correcting for sensor bias using the light reference sample, the dark reference sample, the light sample of the subject, and the dark sample of the subject.


The method may further comprise, prior to using the trained machine learning model to process the measured spectrum, removing outliers from the measured spectrum and generating a mean centered version of the measured spectrum.


The method may further comprise, prior to using the trained machine learning model to process the measured spectrum: applying multiple transforms to the mean centered version of the measured spectrum, wherein the transforms are selected from the group consisting of standard normal variate (SNV), multiplicative scatter correction (MSC), L1 normalization (L1N), L2 normalization (L2N), Savitzky-Golay smoothing (SGS), convolution smoothing (CS), and signal derivative (SD); evaluating performance of each of the multiple transforms to the mean centered version of the measured spectrum; and selecting, from a result of the evaluating, a transformed spectrum, wherein the transformed spectrum is a transformed version of the mean centered version of the measured spectrum.


The method may further comprise selecting at least one range of wavelengths that is a subset of a total wavelength range of the transformed spectrum, and wherein the machine learning model is used to process the transformed spectrum.


The method may further comprise decomposing the transformed spectrum into latent space components, and processing the transformed spectrum using the trained machine learning model may comprise processing the latent space components using respective instances of the machine learning model.


The machine learning model comprise any one or more of a neural additive model, an artificial deep neural network, and a convolutional neural network, for example.


The method may further comprise decomposing the transformed spectrum into latent space components, and processing the measured spectrum using the trained machine learning model may comprise processing the latent space components using respective instances of the machine learning model.


The latent space components may be generated by applying partial least squares or a principal components analysis.


The determining may comprise receiving a sensitivity target and a specificity target, and outputting the physiological state of interest in accordance with the sensitivity and specificity targets.


The physiological state of interest may comprise whether the subject is infected with a virus.


The physiological state of interest may comprise whether the subject has COVID-19.


The physiological state of interest may comprise cannabis or THC impairment. The physiological state of interest may comprise alcohol impairment.


The measuring may be performed using a Fourier Transform Near Infrared spectrometer.


The spectrometer may comprise a platform for receiving a sample container, and wherein the measuring may be performed directly on a finger of an individual.


According to another aspect, there is provided a non-transitory computer readable medium having stored thereon computer program code that is executable by a processor and that, when executed by the processor, causes the processor to perform the method of any of the foregoing aspects and suitable combinations thereof


According to another aspect, there is provided an apparatus comprising: a lamp; at least one spectrometer; an interface comprising a receiver for a body part of a subject; and at least one source fiber and at least one return fiber optically coupling the lamp and the at least one spectrometer to the receiver, wherein the source fiber is positioned to direct light from the lamp to the body part and the return fiber is positioned to receive light transmitted through or reflected by the body part and direct the received light to the at least one spectrometer.


The body part may be a finger and the receiver may comprise a first surface positioned to abut against a pad of the finger and a second surface positioned to abut against a tip of the finger, wherein the at least one source fiber is positioned to direct the light from the lamp to the finger through the second surface and the at least one return fiber is positioned to receive the light transmitted through the body part through the first surface.


The body part may be a finger and the receiver may comprise a first surface positioned to abut against a pad of the finger, wherein the at least one source fiber and the at least one return fiber are positioned to respectively direct the light from the lamp to the finger and to receive the light reflected by the finger.


The at least one spectrometer may comprise a first spectrometer configured to output light in the visible spectrum and a second spectrometer configured to output light in the near infrared spectrum.


The first spectrometer may be further configured to output light in the near infrared spectrum.


The first spectrometer may be configured to output light having wavelengths from approximately 350 nm to 1,000 nm, and the second spectrometer is configured to output light having wavelengths from about 900 nm to about 2,500 nm.


The apparatus may further comprise: a communication port; and a controller communicatively coupled to the at least one spectrometer and the communication port, wherein the controller is configured to: direct light at the body part from the lamp, wherein the light comprises a range of wavelengths from at least one of the near infrared and visible spectra; measure, using the at least one spectrometer, a spectrum of the light after the light has one or both of passed through and been reflected by the body part, wherein the spectrum comprises the range of wavelengths; output the measured spectrum to the communication port.


The apparatus may further comprise: a processor communicatively coupled to the communication port; and a non-transitory computer readable medium having stored thereon computer program code that is executable by the processor and that, when executed by the processor, causes the processor to: 1) compare the measured spectrum against a reference spectrum representative of a known physiological state of interest; and determine whether the subject is in the known physiological state of interest from a similarity between the measured spectrum to the reference spectrum and without direct reference to analytes of the subject, or 2) perform a method comprising determining whether the subject is in a physiological state of interest without direct reference to analytes of the subject, wherein the determining comprises using a trained machine learning model to process the measured spectrum and wherein the trained machine learning model is trained with reference spectra representative of the physiological state of interest.


The at least one spectrometer may comprise a Fourier Transform Near Infrared spectrometer.


According to another aspect, there is provided a method of non-invasively determining a physiological state of interest in a subject comprising: placing a body part in contact with a receptor; directing a source of electromagnetic radiation (EMR) over a range of wavelengths through the receptor and onto body part so that the EMR reaches the blood and interstitial fluid within the body part; measuring the EMR absorbed by, reflected by, or transmitted through, the blood and interstitial fluid of the body part with a detector to obtain a measured spectrum over the range of wavelengths; and comparing the measured spectrum against a reference spectrum to determine the physiological state of interest in the subject.


The measuring may be performed using a Fourier Transform Near Infrared spectrometer.


The spectrometer may comprise a platform for receiving a sample container, and the measuring may be performed directly on a finger of an individual.


According to another aspect, there is provided a method of non-invasively determining a physiological state of interest of a subject comprising, determining one or more than one physiological parameter of the subject; placing a body part in contact with a receptor; directing a source of electromagnetic radiation (EMR) over a range of wavelengths through the receptor and onto body part so that the EMR reaches the blood and interstitial fluid within the body part; measuring the EMR absorbed by, reflected by, or transmitted through, the blood and interstitial fluid of the body part with a detector to obtain a measured spectrum over the range of wavelengths; and comparing the measured spectrum against a reference spectrum to determine the physiological state of interest in the subject.


The measuring may be performed using a Fourier Transform Near Infrared spectrometer.


The spectrometer may comprise a platform for receiving a sample container, and wherein the measuring may be performed directly on a finger of an individual.


According to another aspect, there is provided a device for detecting a physiological state of interest of a subject, comprising: a source of electromagnetic radiation (EMR) that emits a plurality of wavelengths of EMR from about 350 nm to about 2500 nm, the source of EMR being operatively coupled to a power source; a receptor sized to register with, and fit against, a sample, the receptor comprising one or more than one port; one or more than one input radiation guiding element in operable association with the source of EMR, one or more than one output radiation guiding element in operable association with a detector; the one or more than one input radiation guiding element and the one or more than one output radiation guiding element in optical alignment with the one or more than one port located and defining an EMR path within the receptor when the receptor is registered with, and fit against, the sample; the detector for measuring transmitted or reflected EMR received from the sample, the detector operatively coupled to a processing system, the processing system configured to compare a measured spectrum of the transmitted or reflected EMR of the sample against a reference spectrum to derive the physiological state of interest of the sample, wherein, the physiological state of interest is whether the subject has COVID-19.


The detector may comprise a Fourier Transform Near Infrared spectrometer.


This summary does not necessarily describe the entire scope of all aspects. Other aspects, features and advantages will be apparent to those of ordinary skill in the art upon review of the following description of specific embodiments.





BRIEF DESCRIPTION OF THE DRAWINGS

These and other features of the various embodiments will become more apparent from the following description in which reference is made to the appended drawings wherein:



FIG. 1 shows an example of a device in accordance with an example embodiment.



FIG. 2 is prior art (U.S. Pat. No. 6,657,717) and shows the absorbance spectra for globulins, glucose, urea, creatinine, cholesterol and HAS over a range of wavelengths from 500 nm to 1400 nm.



FIG. 3 shows a schematic view of a device according to an example embodiment.



FIG. 4 shows a schematic view of a device according to an example embodiment.



FIG. 5 shows a schematic view of a development pipeline according to an example embodiment.



FIG. 6 shows a schematic of a system for determine state of interest of a subject, according to an example embodiment.



FIG. 7 shows a block diagram illustrating information flow in the system of FIG. 6.



FIGS. 8A and 8B are perspective and sectional views along line 8A-8A, respectively, of a reflectance interface.



FIGS. 9A and 9B are perspective and sectional views along line 9A-9A, respectively, of a transmittance interface.



FIG. 10 is a flow diagram of a method of processing spectral data obtained using the interface of FIGS. 8A and 8B or FIGS. 9A and 9B, according to an example embodiment.



FIG. 11 is a flow diagram of a method of input data processing, comprising part of the method of FIG. 10.



FIG. 12 is a flow diagram of a method of performing outlier detection, comprising part of the method of FIG. 10.



FIG. 13 is a flow diagram of a method of executing an optimization pipeline, comprising part of the method of FIG. 10.



FIG. 14 is a flow diagram depicting application of various models to determine state of interest.



FIG. 15 schematically represents the effect of the preprocessing pipeline of FIG. 10.



FIG. 16 is a flow diagram of a method of executing the optimization pipeline of FIG. 13.



FIG. 17 is a block diagram depicting a configuration block shown in FIG. 16.



FIG. 18 is an example of wavelength reduction performed as part of executing the optimization pipeline of FIG. 10.



FIGS. 19A-C and 20A-C are example models that may be applied to process spectral data to arrive at a state of interest.



FIGS. 21A-C depict a block diagram of an example system that may be used to process spectral data to arrive at a state of interest.



FIG. 22 is an example process diagram in which spectral data is processed using a state model to arrive at a state of interest.



FIG. 23 is an example combined device comprising an interface for obtaining spectroscopic measurements and a spectrometer for processing those measurements.





DETAILED DESCRIPTION

The present disclosure relates to non-invasive methods and a device for determining a physiological state of interest in a subject. The methods are used to determine biochemical profile or fingerprint, that is indicative of the physiological state of interest in the subject.


As used herein, the terms “comprising,” “having,” “including” and “containing,” and grammatical variations thereof, are inclusive or open-ended and do not exclude additional, un-recited elements and/or method steps. The term “consisting essentially of” when used herein in connection with a use or method, denotes that additional elements and/or method steps may be present, but that these additions do not materially affect the manner in which the recited method or use functions. The term “consisting of” when used herein in connection with a use or method, excludes the presence of additional elements and/or method steps. A use or method described herein as comprising certain elements and/or steps may also, in certain embodiments consist essentially of those elements and/or steps, and in other embodiments consist of those elements and/or steps, whether or not these embodiments are specifically referred to. In addition, the use of the singular includes the plural, and “or” means “and/or” unless otherwise stated. The term “plurality” as used herein means more than one, for example, two or more, three or more, four or more, and the like. Unless otherwise defined herein, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art. As used herein, the term “about” refers to an approximately +/−10% variation from a given value. It is to be understood that such a variation is always included in any given value provided herein, whether or not it is specifically referred to. The use of the word “a” or “an” when used herein in conjunction with the term “comprising” may mean “one,” but it is also consistent with the meaning of “one or more,” “at least one” and “one or more than one.”


The expression “body part” or “part of a subject”, as used herein, refers to an element or section of a human to which electromagnetic radiation (EMR) can be directed. The element or section can be, for example, an earlobe, a finger, an arm, a leg, torso, cheek, or a toe.


Described herein are non-invasive methods for determining a physiological state of interest in a subject. For example, the method may involve placing a body part in contact with a receptor of a device. Directing a source of electromagnetic radiation (EMR) over a range of wavelengths, for example from about 350 to about 2500 nm, through the receptor and onto body part so that the EMR reaches the blood and interstitial fluid within the body part, and measuring the EMR absorbed by, reflected by, or transmitted through, the blood and interstitial fluid of the body part with a detector to obtain a spectrum over the range of wavelengths. A quantitative mathematical analysis of the spectrum is performed using an algorithm that determines an amount of two or more than two analytes within the blood and interstitial fluid of the body part. The amount (for example, the concentration or level) of the two or more than two analytes are used to derive a biochemical profile, or fingerprint of the subject.


The amount of the two or more than two analytes (the biochemical fingerprint also termed biochemical profile) may be compared against reference values of the two or more analytes obtained from reference (control) subjects, for example, from a cross section of healthy individuals in order to provide data that may be used to indicate the status of a physiological state of interest of the subj ect. Additionally, the biochemical fingerprint (biochemical profile) may be monitored over time and compared against previous values of the two or more analytes obtained from the same subject in order to obtain data that may be used to manage the physiological state of interest, of the subject.


The biochemical profile may also be combined with physiological parameters to determine the physiological state of interest, of a subject. The “state of being” may include a physiological state of interest for example as a result of intoxication, for example but not limited to intoxication arising from cannabis consumption, alcohol consumption, both cannabis and alcohol consumption, or from the consumption of other intoxicants, for example but not limited to consumption of opiates, fentanyl, amphetamines, phencyclidine, sedatives, anxyolytics, cocaine, caffeine, and nicotine.


The device as described herein may be used to determine if the physiological state of interest of a subj ect is indicative of a state of intoxication. If a state of intoxication is determined, then corrective action may be taken directly, or the result of the test may be forwarded to another party so that corrective action may be taken by the other party. For example, if the device is used for road-side testing, and the operator of the device (and optionally delivering the physiological and behavioral tests) is a law enforcement officer and a result is obtained that is indicative of a state of intoxication for the driver or a car, then the operator of the test, for example the law enforcement officer, may perform corrective action and confiscate the car, suspend the driver's license, press charges and the like. Alternatively, the operator of the test, a health care worker, or the law enforcement officer, may forward the positive result indicating intoxication (impairment) to a third party, for example a justice of the peace, and corrective action may be taken. In some circumstances, for example, where safety is a requirement of the subjects employment, if the subject has been determined to exhibit a positive result indicating intoxication (impairment), then the result may be forwarded to the subject's employer. Examples of situations where safety may be a requirement of the subjects employment, include if the subject is working as an air traffic controller, the subject is a pilot, they operate a commercial vehicle, they operate machinery (large or small) at a work site, they are an operator at a nuclear power facility etc.


By consumption it is meant that the intoxicant, for example, but not limited to, cannabis, alcohol, opiates, fentanyl, amphetamines, phencyclidine, sedatives, anxyolytics, cocaine, caffeine, nicotine, enters, is taken, or is administered to, or by, the subject orally (for example as an edible product), by inhalation (for example, via smoking, an e-cigarette, via a hookah, via an oral spray, snorting, or using an aspirator or inhaler), by transdermal delivery (for example, via a patch, cream, spray, or oil), or intravenously (for example by injection or as a drip solution), or other method.


The physiological state of interest may be detected or defined on the basis of two or more than two metabolites (also termed analytes) and optionally, one or more physiological parameter, or one or more behavioral parameter, or the physiological state of interest may be defined using two or more than two analytes, and both the one or more physiological parameter, and the one or more behavioral parameter. While a state of being may be described using a plurality of metabolites (analytes), these results may be combined with physiological parameter(s), and/or behavioral parameters, to further define the state of being. The metabolites may be determined to have, or not to have, an interdependent role in the state of being. The metabolites, the physiological parameters, or both the metabolites, the physiological parameters may be termed “variables” and these variables may be used to describe the state of being, or these variables considered as markers representative of the state of being. By using multiple metabolites combined with other physiological variables that are correlated in some manner to the physiological state of interest, a more accurate determination of the physiological state of interest may be obtained when compared to determining the physiological state of interest determined using one analyte.


The methods described herein permit management of multiple variables or markers that are correlated with, and that may have an influence on, the physiological state of interest, and therefore can be used to characterize and determine the physiological state of interest of a subject or patient.


A fingerprint indicative of a physiological state resulting from intoxication may be determined using a plurality of metabolites, and one or more physiological parameter, and/or one or more behavioral parameter. These variables (metabolites, and physiological parameters and/or behavioral parameters), may be compared to base-line values that have been determined from healthy individuals, and any deviation from the base-line values is indicative of the physiological state of interest, of the subject. These metabolites, and physiological parameters and/or behavioral parameters may be also compared to values determined overtime from the same subject, and any deviation from the time-course values of these variable, (or markers) may indicates a change in the physiological state of interest in the subject. Therefore, the measured biochemical fingerprint, and physiological parameters and/or behavioral parameters, may be used to monitor, and manage, the physiological state of interest over time.


Non-limiting examples of a physiological state of interest in the subject may include but are not limited to: intoxication, for example as a result of cannabis consumption, alcohol consumption, both cannabis and alcohol consumption, or from the consumption of other intoxicants, for example but not limited to consumption of opiates, fentanyl, amphetamines, phencyclidine, sedatives, anxyolytics, cocaine, caffeine, and nicotine. Non-limiting examples of metabolites that may be determined to obtain a biochemical fingerprint of the corresponding physiological state of interest include (also see Table 1):

    • i) intoxication (cannabis induced): the two or more than two analytes may include: delta-9-tetrahydrocannabinol (THC), THC glucuronide (THCGlu), 11-nor-9-carboxy-THC (THC-COOH), 11-hydroxy THC (11-OH-THC), THC-COOH/11-OH-THC ratio, 11-nor-9-carboxy-THC glucuronide (THC-COOGlu), cannabidol (CBD), cannbinol (CBN), cannabigerol (CBG), delta-9-tetrahydrocannabivarin (THCV), THCV-carboxylic acid, 11-nor-9-carboxy-delta-tetrahydrocannabivarin (THCV-COOH), albumin, apolipoproteins A1 and B (apoA1 and apoB),), total protein, bilirubin, prolactin, triglycerides, creatinine, cortisol, glucose, lactate, Total 4, uric acid, blood urea nitrogen (BUN), blood sugar, calcium, ionized calcium, magnesium, sodium, phosphate, and gamma-aminobutyric acid (GABA);
    • ii) alcohol induced intoxication, then two or more than two analytes may include: alcohol, aldehyde, and lactic acid; or
    • iii) intoxication generally, for example arising from a combination of cannabis and alcohol, or from opiates, fentanyl, amphetamines, phencyclidine, sedatives, anxyolytics, cocaine, caffeine, and nicotine consumption, then two or more than two analytes may include:, then two or more than two analytes may include: delta-9-tetrahydrocannabinol (THC), THC glucuronide (THCGlu), 11-nor-9-carboxy-THC (THC-COOH), 11-hydroxy THC (11-OH-THC), THC-COOH/11-OH-THC ratio, 11-nor-9-carboxy-THC glucuronide (THC-COOGlu), cannabidol (CBD), cannbinol (CBN), cannabigerol (CBG), delta-9-tetrahydrocannabivarin (THCV), THCV-carboxylic acid, 11-nor-9-carboxy-delta-tetrahydrocannabivarin (THCV-COOH), albumin, apolipoproteins A1 and B (apoA1 and apoB), total protein, bilirubin, prolactin, triglycerides, creatinine, cortisol, glucose, lactate, Total 4, uric acid, blood urea nitrogen (BUN), blood sugar, calcium, ionized calcium, magnesium, sodium, phosphate, GABA, alcohol, aldehyde, and lactic acid.









TABLE 1







Non-limiting examples of analytes and physiological parameters


that may be used to determine intoxication as described herein










Description
Influence













C-reactive Protein
C-reactive protein (CRP) is an acute phase protein that
High



increases in the blood with inflammation and infection



such as following a heart attack, surgery, or trauma.



CRP may be used to indicate a change in the level of



stress in a subject.


Creatinine (IDMS)
Measures the level of creatinine in the blood.
Medium



Creatinine is a waste product that forms when



creatine (found in muscle), breaks down. For



example: dehydration. This analyte may be used as a



measurement of kidney function.


Glucose
A blood glucose test measures the amount of glucose
High



in your blood plasma. Cannabis use can reduce the



amount of glucose in blood.


BUN
A blood urea nitrogen (BUN) test measures the amount
Low



of nitrogen in your blood that comes from the waste



product urea. Urea is made when protein is broken



down in your body. May be used to assess kidney



function.


THC
THC
High


Total Protein
A total serum protein test measures the total amount
Low



of protein in the blood. It also measures the amounts



of two major groups of proteins in the blood: albumin



and globulin.


Albumin
A plasma binding protein synthesized by the liver.
High



Albumin helps to maintain osmotic pressure in the



vascular space and also reflects overall nutritional



status. It is also assists in transport of various



substances throughout your body, including hormones,



vitamins, and enzymes.


Prolactin
A prolactin (PRL) test measures how much of a
High



hormone called prolactin you have in your blood. The



hormone is made in your pituitary gland, which is



seated at the base of the brain. It is known to be



affected by physiological stress, which is known to be



a factor in THC impairment.


Potassium
It is mostly found intracellularly, and less so
Low



interstitially and in serum. It can be used to diagnose



and monitor kidney disease, high blood pressure, and



heart disease.


Sodium
Is a cation found mainly in extracellular fluid and may be
High



used to evaluate sate of hydration. Sodium, along with



other electrolytes such as potassium, chloride, and



bicarbonate (or total CO2) may be used to evaluate



metabolic acidosis. Excess alcohol ingestion is known to



cause ketoacidosis. Sensitive but not specific.


Cortisol
A steroid hormone produced by the adrenal cortex. It
High



is used to evaluate pituitary or adrenal function. It is



known to be effected by stress. Sensitive but not



specific.


Lactate
An intermediate breakdown product of glucose
Low



metabolism primarily from anaerobic metabolism in



muscle. Lactic acid levels increases as a result of



strenuous exercise, heart failure, a severe infection



(sepsis), or shock. Lactic acidosis results when there



is oxygen deprivation in the tissues. Lactate may be



used as an indirect estimation of oxygenation, and to



evaluate metabolic acidosis. Acute phase reaction.



Sensitive and potentially specific.


Total T4
A hormone secreted by the thyroid gland. Total T4 is
High



converted into another thyroid hormone (T3;



triiodothyronine). Any change show up in T4 first. T3



and T4 help to control how body stores and uses



energy (metabolism). Sensitive but not specific.


Calcium, ionized
The sum of calcium plus protein bound calcium,
High


calcium
ionized calcium. It is important in cellular transport



mechanisms. The most common cause for low



calcium, ionized calcium is low albumin/protein.



Sensitive but not specific.


Uric Acid
Uric acid is an end product of purine metabolism. High
High



levels can be associated with gout and hypothyroidism.



Cannabis is known to lower the level of uric acid.



Sensitive and specific.


Triglyceride
Triglycerides are a dominant form of fat in the body.
High



They are insoluble in blood. Their levels may be



elevated with high alcohol consumption and immediate



high consumption of carbohydrates, such as junk foods.



Sensitive and potentially specific.


Magnesium
An important cation involved in cellular and bone
High



metabolism. It's needed for proper muscle, nerve, and



enzyme function. It also helps the body make and use



energy. Alcohol abuse lowers the level. Sensitive and



potentially specific.


Creatine Kinase
An enzyme found in the cardiac and skeletal muscle,
High



brain and lung. Levels of CK can rise after a heart



attack, skeletal muscle injury, strenuous exercise, and



from taking certain medicines or supplements.



Sensitive and potentially specific.


GGT
Gamma-glutamyl transferase (GGT) is an enzyme
High



that is found in many organs throughout the body,



with the highest concentrations found in the liver. It



is sensitive to acute alcohol ingestion. Smoking may



cause elevated GGT levels. Acute phase reactor.



Highly sensitive and potentially specific.


AST
The aspartate aminotransferase (AST) is mainly found
Low



in the liver but it is also found in muscle and kidney. It



is involved in amino acid metabolism. It is used to



assess liver damage, such as with alcohol abuse. Slow



responder.


Total Bilirubin
A breakdown product of hemoglobin. It is used to assess liver
Low



function. Sensitive but not specific.


WBC Count
The number of white blood cells (WBC) per ml of blood.
High



WBC's primarily consist of neutrophils, lymphocytes,



monocytes, eosinophils, and basophils. WBC's are



mobilized by inflammation and infection. Sensitive but



not specific.


RBC Count
number of red blood cells (RBC) per cubic ml of whole
Low



blood. This measure is lower with alcohol abuse.


Hemoglobin
Oxygen carrying protein found in RBC's. It is lowered by
High



alcohol abuse. Sensitive but not specific.


Hematocrit
The hematocrit blood test determines the percentage
Medium



of blood volume composed of red blood cells (RBC's).



Used to diagnose anemia.


Neutrophils
A type of granulocytic WBC. They are involved in
High



inflammatory reactions and function as a primary



defence against acute infections. Stress and smoking



can elevate levels of neutrophils. Sensitive and not



specific.


Lymphocytes
Small white blood cells consisting of B lymphocytes that
Medium



control antibody response and T lymphocytes that



control cell mediated response. Sensitive and not



specific.


Eosinophils
Subclass of WBC granulocytes. Involved in allergic
Medium



reactions and in attacking parasites.









Another example of a method for non-invasively determining a physiological state of interest of a subject involves determining the biochemical profile (or fingerprint) as outlined above, along with determining one or more than one physiological parameter of the subject and using the biochemical profile and one or more than one physiological parameter to determine the physiological state of interest, of the subject. In this method, the biochemical profile is determined by placing a body part in contact with a receptor and directing a source of electromagnetic radiation (EMR) over a range of wavelengths, for example from about 350 to about 2500 nm, through the receptor and onto body part so that the EMR reaches the blood and interstitial fluid within the body part. The EMR that is absorbed by, reflected by, or transmitted through, the blood and interstitial fluid of the body part is measured with a detector in order to obtain a spectrum over the range of wavelengths, and a quantitative mathematical analysis of the spectrum is performed using an algorithm to determine an amount of two or more than two analytes within the blood and interstitial fluid of the body part. The amount of the two or more than two analytes are used to derive a biochemical profile, which may be compared against reference values of the two or more analytes. The biochemical profile and the one or more than one physiological parameter may be used to determine the state of being, physiological state of interest for example intoxication as described herein.


In addition to determining a biochemical profile and one or more than one physiological parameter, one or more behavioral parameters may also be determined to characterize the physiological state of interest. Non-limiting examples of behavioral parameters that may be evaluated include, determination of mental acuity (for example but not limited to a name-face test, a fire alarm test, a two delayed recall tests, a misplaced objects test, a shopping list test, a digit symbol test), one or more motor skill test (for example but not limited to, a walk and turn test, a one leg stand test, a horizontal gaze nystagmus test, a divided attention test, a rhomberg balance test), the ability to function at a defined task, for example to operate machinery, drive an automobile (or use a driving simulator), state of physical fitness, standardized field sobriety (Newmeyer, Swortwood, Taylor, et al., 2017, Clin Chem, 63(3), 647-662. doi:10.1373/clinchem.2016.265371), and the like.


For example, a driving simulator, is known to be useful in assessing on-the-road driving tests (Micallef et al., 2018, Fundam Clin Pharmacol. doi:10.1111/fcp.12382) and the simulator may be used as a test in place of driving an automobile). Studies with standardized and objective measures of driving using driving simulators have found impairments in driving following moderate cannabis intake (e.g. 8% THC; in a 500-750 mg cigarette, approximately 40-60 mg of THC). Reported effects include increased weaving, decreased speed, decreased steering control, longer reaction time and, increased headway (Hartman et al., 2015, Drug Alcohol Depend, 154, 25-37. doi :10.1016/j .drugalcdep.2015.06.015; Lenne et al., 2010, Accid Anal Prey, 42(3), 859-866. doi :10.1016/j .aap.2009.04.021; Micallef et al., 2018, Fundam Clin Pharmacol. doi:10.1111/fcp.12382; Anderson, et al., 2010, J Psychoactive Drugs, 42(1), 19-30. doi:10.1080/02791072.2010.10399782; Ronen et al., 2010, Accid Anal Prey, 42(6), 1855-1865. doi:10.1016/j.aap.2010.05.006; Ronen et al., 2008, Accid Anal Prey, 40(3), 926-934. doi :10.1016/j .aap.2007.10.011).


In the case where there is an interrogation of a physiological state of interest, a behavioral parameter data may also be combined with the biochemical profile and physiological parameters to further assist in determining the degree or status of the physiological state of interest. For example, the data (the tested values) may be combined to determine if the subject has surpassed a threshold value, index, ratio, or set of values, and/or, the degree, or the extent to which the subject is exhibiting the physiological state of interest, for example the degree of intoxication. The base-line values that are used to determine the index, ratio, or set of values, against which the tested values are compared, are determined from normalized healthy subjects, analyzed under control conditions.


Non-limiting examples of a physiological state of interest include intoxication arising from cannabis, alcohol, opiates, fentanyl, amphetamines, alcohol, phencyclidine, sedatives, anxyolytics, cocaine; caffeine-induced disorders; and nicotine-induced disorders. For example:

    • i) if the intoxication is cannabis-induced, then the two or more than two analytes may include: delta-9-tetrahydrocannabinol (THC), THC glucuronide (THCGlu), 11-nor-9-carboxy-THC (THC-C 0 OH), 11-hydroxy THC (11-OH-THC), THC-COOH/11-OH-THC ratio, 11-nor-9-carboxy-THC glucuronide (THC-COOGlu), cannabidol (CBD), cannbinol (CBN), cannabigerol (CBG), delta-9-tetrahydrocannabivarin (THCV), THCV-carboxylic acid, 11-nor-9-carboxy-delta-tetrahydrocannabivarin (THCV-COOH), albumin, apolipoproteins A1 and B (apoA1 and apoB), total protein, bilirubin, prolactin, triglycerides, creatinine, cortisol, glucose, lactate, Total 4, uric acid, blood urea nitrogen (BUN), blood sugar, calcium, ionized calcium, magnesium, sodium, phosphate, and GABA, the physiological parameter may include one or more of: heart rate, pulse rate, body temperature, neuropeptide Y, fatty acid amide hydrolase (FAAH), c reactive protein (cRP), creatine kinase (CK), aspartate amino transferase (AAT), asparate aminotransferase (AST), alanine transaminase (ALT), gamma-glutamyl transpeptidase (GGT), white blood cell count (WBC), red blood cell count (RBC), hemoglobin, hematocrit, neutrophils, lymphocytes, eosinophils, hypoactivity; THC in hair, and THC in urine, and the one or more behavioral parameters may include determination of mental acuity (a name-face test, a fire alarm test, a two delayed recall tests, a misplaced objects test, a shopping list test, a digit symbol test), one or more motor skill test (a walk and turn test, a one leg stand test, a horizontal gaze nystagmus test, a divided attention test, a rhomberg balance test), the ability to function at a defined task, to operate machinery, drive an automobile, standardized field sobriety;
    • ii) for alcohol-induced intoxication, then two or more than two analytes may include: alcohol, aldehyde, lactic acid, the physiological parameter may include: heart rate, pulse rate, body temperature, neuropeptide Y, aspartate amino transferase (AAT), alanine transaminase (ALT), gamma-glutamyl transpeptidase (GGT), and the one or more behavioral parameters may include determination of mental acuity (a name-face test, a fire alarm test, a two delayed recall tests, a misplaced objects test, a shopping list test, a digit symbol test), one or more motor skill test (a walk and turn test, a one leg stand test, a horizontal gaze nystagmus test, a divided attention test, a rhomberg balance test), the ability to function at a defined task, to operate machinery, drive an automobile, standardized field sobriety; or
    • iii) for intoxication, for example resulting from cannabis and alcohol, or opiates, fentanyl, amphetamines, phencyclidine, sedatives, anxyolytics, cocaine, caffeine, and nicotine consumption, then two or more than two analytes may include: then two or more than two analytes may include: delta-9-tetrahydrocannabinol (THC), THC glucuronide (THCGlu), 11-nor-9-carboxy-THC (THC-COOH), 11-hydroxy THC (11-OH-THC), THC-COOH/11-OH-THC ratio, 11-nor-9-carboxy-THC glucuronide (THC-COOGlu), cannabidol (CBD), cannbinol (CBN), cannabigerol (CBG), delta-9-tetrahydrocannabivarin (THCV), THCV-carboxylic acid, 11-nor-9-carboxy-delta-tetrahydrocannabivarin (THCV-COOH), albumin, apolipoproteins A1 and B (apoA1 and apoB), total protein, bilirubin, prolactin, triglycerides, creatinine, cortisol, glucose, lactate, Total 4, uric acid, blood urea nitrogen (BUN), blood sugar, calcium, ionized calcium, magnesium, sodium, phosphate, GABA, alcohol, aldehyde, and lactic acid, the physiological parameter may include one or more of: heart rate, pulse rate, body temperature, neuropeptide Y, fatty acid amide hydrolase (FAAH), c reactive protein (cRP), creatine kinase (CK), aspartate amino transferase (AAT), asparate aminotransferase (AST), alanine transaminase (ALT), gamma-glutamyl transpeptidase (GGT), white blood cell count (WBC), red blood cell count (RBC), hemoglobin, hematocrit, neutrophils, lymphocytes, eosinophils, hypoactivity; THC in hair, THC in urine, and the one or more behavioral parameters may include determination of mental acuity (a name-face test, a fire alarm test, a two delayed recall tests, a misplaced objects test, a shopping list test, a digit symbol test), one or more motor skill test (a walk and turn test, a one leg stand test, a horizontal gaze nystagmus test, a divided attention test, a rhomberg balance test), the ability to function at a defined task, to operate machinery, drive an automobile, standardized field sobriety (Newmeyer, Swortwood, Taylor, et al., 2017, Clin Chem, 63(3), 647-662. doi:10.1373/clinchem.2016.265371).


After the data is collected by the device as described herein, the result can be stored, shown on a display, or transmitted to another central CPU for further analysis or display. For example, the data may be transmitted to a central computer that analyses the measured metabolites, in combination with the measured physiological parameters and if desired measured behavioral parameters, in order to determine an overall index of the physiological state for the subject (patient). For example, the collected data may be compared with known data sets previously obtained for the physiological state of interest, and the “index value” determined. The index value may range for example from 0 (the patent is in a sever intoxicated state for the physiological state of interest being analyzed) to 1 (the patient is in an normal or healthy state for the physiological state of interest).


With reference to FIG. 3, the apparatus or device 100 that may be used in the methods described herein, comprises a receptor (also termed primary receptor) 10 shaped so that it can be placed in contact with a region of skin or a body part from a subject 20. A source electromagnetic radiation (EMR; 30) is directed 40 into the receptor 10, and following interaction with two or more than two compounds within the body part 20, the EMR is collected 50 and analyzed. The apparatus 100 may be as shown in FIG. 1, or it may be based on an apparatus as known in the art, for example, but not limited to those disclosed in US 2013/0248695, U.S. Pat. Nos. 5,361758, 5,429,128, WO 93/16629, U.S. Pat. Nos. 6,236,047 6,040,578 or 6,240,306 (all of which are incorporated herein by reference). The EMR 50 that is collected after interaction with compounds within the body part of the subject 20 may be either reflected from, transmitted through, absorbed by, or a combination thereof, the body part of the subject 20 depending upon the apparatus used. The collected EMR signal is directed to a spectrometer 60 and the data processed 70 using one or more than one calibration algorithms to determine the concentration of two, or more than two target compounds within the body part 20, to derive a biochemical profile, and determine the physiological condition or status of the subject. If desired, the data may be transferred wirelessly or by wire 80 to another device 90, for example a cell phone, or an off-site CPU, that comprises a program that can collect the data, display the results or a combination thereof. A non-limiting example of spectra of a range of compounds in blood, and measured non-invasively, is shown in FIG. 2.


The apparatus 100 as shown in FIG. 1 may further comprise a strap attached to the device. The strap may assist in maintaining control of the device while the test is being administered.


In order to minimize the effects of scatter associated with incorrect pressure applied to the receptor during measurement, the device 100 may comprise a pressure sensor 15 positioned around, under, or adjacent to, receptor (primary receptor) 10. If present, the pressure sensor is used to determine if there is too little, or too much, pressure between the skin surface that is placed against the receptor, and the receptor itself. For example, if there is too little or too much pressure against receptor 10 a signal from the pressure sensor 15 may be used to illuminate the bottom section of the device thereby signaling the need for correction and re-adjustment of the body part against the receptor.


The device 100 may comprise more than one receptor. For example, which is not to be considered limiting, when a body part 20, for example an index finger, is being measured by device 100, a second body part, for example the palm of the hand may be positioned over the device so that the second body part may press against the device while the subject is having their first body part measured. In this example, as shown in FIG. 4, device 100 may comprise primary receptor 10 and a secondary receptor 25 located at a location where the second body part may press against the secondary receptor 25.


The secondary receptor 25 may be configured in a similar manner as that of primary receptor 10, so that the EMR 50′ collected after interaction with compounds within the second body part of the subject that are either reflected from, transmitted through, absorbed by, or a combination thereof, is used to determine the presence of the same or another analyte. In the example shown in FIG. 4, the source of EMR 30 is shown as being the same for both the primary 15 and the secondary 25, receptor. In this example, the EMR source may emit a range of wavelengths that may be used to determine the occurrence of one or several analytes within the first and second body parts. For example, if the physiological state of being is intoxication, a first measured analyte, or group of first measured analytes, may be those related to cannabis-induced intoxication as described above, while a second analyte or a second set of analytes may be those related to alcohol-induced intoxication, as described above. If desired, the set of wavelengths emitted by the EMR 30 may also be used to determine the presence of a third analyte or set of analytes, associated with another physiological state of being, for example, opioid-induced intoxication. The third analyte or set of analytes may be interrogated using a third set of wavelengths, the range of these wavelengths may overlap, or be distinct from, the first and/or second set of wavelengths. The third set of wavelengths may be directed to either the first receptor, or the second receptor, or the third set of wavelengths may be directed to a third receptor that is placed within device 100 so that the third receptor contacts a third body part when the body part, for example a hand, is placed on device 100.


Alternatively, separate sources of EMR may be directed to a first receptor and a second receptor in the device 100. For example, one source of EMR, for example a first LED, emitting wavelengths of EMR within a first range of wavelengths, may be directed via optic fiber 40 to primary receptor 15, and a second source of EMR, for example a second LED that emits wavelengths of EMR within a second range of wavelengths, may be directed via optic fiber 40′ to secondary receptor 25. In this configuration, the first set of wavelengths and the second set of wavelengths are different, and the range of wavelengths may overlap, or the range of wavelengths may be distinct, and the device 100 may determine the presence of different analytes within each of the first and second body parts. For example, if the physiological state of being is intoxication, a first measured analyte, or group of first measured analytes, may be those related to cannabis-induced intoxication as described above, while a second analyte or a second set of analytes may be those related to alcohol-induced intoxication, as described above. In a similar manner as noted above, if desired, a third set of wavelengths emitted by the EMR 30 may be used to determine the presence of a third analyte or set of analytes, associated with a third physiological state of being, for example associated with opioid-induced intoxication. In this example, the third set of wavelengths, having a range of wavelengths that may overlap, or be distinct from, the first or second set of wavelengths, may be directed to either the first or second receptors, or the third set of wavelengths may be directed to a third receptor that is also configured within device 100 to touch a third body part.


The primary receptor 10, and the second receptor 25 (and if present the third receptor), of the present embodiment may each be comprised of a single sided probe that can make contact with a skin sample. Such a probe may comprise concentric rings of optic fibers so that each ring is made up by fibers carrying either input or output EMR. If the inner ring of fibers is carrying input EMR, then the outer ring of fibers may carry the output signal, or vice versa. Alternatively, the probe may comprise one or more input optic fibers and a separate set of output optic fibers positioned adjacent the input set of fibers. This type of probe may be used to determine the concentration of two or more than two compounds within the blood and interstitial fluid using reflectance, absorbance, and/or transmittance. During use, the probe may be placed against the skin of the finger, hand, arm, back or elsewhere (see Figurel).


Alternate configurations of an apparatus may also be used for the determination of a compound within a part, as described herein, including, but not limited to those described in US 2013/0248695, U.S. Pat. No. 5,361,758, WO 93/16629, U.S. Pat. Nos. 6,236,047, 5,429,128, 6,040,578 or U.S. Pat. No. 6,240,306 (all of which are incorporated herein by reference), with modification of the calibration algorithms so that they may be used to determine the concentration of two or more than two compounds of interest within each body part, deriving a biochemical profile, and determining a physiological condition of a the subject, as described herein.


The present embodiment provides a method to develop an algorithm that accounts for the differences in concentration of two or more than two compounds within the body part. For example, which is not to be considered limiting, if one of the compounds (or analytes, or metabolites) is glucose, then the concentration of glucose within each of the blood, and the interstitial fluid may be determined. From these values a reference measurement for glucose may be determined, and this reference value used to develop an algorithm. Absorbance values of a body part may be obtained over a set of wavelengths set as a dependent variable, and glucose reference measurement used as an independent variable. These values can then be processed using any suitable statistical procedure, including but not limited to, Partial Least Squares or Multiple Linear Regression to produce an algorithm for blood glucose. This procedure can be repeated for any compound or analyte of interest for which a concentration within blood and/or interstitial fluid is desired.


The concentration of a given compound may be calculated according to the present embodiment by using a calibration equation derived from a statistical analysis, for example but not limited to a least squares best fit, of a plot of the values of concentration of a calibration set of samples of the compound, which are determined using the method of the present embodiment, versus the values of the concentration of the calibration set measured directly by a different method. However, it is to be understood that other statistical tests may be used was known in the art, for example but not limited to multiple linear regression (MLR), partial least squares (PLS), and the like. Any known method for determining the concentration of one, or more than one, compound may be used as would be known to one of skill in the art.


In the case of glucose, as an example, and which is not to be considered limiting, blood glucose levels can be readily determined using well known in vitro techniques as known in the art. The level of glucose in the interstitial compartment may be determined using reverse ionotophoesis. In the case of THC and related metabolites, for example but not limited to, delta-9-tetrahydrocannabinol (THC), THC glucuronide (THCGlu), 11-nor-9-carboxy-THC (THC-COOH), 11-hydroxy THC (11-OH-THC), THC-COOH/11-OH-THC ratio, 11-nor-9-carboxy-THC glucuronide (THC-COOGlu), cannabidol (CBD), cannbinol (CBN), cannabigerol (CBG), delta-9-tetrahydrocannabivarin (THCV), THCV-carboxylic acid, 11 -nor-9-carboxy-delta-tetrahydrocannabivarin (THCV-COOH) levels can be readily determined using well known in vitro techniques as known in the art, for example using liquid chromatography and tandem mass spectrometry (Schwope D., et. al., 2011, Anal. Bioanal Chem. 410:1273-1283), or GC-MS (Marsot, A. et. al. (2016, J. Pharm. Pharm Sci 19:411-422). Detection of other analytes, or compounds, in blood, for example but not limited to, albumin, apolipoproteins A1 and B (apoA1 and apoB), total protein, triglycerides, blood sugar, calcium, ionized calcium, phosphate gamma-aminobutyric acid (GABA), alcohol, aldehyde, lactic acid, hemoglobin, blood urea nitrogen (BUN), albumin, apolipoproteins A1 and B (apoA1 and apoB), total protein, triglycerides, bicarbonate, electrolytes, sodium, potassium, magnesium, calcium, ionized calcium, glycated hemoglobin (A1C), high density lipoprotein (HDL), total cholesterol, omega-3 fatty acid, are well known in the art. These known tests may be used to determine the reference measurement of the compound or analyte in subjects who were exposed to a control treatment, or a range of THC under controlled conditions. These values may then be used as an independent variable in producing an algorithm for the non-invasive determination of corresponding blood analyte. These compound specific algorithms may therefore be used to ensure a proper estimation of the compound within the blood of the body part using non-invasively analyte determination.


By selecting a set of compounds or analytes that are associated with a physiological condition or state of a subject, and using standard measurement techniques, a biochemical profile describing the relative concentrations of these compounds may be determined that is correlated with the physiological condition. For example, if the physiological condition (state of interest) is THC induced intoxication, then two or more analytes, including for example but not limited to, delta-9-tetrahydrocannabinol (THC), THC glucuronide (THCGlu), 11-nor-9-carboxy-THC (THC-COOH), 11-hydroxy THC (11-OH-THC), THC-COOH/11-OH-THC ratio, 11-nor-9-carboxy-THC glucuronide (THC-COOGlu), cannabidol (CBD), cannbinol (CBN), cannabigerol (CBG), delta-9-tetrahydrocannabivarin (THCV), THCV-carboxylic acid, 11-nor-9-carboxy-delta-tetrahydrocannabivarin (THCV-COOH), albumin, apolipoproteins A1 and B (apoA1 and apoB), total protein, bilirubin, prolactin, triglycerides, creatinine, cortisol, glucose, lactate, Total 4, uric acid, blood urea nitrogen (BUN), blood sugar, calcium, ionized calcium, magnesium, sodium, phosphate, and GABA may be determined and the relative amounts of these analytes correlated with the physiological condition of THC-intoxication. As a result, these analytes may be used to obtain a biochemical profile which is an indicator of the physiological condition. By using the methods described herein, whole-blood cannabinoid pharmacokinetics that involves a plurality of analytes directly or indirectly derived from THC metabolism, may be considered in determining the status of the physiological condition (physiological state of being).


By biochemical profile it is meant an output that is derived from a set of values obtained from a set of measured analytes that corresponds to a physiological state of interest.


Additionally, the biochemical profile associated with a target state of being, for example a physiological state of interest, may include a plurality of ghost analytes that are observed to change in response to the state of being (or physiological condition), but whose identify may, or may not, be known. Ghost analytes may be characterized by comparing the absorbance spectra across a range of wavelengths under control or background conditions to obtain baseline values for each of the ghost analytes, with the absorbance spectrum obtained across a range of wavelengths in samples obtained under an induced physiological condition or state, and selecting one or more than one ghost analytes that increase or decrease under the induced physiological condition when compared to the baseline ghost analyte values. For example, a ghost analyte may be identified by analyzing the absorbance pattern across a range of wavelengths and identifying the wavelengths that characterize the ghost analyte (i.e. the ghost analyte displays an increase or decreased absorbance at one or more wavelengths that are specific for the ghost analyte). By characterizing the wavelength pattern associated with a ghost analyte, a baseline ghost analyte value for each ghost analyte may be determined and this value compared with each ghost analyte value determined in response to a physiological condition or state. Ghost analytes may include a plurality of analytes that can be used, along with other known analytes, to obtain a fingerprint or biochemical profile that may be used to define the status of a physiological condition as described herein.


This process for identifying and characterizing ghost analytes, or known analytes, may be repeated and appropriate software and machine learning applied to the acquired data sets to further optimize the predictive accuracy of the set of analytes used to determine the biochemical profile and the status of physiological condition. Using machine learning, complex models based on large data sets may be analyzed to identify acceptable “local minima” within the data sets and enable deep neural networks to be trained on identifying sets of analytes associated with the status of the physiological condition. Support vector machines (SVMs) and/or convolutional neural networks (CNNs) may be used with cross validation to provide insight into the algorithm's ability to generalize learned data representations. Cross validation involves partitioning the data into an arbitrary number of groups and iteratively using one of the groups for testing and the remaining for training (Mirowski, P.W., et. al., 2008 IEEE Workshop on Machine Learning for Signal Processing, Cancun, 2008, pp. 244-249.). Hyperparameter tuning may be used to develop an accurate, and robust final prediction.


In this way, an adaptive machine learning platform may be used to enable multiple machine learning algorithms to be executed for determining a physiological state of interest of a subject. The machine learning platform may include a plurality of machine learning components, each associated with a machine learning algorithm. For example, a machine learning component may be a component that utilizes one or more trained (or otherwise configured) machine learning algorithms that receive empirical data (e.g., data stored in one or more databases, for example including biochemical profile data, physiological parameter data and/or behavioral parameter data as described herein) to determine patterns or predictions that may be features of an underlying mechanism that generated the data and indicative of the physiological state of interest. The machine learning component may be able to utilize observed examples (e.g., from a set of training data) to capture characteristics of interest which may correspond to an unknown underlying probability distribution associated with the physiological state of interest. The adaptive machine learning system may allow a user to update decision-making strategies. Furthermore, the adaptive machine learning system, when operatively linked and communicating with a central CPU, may allow the central CPU to update decision-making strategies that can then be transmitted to the user of the hand held device, comprising one or more processors configured with executable instructions, in real time, or when the database is updated and the updated outputs are uploaded onto the handheld device, comprising one or more processors configured with executable instructions, as described herein.


The procedure for developing a prediction pipeline, for example, but not limited to a THC prediction pipeline, using machine learning may involve three general steps: data pre-processing, model selection and tuning (see FIG. 5):


Data pre-processing: data pre-processing is used to select high-quality training data from the overall data set, and to organize the data. A microprocessor is directed by software to collect and process the data from the device described herein. The collected data is aggregated, which involves appending new samples to existing database tables, formatted, to ensure that datatypes and column headers are consistent across tables in the database, and cleaned (filtering, interpolating, or keep missing values in the data set as needed), and use by a software model pipeline. Organized data is normalized and engineered features are used to predict intoxication. The features used by the final prediction algorithm may include frequency, time, and spatial domain data including all of the analytes being measured by the spectrometer. As required, engineered features may optionally be extracted using a combination of spatial or frequency domain filtering techniques. Software that may be used for data pre-processing includes, but is not limited to NumPY, SciPy, Pandas, Matlotlib or Seaborn. A microprocessor is directed by software to scan the linear array detector and calculate the second derivative of the spectrum computed. The microprocessor can then calculate the concentration of the particular constituents being measured using the absorbance and second derivative values for a number of selected wavelengths


Model Selection: Model architecture is be based on the structure and dimensionality of the data, with known baselines from prior literature on similar datasets being evaluated first before introducing additional complexity. Several algorithms, including but not limited to, support vector machines, deep neural networks, convolution neural networks, and generalized additive models with pairwise interactions may be used in the algorithm development and model selection stage. To ensure that the model complexity (ex. number of nodes, or layers) is appropriately chosen, a bottom up method is used, where an initial logistic regression model acts as a baseline that is compared to each new and increasingly complex model. Machine learning algorithms that may be used to analyze this data include, but are not limited to, support vector machines (SVMs) or convolutional neural networks (CNNs). Machine learning models that are too complex for the problem domain are capable of memorizing individual samples instead of learning the underlying distributions.


Model Evaluation: Cross validation is a data shuffling technique that provides an improved insight into an algorithm's ability to learn an underlying data distribution when limited data is available by testing on a series of training and hold out test splits. For example, data may be split into two subsections (training data and testing data), where the testing data is set aside, and a small percentage of the training data is split again further such that there are several training and testing data sets created from the initial training data subsection. The model is then iteratively trained and validated on these different data sets. The final model is then tested against the original testing data subsection for a final validation. Cross validation makes it easier to observe the generalization ability of the selected model, while also stretching smaller data sets that may otherwise be limited when the data set is split for training (80%) and testing (20%). Software that may be used for model selection includes, but is not limited to NumPY, Pandas, Matlotlib, Seaborn, Tensorflow or Keras.


Tuning stage: the selected model is optimized for performance. For example, hyperparameter search (learning rate, batch size, regularization coefficient, etc.) may be carried out . Hyperparameter tuning is an iterative process that may also be used to develop a robust final prediction. Software that may be used for fine tuning includes, but is not limited to NumPY, Pandas, Tensorflow or Keras.


Continuous Learning: New samples are added to the model using methods including, but not limited to, model weight updating, and champion and challenger. Model weight updating adds another training iteration to the existing model so that the new samples are considered while the model determines the underlying data distribution. Champion and challenger compares the performance of existing models with new models designed on new information with potentially improved insights and or methods. Several continuous learning strategies may be used in tandem.


Therefore, a system is also described herein. The system comprises a computer system that comprises one or more processors programmed with computer program instructions that, when executed, causes the computer system to provide a service platform that enables a developer to obtain training item information for training a machine learning model, for example using the data sets that comprise the biochemical profile, physiological parameters and behavioral parameters as described herein. The training item information indicates inputs and prediction outputs derived from one or more machine learning models' processing of the inputs. Additionally, the computer system may obtain, via the service platform, input/output information derived from one or more machine learning models. The input/output information indicates items provided as input to at least one model of the machine learning models. The service platform may then provide the input/output information derived from the machine learning models to update a first machine learning model, so that the first machine learning model is updated based on the input/output information being provided as input to the first machine learning model. The updated output may then be used to update a device for use within the field, for example, when the device is operatively linked and communicating with a central CPU. In this way, the central CPU (e.g. the service platform) is used update decision-making strategies that are transmitted to the user of the hand held device described herein in real time, or when the database is updated and the updated outputs are uploaded onto the handheld device as described herein.


The biochemical profile may be presented as a ratio of the relative amounts of measured analytes and ghost analytes, or as an index of the measured analytes and ghost analytes. For example, the index may be presented as a proportion of active to inactive analytes, of the various compounds within the blood that are associated, either positively or negatively with the physiological condition.


The above machine leaning process characterizes ghost analytes, known analytes, and other physiological and behavioral parameters, as defined herein, and over time, with repeated data input, the predictive accuracy of the set of analytes, physiological and behavioral parameters is increased.


For example, which is not to be considered limiting, in a THC induced intoxicated sate, a subject may comprise the analyte composition presented in Table 1. It is to be understood that this table presents a simplified data set to exemplify the method described herein. Other analytes, which may include known analytes or ghost analytes, are indicated as analytes A, B . . . C; A′, B′, . . . C′; A″, B″ . . . C″. These other analytes may increase or decrease in response to, or that are correlated with, the physiological condition and they may be included in this analysis. As shown in Table 2, a ratio of the relative abundance of the various selected analytes may be used to obtain a biochemical profile. Alternatively, an index, for example the relative proportion of a group of known “active” analytes (i.e. analytes that are known to be positively correlated with, or produce a THC-intoxicated state, including known and ghost analytes)—to-inactive analytes (i.e. analytes that are known to be correlated with, but not produce the THC-intoxicated state—these analytes may include known and ghost analytes) may be used to derive the biochemical profile.


In the case of smoked marijuana, THC peaks rapidly in the first few minutes after inhaling, and then declines quickly (within hours). THC may then remain at low levels of about 1-2 ng/ml for 8 hours or more. In chronic users, detectable amounts of blood THC can persist for days and therefore this analyte alone is not a reliable indicator of a physiological state of THC-induced intoxication. In chronic users of marijuana, residual THC was detected for 24 to 48 hours or longer at levels of 0.5-3.2 ng/ml in whole blood (1.0-6.4 ng/ml in serum; Skopp G., and Putsch, L., 2004, J. Anal Toxicol. 28: 35-40)









TABLE 2







Analyte composition of a subject before (Control), in an intoxicated state (1 hour after


consumption) and post-intoxicated state (36 hours after consumption). THC: delta-9-


tetrahydrocannabinol (a bio-active analyte); THC-COOH: 11-nor-9-carboxy-THC (a bio-inactive


analyte); 11-OH-THC: 11-hydroxy THC (A bio-active analyte): A (a bio-active analyte), B (a


bio-inactive analyte) . . . C (a bio-inactive analyte), A′, B′, . . . C′, A″, B″ . . . C″:


analytes that increase, decrease, or that have been correlated with, the physiological condition.














Analyte






Ratio


State (ng/ml)
THC*
11-OH-THC
THC-COOH
X
Y . . .
Z
Active:Inactive

















1 hr after
78
5
5
A
B . . .
C
78:5:5:A:B:. . .:C


consumption






(17.7 + A)/(B+ . . . +C)


36 hr after-
8
2
11
A′
B′ . . .
C′
8:11:2:A′:B′:. . . C′


consumption






(1 + A′)/(B′+ . . . +C′)


Control
2
0
3
A″
B″ . . .
C″
2:0:3:A″:B″:. . . C″









(0.7 + A″)/(B″+ . . . +C″)





*THC levels above 3.5-5 ng/ml in blood (or 7-10 ng/ml in serum) indicate likely impairment.






In the example provided in Table 2, in an intoxicated state (1 hr after consumption) the subject exhibits an Active-to-Inactive Index of (17.7+A)/(B+ . . . +C), while the Index associated with a post consumption condition or state, or a control condition or state, are well below this Index (1+A′)/(B′+ . . . +C′), or (0.7+A″)/(B″+ . . . +C″), respectively. In this example, the result may be considered positive for a THC-induced intoxicated state if the Index is greater than a preset value which is determined based on an analysis of the analytes determined in subjects who were exposed to a control treatment, or a range of THC under controlled conditions, for example an Index value of (2+A)/(C+ . . . +C) may be an indication of a positive result for the state of THC-induced intoxication.


Alternatively, the set of ratios for analyte concentrations determined in subjects who were exposed to a control treatment, or a range of THC under controlled conditions, may be compared against the same set of ratio of analytes for the test subject (as presented in Table 2), and these sets of ratios may be used to determine if a threshold value has been achieved indicating a THC-induced intoxicated state.


Other methods of processing the measured analyte concentration to produce a biochemical profile may be used to determine if a threshold value has been obtained and indicating that the subject is positive for the corresponding physiological state may also be used.


A similar approach as described above may be used to determine the biochemical profile for other physiological states or conditions, for example, but not limited to, alcohol-induced intoxication, a combination of cannabis and alcohol induced-intoxication, or from the consumption of opiates, fentanyl, amphetamines, phencyclidine, sedatives, anxyolytics, cocaine, caffeine, and nicotine.


In addition to obtaining a biochemical profile as an indicator of a physiological condition as descried above, additional physiological parameters may be used to further assist in characterization of the physiological condition or state. Furthermore, behavioral parameters may also be considered in combination with the physiological parameters and biochemical profile data that as acquired. Examples of behavioral parameters include a determination of mental acuity (e.g. the name-face test, fire alarm test, two delayed recall tests, misplaced objects test, shopping list test, digit symbol test), one or more motor skill test (walk and turn test, one leg stand test, horizontal gaze nystagmus test, a divided attention test, a rhomberg balance test), the ability to function at a defined task, for example to operate machinery, drive an automobile, standardized field sobriety.


For example if the physiological condition is cannabis or THC-induced intoxication, then physiological parameters may include, for example but not limited to, heart rate, pulse rate, body temperature, neuropeptide Y, fatty acid amide hydrolase (FAAH), C reactive protein (cRP), creatine kinase (CK), aspartate amino transferase (AAT), alanine transaminase (ALT), gamma-glutamyl transpeptidase (GGT), aspartate aminotransferase (AST), white blood cell count (WBC), red blood cell count (RBC), hemoglobin, hematocrit, neutrophils, lymphocytes, eosinophils, hypoactivity; THC concentration in hair, THC concentration in urine, and this physiological parameter data is combined with the biochemical profile data obtained using two or more than two analytes, for example, but not limited to delta-9-tetrahydrocannabinol (THC), THC glucuronide (THCGlu), 11-nor-9-carboxy-THC (THC-COOH), 11-hydroxy THC (11-OH-THC), THC-COOH/11-OH-THC ratio, 11-nor-9-carboxy-THC glucuronide (THC-COOGlu), cannabidol (CBD), cannbinol (CBN), cannabigerol (CBG), delta-9-tetrahydrocannabivarin (THCV), THCV-carboxylic acid, 11-nor-9-carboxy-delta—tetrahydrocannabivarin (THCV-COOH), albumin, apolipoproteins A1 and B (apoA1 and apoB), total protein, bilirubin, prolactin, triglycerides, creatinine, cortisol, glucose, lactate, Total 4, uric acid, blood urea nitrogen (BUN), blood sugar, calcium, ionized calcium, magnesium, sodium, phosphate, and gamma-aminobutyric acid (GABA), and any one or more than one ghost analyte that is associated with THC-induced intoxication, in order to produce an output that defines the status of the physiological condition (THC-induced intoxication). Additionally, behavioral parameters, for example determination of mental acuity (for example but not limited to a name-face test, a fire alarm test, a two delayed recall tests, a misplaced objects test, a shopping list test, a digit symbol test), one or more motor skill test (for example but not limited to, a walk and turn test, a one leg stand test, a horizontal gaze nystagmus test, a divided attention test, a rhomberg balance test), the ability to function at a defined task, for example to operate machinery, drive an automobile, standardized field sobriety may be performed and the data combined with the biochemical profile and physiological parameters to determine the physiological condition.


The results from these methods may be combined to produce a value or index of the biochemical profile, the physiological parameter, the behavioral parameter, or a combination thereof, and these values, or index values, may be used to determine if a threshold value has been obtained or exceeded, by comparing the value or index value against a reference value or reference index value, thereby indicating that the subject is positive for the corresponding physiological state.


Similarly, if the physiological state being evaluated is alcohol-induced intoxication, then the physiological parameters may include, but are not limited to, heart rate, body temperature, neuropeptide Y, aspartate amino transferase (AAT), alanine transaminase (ALT), gamma-glutamyl transpeptidase (GGT). This data is combined with the biochemical profile data obtained using two or more than two analytes including, for example but not limited to, alcohol, aldehyde, lactic acid, and any one or more than one ghost analyte that is associated with alcohol-induced intoxication, to produce an output that defines the status of the corresponding physiological condition that is being tested (alcohol-induced intoxication). Additionally, behavioral parameters, for example determination of mental acuity (for example but not limited to a name-face test, a fire alarm test, a two delayed recall tests, a misplaced objects test, a shopping list test, a digit symbol test), one or more motor skill test (for example but not limited to, a walk and turn test, a one leg stand test, a horizontal gaze nystagmus test, a divided attention test, a rhomberg balance test), the ability to function at a defined task, for example to operate machinery, drive an automobile, standardized field sobriety, may be performed and the data combined with the biochemical profile and physiological parameters to determine the physiological condition.


Furthermore, if the intoxicated state being evaluated is, for example resulting from a combination of cannabis and alcohol, or opiates, fentanyl, amphetamines, phencyclidine, sedatives, anxyolytics, cocaine, caffeine, and nicotine consumption, then the physiological parameter may include one or more of: heart rate, pulse rate, body temperature, neuropeptide Y, fatty acid amide hydrolase (FAAH), c reactive protein (cRP), creatine kinase (CK), aspartate amino transferase (AAT), asparate aminotransferase (AST), alanine transaminase (ALT), gamma-glutamyl transpeptidase (GGT), white blood cell count (WBC), red blood cell count (RBC), hemoglobin, hematocrit, neutrophils, lymphocytes, eosinophils, hypoactivity; THC in hair, THC in urine. This data is combined with the biochemical profile data obtained using two or more than two analytes may include: then two or more than two analytes may include: delta-9-tetrahydrocannabinol (THC), THC glucuronide (THCGlu), 11-nor-9-carboxy-THC (THC-COOH), 11-hydroxy THC (11-OH-THC), THC-COOH/11-OH-THC ratio, 11-nor-9-carboxy-THC glucuronide (THC-COOGlu), cannabidol (CBD), cannbinol (CBN), cannabigerol (CBG), delta-9-tetrahydrocannabivarin (THCV), THCV-carboxylic acid, 11-nor-9-carboxy-delta-tetrahydrocannabivarin (THCV-COOH), albumin, apolipoproteins A1 and B (apoA1 and apoB), total protein, bilirubin, prolactin, triglycerides, creatinine, cortisol, glucose, lactate, Total 4, uric acid, blood urea nitrogen (BUN), blood sugar, calcium, ionized calcium, magnesium, sodium, phosphate, GABA, alcohol, aldehyde, lactic acid, and any one or more than one ghost analyte that is associated with alcohol-induced intoxication, to produce an output that defines the status of the corresponding physiological condition that is being tested (a combination of cannabis and alcohol, or opiates, fentanyl, amphetamines, phencyclidine, sedatives, anxyolytics, cocaine, caffeine, and nicotine consumption). Additionally, behavioral parameters, for example determination of mental acuity (a name-face test, a fire alarm test, a two delayed recall tests, a misplaced objects test, a shopping list test, a digit symbol test), one or more motor skill test (a walk and turn test, a one leg stand test, a horizontal gaze nystagmus test, a divided attention test, a rhomberg balance test), the ability to function at a defined task, to operate machinery, drive an automobile, standardized field sobriety (Newmeyer, Swortwood, Taylor, et al., 2017, Clin Chem, 63(3), 647-662. doi:10.1373/clinchem.2016.265371).


The near infrared region of the electromagnetic spectrum may be used for the measurements of samples as described herein. Measurements may be obtained over a range of wavelengths for example from about 350 nm to about 2500 nm range. Chemical species (analytes and ghost analytes) exhibit characteristic absorption bands within this spectral interval which may be used to characterized each analyte. The near infrared region is well-suited to in vivo diagnostic applications since human tissue is transparent to the incident radiation and therefore sufficient penetration of the radiation is possible to allow accurate quantitative analysis.


The source of EMR used in the present embodiment is preferably near-infrared light, for example but not limited to a polychromatic light source. This type of light source can emit light over a very wide bandwidth including light in the near infrared spectrum. In this case, the light from the light source may pass through a collimator, which is a collection of lenses that concentrate the light into a narrow parallel beam directed at the receptor. The polychromatic light source can be a quartz-halogen or a tungsten-halogen bulb and is powered by a stabilized power source, for example, a DC power supply, or by a battery. This polychromatic light source may be a tungsten-halogen lamp or it may be a collection of LEDs or other light sources selected to emit radiation in the range of about 350-2500 nm, or for example, from about 650 to about 1100 nm.


A receptor is preferably used which is shaped to receive a part of the subject for sampling, for example a clamped part of the skin, or a finger. Alternatively, the receptor could be shaped so that the part of the human, onto which the EMR is to be directed, is placed against the receptor rather than within the receptor. It is preferred that the sampled body part is in close contact with the receptor. Examples of receptors that may be used are provided in US 2013/0248695, U.S. Pat. No. 5,361,758, WO 93/16629, U.S. Pat. Nos. 6,236,047, 6,040,578 or 6,240,306


The EMR is directed onto, and dispersed by, a part of the subject. The dispersed light from the body part, either reflected, transmitted, or both, is collected by using any suitable method, for example, fiber optics, or lenses, and the output signal directed to a diffraction device that separates the wavelengths of light within the output signal into their component parts. Examples of a diffraction device include but are not limited to a diffraction grating or a holographic grating.


The collected signal can comprise EMR that has passed through a part of a subject or has reflected off of a part of the subject, or a combination thereof. The diffracting device may disperses the EMR into its component wavelengths so that the infrared region falls along the length of a detector such as, but not limited to a linear array detector (e.g. a 256 element photo diode array), or a charged couple device (CCD). In the case of an array, the detector has a series of diodes and is preferably electronically scanned by a microprocessor to measure the charge accumulated on each diode, the charge being proportional to the intensity of EMR for each wavelength transmitted through or reflected from the part of the subject in the receptor. The detector is connected to the microprocessor, producing an output spectrum, with the microprocessor analyzing the measurements and ultimately producing a result for each concentration level determined. The result can be stored, shown on a display, or transmitted to another central CPU for further analysis or display. A keyboard may be used to control the device, the central CPU, or both, for example, to specify a particular physiological condition (and corresponding set of analytes) to be measured. The timing and control are activated by the microprocessor to control the device, for example, to determine number and timing of measurements.


After measurements are obtained for the transmittance, reflectance or both, the log of the inverse of these measurements is preferably taken, that is, log 1/T and log 1/R, where T and R represent the transmittance and reflectance respectively. A reference set of measurements is taken of the incident light, being the light generated in the device when no part of the subject is in contact with the receptor. The absorbance is then calculated when a part of the subject is in contact with the receptor as a ratio of measurements compared to the reference set of measurements. If desired, a second derivative of the measurements may be obtained to reduce any variation in the result that may be caused by a change in path length for the light caused by measuring the compound concentration in different thicknesses of the parts of the subject. The second derivative calculation may be used to eliminate base line shifts due to different path lengths or absorbing water bands, and in addition, enhances the separation of overlapping absorption peaks of different constituents of the mixture being analyzed. The microprocessor can collect the plurality of spectra produced and calculate the second derivative of the averaged results.


The results obtained may vary with the temperature of the part of the subject, the device used in the method of the present embodiment may contains a temperature sensor so that the temperature of the analyzed part can be measured rapidly at the time of the spectral sampling. This temperature sensor may comprise a small-mass thermocouple. Computer software can then be used to allow the microprocessor to compensate for spectrum deviations due to the temperature.


The linear array detector is preferably a photo diode array that is positioned to intercept, across its length, the dispersed spectrum from the diffraction grating. The microprocessor is directed by software to scan the linear array detector and calculate the second derivative of the spectrum computed. The microprocessor can then calculate the concentration of the particular constituents being measured using the absorbance and second derivative values for a number of selected wavelengths. A calibration equation is preferably used for each constituent and is determined by the compound being measured.


The measured data may be obtained at the road side, for example by a law enforcement officer, or at a point-of care testing facilities. After the data is collected by the device as described herein, the results may be transmitted to a central computer for further analysis. The measured data may be combined with measured physiological parameters and measured behavioral parameters, and an overall index of the physiological state of interest, for example intoxication, of the subject determined. A plurality of overall indices, each indicative of a physiological state of interest of a subject, may be pooled and delivered for meta-analysis by an interested party, for example health care providers, law enforcement agencies, federal government, or other services that may have an interest in the pooled data.


Also provided herein there is a device for detecting a physiological state of interest of a subject. The device comprises:

    • a source of electromagnetic radiation (EMR; 30) that emits a plurality of wavelengths of EMR from about 350 nm to about 2500 nm, the source of EMR being operatively coupled to a power source;
    • a receptor 10 sized to register with, and fit against, a sample 20, the receptor comprising one or more than one port;
    • one or more than one input radiation guiding element 40 in operable association with the source of EMR, one or more than one output radiation guiding element 50 in operable association with a detector 60,
    • the one or more than one input radiation guiding element and the one or more than one output radiation guiding element in optical alignment with the one or more than one port located and defining an EMR path within the receptor when the receptor is registered with, and fit against, the sample;
    • the detector for measuring transmitted or reflected EMR received from the sample, the detector operatively coupled to a processing system 70;
    • the processing system comprising one or more than one algorithm for determining a concentration for two or more than two analytes in the sample, and using the one or more than one algorithm to derive the physiological state of interest of the sample, wherein, the physiological state of interest is:
    • i) intoxication, then two or more than two analytes may include: delta-9-tetrahydrocannabinol (THC), THC glucuronide (THCGlu), 11-nor-9-carboxy-THC (THC-COOH), 11-hydroxy THC (11-OH-THC), THC-COOH/11-OH-THC ratio, 11-nor-9-carboxy-THC glucuronide (THC-COOGlu), cannabidol (CBD), cannbinol (CBN), cannabigerol (CBG), delta-9-tetrahydrocannabivarin (THCV), THCV-carboxylic acid, 11-nor-9-carboxy-delta-tetrahydrocannabivarin (THCV-COOH), albumin, apolipoproteins A1 and B (apoA1 and apoB), total protein, bilirubin, prolactin, triglycerides, creatinine, cortisol, glucose, lactate, Total 4, uric acid, blood urea nitrogen (BUN), blood sugar, calcium, ionized calcium, magnesium, sodium, phosphate, and gamma-aminobutyric acid (GABA);
    • ii) alcohol induced intoxication, then two or more than two analytes may include: alcohol, aldehyde, and lactic acid; or
    • iii) intoxication generally, for example arising from cannabis, alcohol, opiates, fentanyl, amphetamines, phencyclidine, sedatives, anxyolytics, cocaine, caffeine, and nicotine consumption, then two or more than two analytes may include:, then two or more than two analytes may include: delta-9-tetrahydrocannabinol (THC), THC glucuronide (THCGlu), 11-nor-9-carboxy-THC (THC-COOH), 11-hydroxy THC (11-OH-THC), THC-COOH/11-OH-THC ratio, 11-nor-9-carboxy-THC glucuronide (THC-COOGlu), cannabidol (CBD), cannbinol (CBN), cannabigerol (CBG), delta-9-tetrahydrocannabivarin (THCV), THCV-carboxylic acid, 11-nor-9-carboxy-delta-tetrahydrocannabivarin (THCV-COOH), albumin, apolipoproteins A1 and B (apoA1 and apoB), total protein, bilirubin, prolactin, triglycerides, creatinine, cortisol, glucose, lactate, Total 4, uric acid, blood urea nitrogen (BUN), blood sugar, calcium, ionized calcium, magnesium, sodium, phosphate, GABA, alcohol, aldehyde, and lactic acid.


For example, the device comprising the receptor 10 may be housed in a ‘computer mouse-like’ housing. To operate the device, an operator activates the program, for example, an App on their cell phone 90 or a program on a remote lap top and turns on the device. The device may comprise an outer layer which is translucent in color and when the device is turned on, the outer layer may turn a yellow hue indicating a standby state. The subject being tested places their body part, for example a finger 20 into a small cavity on the top of the mouse-like receptor 10. Once the finger of the individual is inserted into the cavity and the cavity —body part interface is dark, the light source 30, for example LED's, within the device will be triggered to scan the finger and take a measurement. When the cavity is dark, and the scanning begins the outer layer of the device may change from yellow to green indicating that a sample is being obtained. After a period of time the scanning will be complete, the outer layer of the device may turn red signaling that the test is complete, and the finger can be removed. Data obtained from the test may be processed in the device, or sent to a remote CPU for further processing, for example via blue tooth for further processing.


The updated may then be used to update a device for use within the field, for example, when the device is operatively linked and communicating with a remote or central CPU. In this way, the central CPU may be used update decision-making strategies that are transmitted to the user of the hand held device described herein in real time, or when the database is updated and the updated outputs are uploaded onto the handheld device as described herein.


When used, the device as described above, and based on the biochemical profile, or the biochemical profile in combination with physiological parameter(s), or the biochemical profile in combination with physiological parameter(s) and behavioral parameter(s), may determine that the physiological state of interest of a subject indicative of a state of intoxication has been realized and corrective action may be required. For example, if the device (and optionally physiological parameters and behavioral parameters) is used for road-side testing by a law enforcement officer (a non-limiting example of an operator of the test) and the driver of the car is determined to be in an intoxicated state, then the law enforcement officer may perform corrective action and confiscate the car, suspend the driver's license, press charges and the like.


In other circumstances, the result derived from the device (and optionally physiological parameters and behavioral parameters) may be forwarded to a third party so that corrective action may be taken by the third party. For example, the operator of the test, for example but not limited to, a health care practitioner or the law enforcement officer, may forward the positive result indicating intoxication (impairment) to a third party, for example, a justice of the peace, and corrective action may be taken. Alternatively, if safety is a requirement of the subjects employment, and the subject has been determined to exhibit a positive result indicating intoxication (impairment), then the result may be forwarded to the subject's employer. Examples of situations where safety may be a requirement of the subjects employment, include if the subject is working as an air traffic controller, the subject is a pilot, they operate a commercial vehicle, they operate machinery (large or small) at a work site, they are an operator at a nuclear power facility etc.


To test or calibrate the device a synthetic sample or ‘phantom finger’ (U.S. Pat. No. 6,657,717, which is incorporated herein by reference) comprised of pre-defined materials may be applied against the receptor. Alternatively, the operator of the device may use their own corresponding body part or finger.


Example: Trial to Assess Physiological State of Interest of a Subject

Participants: approximately 500 patients (3 different people per day, every 5 days a week, every 4 weeks a month, for 8 months). Status of patients regarding usage determined and indexed as heavy users to minimal users. Participants are screened for psychiatric disorders using the Structured Clinical Interview for DSM-IV axis I Disorders (SCID-I). All subjects with a psychiatric disorder warranting treatment are excluded from the study. Females use an approved method of birth control for the duration of the study. Participants refrain from use of cannabis for 72 prior to the text sessions. To determine the participants baseline (pre-test) level of THC, a pre-test saliva and urine sample are obtained and tested for THC. Additionally, a breathalyzer test is performed to detect recent alcohol use.


Demographics: ages 18 to 70, multiple races, multiple skin colors, nationalities, genders, various weights, night time and day time testing, heavy users and novices.


Tests: Participants are asked to eat a light breakfast (e.g. a muffin or bagel) before each test session. Tests are performed on sober patients (background), followed by inducing intoxication and testing from before consumption of the intoxicant to four to six hours after consumption of the intoxicant.


Participants complete subjective effects questionnaires and a baseline driving trial. Once these procedures are complete, participants are given one oral dose of cannabis. Blood, subjective tests, cognitive tests and driving trials, vitals and visual analog scale (VAS) are conducted at regular intervals over 7 hours after dosing.


Intoxication includes a) administering repeated loading of cannabis (as an oil; of from 0 to 75 mg), and b) administering repeated loading of known amount of cannabis (of from 0 to 75 mg), and known amounts of alcohol. Participants drive the driving simulator before and after ingesting oral THC in a single session. Blood is drawn before and during the treatment in order to determine the status of analytes (see below) brought about as a result of a state of impairment. During the test period, each participant reaches an intoxicated state brought about by THC as measured by motor skill and mental acuity impairment.















Study Session



Cannabis smoking marks Time 0





















Base-
5
15
30
60
90
2
3
4
5
6



Eligibility
line −120 m
m
m
m
m
m
h
h
h
h
h























Driving Trial
X
X


X

X

X





Breath tests (alcohol.)
X
X


Infrared detections

X
X
X
X
X
X
X
X
X
X
X


Physical Exam,
X


Psychiatric Exam


(SCID)


Vital Signs
X
X
X
X
X
X
X
X
X
X
X
X


Urine: Point-of-
X
X


care drug screen


Urine: BSS
X
X


Urine: Point-of-
X
X


care pregnancy test


Blood: Biochemistry,
X


Hematology


Blood: THC and

X
X
X
X
X
X
X


metabolites quantification


Blood: other analytes

X


X

X
X



X


Saliva: THC detections
X
X


X

X
X



X














[PDCB]


VAS

X
X
X
X
X
X
X
X
X
X
X


Verbal free recall

X


X

X

X


Demographics and self-

X


report questionnaire


Timeline follow back
X
X









Physical Measurements: 1) Breath sample for alcohol; 2) Physical examination 3); Vital signs (temperature, pulse, blood pressure, respiration rate), height and weight; 4) Blood samples to measure biochemistry and haematology, THC, CBD and metabolites quantification (see below); 5) Urine sample for toxicology screening for drugs of abuse (point-of-care testing); 6) Broad spectrum urine screen; 7) Urine sample for pregnancy testing (point-of-care testing); 8) Saliva sample for determination of saliva THC.


Behavioural Information: 1) Psychiatric examination: Structured Clinical Interview for DSM-IV Axis I Disorders (SCID-I); 2) Timeline Follow Back (TLFB) for 3 months at eligibility assessment and for 7 days at the test session. In the TLFB participants report use of a substance each day for a number days prior to the assessment.


Cognitive/Psychomotor Test: Verbal Free Recall Task for verbal learning and memory.


Subjective Assessment for Cannabis Effects: Self-reports of cannabis effects using Visual Analog Scales (VAS)


The tests are repeated per participant. In addition to a physical examination, the following data is obtained from each participant:

    • i) a blood draw, and a urine sample, are timed to measure THC in the blood (inhalation of, and/or edible consumption of, cannabis) and analytes, and other physiological parameters are determined, including body temperature, pulse, blood pressure, rate of respiration, C-reactive protein, creatine (IDMS), glucose, blood urea nitrogen (BUN), THC, total protein, albumin, prolactin, potassium, sodium, cortisol, lactate, Total T4, calcium, ionized calcium, uric acid, triglyceride, magnesium, creatine kinase, gamma-glutamyl transferase (GGT), aspartate aminotransferase (AST), total bilirubin, WBC count, RBD count, hemoglobin, hematocrit, neutrophils, lymphocytes, eosinophils.


For patients receiving cannabis, two or more of the following analytes may be determined: delta-9-tetrahydrocannabinol (THC), THC glucuronide (THCGlu), 11-nor-9-carboxy-THC (THC-COOH), 11-hydroxy THC (11-OH-THC), THC-COOH/11-OH-THC ratio, 11-nor-9-carboxy-THC glucuronide (THC-COOGlu), cannabidol (CBD), cannbinol (CBN), cannabigerol (CBG), delta-9-tetrahydrocannabivarin (THCV), THCV-carboxylic acid, 11-nor-9-carboxy-delta-tetrahydrocannabivarin (THCV-COOH), albumin, apolipoproteins A1 and B (apoA1 and apoB), total protein, bilirubin, prolactin, triglycerides, creatinine, cortisol, glucose, lactate, Total 4, uric acid, blood urea nitrogen (BUN), blood sugar, calcium, ionized calcium, magnesium, sodium, phosphate, and GABA. Ghost analytes are also tracked to determine which my correlate with, and can be used to determine, the state of cannabis-induced intoxication. The physiological parameter may include one or more of: heart rate, pulse rate, body temperature, neuropeptide Y, fatty acid amide hydrolase (FAAH), c reactive protein (cRP), creatine kinase (CK), aspartate amino transferase (AAT), asparate aminotransferase (AST), alanine transaminase (ALT), gamma-glutamyl transpeptidase (GGT), white blood cell count (WBC), red blood cell count (RBC), hemoglobin, hematocrit, neutrophils, lymphocytes, eosinophils, hypoactivity; THC in hair, THC in urine.


For patients receiving alcohol in addition to cannabis, additional analytes for testing may include: alcohol, aldehyde and lactic acid. Furthermore, ghost analytes are also tracked to determine which my correlate with, and can be used to determine, the state of cannabis and alcohol-induced intoxication. The physiological parameter may include measurement of: heart rate, body temperature, neuropeptide Y, aspartate amino transferase (AAT), alanine transaminase (ALT), gamma-glutamyl transpeptidase (GGT).

    • ii) two scans using the non-invasive device described herein, are obtained from the patients finger at the same time as each blood draw (step i) is obtained;
    • iii) a measure of mental acuity and motor skill function is determined following each cannabis, or cannabis and alcohol loading. The behavioral parameters may include determination of mental acuity (a name-face test, a fire alarm test, a two delayed recall tests, a misplaced objects test, a shopping list test, a digit symbol test), one or more motor skill test (a walk and turn test, a one leg stand test, a horizontal gaze nystagmus test, a divided attention test, a rhomberg balance test), the ability to function at a defined task, to operate machinery (simulated), drive an automobile (simulation, to determine impaired driving skills, including reaction time, collisions, mean speed, mean speed while distracted, lateral control, lateral control while distracted), for each test group (i.e. patients receiving cannabis, or patients receiving both cannabis and alcohol).


Driver Simulator Testing: the simulator consists of a driver's side instrument cluster, steering wheel, controls, and center console as in a GM compact car. The steering wheel, brake and accelerator pedals provide dynamic force feedback. The visual system comprises three 50-inch screens providing a 180° field of view in the front, and two 17-inch side displays providing visual feedback for the left and right blind zones.


Participants receive a series of simulator training trials at the start of the study session to become familiar with the simulated vehicle's steering, accelerator, and braking controls. The driving simulations used for main effects testing consist of a series of driving events designed to assess mechanisms by which cannabis consumption may impact driver performance. Dependent measures of driver performance (e.g., standard deviation of lateral position, mean speed) are based on global performance throughout the entire simulation as well as event-specific performance and are measured using the simulator software. Risk-taking behaviour is assessed by measuring average speed throughout the simulation. A divided attention task is included in some of the driving scenarios to increase cognitive load and to better simulate real-world conditions.


Analysis: Deep neural network (DNN) architecture supplemented by use of generalized additive models (GAMs; which provide the interoperability of logistic regression, with the capability to solve non-linear problems) are used to provide enhanced interpretability of model decisions.


Safety of the patients pre, during and post intoxication is ensured.


In at least some embodiments, the state of interest may correspond to whether a subject has COVID-19. The analytes relevant for this determination in at least some example embodiments follow:


















Increased
Decreased




in response
in response




to COVID-19
to COVID-19



Analyte
infection
infection
















BLOOD CELLS











WBC, total

X*





X



Lymphocytes

X





X





X



T cells

X



Leukocytes
X*



Monocyte %

X*



Eosinophil %

X





X



Basophil %

X*



Neutrophil-
X



Lymphocyte



Ratio (NLR)



Platelets
X



Platelet-
X



Lymphocyte



Ratio (PLR)







COAGULATION FACTORS











D-dimer
X




Fibrin/Fibrinogen
X



degradation



products (FDP)



Fibrinogen
X



Antithrombin

X



Prothrombin

X



time activity



Thrombin time

X







INFLAMMATORY MARKERS











C-reactive
X




protein
X



Erythrocyte
X



sedimentation rate



Procalcitonin






X



LDH
X



Albumin

X



Serum ferritin
X



IL-2R
X




X



IL-6
X




X



IL-8






X*



IL-10
X*



TNFa
X




X



Immuno-





globulins



(IgA, IgG, IgM)



Complement





proteins (C3, C4)











Determining State of Interest without Direct Reference to Analytes


In at least some example embodiments, determining state of interest may be performed without direct reference to analytes of a subject and, consequently, without sampling blood of the subject. Instead, one or more reference spectra are empirically determined to correspond to a particular reference state of interest of a subject, and one or more measured spectra are compared to those one or more reference spectra. Based on that comparison, which in certain embodiments leverages machine learning, a processor determines whether the state of interest of the subject corresponds to the reference state of interest. For example, the reference state of interest may be indicative of a particular disease, such as COVID-19. The relationship between the reference spectra and the state of interest is established without having to sample and analyze blood. Determining state of interest is done without direct reference to analytes of the subject since the relationship of particular analytes to wavelengths represented in the measured spectra, and the influence of particular components of the spectra to state of interest, are unknown.


In at least some example embodiments, a method for determining a subject's state of interest starts with directing light at a body part of a subject such that the light passes through or is reflected by blood and interstitial fluid of the body part. The light that is incident on the body part comprises a range of wavelengths from at least one of the near infrared and visible spectra. A spectrum of the light is then measured using a spectrometer after the light has one or both of passed through and been reflected by the body part. The measured spectrum comprises the range of wavelengths incident on the body part. A processor then compares the measured spectrum against a reference spectrum representative of a known physiological state of interest, such as whether a subject has a disease such as COVID-19. After that comparison, the processor determines whether the subject is in the known physiological state of interest from a similarity between the measured spectrum to the reference spectrum, and without having to directly reference analytes of the subj ect.


Referring now to FIG. 6, there is shown a system 600 for determining a physiological state of interest of a subject. The system 600 comprises a computer 602, an enclosure 604, and an interface 606. The computer 602 is electrically coupled to the enclosure 604 via an electrical cable 610 to permit communication between the computer 602 and enclosure 604, and the enclosure 604 is optically coupled to the interface 606 via a fiber optic cable 608 to permit one or more spectrometers within the enclosure 604 to obtain one or more spectra via the interface 608. The computer 602 comprises a first processor communicatively coupled to a first computer readable medium and to a communications port; the enclosure 604 comprises a second processor communicatively coupled to one or more spectrometers, a lamp, a display, a second computer readable medium, an electrical communications port, and an optical communications port; and the interface comprises an optical communications port to receive and transmit an optical signal. As discussed in further detail below, the computer 602 is responsible for initiating the measurement sequence, configuring the one or more spectrometers in the enclosure 604, and recording raw spectral data measured by the one or more spectrometers.


More particularly, the one or more spectrometers within the enclosure 604 comprise a first spectrometer configured to measure spectra in the visible and near infrared ranges (e.g., from approximately 350 nm to 1,000 nm) (the “VIS-NIR spectrometer”), and a second spectrometer configured to measure spectra in the near infrared range (e.g., from approximately 900 nm to 2,500 nm) (the “NIR spectrometer”). A single lamp emits the light used for spectral readings. The enclosure's 604 display shows the lamp's total hours of use.


The fiber optic cable 608 comprises multiple optical fibers 806. Depending on whether the interface 606 relies on transmission or reflectance, as discussed further below in respect of FIGS. 8A, 8B, 9A, and 9B, one or more of the optical fibers 806 is used to transmit light from the lamp to the interface 606, and one or more optical fibers 806 is used to return light from the interface 606 to the spectrometers for spectral analysis. In at least some example embodiments, the cable 608 comprises eight optical fibers 806: a source fiber 806a transmits light from the lamp to the interface 606, three return fibers 806b return light from the interface 606 after it has interacted with the subject to the VIS-NIR spectrometer, and four return fibers 806b return light from the interface 606 after it has interacted with the subject to the NIR spectrometer.


While in the depicted embodiment the spectrometers are in the enclosure 604 and the interface 606 is distinct from the enclosure 604, in alternative embodiments (not depicted) a single housing contains a processor, computer readable medium, display, spectrometers, lamp, and also acts as the interface 606. For example, in at least some example embodiments the enclosure 604 and interface 606 may be a combined device that applies Fourier Transform Near Infrared Spectroscopy (“FT-NIR”) to obtain transmission and/or reflectance measurements. An example combined device is a Bruker™ Tango™ FT-NIR spectrometer from Bruker Optics Inc. marketed for use with non-human specimens, such as inorganic materials. Despite not being designed to directly perform spectroscopy on humans, the Bruker™ Tango™ FT-NIR spectrometer, and FT-NIR spectroscopy more generally, may in at least some embodiments be used to capture the spectroscopic readings analyzed herein. While the combined device is described as applying FT-NIR spectroscopy, in other embodiments the FT-NIR spectrometer may be used as a spectrometer within the enclosure 604 in conjunction with the interfaces 606 as described above.


Referring now to FIG. 7, there is shown a block diagram illustrating information flow in the system 600 of FIG. 6. FIG. 7 shows how the laptop 602, enclosure 604, and interface 606 communicate at different times t1 through t4. To start, the computer 602 initiates a program and sends an initiation signal to the enclosure 604 at time t1. In response to the initiation signal the one or more spectrometers within the enclosure 604 begin a spectral acquisition process, which comprises sending near infrared and/or visible light through the fiber optic cable 608 to the interface 606 at time t2. More particularly, each of the spectrometers takes four measurements: a light reference sample, a dark reference sample, a light sample, and a dark sample. As discussed further below, the light and dark reference samples are taken to mitigate or eliminate sensor variance across various devices, including the interfaces 606 and spectrometers within the enclosure 604. The light reference sample and light sample are taken with the lamp shutter open, while the dark reference sample and dark sample are taken with the lamp shutter closed. The light and dark samples are taken using a body part of a subject, such as the subject's finger.


The light is scattered/reflected/transmitted and then returned to the one or more spectrometers within the enclosure 604 via the cable 608 at time t3. Spectra are measured within the enclosure 604 by the one or more spectrometers and converted to an electrical signal returned to the computer 602 at time t4 for storage. The spectral data is stored as eight separate files: for each of the VIS-NIR and NIR spectrometers, one file for the light reference sample, one file for the light sample, one file for the dark reference sample, and one file for the dark sample. One or both of the processors in the computer 602 and enclosure 604 perform the processing on the eight files that compares one or more measured spectra to one or more reference spectra, and makes a determination as to the subject's state of interest. The computer's 602 display displays the measured spectra being recorded by each of the spectrometers, and also the final determination of the subject's state of interest.


Referring now to FIGS. 8A and 8B, there are respectively shown perspective and sectional views (along line 8A-8A) of an embodiment of the interface 606 that relies on reflectance (“reflectance interface”) to capture the light to be measured by the one or more spectrometers in the enclosure 604. The reflectance interface 606 is generally shaped to receive the fingers of a subject. Extending through a rear side thereof is the fiber optic cable 608, which comprises multiple fiber optic fibers. At least one of these fibers is the source fiber 806a, which carries light from the one or more spectrometers in the enclosure 604 to the interface 606, and at least one of these fibers is the return fiber 806b, which carries light after it has interacted with the subject's finger back to the spectrometers in the enclosure 604.


Sitting on a receiver for the subject's finger on the top side of the reflectance interface 606 is a reference puck 802. An example reference puck is manufactured from Spectralon TM reflectance material from Labsphere, Inc. The reference puck is used in place of the subject's finger when obtaining the dark reference sample and light reference sample, mentioned above, each time a subject's body part is measured.


In FIG. 8B, a light path 808 is shown that indicates how light travels when interacting with the subject's finger. Namely, the reflectance interface 606 comprises a first surface on which the reference puck 802 rests that is positioned to abut against a pad of the finger. The source and return fibers 806a,b are positioned from beneath the first surface and respectively transmit and receive light through the first surface. More particularly, as indicated by the light path 808, the light exits the source fiber, reflects off the finger, and returns to the return fibers 806b for transmission back to the spectrometers in the enclosure 604.


Referring now to FIGS. 9A and 9B, there are respectively shown perspective and sectional views (along line 9A-9A) of an embodiment of the interface 606 that relies on transmittance (“transmittance interface”) to capture the light to be measured by the one or more spectrometers in the enclosure 604. The transmittance interface 606 is generally shaped to receive the fingers of a subject. Extending through a front side thereof are the source and return fibers 806a,b. As with the reflectance interface 606, the reference puck 802 sits on the top side of the transmittance interface 606 and is used when obtaining the dark reference sample and the light reference sample.


In contrast to the reflectance interface 606, the finger receiver of the transmittance interface 606 comprises a first surface positioned to abut against a pad of the finger and a second surface positioned to abut against a tip of the finger. The source fiber 806a is positioned to direct the light from the lamp to the finger through the second surface and the return fibers 806b are positioned to receive the light transmitted through the finger through the first surface. The light path 808 indicates the direction in which light travels from the source fiber 806a through the finger and to the return fibers 806b.


Referring now to FIG. 23, there is depicted a top plan view of a combined device 2300 that collectively performs the functionality of the interface 606 and enclosure 604. More particularly, the device 2300 comprises the housing 604 within which are a processor, a computer readable medium, an FT-NIR spectrometer, and a lamp, as described above in respect of FIG. 6. A display 2302 is communicatively coupled to the processor and mounted to a top side of the housing 604; the display's 2302 top edge is visible in FIG. 23. Also mounted to a top side of the housing 604 is the device's 2300 interface 606. The interface 606, analogous to that used in a Bruker™ Tango™ FT-NIR spectrometer, is designed for use with samples held in container such as a cup, a Petri dish, or a vial. More particularly, the interface 606 comprises a platform 2304 within which is an opening 2306. In conventional usage, a container containing the sample is placed on the platform 2304, and optical fibers (not depicted) within the opening 2306 transmit light to and collect light from the sample through the container. The FT-NIR spectrometer within the housing 604 receives and processes the optical data. In contrast to this conventional usage, in at least some example embodiments an individual may place their finger directly on the opening 2306 and FT-NIR measurements may be obtained by directly measuring the individual in a manner analogous to that described above in FIGS. 8A and 8B and/or 9A and 9B.


Referring now to FIG. 10, there is shown a flow diagram 1000 of a method of processing spectral data using one or both of the transmittance and reflectance interfaces 606. In the present embodiment, the method is performed by the processor in the computer 602. In at least some other embodiments, the method may be alternatively performed by the processor in the enclosure 604, or collectively by the processors in the computer 602 and enclosure 604. More particularly, the processor performs input (spectral) data processing at block 1002. Once the input data is processed at block 1002, the processor performs outlier detection at block 1003 and executes a preprocessing pipeline at block 1004. During system optimization, executing the preprocessing pipeline comprises executing an iterative optimization pipeline at block 1010. FIG. 15 schematically represents the effect of the preprocessing pipeline executed at block 1004 in which input data in an array of approximately 2.5 k×1 k is reduced by virtue of preprocessing to approximate 100×1 k. After preprocessing, the processor applies a model 1006 to determine state of interest from the processed spectral data. Example implementations of blocks 1002, 1003, 1004, 1006, and 1010 are discussed in further detail in respect of FIGS. 11-21C, below.


In FIG. 11, the processor at block 1102 joins raw spectral data files together. Table 3, below, describes the spectral data received from the VIS-NIR and NIR spectrometers:









TABLE 3







Data from VIS-NIR and NIR Spectrometers












#
#


Name
Description
Scans
Columns





Light Reference
Light reference measurement used
~1-
  ~2k


(VIS-NIR
to calibrate the VIS-NIR
100


spectrometer)
spectrometer (~350 nm-1000 nm)


Dark Reference
Dark reference measurement used
~1-
  ~2k


(VIS-NIR
to calibrate the VIS-NIR
100


spectrometer)
spectrometer (~350 nm-1000 nm)


Light Sample
Light sample measurement using
~1-
  ~2k


(VIS-NIR
the VIS-NIR spectrometer
10k


spectrometer)
(~350 nm-1000 nm)


Dark Sample
Dark sample measurement using
~1-
  ~2k


(VIS-NIR
the VIS-NIR spectrometer
10k


spectrometer)
(~350 nm-1000 nm)


Light Reference
Light reference measurement used
~1-
~512


(NIR
to calibrate the NIR spectrometer
100


spectrometer)
(~900 nm-2500 nm)


Dark Reference
Dark reference measurement used
~1-
~512


(NIR
to calibrate the NIR spectrometer
100


spectrometer)
(~900 nm-2500 nm)


Light Sample
Light sample measurement using
~1-
~512


(NIR
the NIR spectrometer (~900 nm-
10k


spectrometer)
2500 nm)


Dark Sample
Dark sample measurement using
~1-
~512


(NIR
the NIR spectrometer (~900 nm-
10k


spectrometer)
2500 nm)









For example, in Table 3 the VIS-NIR spectrometer has a measurement range of ˜300 nn to 1000 nm and the NIR spectrometer has a measurement range of ˜900 nm to 2500 nm. However, accuracy and/or precision may decrease for one or both spectrometers near the end of the ranges. Consequently, for a spectrum spanning 400 nm to 2500 nm for example, the VIS-NIR spectrometer may be used to acquire measurements from 400 nm to up to 900 nm, while the NIR spectrometer may be used to acquire measurements from 900 nm up to 2500 nm.


Part of joining the spectral files comprises performing block 1104, in which the processor combines raw spectra data for the light reference sample and dark reference sample measurements. Subsequent to this the processor at block 1106 downsamples spectra data for the light sample and dark sample measurements. Following blocks 1104 and 1106, the processor at block 1108 merges spectral data and corrects for sensor bias. More particularly, the reference data is used to mitigate sensor variance across the different spectrometers. Each of the reference samples (both light and dark) are averaged together to reduce the effect of noise. More particularly, the processor applies the following when merging spectral data to correct for sensor bias:











log
10





R
L

-

R
D




S
L

-

S
D




+


log
10




IT
S


IT
R







(
1
)







where RL and RD are respectively the light reference sample and dark reference sample, SL and SD are respectively the light sample and dark sample, and ITS and ITR are respectively the integration time for the subject and reference.


Spectral data for the non-reference samples (both light and dark) contain temporal data that, in at least some embodiments, is retained. For example, in at least some embodiments the processor averages samples within a certain percentage of the maximum, median, and minimum values for regions that vary with time (e.g., a pulse). The end result is a matrix containing multiple entries for each wavelength across both spectrometers.



FIG. 12 is a flow diagram of a method of performing outlier detection at block 1003 of FIG. 10, which helps to correct for light scattering. At block 1202, the processor standardizes input data values to have zero mean and unit variance for each subject:











X
i

-

X

i

_

mean




X
std





(
2
)







and proceeds to block 1204 where a trained deep Auto-Encoder is used to identify good vs. bad samples (“good” is defined heuristically based on the fact that the visual representation of the spectra matches the expected shape based on all samples acquired). The processor compares the mean squared error [Error(x,y)] between the input data and the reconstructed Auto-Encoder output. The processor at block 1206 determines outliers as the samples where the mean squared error is greater than 2 times the total population's standard deviation; in at least some other embodiments, outliers may be determined to be samples a different number of standard deviations from the mean (e.g., lx). The deep Auto-Encoder is applied directly to spectral data and is accordingly independent of changes to state of interest; therefore, the Auto-Encoder only needs to be trained once for a particular configuration. After identifying the outliers, the processor at block 1208 returns a mean centered version (Xi−Xi_mean) of the original input spectral data with the outliers removed. While in at least some example embodiments the Auto-Encoder is applied to the raw data, in at least some other embodiments the Auto-Encoder is applied to mean centered or scaled input data (without the outliers having yet been removed) depending on performance and robustness needs.



FIG. 13 is a flow diagram of a method of executing an optimization pipeline, as described above at block 1010 of FIG. 10. The optimization pipeline iteratively tests various combinations of spectroscopic transformations to arrive at a particular, and optimally ideal, transformation or sequence of transformations for a particular target, exemplified by the subject's state of interest. FIG. 13 is performed iteratively during initial training of the system 600. The processor at block 1302 smooths/filters the data once outliers have been removed; transforms smoothed/filtered data at block 1304; and reduces the number of wavelengths represented in the transformed data at block 1306. These operations are described further below in respect of FIGS. 16-18.



FIG. 16 depicts a particular manner of executing the optimization pipeline of FIG. 13. In FIG. 16, the mean centered input data output from the outlier detection block 1003 is input to a configuration block 1604; the output of the configuration block 1604 is the input to a wavelength reduction block 1606; and the output of the wavelength reduction block 1606 is the input to block 1608 in which the processor applies partial least squares. Applying partial least squares is one example of performing a decomposition into latent space components; in other embodiments different latent space decompositions are possible, such as by applying a principal components analysis.



FIG. 17 is a block diagram depicting the configuration block 1604 of FIG. 16. The configuration block 1604 comprises first through fifth configuration sub-blocks 1702a-e, a standardization and scaling sub-block 1704, and a PLS sub-block 1402. The configuration sub-blocks 1702a-e are used by the processor to iteratively test different sequences of transforms to determine a preferred, and in some cases optimal, configuration of transforms for determining state of interest. Example transforms comprise standard normal variate (SNV), multiplicative scatter correction (MSC), L1 normalization (L1N), L2 normalization (L2N), Savitzky-Golay smoothing (SGS), convolutional smoothing (CS), and signal derivative (SD), which are spectra and row specific. Following application of the transforms by the sub-blocks 1702a-e, the processor applies standardization and scaling at the standardization and scaling sub-block 1704; and after standardization and scaling, the processor applies PLS to the standardized and scaled data using the PLS sub-block 1402. The output of the PLS sub-block 1402 is used when assessing the performance of the particular combination of transforms applied by the sub-blocks 1702a-e.


While five configuration sub-blocks 1702a-e are shown in FIG. 17, in at least some other embodiments (not depicted) a different number of sub-blocks may be used so that the processor can test more or fewer combinations of transforms.



FIG. 18 is an example of wavelength reduction performed by the wavelength reduction block 1606 of FIG. 16. The output of the configuration block 1604 is used as input to the wavelength reduction block 1606. In the wavelength reduction block 1606, the processor applies a genetic algorithm to reduce the number of wavelengths for further processing. As shown in FIG. 18, in this example the processor has deleted four continuous ranges of wavelengths from the spectral data.


After wavelength reduction, the processor again applies PLS to the spectral data to find an ideally optimal fit for the reduced wavelength data. Once the processor arrives at a suitable fit, the components generated as a result of applying PLS are extracted to be used in the model selection process.



FIG. 14 is a flow diagram depicting application of various models to determine state of interest. In FIG. 14, and as discussed further below, the processor may apply linear regression (LR) at block 1403 based on the PLS components output by block 1402 to arrive at the target 1408 directly; apply PLS in conjunction with a neural additive model (NAM) at block 1404; or apply a neural network (NN) at block 1406 to arrive at the target 1408.


Referring now to FIGS. 19A-C, there are respectively depicted three example models that leverage the preprocessing pipeline described above to arrive at the target 1408. In FIG. 19A, PLS and a logistic regression model is applied to classify the subject's state of interest. In FIG. 19B, PLS and NAM are applied to classify the subject's state of interest. NAMs extract non-linear relationships between 1, 2, or 3 input variables and the target 1408. In FIG. 19C, NAM is applied to automatically determine key wavelength contributions and transforms. NAM may be applied to standardized spectra, by treating each wavelength as an independent variable as opposed to using PLS-derived components, or to PLS-derived components.


Referring now to FIGS. 20A-C, there are respectively depicted three example models that may be used at block 1006. In FIG. 20A, PLS and an artificial deep neural network (DNN) are applied to the standardized spectral data to arrive at state of interest. In FIG. 20B, a DNN alone is applied to the standardized spectral data to arrive at the state of interest. And in FIG. 20C, to leverage application specific temporal data, a convolutional neural network (CNN) is applied to process the standardized spectral data. In FIGS. 19B, 19C, and 20A-C, “x” as an input represents standardized wavelengths while “c” represents PLS-derived components.


Referring now to FIGS. 21A-C, there is depicted a block diagram of an example system 2100 that may be used to process spectral data to arrive at a state of interest. In FIG. 21A, an example of the downsampling at block 1106 is the Mean block annotated with “(Reduce size)” in FIG. 21A. An example of the output of the input data processing block 1002 is the 1×2560 vector at the bottom of FIG. 21B. In FIG. 21C, the “PLS components” box shrinks vector size from 1×2560 to 1×10, with 10 being the number of PLS-derived components.


More particularly, in FIGS. 21A-C spectral data in the form of eight input files 2102. Starting with FIG. 21A, as discussed above the input files 2102 comprise four files 2104 from the VIS-NIR spectrometer (“VIS files 2104”) and four files 2106 from the NIR spectrometer (“NIR files 2106”). Also as discussed above, the VIS files 2104 comprise a dark and a light reference sample 2108a,b and a dark and a light sample of the subject 2108c,d representing visible light measurements; and the NIR files 2106 analogously comprise a dark and a light reference sample 2108e,f and a dark and a light sample of the subject 2108g,h representing NIR measurements. The mean of the dark and light reference samples 2108a,b,e,f from the VIR and NIR files 2104,2106 is obtained at blocks 2110a and 2110b respectively in order to reduce their size via downsampling as described above in respect of block 1106; the results of the downsampling are respectively two 1×2048 files and two 1×512 files. These two downsampled dark reference samples 2108a,b, the two dark samples 2108c,d, the two downsampled light reference samples 2108e,f, and the two light samples 2108g,h are then smoothed via convolution with 2D kernels at convolution blocks 2112a-d, respectively. Each of files in the pairs is subtracted from the other file of the pair, resulting in a VIS reference file 2115a, a VIS subject file 2115b, a NIR reference file 2115c, and a VIS reference file 2115d.


Moving from FIG. 21A to FIG. 21B, the VIS reference file 2115a, the VIS subject file 2115b, the NIR reference file 2115c, and the dark VIS reference file 2115d are processed using Equation (1) at blocks 2116a and 2118a (for the VIS files 2115a,b) and blocks 2116band 2118b (for the NIR files 2115c,d) to result in a 1500×2048 file generated from the VIS files 2115a,b and a 1500×512 file generated from the NIR files 2115c,d.


The 1500×2048 and 1500×512 files are downsampled at blocks 2120a and 2122a (for the 1500×2048 file) and at blocks 2120b and 2122b (for the 1500×512 file) to result in a 1×2048 VIS file and a 1×512 NIR file. This downsampling comprises identifying a number of local maxima within the data at blocks 2120a,b to reduce variance and then determining the mean of those maxima at blocks 2122a,b. The resulting 1×2048 VIS file and 1×512 NIR file are processed according to Equation (2) at blocks 2124a and 2126a (for the 1×2048 VIS file) and blocks 2124b and 2126b(for the 1×512 NIR file). The resulting 1×2048 and 1×512 files are concatenated at block 2128 to result in a single 1×2560 file 2129.


Moving from FIG. 21B to FIG. 21C, at block 2130 the PLS components of the 1×2560 file 2129 are determined as described in respect of block 1402. Ten components are identified as represented in a 1×10 file; while ten components are selected in the depicted embodiment, any number of different components may be selected in an alternative embodiment. This 1×10 file is normalized (e.g., between 0 and 1) at block 2134. The elements of the normalized 1×10 file are respectively input to first through tenth trained neural networks 2136a-j, such as those depicted in FIGS. 20A-C. As one specific example, each of the first through tenth trained neural networks 2136a-j may comprise a NAM. The outputs of the trained neural networks 2136a-1 are summed together to result in a 1×1 file that is used as input to an activation function, such as a sigmoid function, at block 2138. The output of the activation function is used to determine state of interest such as described in respect of FIG. 22 and COVID-19 below.


While FIGS. 21A-C depict an example embodiment, variations to this depicted embodiment are possible. For example, FIGS. 21A-C depict files storing data as vectors of various dimensions. In alternative embodiments, the dimensions of these vectors may vary as desired. For example, more or fewer than ten PLS components may be selected, in which case the size of the 1x10 file output from block 2130 would correspondingly vary, as would the number of neural networks used to respectively process those components.


As another example variation, FIG. 21A shows eight input files 2102 generated from eight measurements corresponding to two dark and two light samples from each of the VIS-NIR and NIR spectrometer. However, alternative embodiments may feature more or fewer than eight files 2102, and more or fewer than eight measurements. For example, data from multiple measurements may be stored in the form of a single file. As another example, data may be acquired using more or fewer measurements. In at least some alternative embodiments, either NIR or visible light is used, and consequently only four measurements would be made. And, a spectrometer that performs pre-processing itself such as an FT-NIR spectrometer like the Bruker™ Tango™ FT-NIR spectrometer, outputs spectral data analogous to the 1×2560 file 2129 of FIGS. 21B and 21C. Consequently, this spectral data may be immediately processed by a single neural network without being used to generate PCA or PLS components.


As another example variation, a single trained neural network may directly receive as input the processed spectral data as represented by the 1×2560 file 2129 in place of determining PLS components at block 2130 and subsequently processing those components using the first through tenth neural networks 2134a-j. Rather, that single trained neural network analyzes all of that processed spectral data as opposed to specific PCS components; a NAM with CNN activations may be used for this purpose, for example, or any of the networks shown in FIGS. 19A-C. While one advantage of processing only specific PCS components is the ability to determine which PCS components are responsible for influencing a state of interest, bypassing the use of PCS components may allow for state of interest to be determined based on more complex relationships not traceable to one or more specific PCS components.


Training the neural networks referred to above in respect of FIGS. 21A-C may be done using as training data the input files 2102 representing one or more spectral measurements of a subject paired with the subject's state of interest. For example, where state of interest is whether the subject has a disease such as COVID-19, the each pair of training data comprises one or more spectral readings of the subject, together with state of interest in the form of whether the subject has COVID-19. At inference, spectral measurements as described herein may be input as the input files 2102, and the trained neural networks 2134a-1 accordingly output whether the subject has COVID-19.


Referring now to FIG. 22, there is depicted an example process diagram 2200 in which a state model 2202 trained to perform the method depicted in FIG. 10 is used to directly process spectral data to determine a state of interest without directly referencing analytes and through a blood sample. In the example of FIG. 22, the state of interest is whether the subject has COVID-19, a disease caused by a coronavirus. Coronaviruses, however, range in severity from a common cold (mild) to COVID-19 (potentially lethal). Very high sensitivity to a state of interest corresponding to infection may lead to higher cases of the common cold being classified as COVID-19, while optimizing for specificity increases the likelihood of only true-positive COVID-19 cases being identified at the potential cost of more false negatives. The state model 2202 is adjustable in terms of sensitivity and specificity to adjust to situations in which it may be desired to identify anyone who has a coronavirus infection regardless of strain/lethality, and other situations in which the desire is solely to identify those with COVID-19.


For example, the target 1408 may be a binary state of interest (e.g., COVID-19 positive or negative), which can be represented by either a 0 or a 1. To adjust sensitivity and specificity, the threshold at which the state of interest is determined to be 0 or 1 is adjustable. For example, increasing sensitivity can be done by instructing the processor to lower the threshold at which a positive state of interest is identified (e.g., from 0.5 to 0.25), while increasing specificity can be done be instructing the processor to increase the threshold at which a positive state of interest is identified (e.g., from 0.5 to 0.75).


All citations are hereby incorporated by reference.


The embodiments have been described above with reference to flow, sequence, and block diagrams of methods, apparatuses, systems, and computer program products. In this regard, the depicted flow, sequence, and block diagrams illustrate the architecture, functionality, and operation of implementations of various embodiments. For instance, each block of the flow and block diagrams and operation in the sequence diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified action(s). In some alternative embodiments, the action(s) noted in that block or operation may occur out of the order noted in those figures. For example, two blocks or operations shown in succession may, in some embodiments, be executed substantially concurrently, or the blocks or operations may sometimes be executed in the reverse order, depending upon the functionality involved. Some specific examples of the foregoing have been noted above but those noted examples are not necessarily the only examples. Each block of the flow and block diagrams and operation of the sequence diagrams, and combinations of those blocks and operations, may be implemented by special purpose hardware-based systems that perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.


The controller(s) and processor(s) used in the foregoing embodiments may comprise, for example, a processing unit (such as a processor, microprocessor, or programmable logic controller) communicatively coupled to a non-transitory computer readable medium having stored on it program code for execution by the processing unit, microcontroller (which comprises both a processing unit and a non-transitory computer readable medium), field programmable gate array (FPGA), system-on-a-chip (SoC), an application-specific integrated circuit (ASIC), or an artificial intelligence accelerator. Examples of computer readable media are non-transitory and include disc-based media such as CD-ROMs and DVDs, magnetic media such as hard drives and other forms of magnetic disk storage, semiconductor based media such as flash media, random access memory (including DRAM and SRAM), and read only memory.


It is contemplated that any part of any aspect or embodiment discussed in this specification can be implemented or combined with any part of any other aspect or embodiment discussed in this specification.


In construing the claims, it is to be understood that the use of computer equipment, such as a processor, to implement the embodiments described herein is essential at least where the presence or use of that computer equipment is positively recited in the claims.


One or more example embodiments have been described by way of illustration only. This description is being presented for purposes of illustration and description, but is not intended to be exhaustive or limited to the form disclosed. It will be apparent to persons skilled in the art that a number of variations and modifications can be made without departing from the scope of the claims.

Claims
  • 1. A method comprising: (a) directing light at a body part of a subject such that the light passes through or is reflected by blood and interstitial fluid of the body part, wherein the light incident on the body part comprises a range of wavelengths from at least one of the near infrared and visible spectra;(b) measuring a spectrum of the light after the light has one or both of passed through and been reflected by the body part, wherein the spectrum comprises the range of wavelengths; and(c) determining whether the subj ect is in a physiological state of interest without direct reference to analytes of the subject, wherein the determining comprises using a trained machine learning model to process the measured spectrum and wherein the trained machine learning model is trained with reference spectra representative of the physiological state of interest.
  • 2. The method of claim 1, wherein the light incident on the body part comprises a range of wavelengths from both of the near infrared and visible spectra.
  • 3. The method of claim 1, wherein the spectrum is measured on the light that has passed through the body part.
  • 4. The method of claim 1, wherein the spectrum is measured on the light that has been through the body part and that has been reflected by the body part.
  • 5. The method of claim 4, wherein the measured spectrum comprises a light reference sample, a dark reference sample, a light sample of the subject, and a dark sample of the subject, and wherein the comparing comprises correcting for sensor bias using the light reference sample, the dark reference sample, the light sample of the subject, and the dark sample of the subject.
  • 6. The method of claim 1, further comprising, prior to using the trained machine learning model to process the measured spectrum, removing outliers from the measured spectrum and generating a mean centered version of the measured spectrum.
  • 7. The method of claim 6, further comprising, prior to using the trained machine learning model to process the measured spectrum: (a) applying multiple transforms to the mean centered version of the measured spectrum, wherein the transforms are selected from the group consisting of standard normal variate (SNV), multiplicative scatter correction (MSC), L1 normalization (L1N), L2 normalization (L2N), Savitzky-Golay smoothing (SGS), convolution smoothing (CS), and signal derivative (SD);(b) evaluating performance of each of the multiple transforms to the mean centered version of the measured spectrum; and(c) selecting, from a result of the evaluating, a transformed spectrum, wherein the transformed spectrum is a transformed version of the mean centered version of the measured spectrum.
  • 8. The method of claim 7, further comprising selecting at least one range of wavelengths that is a subset of a total wavelength range of the transformed spectrum, and wherein the machine learning model is used to process the transformed spectrum.
  • 9. The method of claim 8, further comprising decomposing the transformed spectrum into latent space components, and wherein processing the transformed spectrum using the trained machine learning model comprises processing the latent space components using respective instances of the machine learning model.
  • 10. The method of claim 1, wherein the machine learning model comprises a neural additive model.
  • 11. The method of claim 1, wherein the machine learning model comprises an artificial deep neural network.
  • 12. The method of claim 1, wherein the machine learning model comprises a convolutional neural network.
  • 13. The method of claim 1, , further comprising decomposing the transformed spectrum into latent space components, and wherein processing the measured spectrum using the trained machine learning model comprises processing the latent space components using respective instances of the machine learning model.
  • 14. The method of claim 9, wherein the latent space components are generated by applying partial least squares or a principal components analysis.
  • 15. The method of claim 1, wherein the determining comprises receiving a sensitivity target and a specificity target, and outputting the physiological state of interest in accordance with the sensitivity and specificity targets.
  • 16. The method of claim 1, wherein the physiological state of interest comprises whether the subject is infected with a virus.
  • 17. The method of claim 1, wherein the physiological state of interest comprises whether the subject has COVID-19.
  • 18. The method of claim 1, wherein the physiological state of interest comprises THC impairment.
  • 19. The method of claim 1, wherein the physiological state of interest comprises alcohol impairment.
  • 20. The method of claim 1, wherein the measuring is performed using a Fourier Transform Near Infrared spectrometer.
  • 21. The method of claim 20, wherein the spectrometer comprises a platform for receiving a sample container, and wherein the measuring is performed directly on a finger of an individual.
  • 22. A non-transitory computer readable medium having stored thereon computer program code that is executable by a processor and that, when executed by the processor, causes the processor to perform a method comprising: (a) directing light at a body part of a subject such that the light passes through or is reflected by blood and interstitial fluid of the body part, wherein the light incident on the body part comprises a range of wavelengths from at least one of the near infrared and visible spectra;(b) measuring a spectrum of the light after the light has one or both of passed through and been reflected by the body part, wherein the spectrum comprises the range of wavelengths; and(c) determining whether the subject is in a physiological state of interest without direct reference to analytes of the subject, wherein the determining comprises using a trained machine learning model to process the measured spectrum and wherein the trained machine learning model is trained with reference spectra representative of the physiological state of interest.
  • 23.-31. (canceled)
CROSS-REFERENCE TO RELATED APPLICATIONS

The present application claims priority to U.S. provisional patent application no. 63/149,199 filed on Feb. 12, 2021, and entitled, “Non-Invasive Determination of a Physiological State of Interest of a Subject” and to United States provisional patent application no. 63/218,477 filed on Jul. 5, 2021, and entitled, “Non-Invasive Determination of a Physiological State of Interest of a Subject”, the entireties of both of which are hereby incorporated by reference.

PCT Information
Filing Document Filing Date Country Kind
PCT/CA2022/050208 2/11/2022 WO
Provisional Applications (2)
Number Date Country
63149199 Feb 2021 US
63218477 Jul 2021 US