US20070258083A1 - Noise reduction for analyte detection systems - Google Patents
Noise reduction for analyte detection systems Download PDFInfo
- Publication number
- US20070258083A1 US20070258083A1 US11/734,261 US73426107A US2007258083A1 US 20070258083 A1 US20070258083 A1 US 20070258083A1 US 73426107 A US73426107 A US 73426107A US 2007258083 A1 US2007258083 A1 US 2007258083A1
- Authority
- US
- United States
- Prior art keywords
- sample
- signal
- detector
- fluid
- optical
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
- 239000012491 analyte Substances 0.000 title claims abstract description 120
- 238000001514 detection method Methods 0.000 title claims description 7
- 230000009467 reduction Effects 0.000 title description 4
- 238000005259 measurement Methods 0.000 claims abstract description 161
- 238000000034 method Methods 0.000 claims abstract description 81
- 238000012544 monitoring process Methods 0.000 claims abstract description 32
- 230000001360 synchronised effect Effects 0.000 claims abstract description 7
- 230000003287 optical effect Effects 0.000 claims description 104
- 230000000694 effects Effects 0.000 claims description 15
- 238000012545 processing Methods 0.000 claims description 7
- 230000000116 mitigating effect Effects 0.000 claims description 2
- 238000004422 calculation algorithm Methods 0.000 abstract description 38
- 210000001124 body fluid Anatomy 0.000 abstract description 5
- 239000000523 sample Substances 0.000 description 235
- 239000012530 fluid Substances 0.000 description 157
- 238000001228 spectrum Methods 0.000 description 94
- 239000000306 component Substances 0.000 description 70
- 210000004027 cell Anatomy 0.000 description 59
- 238000012806 monitoring device Methods 0.000 description 58
- 210000004369 blood Anatomy 0.000 description 57
- 239000008280 blood Substances 0.000 description 57
- WQZGKKKJIJFFOK-GASJEMHNSA-N Glucose Natural products OC[C@H]1OC(O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-GASJEMHNSA-N 0.000 description 52
- 239000008103 glucose Substances 0.000 description 51
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 49
- 239000011780 sodium chloride Substances 0.000 description 48
- 238000012360 testing method Methods 0.000 description 37
- 238000004458 analytical method Methods 0.000 description 23
- NOESYZHRGYRDHS-UHFFFAOYSA-N insulin Chemical compound N1C(=O)C(NC(=O)C(CCC(N)=O)NC(=O)C(CCC(O)=O)NC(=O)C(C(C)C)NC(=O)C(NC(=O)CN)C(C)CC)CSSCC(C(NC(CO)C(=O)NC(CC(C)C)C(=O)NC(CC=2C=CC(O)=CC=2)C(=O)NC(CCC(N)=O)C(=O)NC(CC(C)C)C(=O)NC(CCC(O)=O)C(=O)NC(CC(N)=O)C(=O)NC(CC=2C=CC(O)=CC=2)C(=O)NC(CSSCC(NC(=O)C(C(C)C)NC(=O)C(CC(C)C)NC(=O)C(CC=2C=CC(O)=CC=2)NC(=O)C(CC(C)C)NC(=O)C(C)NC(=O)C(CCC(O)=O)NC(=O)C(C(C)C)NC(=O)C(CC(C)C)NC(=O)C(CC=2NC=NC=2)NC(=O)C(CO)NC(=O)CNC2=O)C(=O)NCC(=O)NC(CCC(O)=O)C(=O)NC(CCCNC(N)=N)C(=O)NCC(=O)NC(CC=3C=CC=CC=3)C(=O)NC(CC=3C=CC=CC=3)C(=O)NC(CC=3C=CC(O)=CC=3)C(=O)NC(C(C)O)C(=O)N3C(CCC3)C(=O)NC(CCCCN)C(=O)NC(C)C(O)=O)C(=O)NC(CC(N)=O)C(O)=O)=O)NC(=O)C(C(C)CC)NC(=O)C(CO)NC(=O)C(C(C)O)NC(=O)C1CSSCC2NC(=O)C(CC(C)C)NC(=O)C(NC(=O)C(CCC(N)=O)NC(=O)C(CC(N)=O)NC(=O)C(NC(=O)C(N)CC=1C=CC=CC=1)C(C)C)CC1=CN=CN1 NOESYZHRGYRDHS-UHFFFAOYSA-N 0.000 description 20
- 210000002381 plasma Anatomy 0.000 description 19
- 230000006870 function Effects 0.000 description 18
- 230000003595 spectral effect Effects 0.000 description 18
- 238000001802 infusion Methods 0.000 description 17
- 239000000463 material Substances 0.000 description 17
- 238000004891 communication Methods 0.000 description 15
- 239000002699 waste material Substances 0.000 description 14
- 102000001554 Hemoglobins Human genes 0.000 description 13
- 108010054147 Hemoglobins Proteins 0.000 description 13
- 239000003599 detergent Substances 0.000 description 13
- 239000000126 substance Substances 0.000 description 13
- 239000003814 drug Substances 0.000 description 12
- 241000237858 Gastropoda Species 0.000 description 11
- 238000010521 absorption reaction Methods 0.000 description 11
- 238000000862 absorption spectrum Methods 0.000 description 11
- 230000005855 radiation Effects 0.000 description 11
- 239000000243 solution Substances 0.000 description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 11
- 102000004877 Insulin Human genes 0.000 description 10
- 108090001061 Insulin Proteins 0.000 description 10
- 239000003146 anticoagulant agent Substances 0.000 description 10
- 230000008901 benefit Effects 0.000 description 10
- 229940125396 insulin Drugs 0.000 description 10
- 238000000528 statistical test Methods 0.000 description 10
- HTTJABKRGRZYRN-UHFFFAOYSA-N Heparin Chemical compound OC1C(NC(=O)C)C(O)OC(COS(O)(=O)=O)C1OC1C(OS(O)(=O)=O)C(O)C(OC2C(C(OS(O)(=O)=O)C(OC3C(C(O)C(O)C(O3)C(O)=O)OS(O)(=O)=O)C(CO)O2)NS(O)(=O)=O)C(C(O)=O)O1 HTTJABKRGRZYRN-UHFFFAOYSA-N 0.000 description 9
- 229940079593 drug Drugs 0.000 description 9
- 229960002897 heparin Drugs 0.000 description 9
- 229920000669 heparin Polymers 0.000 description 9
- 239000011159 matrix material Substances 0.000 description 9
- 230000008569 process Effects 0.000 description 9
- 229940127219 anticoagulant drug Drugs 0.000 description 8
- 239000003978 infusion fluid Substances 0.000 description 8
- 230000035945 sensitivity Effects 0.000 description 8
- 230000001133 acceleration Effects 0.000 description 7
- 239000012503 blood component Substances 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 230000009471 action Effects 0.000 description 6
- 238000002834 transmittance Methods 0.000 description 6
- JVTAAEKCZFNVCJ-UHFFFAOYSA-M Lactate Chemical compound CC(O)C([O-])=O JVTAAEKCZFNVCJ-UHFFFAOYSA-M 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 5
- 238000004140 cleaning Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 5
- 238000005096 rolling process Methods 0.000 description 5
- 238000005070 sampling Methods 0.000 description 5
- 238000001914 filtration Methods 0.000 description 4
- 238000011010 flushing procedure Methods 0.000 description 4
- 208000015181 infectious disease Diseases 0.000 description 4
- 230000033001 locomotion Effects 0.000 description 4
- 238000011002 quantification Methods 0.000 description 4
- 230000002829 reductive effect Effects 0.000 description 4
- 238000004611 spectroscopical analysis Methods 0.000 description 4
- 238000012549 training Methods 0.000 description 4
- 238000012546 transfer Methods 0.000 description 4
- 206010053567 Coagulopathies Diseases 0.000 description 3
- 238000002835 absorbance Methods 0.000 description 3
- 230000035602 clotting Effects 0.000 description 3
- 230000003750 conditioning effect Effects 0.000 description 3
- 238000012790 confirmation Methods 0.000 description 3
- 238000011109 contamination Methods 0.000 description 3
- 230000001419 dependent effect Effects 0.000 description 3
- 238000010790 dilution Methods 0.000 description 3
- 239000012895 dilution Substances 0.000 description 3
- 238000011049 filling Methods 0.000 description 3
- 239000007789 gas Substances 0.000 description 3
- 230000036541 health Effects 0.000 description 3
- 238000003780 insertion Methods 0.000 description 3
- 230000037431 insertion Effects 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 230000010076 replication Effects 0.000 description 3
- 238000013179 statistical model Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 239000013598 vector Substances 0.000 description 3
- -1 “bubble”) sensors Substances 0.000 description 3
- RZVAJINKPMORJF-UHFFFAOYSA-N Acetaminophen Chemical compound CC(=O)NC1=CC=C(O)C=C1 RZVAJINKPMORJF-UHFFFAOYSA-N 0.000 description 2
- 229940127090 anticoagulant agent Drugs 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 2
- WQZGKKKJIJFFOK-VFUOTHLCSA-N beta-D-glucose Chemical compound OC[C@H]1O[C@@H](O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-VFUOTHLCSA-N 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- WUKWITHWXAAZEY-UHFFFAOYSA-L calcium difluoride Chemical compound [F-].[F-].[Ca+2] WUKWITHWXAAZEY-UHFFFAOYSA-L 0.000 description 2
- 229910001634 calcium fluoride Inorganic materials 0.000 description 2
- 238000005119 centrifugation Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 230000001276 controlling effect Effects 0.000 description 2
- 238000013500 data storage Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 230000002526 effect on cardiovascular system Effects 0.000 description 2
- 235000013305 food Nutrition 0.000 description 2
- 230000002641 glycemic effect Effects 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 238000001990 intravenous administration Methods 0.000 description 2
- 239000006101 laboratory sample Substances 0.000 description 2
- 230000000670 limiting effect Effects 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 238000002483 medication Methods 0.000 description 2
- 230000015654 memory Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000000843 powder Substances 0.000 description 2
- 102000004169 proteins and genes Human genes 0.000 description 2
- 108090000623 proteins and genes Proteins 0.000 description 2
- 239000013558 reference substance Substances 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 230000000717 retained effect Effects 0.000 description 2
- 238000007789 sealing Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 238000009987 spinning Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- REAWNMHCBIUKLZ-ZYHUDNBSSA-N (3R,4R)-4-(hydroxymethyl)-3-(6-methylheptanoyl)oxolan-2-one Chemical compound CC(C)CCCCC(=O)[C@H]1[C@H](CO)COC1=O REAWNMHCBIUKLZ-ZYHUDNBSSA-N 0.000 description 1
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 1
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 1
- 206010034701 Peroneal nerve palsy Diseases 0.000 description 1
- 208000001647 Renal Insufficiency Diseases 0.000 description 1
- XSQUKJJJFZCRTK-UHFFFAOYSA-N Urea Chemical compound NC(N)=O XSQUKJJJFZCRTK-UHFFFAOYSA-N 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 239000003570 air Substances 0.000 description 1
- 230000000844 anti-bacterial effect Effects 0.000 description 1
- 230000004872 arterial blood pressure Effects 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 239000003181 biological factor Substances 0.000 description 1
- 230000031018 biological processes and functions Effects 0.000 description 1
- BJQHLKABXJIVAM-UHFFFAOYSA-N bis(2-ethylhexyl) phthalate Chemical compound CCCCC(CC)COC(=O)C1=CC=CC=C1C(=O)OCC(CC)CCCC BJQHLKABXJIVAM-UHFFFAOYSA-N 0.000 description 1
- 210000000601 blood cell Anatomy 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- BPKIGYQJPYCAOW-FFJTTWKXSA-I calcium;potassium;disodium;(2s)-2-hydroxypropanoate;dichloride;dihydroxide;hydrate Chemical compound O.[OH-].[OH-].[Na+].[Na+].[Cl-].[Cl-].[K+].[Ca+2].C[C@H](O)C([O-])=O BPKIGYQJPYCAOW-FFJTTWKXSA-I 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000004202 carbamide Substances 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 238000004040 coloring Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 239000013256 coordination polymer Substances 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 239000008121 dextrose Substances 0.000 description 1
- 229910001873 dinitrogen Inorganic materials 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000037406 food intake Effects 0.000 description 1
- 238000005534 hematocrit Methods 0.000 description 1
- 235000008216 herbs Nutrition 0.000 description 1
- 239000002117 illicit drug Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 201000006370 kidney failure Diseases 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 239000004816 latex Substances 0.000 description 1
- 229920000126 latex Polymers 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000007257 malfunction Effects 0.000 description 1
- 229940050561 matrix product Drugs 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 239000005416 organic matter Substances 0.000 description 1
- 229960005489 paracetamol Drugs 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000036470 plasma concentration Effects 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 230000001698 pyrogenic effect Effects 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 239000010409 thin film Substances 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/14532—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue for measuring glucose, e.g. by tissue impedance measurement
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/1455—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
- A61B5/1459—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters invasive, e.g. introduced into the body by a catheter
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/27—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands using photo-electric detection ; circuits for computing concentration
- G01N21/274—Calibration, base line adjustment, drift correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/01—Arrangements or apparatus for facilitating the optical investigation
- G01N21/03—Cuvette constructions
- G01N21/07—Centrifugal type cuvettes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
Definitions
- Some embodiments disclosed herein relate to methods and apparatus for determining the concentration of an analyte in a sample, such as an analyte in a sample of bodily fluid, as well as methods and apparatus which can be used to support the making of such determinations.
- Some embodiments use a synchronous demodulator and digital filter to reduce microphonic signal content. Some embodiments monitor the microphonic signal content and “holds off” on making a measurement until vibrations subside. In some embodiments, monitoring is performed using an accelerometer or other vibration sensor. In another embodiment, an algorithm is used to examine the detector output signal and detect excessive microphonic components thereby eliminating the need for the accelerometers.
- a method for determining the concentration of an analyte in a sample comprises: providing an optical detector signal from an optical detector, the signal having information relating to the concentration of an analyte in a sample; passing the optical detector signal through a first demodulator; providing a vibration sensor signal, the vibration signal having information relating to the vibration of the optical detector; passing the vibration sensor signal through a second demodulator; providing a threshold value for the vibration sensor signal that is calculated to correspond to an accuracy parameter of a system concentration output; and controlling measurement by the optical detector based on whether or not the threshold value is exceeded.
- Some embodiments comprise a system for improving accuracy of a mobile analyte concentration measurement apparatus.
- the system can comprise: a sample detector configured to provide a detector output signal; a first signal conditioner configured to receive the detector output signal; a vibration sensor mounted to detect vibration of the sample detector and configured to provide a sensor output signal; a second signal conditioner configured to receive the sensor output signal; and a controller configured to communicate with the sample detector and the vibration sensor and prevent the sample detector from detecting when the sensor output signal exceeds a threshold value.
- Some embodiments comprise a method of mitigating the effects of vibration on an optical analyte detection system measurements.
- the method can comprise: monitoring microphonic effects to determine when a threshold is exceeded; and automatically delaying optical measurement until after the threshold is no longer exceeded.
- monitoring microphonic effects comprises processing a detector output signal.
- monitoring microphonic effects comprises monitoring output from a vibration sensor that is physically associated with a detector.
- monitoring output from a vibration sensor further comprises monitoring output from an accelerometer.
- FIG. 1 shows an embodiment of an apparatus for withdrawing and analyzing fluid samples
- FIG. 2 illustrates how various other devices can be supported on or near an embodiment of apparatus illustrated in FIG. 1 ;
- FIG. 3 illustrates an embodiment of the apparatus in FIG. 1 connected to a patient
- FIG. 4 is a block diagram of an embodiment of a system for withdrawing and analyzing fluid samples
- FIG. 5 schematically illustrates an embodiment of a fluid system within a system for withdrawing and analyzing fluid samples
- FIG. 6 is an oblique schematic depiction of an embodiment of a modular monitoring device
- FIG. 7 shows a cut-away side view of an embodiment of a monitoring device
- FIG. 8 illustrates an embodiment of a disposable cartridge that can interface with a fluid system
- FIG. 9 schematically illustrates an embodiment of an optical system that comprises a spectroscopic analyzer adapted to measure spectra of a fluid sample
- FIG. 10 is a flowchart that schematically illustrates an embodiment of a spectroscopic method for determining the concentration of an analyte of interest in a fluid sample
- FIG. 11 is a flowchart that schematically illustrates an embodiment of a method for estimating the concentration of an analyte in the presence of interferents
- FIG. 12 is a flowchart that schematically illustrates an embodiment of a method for performing a statistical comparison of the absorption spectrum of a sample with the spectrum of a sample population and combinations of individual library interferent spectra;
- FIG. 13 is a flowchart that schematically illustrates an example embodiment of a method for estimating analyte concentration in the presence of the possible interferents
- FIGS. 14A and 14B schematically illustrate the visual appearance of embodiments of a user interface for a system for withdrawing and analyzing fluid samples
- FIG. 15 schematically depicts various components and/or aspects of a patient monitoring system and the relationships among the components and/or aspects;
- FIG. 16 is a block diagram of a system for reducing noise.
- FIG. 17 is a diagram of a lock-in amplifier system that can help reduce noise.
- FIG. 18 shows example signal levels of a signal from an input device.
- FIG. 19 shows a signal output from a detector before that signal reaches a lock-in amplifier.
- FIG. 20 shows signals after they have passed through a lock-in amplifier system such as the one depicted in FIG. 17 .
- FIG. 21 shows a block diagram of a system for dealing with noise in an analyte detection environment.
- FIG. 22 illustrates an output can be taken from a detector channel and fed directly into an accelerometer demodulation multiplier.
- FIG. 1 shows an embodiment of an apparatus 100 for withdrawing and analyzing fluid samples.
- the apparatus 100 includes a monitoring device 102 .
- the monitoring device 102 can be an “OptiScanner®,” available from OptiScan Biomedical Corporation of Hayward, Calif.
- the apparatus 100 can measure one or more physiological parameters, such as the concentration of one or more substance(s) in a sample fluid.
- the sample fluid can be, for example, whole blood from a patient 302 (see, e.g., FIG. 3 ).
- the apparatus 100 can also deliver an infusion fluid to the patient 302 .
- the monitoring device 102 includes a display 104 such as, for example, a touch-sensitive liquid crystal display.
- the display 104 can provide an interface that includes alerts, indicators, charts, and/or soft buttons.
- the device 102 also can include one or more inputs and/or outputs 106 that provide connectivity.
- the device 102 is mounted on a stand 108 .
- the stand 108 can be easily moved and includes one or more poles 110 and/or hooks 112 .
- the poles 110 and hooks 112 can be configured to accommodate other medical implements, including, for example, infusion pumps, saline bags, arterial pressure sensors, other monitors and medical devices, and so forth.
- FIG. 2 illustrates how various other devices can be supported on or near the apparatus 100 illustrated in FIG. 1 .
- the poles 110 of the stand 108 can be configured (e.g., of sufficient size and strength) to accommodate multiple devices 202 , 204 , 206 .
- one or more COLLEAGUE® volumetric infusion pumps available from Baxter International Inc. of Deerfield, Ill. can be accommodated.
- one or more Alaris® PC units available from Cardinal Health, Inc. of Dublin, Ohio can be accommodated.
- various other medical devices can be integrated with the disclosed monitoring device 102 such that multiple devices function in concert for the benefit of one or multiple patients without the devices interfering with each other.
- FIG. 3 illustrates the apparatus 100 of FIG. 1 as it can be connected to a patient 302 .
- the monitoring device 102 can be used to determine the concentration of one or more substances in a sample fluid.
- the sample fluid can come from a fluid container in a laboratory setting, or it can come from a patient 302 , as illustrated here. In some preferred embodiments, the sample fluid is whole blood.
- the monitoring device 102 can also deliver an infusion fluid to the patient 302 .
- An infusion fluid container 304 e.g., a saline bag
- infusion fluid e.g., saline and/or medication
- the monitoring device 102 can be in fluid communication with both the container 304 and the sample fluid source (e.g., the patient 302 ), through tubes 306 .
- the infusion fluid can comprise any combination of fluids and/or chemicals. Some advantageous examples include (but are not limited to): water, saline, dextrose, lactated Ringer's solution, drugs, and insulin.
- the illustrated monitoring device 102 allows the infusion fluid to pass to the patient 302 and/or uses the infusion fluid itself (e.g., as a flushing fluid or a standard with known optical properties, as discussed further below). In some embodiments, the monitoring device 102 may not employ infusion fluid. The monitoring device 102 may thus draw samples without delivering any additional fluid to the patient 302 .
- the monitoring device 102 can include, but is not limited to, fluid handling and analysis apparatuses, connectors, passageways, catheters, tubing, fluid control elements, valves, pumps, fluid sensors, pressure sensors, temperature sensors, hematocrit sensors, hemoglobin sensors, colorimetric sensors, gas (e.g., “bubble”) sensors, fluid conditioning elements, gas injectors, gas filters, blood plasma separators, and/or communication devices (e.g., wireless devices) to permit the transfer of information within the monitoring device 102 or between the monitoring device 102 and a network.
- fluid handling and analysis apparatuses e.g., connectors, passageways, catheters, tubing, fluid control elements, valves, pumps, fluid sensors, pressure sensors, temperature sensors, hematocrit sensors, hemoglobin sensors, colorimetric sensors, gas (e.g., “bubble”) sensors, fluid conditioning elements, gas injectors, gas filters, blood plasma separators, and/or communication devices (e.g., wireless devices) to permit the transfer of information within the monitoring
- one or more components of the apparatus 100 can be located at another facility, room, or other suitable remote location.
- One or more components of the monitoring device 102 can communicate with one or more other components of the monitoring device 102 (or with other devices) by communication interface(s) such as, but not limited to, optical interfaces, electrical interfaces, and/or wireless interfaces. These interfaces can be part of a local network, internet, wireless network, or other suitable networks.
- FIG. 4 is a block diagram of a system 400 for withdrawing and analyzing fluid samples.
- the monitoring device 102 can comprise such a system.
- the system 400 includes a fluid source 402 connected to a fluid system 404 .
- the fluid system 404 prepares fluid samples that are analyzed by an optical system 412 .
- the system 400 includes a display controller 414 and an algorithm processor 416 that assist in fluid sample analysis and presentation of data.
- the sampling and analysis system 400 is a mobile point of care apparatus that monitors physiological parameters such as, for example, blood glucose concentration.
- Tubes and connectors within the system 400 can be coated with an antibacterial coating to reduce the risk of infection.
- Connectors between at least some components of the system 400 can include a self-sealing valve, such as a spring valve, in order to reduce the risk of contact between port openings and fluids, and to guard against fluid escaping from the system.
- the sampling and analysis system 400 includes a fluid source 402 that contains fluid to be sampled.
- the fluid system 404 of the sampling and analysis system 400 is connected to a fluid source 402 from which fluid samples can be drawn.
- the fluid source 402 can be, for example, a patient's blood vessel such as a vein or an artery, a container such as a decanter or a tube, or any other corporeal or extracorporeal fluid source.
- the fluid to be sampled can be, for example, blood, plasma, or another bodily fluid.
- the fluid system 404 withdraws a sample of fluid from the fluid source 402 for analysis, centrifuges at least a portion of the sample, and prepares at least a portion of the sample for analysis by an optical sensor such as a spectrophotometer. In some embodiments, at least a portion of the sample is returned to the fluid source 402 . At least some of the sample, such as portions of the sample that are mixed with other materials or portions that are otherwise altered during the sampling and analysis process, can also be placed in a waste bladder. The waste bladder can be integrated within the fluid system 404 or supplied by a user of system 400 .
- the fluid system 404 can also be connected to a saline source, a detergent source, and/or an anticoagulant source, each of which can be supplied by a user or integrated within fluid system 404 .
- Components of the fluid system 404 can be modularized into one or more non-disposable, disposable, and/or replaceable subsystems. In the embodiment shown in FIG. 4 , components of the fluid system 404 are separated into a non-disposable subsystem 406 , a first disposable subsystem 408 , and a second disposable subsystem 410 .
- the non-disposable subsystem 406 can include components that do not generally require regular replacement during the useful lifetime of the system 400 .
- the non-disposable subsystem 406 of the fluid system 404 includes one or more reusable valves and sensors.
- the non-disposable subsystem 406 can include one or more pinch valves (or non-disposable portions thereof), ultrasonic bubble sensors, non-contact pressure sensors, and optical blood dilution sensors.
- the non-disposable subsystem 406 can also include one or more pumps (or non-disposable portions thereof).
- the components of the non-disposable subsystem 406 are not directly exposed to fluids and/or are not readily susceptible to contamination.
- First and second disposable subsystems 408 , 410 can include components that are regularly replaced under certain circumstances in order to facilitate the operation of the system 400 .
- the first disposable subsystem 408 can be replaced after a certain period of use, such as a few days, has elapsed. Replacement may be necessary, for example, when a bladder within the first disposable subsystem 408 is filled to capacity. Such replacement may mitigate fluid system performance degradation associated with and/or contamination wear on system components.
- the first disposable subsystem 408 includes components that may contact fluids such as patient blood, saline, flushing solutions, anticoagulants, and/or detergent solutions.
- the first disposable subsystem 408 can include one or more tubes, fittings, cleaner pouches and/or waste bladders.
- the components of the first disposable subsystem 408 can be sterilized in order to decrease the risk of infection and can be configured to be easily replaceable.
- the second disposable subsystem 410 can be designed to be replaced under certain circumstances.
- the second disposable subsystem 410 can be replaced when the patient being monitored by the system 400 is changed.
- the components of the second disposable subsystem 410 may not need replacement at the same intervals as the components of the first disposable subsystem 408 .
- the second disposable subsystem 410 can include a flow cell and/or at least some components of a centrifuge, components that may not become filled or quickly worn during operation of the system 400 .
- Replacement of the second disposable subsystem 410 can decrease or eliminate the risk of transferring fluids from one patient to another during operation of the system 400 , enhance the measurement performance of system 400 , and/or reduce the risk of contamination or infection.
- the flow cell of the second disposable subsystem 410 receives the sample obtained from the fluid source 402 via the fluidics of the first disposable subsystem 408 .
- the flow cell is a container that can hold fluid for the centrifuge and provide a window to the sample for analysis by a spectrometer.
- the flow cell includes windows that are made of a material that is substantially transparent to electromagnetic radiation in the mid-infrared range of the spectrum.
- the flow cell windows can be made of calcium fluoride.
- An injector can provide a fluidic connection between the first disposable subsystem 408 and the flow cell.
- the injector can be removed from the flow cell to allow for free spinning of the flow cell during centrifugation.
- the components of the sample are separated by centrifuging at a high speed for a period of time before measurements are performed by the optical system 412 .
- a blood sample can be centrifuged at 7200 RPM for 2 minutes in order to separate plasma from other blood components for analysis. Separation of a sample into the components can permit measurement of solute (e.g., glucose) concentration in plasma, for example, without interference from other blood components.
- solute e.g., glucose
- This kind of post-separation measurement (sometimes referred to as a “direct measurement”) has advantages over a solute measurement taken from whole blood because the proportions of plasma to other components need not be known or estimated in order to infer plasma glucose concentration.
- An anticoagulant such as, for example, heparin can be added to the sample before centrifugation to prevent clotting.
- the fluid system 404 can be used with a variety of anticoagulants, including anticoagulants supplied by a hospital or other user of the monitoring system 400 .
- a detergent solution formed by mixing detergent powder from a pouch connected to the fluid system 404 with saline can be used to periodically clean residual protein and other sample remnants from one or more components of the fluid system 404 , such as the flow cell. Sample fluid to which anticoagulant has been added and used detergent solution can be transferred into the waste bladder.
- the system 400 shown in FIG. 4 includes an optical system 412 that can measure optical properties (e.g., transmission) of a fluid sample (or a portion thereof).
- the optical system 412 measures transmission in the mid-infrared range of the spectrum.
- the optical system 412 includes a spectrometer that measures the transmission of broadband infrared light through a portion of a flow cell filled with fluid. The spectrometer need not come in direct contact with the sample.
- the term “flow cell” is a broad term that carries its ordinary meaning as an object that can provide a place for fluid. The fluid can enter the flow cell by flowing.
- the optical system 412 includes a filter wheel that contains one or more filters. In some embodiments, twenty-five filters are mounted on the filter wheel.
- the optical system 412 includes a light source that passes light through a filter and the flow cell to a detector. In some embodiments, a stepper motor moves the filter wheel in order to position a selected filter in the path of the light.
- An optical encoder can also be used to finely position one or more filters.
- the system 400 shown in FIG. 4 includes a display controller 414 that provides for communication of information to a user of the system 400 .
- the display controller 414 can include a display processor that controls or produces an interface to communicate information to the user.
- the display controller 414 can include a display screen.
- One or more parameters such as, for example, blood glucose concentration, system 400 operating parameters, and/or other operating parameters can be displayed on a monitor (not shown) associated with the system 400 .
- An example of one way such information can be displayed is shown in FIGS. 14A and 14B .
- the display controller 414 can communicate measured physiological parameters and/or operating parameters to a computer system over a communications connection.
- the system 400 shown in FIG. 4 includes an algorithm processor 416 that can receive optical density (OD) values (or other analog or digital optical data) from the optical system 412 .
- the algorithm processor 416 calculates one or more physiological parameters by adjusting the coefficients of a model, if necessary, and computing the physiological parameters using an equation having the adjusted coefficients.
- the algorithm processor 416 , the display controller 414 , and any embedded controllers within system 400 can be connected to one another with a communications bus.
- FIG. 5 schematically illustrates a fluid system 510 .
- T# indicates a section of tubing.
- C# indicates a connector that joins multiple tubing sections.
- V# indicates a valve.
- BS # indicates a bubble sensor or ultrasonic air detector.
- N# indicates a needle (e.g., a needle that injects sample into a flow cell).
- PS# indicates a pressure sensor (e.g., a reusable pressure sensor).
- Pump# indicates a fluid pump (e.g., a syringe pump with a disposable body and reusable drive).
- Hb 12 indicates a sensor for hemoglobin (e.g., a dilution sensor that can detect hemoglobin optically).
- most lines can be filled with saline that can be introduced into the system through the tubes 514 and 516 , and which can come from an infusion pump 518 and/or a saline bag 520 .
- the infusion pump 518 and the saline bag 520 can be provided separately from the system 510 .
- a hospital can use existing saline bags and infusion pumps to interface with the described system.
- the valve 521 can be open to allow saline to flow into the tube 512 (T 1 ).
- a first pump 522 draws sample fluid to be analyzed (e.g. blood) from a fluid source (e.g., a laboratory sample container, a living patient, etc.) up into the patient tube 512 (T 1 ), through the open valve F 23 (V 0 ), through the first connector 524 (C 1 ), past the hemoglobin sensor 526 (Hb 12 ), and into the looped tube 528 (T 4 ).
- a fluid source e.g., a laboratory sample container, a living patient, etc.
- the valve 529 (V 7 a ) is open to fluid flow, but the valves 531 (V 1 a ) and 533 (V 3 a ) can be closed and therefore block (or substantially block) fluid flow.
- the lines are filled with saline and the hemoglobin (Hb) level is zero.
- the tubes that are filled with saline are in fluid communication with the a sample source (not shown).
- the sample source can be the vessels of a living human or a pool of liquid in a laboratory sample container, for example.
- the first pump 522 draws a relatively continuous column of fluid that first comprises generally nondiluted saline, then a mixture of saline and sample fluid (e.g., blood), and then eventually nondiluted sample fluid.
- the sample fluid is blood.
- the hemoglobin sensor 526 detects the level of Hemoglobin in the sample fluid. As blood starts to arrive at the hemoglobin sensor 526 (Hb 12 ), the hemoglobin level rises. When the hemoglobin level reaches a preset value (e.g., which can occur after a draw of approximately 2 mL depending on the size of the catheter used) there is a nondiluted sample present at the first connector 524 (C 1 ).
- a nondiluted sample can be, for example, a blood sample that is not diluted with saline solution, but instead has the characteristics of the rest of the blood flowing through a patient's body.
- a loop of tubing 530 (e.g., a 1-mL loop) can be advantageously positioned as illustrated to help insure that undiluted fluid (e.g., undiluted blood) is present at the first connector 524 (C 1 ) when the hemoglobin sensor 526 registers that the preset Hb threshold is crossed.
- the loop of tubing 530 provides additional length to the tube 528 (T 4 ) to make it less likely that the portion of the fluid column in the tubing at the first connector 524 (C 1 ) has advanced all the way past the mixture of saline and sample fluid, and the nondiluted blood portion of that fluid has reached the first connector 524 (C 1 ).
- a second pump 532 draws four “slugs” of blood into the tubing 534 (T 3 ).
- the term “slug” refers to a continuous column of fluid. Slugs can be separated from one another by injecting (or sucking in) small amounts of air to create bubbles at intervals in the tube. In the illustrated embodiment, blood slugs are alternated with air bubbles by maintaining the valve 523 (V 0 ) closed, maintaining the valve 533 (V 3 a ) open, and alternately closing and opening the valves 529 (V 7 a ) and 531 (V 1 a ) such that one is closed while the other one is open.
- four blood slugs are created.
- the first three blood slugs are approximately 15 ⁇ L and the fourth is approximately 35 ⁇ L.
- heparin can be inserted into each slug.
- a heparin vial 538 e.g., an insertable vial provided independently by the user of the system 510
- a shuttle valve 541 can connect to both the tube 540 and the tube 534 (T 3 ).
- the valve can open the tube 540 to a suction force (e.g., created by the pump 532 ), allowing heparin to be drawn from the vial 538 into the valve 541 .
- the shuttle valve 541 can slide the heparin over into fluid communication with the tube 534 .
- the shuttle valve 541 can then return to its previous position.
- heparin can be shuttled from the tube 540 into the tube 534 (T 3 ) such that each blood slug contains a precisely controlled amount of heparin.
- the first pump 522 pushes the blood out of the tube 528 (T 4 ) and back to the patient by opening the valve 523 (V 0 ), closing the valves 531 (V 1 a ) and 533 (V 3 a ), and keeping the valve 529 (V 7 a ) open.
- the tube 528 (T 4 ) is preferably flushed with approximately 2 mL of saline.
- less than two minutes elapses between the time that blood is drawn from the patient and the time that the blood is returned to the patient after formation of the blood slugs.
- the four slugs are pushed up the tube 534 (T 3 ), through the second connector 546 (C 2 ), and into the flow cell 548 , which can be located on the centrifuge wheel 550 .
- the bubble sensor 552 (BS 14 ) can identify the fourth slug by identifying and counting how many air bubbles (or inter-slug spaces) pass by the sensor.
- the fourth slug can be identified, and the pump 522 can stop forcing the fluid column through the tube 534 so that the fourth slug remains within the flow cell 548 .
- the first three blood slugs can serve to flush any residual saline out the flow cell 548 .
- the three leading slugs can be deposited in the waste bladder 554 by passing through the tube F 56 (T 6 ) and through the valve 557 (V 4 a ).
- the fourth blood slug is centrifuged for two minutes at 7200 RPM. This separates the whole blood into its components, isolates the plasma, and positions the plasma in the flow cell 548 for measurement.
- the centrifuge 550 is stopped with the flow cell 548 in a beam of radiation (not shown) for analysis.
- the radiation, a detector, and logic can be used to analyze the plasma spectroscopically (e.g., for glucose and/or lactate concentration).
- the second pump 532 flushes the flow cell 548 and sends the flushed contents to the waste bladder 554 .
- This flush can be done with a cleaning solution from the terg tank 558 .
- the second pump 532 is in fluid communication with the tube 560 (T 9 ) and the terg tank 558 because the valve 559 (V 7 b ) is open.
- the second pump 532 forces cleaning solution from the terg tank 558 through the open valve 561 and the tube 562 (T 7 ) when the valve 559 is open.
- the cleaning flush can pass through the flow cell 548 , through the second connector 546 , through the tube 564 (T 5 ) and the open valve 563 (V 2 b ), and into the waste bladder 554 . Following this flush,
- the first pump 522 (pump # 1 ) can flush the cleaning solution out of the flow cell 548 using saline in drawn from the saline bag 520 .
- This flush pushes saline through the tube 528 (T 4 ), the tube 534 (T 3 ), the flow cell 548 , and the tube 556 (T 6 ).
- valves are open for this flush: 529 (V 7 a ), 533 (V 3 a ), 557 (V 4 a ), and the following valves are closed: 542 (PV 1 ), 523 (V 0 ), 531 (V 1 a ), 566 (V 3 b ), 563 (V 2 b ), and 561 (V 4 b ).
- a low flow of saline (e.g., 1-5 mL/hr) is preferably moved through the patient tube 512 (T 1 ) and into the patient to keep the patient's vessel open (e.g., to establish a keep vessel open, or “KVO” flow).
- the source of this KVO flow can be the infusion pump 518 , the third pump 568 (pump # 3 ), or the first pump 522 (pump # 1 ).
- the infusion pump 518 can run continuously throughout the measurement cycle described above. This continuous flow can advantageously avoid any alarms that may be triggered if the infusion pump 518 senses that the flow has stopped or changed in some other way.
- the third pump 568 when the valve 521 closes to allow pump 522 (pump # 1 ) to withdraw fluid from a fluid source (e.g., a patient), the third pump 568 (pump # 3 ) can withdraw fluid through the connector 570 , thus allowing the infusion pump 518 to continue pumping normally as if the fluid path was not blocked by the valve 521 . If the measurement cycle is about two minutes long, this withdrawal by the third pump 568 can continue for approximately two minutes.
- the third pump 568 can reverse and insert the saline back into the system at a low flow rate.
- the time between measurement cycles is longer than the measurement cycle itself (e.g., longer than two minutes). Accordingly, the third pump 568 can insert fluid back into the system at a lower rate than it withdrew that fluid. This can help prevent an alarm by the infusion pump.
- FIG. 6 is an oblique schematic depiction of a modular monitoring device 600 .
- the modular monitoring device 600 includes a body portion 602 having receptacles 604 , 606 .
- the receptacles 604 , 606 include connectors with which disposable cassettes 610 , 612 can interface.
- portions of the fluidic system that directly contact fluid are incorporated into one or more removable cassettes.
- a first cassette 610 can be used to store at least a portion of the fluid system 510 described previously, including portions that contact sample fluids, saline, detergent solution, and/or anticoagulant.
- a non-disposable fluidics subsystem 608 is disposed within the body portion 602 of the monitoring device 600 .
- the first cassette 610 can include one or more openings that allow portions of the non-disposable fluidics subsystem 608 to interface with the cassette 610 .
- the non-disposable fluidics subsystem 608 can include one or more pinch valves that are designed to extend through such openings to engage one or more sections of tubing. When the first cassette 610 is inserted into a corresponding first receptacle 604 , actuation of the pinch valves can selectively close sections of tubing within the cassette.
- the non-disposable fluidics subsystem 608 can also include one or more sensors that interface with connectors, tubing sections, or pumps located within the first cassette 610 .
- the monitoring device 600 includes an optical system 614 disposed within the body portion 602 .
- the optical system 614 can include a light source and a detector that are adapted to perform measurements on fluids within a flow cell.
- the flow cell is disposed within a second cassette 612 .
- the second cassette 612 can include an optical window through which the optical system 614 can shine radiation for measuring properties of a fluid in the flow cell when the cassette is inserted into a corresponding second receptacle 606 .
- the optical system 614 can include other components (some of which may interface with the second cassette 612 ) such as, for example, a power supply, a centrifuge motor, a filter wheel, and/or a beam splitter.
- the first cassette 610 and the second cassette 612 are adapted to be in fluid communication with each other.
- the first cassette 610 can include a retractable injector that injects fluids into a flow cell disposed in the second cassette 612 .
- the injector can be retracted to allow the centrifuge to rotate the flow cell freely.
- a fluid communication path can be provided by components disposed within the body portion 602 of the monitoring device 600 .
- the body portion 602 of the monitoring device 600 can also include one or more connectors for an external battery (not shown).
- the external battery can serve as a backup emergency power source in the event that a primary emergency power source such as, for example, an internal battery (not shown) is exhausted.
- FIG. 7 shows a cut-away side view of a monitoring device 700 (which can correspond, for example, to the device 102 shown in FIG. 1 ).
- the device 700 includes a casing 702 that can include one or more receptacles. Depicted in FIG. 7 are examples of ways in which components of the device 700 mounted within the casing 702 can interact with components of the device 700 disposed within cassettes inserted into the receptacles. Not all components of the device 700 are shown in FIG. 7 .
- a first cassette 704 having a variety of components is shown inserted into a receptacle formed in the casing 702 .
- a second cassette 706 is also inserted into a receptacle. Components mounted within the cassettes are indicated with dashed lines in FIG. 7 , while components mounted within the casing 702 are depicted with solid lines.
- one or more actuators 708 housed within the casing 702 operate syringe pumps 710 located within the first cassette 704 .
- the pumps 710 are connected to sections of tubing 716 that move fluid among various components of the system.
- the movement of fluid is at least partially controlled by the action of one or more pinch valves 712 positioned within the casing 702 .
- the pinch valves 712 have arms 714 that extend within the first cassette 704 . Movement of the arms 714 can constrict a section of tubing 716 in order to create an effective seal.
- the second cassette 706 includes a flow cell 720 that engages a centrifuge motor 718 mounted within the casing 702 of the device 700 when the cassette is inserted into a receptacle.
- a filter wheel motor 722 disposed within the housing 702 rotates a filter wheel 724 in order to align a filter with a window of the flow cell 720 .
- An optical light path including a light source 726 within the housing 702 routes a beam of infrared light through the filter and the flow cell 720 .
- a detector 728 measures the optical density of the light transmitted through the filter and flow cell 720 .
- FIG. 8 illustrates a disposable cartridge 800 that can interface with a fluid system such as the fluid system 510 of FIG. 5 .
- the disposable cartridge 800 can be configured for insertion into a receptacle of the device 700 shown in FIG. 7 .
- the cartridge 800 includes one or more features that ease insertion of the cartridge 800 into a corresponding receptacle.
- the cartridge 800 can be shaped so as to promote insertion of the cartridge 800 in the correct orientation.
- the cartridge 800 can also include labeling or coloring affixed to or integrated with the cartridge's exterior casing that help a handler insert the cartridge 800 into a receptacle properly.
- the cartridge 800 can include one or more ports for connecting to material sources.
- one port 802 can be configured to attach to an anticoagulant source 804 .
- Other ports can be provided to connect to, for example, a saline source, an infusion pump, a sample source, and/or a source of nitrogen gas.
- the ports can be connected to sections of tubing within the cartridge 800 . In some embodiments, the sections of tubing are opaque or covered so that fluids within the tubing cannot be seen.
- the cartridge 800 shown in FIG. 8 includes one or more injector needles 806 .
- the injector needles 806 can be configured to inject at least a portion of a sample into a flow cell (not shown).
- the housing of the cartridge 800 can include a tubing space 808 for one or more sections of tubing.
- the body of the cartridge 800 includes one or more apertures 809 through which various components, such as, for example, pinch valves and sensors, can interface with the fluidics contained in the cartridge 800 .
- the sections of tubing found in the tubing space 808 can be aligned with the apertures 809 in order to implement at least some of the functionality shown in the fluid system 510 of FIG. 5 .
- the cartridge 800 can include a pouch space 810 for storing one or more components of the fluid system 510 .
- one or more pouches and/or bladders can be disposed in the pouch space 810 .
- a cleaner pouch and a waste bladder are housed in the pouch space 810 .
- the waste bladder can be placed under the cleaner pouch such that, as detergent is removed from the cleaner pouch, the waste bladder has more room to fill.
- the components placed in the pouch space 810 can also be placed side-by-side or in any other suitable configuration.
- the pouch space 810 can be isolated from the rest of the cartridge 800 by one or more walls 811 .
- One or more connectors 812 , 814 can be formed adjacent to the pouch space 810 to provide communication between components housed in the pouch space 810 and other components of the fluid system 510 .
- the cartridge 800 can include one or more pumps 816 that facilitate movement of fluid within the fluid system 510 .
- Each of the pumps 816 can be, for example, a syringe pump having a plunger.
- the plunger can include a portion 818 configured to interface with an actuator housed outside the cartridge 800 .
- the portion 818 of the pump that interfaces with an actuator can be exposed to the exterior of the cartridge 800 housing by one or more apertures in the housing.
- the disposable cartridge 800 is designed for single patient use.
- the cartridge 800 may also be designed for replacement after a period of operation. For example, in some embodiments, if the cartridge 800 is installed in a continuously operating monitoring device that performs four measurements per hour, the waste bladder may become filled or the detergent in the cleaner pouch depleted after about three days. The cartridge 800 can be replaced before the detergent and waste bladder are exhausted.
- the cartridge 800 can be configured for easy replacement.
- the cartridge 800 is designed to have an installation time of only several minutes.
- the cartridge can be designed to be installed in less than about five minutes.
- various portions of the fluidics contained in the cartridge 800 can be primed by automatically filling the fluidics with saline.
- the saline can be mixed with detergent powder from the cleaner pouch in order to create a cleaning solution.
- the cartridge 800 can also be designed to have a relatively brief shut down time.
- the shut down process can be configured to take less than about five minutes.
- the shut down process can include flushing the patient line; sealing off the insulin pump connection, the saline source connection, and the sample source connection; and taking other steps to decrease the risk that fluids within the used cartridge 800 will leak after disconnection from the monitoring device.
- the cartridge 800 is designed to fit within standard waste containers found in a hospital, such as a standard biohazard container.
- the cartridge 800 can be less than one foot long, less than one foot wide, and less than two inches thick.
- the cartridge 800 is designed to withstand a substantial impact, such as that caused by hitting the ground after a four foot drop, without damage to the housing or internal components.
- the cartridge 800 is designed to withstand significant clamping force applied to its casing.
- the cartridge 800 can be built to withstand five pounds per square inch of force without damage.
- the cartridge 800 is non pyrogenic and/or latex free.
- the system 400 comprises the optical system 412 for analysis of a fluid sample.
- the optical system 412 comprises one or more optical components including, for example, a spectrometer, a photometer, a reflectometer, or any other suitable device for measuring optical properties of the fluid sample.
- the optical system 412 may perform one or more optical measurements on the fluid sample including, for example, measurements of transmittance, absorbance, reflectance, scattering, and/or polarization.
- the optical measurements may be performed in one or more wavelength ranges including, for example, infrared (IR) and/or optical wavelengths.
- IR infrared
- the measurements from the optical system 412 are communicated to the algorithm processor 416 for analysis.
- the algorithm processor 416 computes concentration of analyte(s) (and/or interferent(s)) of interest in the fluid sample.
- Analytes of interest include, e.g., glucose and lactate in whole blood or blood plasma.
- FIG. 9 schematically illustrates an embodiment of the optical system 412 that comprises a spectroscopic analyzer 910 adapted to measure spectra of a fluid sample such as, for example, blood or blood plasma.
- the analyzer 910 comprises an energy source 912 disposed along an optical axis X of the analyzer 910 .
- the energy source 912 When activated, the energy source 912 generates an electromagnetic energy beam E, which advances from the energy source 912 along the optical axis X.
- the energy source 912 comprises an infrared energy source
- the energy beam E comprises an infrared beam.
- the infrared energy beam E comprises a mid-infrared energy beam or a near-infrared energy beam.
- the energy beam E may include optical and/or radio frequency wavelengths.
- the energy source 912 may comprise a broad-band and/or a narrow-band source of electromagnetic energy.
- the energy source 912 comprises optical elements such as, e.g., filters, collimators, lenses, mirrors, etc., that are adapted to produce a desired energy beam E.
- the energy beam E is an infrared beam in a wavelength range between about 2 ⁇ m and 20 ⁇ m.
- the energy beam E comprises an infrared beam in a wavelength range between about 4 ⁇ m and 10 ⁇ m.
- water generally is the main contributor to the total absorption together with features from absorption of other blood components, particularly in the 6 ⁇ m-10 ⁇ m range.
- the 4 ⁇ m to 10 ⁇ m wavelength band has been found to be advantageous for determining glucose concentration, because glucose has a strong absorption peak structure from about 8.5 ⁇ m to 10 ⁇ m, whereas most other blood components have a relatively low and flat absorption spectrum in the 8.5 ⁇ m to 10 ⁇ m range. Two exceptions are water and hemoglobin, which are interferents in this range.
- the energy beam E may be temporally modulated to provide increased signal-to-noise ratio (S/N) of the measurements provided by the analyzer 910 as further described below.
- the beam E is modulated at a frequency of about 10 Hz or in a range from about 1 Hz to about 30 Hz.
- a suitable energy source 912 may be an electrically modulated thin-film thermoresistive element such as the HawkEye IR-50 available from Hawkeye Technologies of Milford, Conn.
- the energy beam E propagates along the optical axis X and passes through an aperture 914 and a filter 915 thereby providing a filtered energy beam E f .
- the aperture 914 helps collimate the energy beam E and may include one or more filters adapted to reduce the filtering burden of the filter 915 .
- the aperture 914 may comprise a broadband filter that substantially attenuates beam energy outside a wavelength band between about 4 ⁇ m to about 10 ⁇ m.
- the filter 915 may comprise a narrow-band filter that substantially attenuates beam energy having wavelengths outside of a filter passband (which may be tunable or user-selectable in some embodiments).
- the filter passband may be specified by a half-power bandwidth (“HPBW”).
- the filter 915 may have an HPBW in a range from about 0.01 ⁇ m to about 1 ⁇ m. In one embodiment, the bandwidths are in a range from about 0.1 ⁇ m to 0.35 ⁇ m. Other filter bandwidths may be used.
- the filter 915 may comprise a varying-passband filter, an electronically tunable filter, a liquid crystal filter, an interference filter, and/or a gradient filter. In some embodiments, the filter 915 comprises one or a combination of a grating, a prism, a monochrometer, a Fabry-Perot etalon, and/or a polarizer. Other optical elements as known in the art may be utilized as well.
- the analyzer 910 comprises a filter wheel assembly 921 configured to dispose one or more filters 915 along the optical axis X.
- the filter wheel assembly 921 comprises a filter wheel 918 , a filter wheel motor 916 , and a position sensor 920 .
- the filter wheel 918 may be substantially circular and have one or more filters 915 or other optical elements (e.g., apertures, gratings, polarizers, etc.) disposed around the circumference of the wheel 918 .
- the number of filters 915 in the filter wheel 916 may be, for example, 1, 2, 5, 10, 15, 20, 25, or more.
- the motor 916 is configured to rotate the filter wheel 918 to dispose a desired filter 915 (or other optical element) in the energy beam E so as to produce the filtered beam E f .
- the motor 916 comprises a stepper motor.
- the position sensor 920 determines the angular position of the filter wheel 916 , and communicates a corresponding filter wheel position signal to the algorithm processor 416 , thereby indicating which filter 915 is in position on the optical axis X.
- the position sensor 920 may be a mechanical, optical, and/or magnetic encoder.
- An alternative to the filter wheel 918 is a linear filter translated by a motor.
- the linear filter may include an array of separate filters or a single filter with properties that change along a linear dimension.
- the filter wheel motor 916 rotates the filter wheel 918 to position the filters 915 in the energy beam E to sequentially vary the wavelengths or the wavelength bands used to analyze the fluid sample.
- each individual filter 915 is disposed in the energy beam E for a dwell time during which optical properties in the passband of the filter are measured for the sample.
- the filter wheel motor 916 then rotates the filter wheel 918 to position another filter 915 in the beam E.
- 25 narrow-band filters are used in the filter wheel 918 , and the dwell time is about 2 seconds for each filter 915 .
- a set of optical measurements for all the filters can be taken in about 2 minutes, including sampling time and filter wheel movement.
- the dwell time may be different for different filters 915 , for example, to provide a substantially similar S/N ratio for each filter measurement.
- the filter wheel assembly 921 functions as a varying-passband filter that allows optical properties of the sample to be analyzed at a number of wavelengths or wavelength bands in a sequential manner.
- the filter wheel 918 includes 25 finite-bandwidth infrared filters having a Gaussian transmission profile and full-width half-maximum (FWHM) bandwidth of 28 cm ⁇ 1 corresponding to a bandwidth that varies from 0.14 ⁇ m at 7.08 ⁇ m to 0.28 ⁇ m at 10 ⁇ m.
- the central wavelength of the filters are, in microns: 7.082, 7.158, 7.241, 7.331, 7.424, 7.513, 7.605, 7.704, 7.800, 7.905, 8.019, 8.150, 8.271, 8.598, 8.718, 8.834, 8.969, 9.099, 9.217, 9.346, 9.461, 9.579, 9.718, 9.862, and 9.990.
- the filtered energy beam E f propagates to a beamsplitter 922 disposed along the optical axis X.
- the beamsplitter 922 separates the filtered energy beam E f into a sample beam E s and a reference beam E r .
- the reference beam E r propagates along a minor optical axis Y, which in this embodiment is substantially orthogonal to the optical axis X.
- the energies in the sample beam E s and the reference beam E r may comprise any suitable fraction of the energy in the filtered beam E f .
- the sample beam E s comprises about 80%
- the reference beam E r comprises about 20%, of the filtered beam energy E f .
- a reference detector 936 is positioned along the minor optical axis Y.
- An optical element 934 such as a lens, may be used to focus or collimate the reference beam E r onto the reference detector 936 .
- the reference detector 936 provides a reference signal, which can be used to monitor fluctuations in the intensity of the energy beam E emitted by the source 912 . Such fluctuations may be due to drift effects, aging, wear, or other imperfections in the source 912 .
- the algorithm processor 416 may utilize the reference signal to identify changes in properties of the sample beam E s that are attributable to changes in the emission from the source 912 and not to the properties of the fluid sample.
- the analyzer 910 may advantageously reduce possible sources of error in the calculated properties of the fluid sample (e.g., concentration).
- the beamsplitter 922 is not used, and substantially all of the filtered energy beam E f propagates to the fluid sample.
- the sample beam E s propagates along the optical axis X, and a relay lens 924 transmits the sample beam E s into a sample cell 948 so that at least a fraction of the sample beam E s is transmitted through at least a portion of the fluid sample in the sample cell 948 .
- a sample detector 930 is positioned along the optical axis X to measure the sample beam E s that has passed through the portion of the fluid sample.
- An optical element 928 such as a lens, may be used to focus or collimate the sample beam E s onto the sample detector 930 .
- the sample detector 930 provides a sample signal that can be used by the algorithm processor 416 as part of the sample analysis.
- the sample cell 948 comprises the flow cell 648 located toward the circumference of the centrifuge wheel 650 .
- the flow cell 648 comprises windows that are substantially transmissive to energy in the sample beam E s .
- the windows may comprise calcium fluoride.
- the flow cell 648 is in fluid communication with an injector system that permits filling the flow cell 648 with a fluid sample (e.g., whole blood) and flushing the flow cell 648 (e.g., with saline or a detergent).
- the injector system may disconnect after filling the flow cell 648 with the fluid sample to permit free spinning of the centrifuge wheel 650 by centrifuge motor 926 .
- the fluid sample e.g., a whole blood sample
- the fluid sample is spun at about 7200 rpm for about 2 minutes to separate blood plasma for spectral analysis.
- an anti-clotting agent such as heparin may be added to the fluid sample before centrifuging to reduce clotting.
- the embodiment of the analyzer 910 illustrated in FIG. 9 advantageously permits direct measurement of the concentration of analytes in the plasma sample rather than by inference of the concentration from measurements of a whole blood sample.
- An additional advantage is that relatively small volumes of fluid may be spectroscopically analyzed.
- the fluid sample volume is between about 1 ⁇ L and 80 ⁇ L and is about 25 ⁇ L in some embodiments.
- the flow cell 648 is disposable and is intended for use with a single patient or for a single measurement.
- the reference detector 936 and the sample detector 930 comprise broadband pyroelectric detectors.
- some pyroelectric detectors are sensitive to vibrations.
- the output of a pyroelectric infrared detector is the sum of the exposure to infrared radiation and to vibrations of the detector.
- the sensitivity to vibrations also known as “microphonics,” can introduce a noise component to the measurement of the reference and sample energy beams E r , E s using some pyroelectric infrared detectors.
- the analyzer 910 may be desirable for the analyzer 910 to provide high signal-to-noise ratio measurements, such as, e.g., S/N in excess of 100 dB, some embodiments of the analyzer 910 utilize one or more vibrational noise reduction apparatus or methods.
- the analyzer 910 may be mechanically isolated so that high S/N spectroscopic measurements can be obtained for vibrations below an acceleration of about 1.5 G.
- vibrational noise can be reduced by using a temporally modulated energy source 912 combined with an output filter.
- the energy source 912 is modulated at a known source frequency, and measurements made by the detectors 936 and 930 are filtered using a narrowband filter centered at the source frequency.
- the energy output of the source 912 is sinusoidally modulated at 10 Hz, and outputs of the detectors 936 and 930 are filtered using a narrow bandpass filter of less than about 1 Hz centered at 10 Hz. Accordingly, microphonic signals that are not at 10 Hz are significantly attenuated.
- the modulation depth of the energy beam E may be greater than 50% such as, for example, 80%.
- the duty cycle of the beam may be between about 30% and 70%.
- the temporal modulation may be sinusoidal or any other waveform.
- detector output may be filtered using a synchronous demodulator and digital filter.
- the demodulator and filter are software components that may be digitally implemented in a processor such as the algorithm processor 416 . Synchronous demodulators, coupled with low pass filters, are often referred to as “lock in amplifiers.”
- the analyzer 910 may also include a vibration sensor 932 (e.g., one or more accelerometers) disposed near one (or both) of the detectors 936 and 930 .
- the output of the vibration sensor 932 is monitored, and suitable actions are taken if the measured vibration exceeds a vibration threshold. For example, in some embodiments, if the vibration sensor 932 detects above-threshold vibrations, the system discards any ongoing measurement and “holds off” on performing further measurements until the vibrations drop below the threshold. Discarded measurements may be repeated after the vibrations drop below the vibration threshold. In some embodiments, if the duration of the “hold off” is sufficiently long, the fluid in the sample cell 930 is flushed, and a new fluid sample is delivered to the cell 930 for measurement.
- a vibration sensor 932 e.g., one or more accelerometers
- the vibration threshold may be selected so that the error in analyte measurement is at an acceptable level for vibrations below the threshold.
- the threshold corresponds to an error in glucose concentration of 5 mg/dL.
- the vibration threshold may be determined individually for each filter 915 .
- the analyzer 910 include a temperature system (not shown in FIG. 9 ) for monitoring and/or regulating the temperature of system components (such as the detectors 936 , 930 ) and/or the fluid sample.
- a temperature system may include temperature sensors, thermoelectrical heat pumps (e.g., a Peltier device), and/or thermistors, as well as a control system for monitoring and/or regulating temperature.
- the control system comprises a proportional-plus-integral-plus-derivative (PID) control.
- PID proportional-plus-integral-plus-derivative
- the temperature system is used to regulate the temperature of the detectors 930 , 936 to a desired operating temperature, such as 35 degrees Celsius.
- the analyzer 910 illustrated in FIG. 9 can be used to determine optical properties of a substance in the sample cell 948 .
- the substance may include whole blood, plasma, saline, water, air or other substances.
- the optical properties include measurements of an absorbance, transmittance, and/or optical density in the wavelength passbands of some or all of the filters 915 disposed in the filter wheel 918 .
- a measurement cycle comprises disposing one or more filters 915 in the energy beam E for a dwell time and measuring a reference signal with the reference detector 936 and a sample signal with the sample detector 930 .
- the number of filters 915 used in the measurement cycle will be denoted by N, and each filter 915 passes energy in a passband around a center wavelength ⁇ i , where i is an index ranging over the number of filters (e.g., from 1 to N).
- the set of optical measurements from the sample detector 936 in the passbands of the N filters 915 provide a wavelength-dependent spectrum of the substance in the sample cell 948 .
- the spectrum will be denoted by C s ( ⁇ i ), where C s may be a transmittance, absorbance, optical density, or some other measure of an optical property of the substance.
- the spectrum is normalized with respect to one or more of the reference signals measured by the reference detector 930 and/or with respect to spectra of a reference substance (e.g., air or saline).
- the measured spectra are communicated to the algorithm processor 416 for calculation of the concentration of the analyte(s) of interest in the fluid sample.
- the analyzer 910 performs spectroscopic measurements on the fluid sample (known as a “wet” reading) and on one or more reference samples. For example, an “air” reading occurs when the sample detector 936 measures the sample signal without the sample cell 948 in place along the optical axis X. A “water” or “saline” reading occurs when the sample cell 948 is filled with water or saline, respectively.
- the algorithm processor 416 may be programmed to calculate analyte concentration using a combination of these spectral measurements.
- a pathlength corrected spectrum is calculated using wet, air, and reference readings.
- the transmittance at wavelength ⁇ i denoted by T i
- T i (S i (wet)/R i (wet))/(S i (air)/R i (air)), where S i denotes the sample signal from the sample detector 936 and R i denotes the corresponding reference signal from the reference detector 930 .
- the analyzer 910 takes a set of wet readings in each of the N filter passbands and then takes a set of air readings in each of the N filter passbands. In other embodiments, the analyzer 910 may take an air reading before (or after) the corresponding wet reading.
- the absorption coefficient ⁇ i of a substance may be written as the product of an absorptivity per mole times a molar concentration of the substance.
- FIG. 9 schematically illustrates the pathlength L of the sample cell 948 .
- the pathlength L may be determined from spectral measurements made when the sample cell 948 is filled with a reference substance.
- one or more water (or saline) readings can be used to determine the pathlength L from measurements of the transmittance (or optical density) through the cell 948 .
- several readings are taken in different wavelength passbands, and a curve-fitting procedure is used to estimate a best-fit pathlength L.
- the pathlength L may be estimated using other methods including, for example, measuring interference fringes of light passing through an empty sample cell 948 .
- the pathlength L may be used to determine the absorption coefficients of the fluid sample at each wavelength.
- Molar concentration of an analyte of interest can be determined from the absorption coefficient and the known molar absorptivity of the analyte.
- a sample measurement cycle comprises a saline reading (at one or more wavelengths), a set of N wet readings, followed by a set of N air readings. As discussed above, the sample measurement cycle can be performed in about 2 minutes when the filter dwell times are about 2 seconds.
- a detergent cleaner may be flushed through the flow cell 648 to reduce buildup of organic matter (e.g., proteins) on the windows of the flow cell 648 . The detergent is then flushed to a waste bladder.
- the system stores information related to the spectral measurements so that the information is readily available for recall by a user.
- the stored information may include wavelength-dependent spectral measurements (including fluid sample, air, and/or saline readings), computed analyte values, system temperatures and electrical properties (e.g., voltages and currents), and any other data related to use of the system (e.g., system alerts, vibration readings, S/N ratios, etc.).
- the stored information may be retained in the system for a time period such as, for example, 30 days. After this time period, the stored information may be communicated to an archival data storage system and then deleted from the system. In certain embodiments, the stored information is communicated to the archival data storage system via wired or wireless methods, e.g., over a hospital information system (HIS).
- HIS hospital information system
- the algorithm processor 416 ( FIG. 4 ) (or any other suitable processor) may be configured to receive from the analyzer 910 the wavelength-dependent optical measurements of the fluid sample.
- the optical densities OD i in each of the N filter passbands centered around wavelengths ⁇ i are communicated to the processor 416 , which analyzes the optical densities to measure and quantify one or more analytes in the presence of interferents.
- Interferents can comprise components of a material sample being analyzed for an analyte, where the presence of the interferent affects the quantification of the analyte.
- an interferent in the spectroscopic analysis of a sample to determine an analyte concentration, could be a compound having spectroscopic features that overlap with those of the analyte.
- the presence of such an interferent can introduce errors in the quantification of the analyte. More specifically, the presence of interferents can affect the sensitivity of a measurement technique to the concentration of analytes of interest in a material sample, especially when the system is calibrated in the absence of, or with an unknown amount of, the interferent.
- interferents can be classified as being endogenous (i.e., originating within the body) or exogenous (i.e., introduced from or produced outside the body).
- endogenous interferents include those blood components having origins within the body that affect the quantification of glucose, and may include water, hemoglobin, blood cells, and any other component that naturally occurs in blood.
- Exogenous interferents include those blood components having origins outside of the body that affect the quantification of glucose, and can include items administered to a person, such as medicaments, drugs, foods or herbs, whether administered orally, intravenously, topically, etc.
- interferents can comprise components which are possibly, but not necessarily, present in the sample type under analysis.
- a medicament such as acetaminophen is possibly, but not necessarily, present in this sample type.
- water is necessarily present in such blood or plasma samples.
- FIG. 10 is a flowchart that schematically illustrates an embodiment of a spectroscopic method 1010 for determining the concentration of an analyte of interest in a fluid sample in the presence of one or more possible interferents.
- spectral measurements of the fluid sample are obtained.
- the analyzer 910 may be used to obtain optical measurements Cs( ⁇ i ) of the fluid sample in a number N of filter passbands centered around wavelengths ⁇ i .
- quality of the spectral measurements is determined regardless of the concentration of the analyte of interest of the presence of possible interferents.
- one or more of poor quality spectral measurements Cs( ⁇ i ) may be rejected (e.g., as having a S/N ratio that is too low), and the method 1010 performed on the remaining, sufficiently high-quality measurements.
- additional spectral measurements of the fluid sample are obtained to replace one or more of the poor quality measurements.
- the spectral measurements are tested to determine the possible presence of interferents.
- the system may utilize spectroscopic signatures of possible interferents to test for their presence.
- the method 1010 proceeds to block 1022 in which analyte concentration is determined.
- analyte concentration is determined using a hybrid linear algorithm (HLA) in which analyte concentration is estimated from measured spectra using one or more calibration coefficients and an offset.
- HLA hybrid linear algorithm
- the test determines that one or more interferents are present at concentrations above threshold, then, in block 1018 , the above-threshold interferents are identified.
- the method 1010 proceeds to block 1020 in which the analyte concentration algorithm is adapted to account for the presence of one or more of the identified interferents. For example, in embodiments using HLA, the calibration coefficients may be adjusted to compensate for the presence of some or all of the identified interferents.
- the method 1010 proceeds to block 1022 in which analyte concentration is determined as further described below.
- Certain disclosed analysis methods are particularly effective if each analyte and interferent has a characteristic signature in the measurement (e.g., a characteristic spectroscopic feature), and if the measurement is approximately affine (e.g., includes a linear term and an offset) with respect to the concentration of each analyte and interferent.
- a calibration process is used to determine a set of one or more calibration coefficients and one or more optional offset values that permits the quantitative estimation of an analyte.
- the calibration coefficients and the offsets may be used to calculate an analyte concentration from spectroscopic measurements of a material sample (e.g., the concentration of glucose in blood plasma).
- the concentration of the analyte is estimated by multiplying the calibration coefficient by a measurement value (e.g., an optical density) to estimate the concentration of the analyte.
- a measurement value e.g., an optical density
- Both the calibration coefficient and measurement can comprise arrays of numbers.
- the measurement comprises the spectra C s ( ⁇ i ) measured at the wavelengths ⁇ i
- the calibration coefficient and optional offset comprise an array of values corresponding to each wavelength ⁇ i .
- a hybrid linear algorithm is used to estimate analyte concentration in the presence of a set of interferents, while retaining a high degree of sensitivity to the desired analyte.
- the data used to accommodate the random set of interferents may include (a) signatures of each of the members of the family of potential additional substances and (b) the typical quantitative level at which each additional substance, if present, is likely to appear.
- the calibration constant (and optional offset) are adjusted to minimize or reduce the sensitivity of the calibration to the presence of interferents that are identified as possibly being present in the fluid sample.
- the analyte analysis method uses a set of training spectra each having known analyte concentration(s) and produces a calibration that minimizes the variation in estimated analyte concentration with interferent concentration.
- the resulting calibration coefficient measures sensitivity of the measurement to analyte concentration(s) and, on average, is not sensitive to interferent concentrations.
- the training spectra need not include a spectrum from the individual whose analyte concentration is to be determined. That is, the term “training” when used in reference to the disclosed methods does not require training using measurements from the individual whose analyte concentration will be estimated (e.g., by analyzing a bodily fluid sample drawn from the individual).
- sample Population is a broad term and includes, without limitation, a large number of samples having measurements that are used in the computation of a calibration—in other words, used to train the method of generating a calibration.
- the Sample Population measurements can each include a spectrum (analysis measurement) and a glucose concentration (analyte measurement).
- the Sample Population measurements are stored in a database, referred to herein as a “Population Database.”
- the Sample Population may or may not be derived from measurements of material samples that contain interferents to the measurement of the analyte(s) of interest.
- One distinction made herein between different interferents is based on whether the interferent is present in both the Sample Population and the sample being measured, or only in the sample.
- the term “Type-A interferent” refers to an interferent that is present in both the Sample Population and in the material sample being measured to determine an analyte concentration. In certain methods it is assumed that the Sample Population includes only interferents that are endogenous, and does not include any exogenous interferents, and thus Type-A interferents are endogenous.
- Type-A interferents depends on the measurement and analyte(s) of interest, and may number, in general, from zero to a very large number (e.g., greater than 300).
- the material sample being measured for example a fluid sample in the sample cell 948 , may also include interferents that are not present in the Sample Population.
- Type-B interferent refers to an interferent that is either: 1) not found in the Sample Population but that is found in the material sample being measured (e.g., an exogenous interferent), or 2) is found naturally in the Sample Population, but is at abnormally high concentrations in the material sample (e.g., an endogenous interferent).
- Examples of a Type-B exogenous interferent may include medications, and examples of Type-B endogenous interferents may include urea in persons suffering from renal failure. For example, in mid-infrared spectroscopic absorption measurements of glucose in blood (or blood plasma), water is present in all fluid samples, and is thus a Type-A interferent.
- the selected drug is a Type-B interferent.
- the ingestion or injection of some medicines or illicit drugs can result in very high and rapidly changing concentrations of exogenous interferents.
- a list of one or more possible Type-B Interferents is referred to herein as forming a “Library of Interferents,” and each interferent in the library is referred to as a “Library Interferent.”
- the Library Interferents include exogenous interferents and endogenous interferents that may be present in a material sample due, for example, to a medical condition causing abnormally high concentrations of the endogenous interferent.
- FIG. 11 is a flowchart that schematically illustrates an embodiment of a method 1100 for estimating the concentration of an analyte in the presence of interferents.
- a measurement of a sample is obtained, and in block 1120 data relating to the obtained measurement is analyzed to identify possible interferents to the analyte.
- a model is generated for predicting the analyte concentration in the presence of the identified possible interferents, and in block 1140 the model is used to estimate the analyte concentration in the sample from the measurement.
- the model generated in block 1130 is selected to reduce or minimize the effect of identified interferents that are not present in a general population of which the sample is a member.
- the sample measurement in block 1110 is an absorption spectrum, Cs( ⁇ i ), of a measurement sample S that has, in general, one analyte of interest, glucose, and one or more interferents.
- the absorption spectrum may comprise the set of optical densities OD i measured by the analyzer 910 .
- the sample S includes Type-A interferents, at concentrations preferably within the range of those found in the Sample Population.
- a statistical comparison of the absorption spectrum of the sample S with a spectrum of the Sample Population and combinations of individual Library Interferent spectra is performed.
- the statistical comparison provides a list of Library Interferents that are possibly contained in sample S and may include either no Library Interferents or one or more Library Interferents.
- a set of spectra are generated using the spectra of the Sample Population and their respective known analyte concentrations and known spectra of the Library Interferents identified in block 1120 .
- the generated spectra are used to calculate a calibration coefficient ⁇ ( ⁇ i ) that can be used with the sample measurements Cs( ⁇ i ) to provide an estimate of the analyte concentration, g est .
- block 1120 includes forming a statistical Sample Population model (block 1210 ), assembling a library of interferent data (block 1220 ), comparing the obtained measurement and statistical Sample Population model with data for each interferent from an interferent library (block 1230 ), performing a statistical test for the presence of each interferent from the interferent library (block 1240 ), and identifying possible interferents that pass the statistical test (block 1250 ).
- the acts of block 1220 can be performed once or can be updated as necessary.
- the acts of blocks 1230 , 1240 , and 1250 can either be performed sequentially for all Library Interferents or can be repeated sequentially for each interferent.
- a Sample Population Database is formed that includes a statistically large Sample Population of individual spectra taken over the same wavelength range as the sample spectrum, C s ( ⁇ i ).
- the Sample Population does not have any of the Library Interferents present, and the material sample has interferents contained in the Sample Population and one or more of the Library Interferents. Stated in terms of Type-A and Type-B interferents, the Sample Population has Type-A interferents, and the material sample has Type-A and may have Type-B interferents.
- the statistical sample model comprises a mean spectrum and a covariance matrix calculated for the Sample Population. For example, if each spectrum measured at N wavelengths ⁇ i is represented by an N ⁇ 1 array, C, then the mean spectrum, ⁇ , is an N ⁇ 1 array having values at each wavelength averaged over the range of spectra in the Sample Population.
- additional statistical parameters may be included in the statistical model of the Sample Population spectra.
- a Library of Interferents may be assembled in block 1220 .
- a number of possible interferents can be identified, for example, as a list of possible medications or foods that might be ingested by the population of patients at issue. Spectra of these interferents can be obtained, and a range of expected interferent concentrations in the blood, or other expected sample material, can be estimated.
- Information in the Library may be assembled once and accessed as needed. For example, the Library and the statistical model of the Sample Population may be stored in a storage device associated with the algorithm processor 416 ( FIG. 4 ).
- the obtained measurement data e.g., the sample spectrum
- the statistical Sample Population model e.g., the mean spectrum and the covariance matrix
- the statistical test for the presence of an interferent comprises the following actions.
- the measured spectrum of the fluid sample, C s is modified for each interferent of the library by analytically subtracting, wavelength-by-wavelength, the spectrum of the interferent.
- the interferent concentration is assumed to be in a range from a minimum value, Tmin, to a maximum value, Tmax.
- Tmin may be zero or, alternatively, be a value between zero and Tmax, such as some fraction of Tmax.
- the statistical test for determining the presence of possible interferents in block 1240 further comprises determining a Mahalanobis distance (MD) between the modified spectrum C′ s (T) and the statistical model ( ⁇ , V) of the Sample Population.
- the value of MD 2 found from Eq. (1) is referred to herein as the “squared Mahalanobis distance” or the “MD 2 score.”
- the MD 2 score is used in various embodiments of the statistical test for determining the presence of an interferent.
- a list of possible interferents may be identified as the particular Library Interferents that pass one or more statistical tests for being present in the sample.
- One or more tests may be used, alone or in combination, to identify the possible interferents. For example, if a statistical test indicates that the interferent is present in negative concentrations, then this non-physical result is used to exclude the possible interferent from the list of possible interferents. In some embodiments, only the single most probable interferent is included on the list.
- the concentration T is varied from Tmin to Tmax (e.g., evaluate C′ s (T) over a range of values of T in Eq. (1)). If the minimum value of MD (or MD 2 ) in this interval is below a minimum threshold, then the test indicates the probable presence of the interferent in the sample.
- the minimum threshold MD 2 is chosen relative to quantiles of a ⁇ 2 random variable having N degrees of freedom, where N is the number of wavelengths in the spectrum C s . In some embodiments, the 95% quantile is used as the minimum threshold.
- the maximum threshold generally is empirically determined. In one embodiment, it is found that a maximum threshold value is in a range from about 50 to about 200.
- interferent concentration is assumed to have a log-normal distribution with a value of 95% at the assumed maximum interferent concentration in the fluid and a standard deviation of one half the mean.
- the RPS is then compared to a minimum value (P min ).
- Possible interferents are identified as interferents having an RPS greater than P min .
- Possible interferents are denoted herein with the variable ⁇ .
- the value of P min is empirically determined from an analysis of the measurements. For example, a value of 0.70 may result in a single possible interferent (a “single interferent identification”) and a value of 0.3 may result in three possible interferents (a multiple interferent identification).
- one or more of the above statistical tests are used to determine a list of possible interferents ⁇ that may be present in the fluid sample.
- the list of possible interferents includes only the single most probable interferent.
- the list of possible interferents ⁇ may include each of the interferents in the Library of Interferents.
- FIG. 13 is a flowchart that schematically illustrates an example embodiment of the acts of block 1130 .
- synthesized Sample Population measurements are generated to form an Interferent Enhanced Spectral Database (IESD).
- IESD Interferent Enhanced Spectral Database
- the spectra in the IESD are partitioned into a calibration set and a test.
- the calibration set is used to generate a calibration coefficient
- the calibration coefficient is used to estimate the analyte concentration of the test set.
- errors in the estimated analyte concentration of the test set are calculated, and in block 1360 an average calibration coefficient is calculated based on errors in the test set(s).
- the average calibration coefficient is applied to the measured spectra to determine an estimated single-interferent analyte concentration.
- the blocks 1310 - 1360 are performed for each possible interferent ⁇ to provide a corresponding “single-interferent” average calibration coefficient for each particular interferent. In other embodiments, the blocks 1310 - 1360 are performed only for the single most probable interferent in the list identified in block 1250 .
- synthesized Sample Population spectra are generated by adding a random concentration of one of the possible interferents ⁇ to each Sample Population spectrum. These spectra are referred to herein as an Interferent-Enhanced Spectral Database or IESD.
- the IESD is formed as follows. A plurality of Randomly-Scaled Single Interferent Spectra (RSIS) are formed by combinations of the interferent ⁇ having spectrum IF ⁇ multiplied by the maximum concentration Tmax m , which is scaled by a random factor between zero and one.
- RSIS Randomly-Scaled Single Interferent Spectra
- the scaling places the maximum concentration at the 95 th percentile of a log-normal distribution in order to generate a wide range of concentrations.
- the log-normal distribution has a standard deviation equal to half of its mean value.
- CIS Combination Interferent Spectra
- the CIS and replicates of the Sample Population spectra are combined to form the IESD. Since the interferent spectra and the Sample Population spectra may have been obtained from measurements having different optical pathlengths, the CIS may be scaled to the same pathlength as the Sample Population spectra.
- the Sample Population Database is then replicated R times, where R depends on factors including the size of the Database and the number of interferents.
- the IESD includes R copies of each of the Sample Population spectra, where one copy is the original Sample Population Data, and the remaining R ⁇ 1 copies each have one randomly chosen CIS spectra added. Accordingly, each of the IESD spectra has an associated analyte concentration from the Sample Population spectra used to form the particular IESD spectrum.
- a 10-fold replication of the Sample Population Database is used for 130 Sample Population spectra obtained from 58 different individuals and 18 Library Interferents.
- a smaller replication factor may be used if there is greater spectral variety among the Library Interferent spectra, and a larger replication factor may be used if there is a greater number of Library Interferents.
- the blocks 1320 - 1350 may be executed to repeatedly combine different spectra of the IESD to statistically average out effects of the interferent ⁇ .
- the IESD may be partitioned into two subsets: a calibration set and a test set. Repeated partitioning of the IESD into different calibration and test sets improves the statistical significance of the calibration constant.
- the calibration set includes a random selection of some of the IESD spectra, and the test set includes the remaining unselected IESD spectra. In one embodiment, the calibration set includes approximately two-thirds of the IESD spectra.
- the calibration set is used to generate a calibration coefficient for estimating the analyte concentration from a sample measurement.
- a glucose absorption spectrum is obtained and indicated as ⁇ G .
- the calibration coefficient, ⁇ is calculated from C′ and ⁇ G , according to the following 5 actions:
- the calibration coefficient is used to estimate the analyte concentration for the spectra in the test set. For example, each spectrum of the test set is multiplied by the calibration coefficient ⁇ to calculate an estimated glucose concentration. Since each spectrum in the test set has a known glucose concentration, the error between the calculated and known glucose concentration may be calculated, in block 1350 .
- Blocks 1320 - 1350 may be repeated for a number of different random combinations of calibration sets.
- the number of combinations may be in a range from hundreds to thousands.
- an average calibration constant is calculated from the calibration coefficient and the error from the many calibration and test sets.
- the average calibration coefficient may be computed as a weighted average of the individual calibration coefficients from the combinations.
- one embodiment of a method of computing a calibration constant based on an identified interferent ⁇ can be summarized as follows:
- the average calibration coefficient determined in block 1360 corresponds to a single interferent ⁇ from the list of possible interferents and is denoted herein as a single-interferent calibration coefficient ⁇ avg ( ⁇ ).
- the blocks 1310 - 1370 can be repeated for each of the interferents on the list of possible interferents, thereby providing an array of estimated, single-interferent analyte concentrations.
- the blocks 1310 - 1360 are performed only once for the single most probable interferent on the list (e.g., the array of single-interferent analyte concentrations has a single member).
- the array of single-interferent concentrations, g( ⁇ ) are combined to determine an estimated analyte concentration, g est , for the fluid sample.
- a weighting function p( ⁇ ) is determined for each of the interferents on the list of possible interferents.
- the Raw Probability Score (RPS) (described above following Eq. (2)) is used in determining the weighting function.
- the RPS's determined for the interferents on the list of possible interferents are rescaled to unit probability.
- the estimated analyte concentration, g est is determined by combining the single-interferent estimates, g( ⁇ ), and the weighting functions, p( ⁇ ), to generate a likelihood-weighted average analyte concentration.
- the algorithm processor 416 may be configured, additionally or alternatively, to implement other methods for determining analyte concentration. For example, in certain embodiments, a parameter-free interferent rejection algorithm is implemented. In certain such embodiments, a sample measurement is obtained, substantially as described above in reference to block 1110 of FIG. 11 . The algorithm processor 416 , in block 1120 , analyzes the obtained measurement to identify possible interferents. For example, the algorithm processor 416 may form a statistical sample population model and calculate statistical sample population parameters including mean spectra and covariance matrix (e.g., as described above with reference to block 1210 of FIG. 12 ). The processor 416 may then assemble a library of interferent data (e.g., as described above with reference to block 1220 of FIG. 12 ). The library may include interferent spectra, maximum plasma concentration, and a common random concentration distribution function for each interferent. In some embodiments, the processor 416 calculates a common variance (denoted by v) of the common random concentration distribution function.
- the library may be divided into groups comprising some or all combinations of a number K of the library interferents.
- the number K may be an integer such as 1, 2, 3, 4, 5, 6, 7, 8, 15, 20, or more.
- a statistical test may then be performed to determine how well some or all of the groups of K library interferents fits the statistical population model. For example, the statistical test may provide a value for the Mahalanobis distance (of distance squared) for each group and/or an estimate of the concentration of some or all of the library interferents. In some embodiments, groups in which one or more estimated concentrations are negative are eliminated as being unphysical.
- some or all groups having negative estimated concentrations may be retained, because they may indicate that the estimated concentration is lower than a standard or reference concentration (e.g., due to dilution of the sample measurement by saline or another fluid).
- a subset of the remaining groups may be selected, which provide the most likely interferents.
- the subset may comprise the groups having a number N of the smallest values of the Mahalanobis distance (or distance squared). In various embodiments, the number N may be 1, 2, 5, 10, 20, 100, 200, or more.
- the subset is used to form a model group comprising some or all combinations of a number L of the groups in the subset.
- the algorithm processor 416 may then, in block 1130 of FIG. 11 , generate a model for predicting the analyte concentration from the obtained sample measurement.
- a model for predicting the analyte concentration For example, in some implementations, for some or all of the model groups, an average group interferent calibration coefficient is calculated, which accounts for the presence of the distinct interferents in any particular model group.
- the group interferent calibration coefficient may be calculated according to blocks 1310 - 1360 of FIG. 13 in some embodiments. In these embodiments, the group interferents are used, rather than a single interferent, in block 1310 to generate synthesized sample population spectra by adding random concentrations of each interferent present in the group to form an Interferent Enhanced Spectral Database (IESD).
- IESD Interferent Enhanced Spectral Database
- the IESD is partitioned into a calibration set and a test set.
- the calibration set is used to generate a calibration coefficient for estimating the analyte concentration in the presence of the interferents in the group.
- the calibration coefficient is used to estimate the analyte concentration of the test set, assuming the presence of that interferent group's interferents.
- the error is calculated in the estimated analyte concentration for the test set. Blocks 1320 - 1350 may be repeated one or more times to obtain group interferent calibration coefficients and errors for different combinations of calibration and test sets.
- an average group interferent calibration coefficient for each group is calculated from the results determined from blocks 1320 - 1350 .
- the algorithm processor 416 may then use an average calibration coefficient to estimate analyte concentration from the obtained sample measurement.
- the average calibration coefficient is determined from an average of the group interferent calibration coefficients determined in block 1360 .
- the average may be a straight average or a weighted average in various embodiments.
- the analyte concentration is determined by multiplying this average calibration coefficient by the measured spectra.
- the algorithm processor 416 uses different algorithms in block 1130 of FIG. 11 to determine an average calibration coefficient. For example, in some embodiments, every IESD is used as a calibration set, and there is no partition of the IESD into a calibration set and a test set and no error estimate is calculated. Accordingly, in some of these embodiments, the algorithm processor 416 may not perform blocks 1320 , 1340 , and 1350 . The average calibration coefficient is determined, in block 1360 (or block 1370 ) from all the groups in the IESD.
- the average group calibration coefficient may be determined from the following actions.
- a possible advantage of the foregoing algorithms is more rapid execution time by the algorithm processor 416 , because the calibration coefficient is computed directly, without synthesizing spectra or breaking the data into calibration and test sets.
- regression, partial least squares, and/or principal component resolution techniques may be used to determine the average group calibration coefficient.
- the system 400 may include a display controller 414 , for example, as depicted in FIG. 4 .
- the display controller 414 may comprise an input device including, for example, a keypad or a keyboard, a mouse, a touchscreen display, and/or any other suitable device for inputting commands and/or information.
- the display controller 414 may also include an output device including, for example, an LCD monitor, a CRT monitor, a touchscreen display, a printer, and/or any other suitable device for outputting text, graphics, images, videos, etc.
- a touchscreen display is advantageously used for both input and output.
- the display controller 414 may include a user interface 1400 by which users can conveniently and efficiently interact with the system 400 .
- the user interface 1400 may be displayed on the output device of the system 400 (e.g., the touchscreen display).
- FIGS. 14A and 14B schematically illustrate the visual appearance of embodiments of the user interface 1400 .
- the user interface 1400 may show patient identification information 1402 , which may include patient name and/or a patient ID number.
- the user interface 1400 also may include the current date and time 1404 .
- An operating graphic 1406 shows the operating status of the system 400 .
- the operating status is “Running,” which indicates that the system 400 is fluidly connected to the patient (“Jill Doe”) and performing normal system functions such as infusing fluid and/or drawing blood.
- the user interface 1400 can include one or more analyte concentration graphics 1408 , 1412 , which may show the name of the analyte and its last measured concentration.
- the graphic 1408 in FIG. 14A shows “Glucose” concentration of 150 mg/dl, while the graphic 1412 shows “Lactate” concentration of 0.5 mmol/L.
- the particular analytes displayed and their measurement units may be selected by the user.
- the size of the graphics 1408 , 1412 may be selected to be easily readable out to a distance such as, e.g., 30 feet.
- the user interface 1400 may also include a next-reading graphic 1410 that indicates the time until the next analyte measurement is to be taken. In FIG. 14A , the time until next reading is 3 minutes, whereas in FIG. 14B , the time is 6 minutes, 13 seconds.
- the user interface 1400 may include an analyte concentration status graphic 1414 that indicates status of the patient's current analyte concentration compared with a reference standard.
- the analyte may be glucose
- the reference standard may be a hospital ICU's tight glycemic control (TGC).
- TGC tight glycemic control
- the status graphic 1414 displays “High Glucose,” because the glucose concentration (150 mg/dl) exceeds the maximum value of the reference standard.
- the status graphic 1414 displays “Low Glucose,” because the current glucose concentration (79 mg/dl) is below the minimum reference standard.
- the status graphic 1414 may indicate normal (e.g., “Normal Glucose”), or it may not be displayed at all.
- the status graphic 1414 may have a background color (e.g., red) when the analyte concentration exceeds the acceptable bounds of the reference standard.
- the user interface 1400 may include one or more trend indicators 1416 that provide a graphic indicating the time history of the concentration of an analyte of interest.
- the trend indicator 1416 comprises a graph of the glucose concentration (in mg/dl) versus elapsed time (in hours) since the measurements started.
- the graph includes a trend line 1418 indicating the time-dependent glucose concentration.
- the trend line 1418 may include measurement error bars and may be displayed as a series of individual data points.
- the glucose trend indicator 1416 is shown as well as a trend indicator 1430 and trend line 1432 for the lactate concentration.
- a user may select whether none, one, or both trend indicators 1416 , 1418 are displayed.
- one or both of the trend indicators 1416 , 1418 may appear only when the corresponding analyte is in a range of interest such as, for example, above or below the bounds of a reference standard.
- the user interface 1400 may include one or more buttons 1420 - 1426 that can be actuated by a user to provide additional functionality or to bring up suitable context-sensitive menus and/or screens. For example, in the embodiments shown in FIGS. 14A and 14B , four buttons 1420 - 1426 are shown, although fewer or more buttons are used in other embodiments.
- the button 1420 (“End Monitoring”) may be pressed when one or both of the disposable cassettes 610 , 612 (see FIG. 6 ) are to be removed. In many embodiments, because the cassettes 610 , 612 are not reusable, a confirmation window appears when the button 1420 is pressed. If the user is certain that monitoring should stop, the user can confirm this by actuating an affirmative button in the confirmation window.
- the system 400 performs appropriate actions to cease fluid infusion and blood draw and to permit ejection of one (or both) cassettes 610 , 612 .
- the button 1422 (“Pause”) may be actuated by the user if patient monitoring is to be interrupted but is not intended to end.
- the “Pause” button 1422 may be actuated if the patient is to be temporarily disconnected from the system 400 (e.g., by disconnecting the tubes 306 ). After the patient is reconnected, the button 1422 may be pressed again to resume monitoring. In some embodiments, after the “Pause” button 1422 has been pressed, the button 1422 displays “Resume.”
- the button 1424 (“Delay 5 Minutes”) causes the system 400 to delay the next measurement by a delay time period (e.g., 5 minutes in the depicted embodiments). Actuating the delay button 1424 may be advantageous if taking a reading would be temporarily inconvenient, for example, because a health care professional is attending to other needs of the patient.
- the delay button 1424 may be pressed repeatedly to provide longer delays. In some embodiments, pressing the delay button 1424 is ineffective if the accumulated delay exceeds a maximum threshold.
- the next-reading graphic 1410 automatically increases the displayed time until the next reading for every actuation of the delay button 1424 (up to the maximum delay).
- the button 1426 (“Dose History”) may be actuated to bring up a dosing history window that displays patient dosing history for an analyte or medicament of interest.
- the dosing history window displays insulin dosing history of the patient and/or appropriate hospital dosing protocols.
- a nurse attending the patient can actuate the dosing history button 1426 to determine the time when the patient last received an insulin dose, the last dosage amount, and/or the time and amount of the next dosage.
- the system 400 may receive the patient dosing history via wired or wireless communications from a hospital information system.
- the user interface 1400 may include additional and/or different buttons, menus, screens, graphics, etc. that are used to implement additional and/or different functionalities.
- FIG. 15 schematically depicts various components and/or aspects of a patient monitoring system 15130 and how those components and/or aspects relate to each other. Some of the depicted components can be included in a kit containing a plurality of components. Some of the depicted components, including, for example, the components represented within the dashed rounded rectangle 15140 of FIG. 15 , are optional and/or can be sold separately from other components.
- the patient monitoring system 15130 shown in FIG. 15 includes a monitoring device 15132 .
- the monitoring device 15132 can provide monitoring of physiological parameters of a patient.
- the monitoring device 15132 measures glucose and/or lactate concentrations in the patient's blood.
- the measurement of such physiological parameters is substantially continuous.
- the monitoring device 15132 may also measure other physiological parameters of the patient.
- the monitoring device 15132 is used in an intensive care unit (ICU) environment. In some embodiments, one monitoring device 15132 is allocated to each patient room in an ICU.
- ICU intensive care unit
- the patient monitoring system 15130 can include an optional interface cable 15142 .
- the interface cable 15142 connects the monitoring device 15132 to a patient monitor (not shown).
- the interface cable 15142 can be used to transfer data from the monitoring device 15132 to the patient monitor for display.
- the patient monitor is a bedside cardiac monitor having a display that is located in the patient room.
- the interface cable 15142 transfers data from the monitoring device 15132 to a central station monitor and/or to a hospital information system (HIS).
- HIS hospital information system
- the ability to transfer data to a central station monitor and/or to a HIS may depend on the capabilities of the patient monitor system.
- an optional bar code scanner 15144 is connected to the monitoring device 15132 .
- the bar code scanner 15144 is used to enter patient identification codes, nurse identification codes, and/or other identifiers into the monitoring device 15132 .
- the bar code scanner 15144 contains no moving parts.
- the bar code scanner 15144 can be operated by manually sweeping the scanner 15144 across a printed bar code or by any other suitable means.
- the bar code scanner 15144 includes an elongated housing in the shape of a wand.
- the patient monitoring system 15130 includes a fluidic system kit 15134 connected to the monitoring device 15132 .
- the fluidic system kit 15134 includes fluidic tubes that connect a fluid source to an analytic subsystem.
- the fluidic tubes can facilitate fluid communication between a blood source or a saline source and an assembly including a flow cell and/or a centrifuge.
- the fluidic system kit 15134 includes many of the components that enable operation of the monitoring device 15132 .
- the fluidic system kit 15134 can be used with anti-clotting agents (such as heparin), saline, a saline infusion set, a patient catheter, a port sharing IV infusion pump, and/or an infusion set for an IV infusion pump, any or all of which may be made by a variety of manufacturers.
- the fluidic system kit 15134 includes a monolithic housing that is sterile and disposable.
- at least a portion of the fluidic system kit 15134 is designed for single patient use.
- the fluidic system kit 15134 can be constructed such that it can be economically discarded and replaced with a new fluidic system kit 15134 for every new patient to which the patient monitoring system 15130 is connected.
- the fluidic system kit 15134 can be designed to be discarded after a certain period of use, such as a day, several days, several hours, three days, a combination of hours and days such as, for example, three days and two hours, or some other period of time. Limiting the period of use of the fluidic system kit 15134 may decrease the risk of malfunction, infection, or other conditions that can result from use of a medical apparatus for an extended period of time.
- the fluidic system kit 15134 includes a connector with a luer fitting for connection to a saline source.
- the connector may be, for example, a three-inch pigtail connector.
- the fluidic system kit 15134 can be used with a variety of spikes and/or IV sets used to connect to a saline bag.
- the fluidic system kit 15134 also includes a three-inch pigtail connector with a luer fitting for connection to one or more IV pumps.
- the fluidic system kit 15134 can be used with one or more IV sets made by a variety of manufacturers, including IV sets obtained by a user of the fluidic system kit 15134 for use with an infusion pump.
- the fluidic system kit 15134 includes a tube with a low dead volume luer connector for attachment to a patient vascular access point.
- the tube can be approximately seven feet in length and can be configured to connect to a proximal port of a cardiovascular catheter.
- the fluidic system kit 15134 can be used with a variety of cardiovascular catheters, which can be supplied, for example, by a user of the fluidic system kit 15134 .
- the monitoring device 15132 is connected to a support apparatus 15136 , such as an IV pole.
- the support apparatus 15136 can be customized for use with the monitoring device 15132 .
- a vendor of the monitoring device 15132 may choose to bundle the monitoring device 15132 with a custom support apparatus 15136 .
- the support apparatus 15136 includes a mounting platform for the monitoring device 15132 .
- the mounting platform can include mounts that are adapted to engage threaded inserts in the monitoring device 15132 .
- the support apparatus 15136 can also include one or more cylindrical sections having a diameter of a standard IV pole, for example, so that other medical devices, such as IV pumps, can be mounted to the support apparatus.
- the support apparatus 15136 can also include a clamp adapted to secure the apparatus to a hospital bed, an ICU bed, or another variety of patient conveyance device.
- the monitoring device 15132 is electrically connected to an optional computer system 15146 .
- the computer system 15146 can be used to communicate with one or more monitoring devices. In an ICU environment, the computer system 15146 can be connected to at least some of the monitoring devices in the ICU.
- the computer system 15146 can be used to control configurations and settings for multiple monitoring devices (for example, the system can be used to keep configurations and settings of a group of monitoring devices common).
- the computer system 15146 can also run optional software, such as data analysis software 15148 , HIS interface software 15150 , and insulin dosing software 15152 .
- the computer system 15146 runs optional data analysis software 15148 that organizes and presents information obtained from one or more monitoring devices.
- the data analysis software 15148 collects and analyzes data from the monitoring devices in an ICU.
- the data analysis software 15148 can also present charts, graphs, and statistics to a user of the computer system 15146 .
- the computer system 15146 runs optional hospital information system (HIS) interface software 15150 that provides an interface point between one or more monitoring devices and an HIS.
- HIS interface software 15150 may also be capable of communicating data between one or more monitoring devices and a laboratory information system (LIS).
- LIS laboratory information system
- the computer system 15146 runs optional insulin dosing software 15152 that provides a platform for implementation of an insulin dosing regimen.
- the hospital tight glycemic control protocol is included in the software. The protocol allows computation of proper insulin doses for a patient connected to a monitoring device 15146 .
- the insulin dosing software 15152 can communicate with the monitoring device 15146 to ensure that proper insulin doses are calculated.
- FIG. 16 is a block diagram of an embodiment of a system 1600 for reducing noise and/or unwanted elements in a signal.
- the system 1600 includes a detector 1602 , such as, for example, a sample detector 930 or a reference detector 936 .
- a detector 1602 such as, for example, a sample detector 930 or a reference detector 936 .
- One type of detector that can serve as a sample detector 930 or a reference detector 936 is a pyroelectric infrared detector. It is typical for pyroelectric detectors to also be sensitive to vibrations. Thus, for example, the output of a pyroelectric infrared detector is the sum of the exposure to infrared radiation and to vibrations of the detector.
- the sensitivity to vibrations can introduce a large noise component to the measurement of radiation using a pyroelectric infrared detector. It is desirable for a spectrometer, such as the analyte detection system 910 , to have a high signal-to-noise ratio, such as a S/N in excess of 100 dB. It can be difficult to achieve this low noise level in the presence of vibrations, even with good mechanical isolation.
- the presently disclosed system includes one or more of the following techniques for reducing the vibrational noise component of measurements using detectors (e.g., the sample detector 930 or the reference detector 936 ).
- Some embodiments for reducing vibrational noise include the use of a modulated infrared source combined with an output filter.
- the analyzer 910 is the infrared source.
- the infrared source of the analyzer 910 can be modulated at a known frequency, and the detector output can be filtered using a narrow band filter centered about the known source frequency.
- the energy source 912 has an energy output that is sine-wave modulated at 10 Hz, and the output of the detector(s) (e.g., the detectors 930 , 936 ) is filtered using a narrow-band pass filter.
- the narrow-band pass filter can have a frequency of less than 1 Hz and be centered about 10 Hz, for example.
- Microphonic (also referred to as vibration-induced) signals that are not exactly 10 Hz can be significantly attenuated with this arrangement.
- the detector output is filtered using signal conditioning 1604 , including a synchronous demodulator and digital filter.
- the demodulator can be a software component implemented in the signal processing computer. Synchronous demodulators, coupled with low pass filters, can be referred to as “lock-in amplifiers.”
- FIG. 17 illustrates an example of signal conditioning 1604 including a lock-in amplifier system 1710 .
- An input device 1712 such as, for example, the analyzer 910 of FIG. 9 or the detector 1602 of FIG. 16 , is shown at the left.
- the input device 1712 is part of a circuit 1713 that provides one or more inputs into the lock-in amplifier system 1710 .
- the circuit may contain one or more resistors R 1 , R 2 and may also have a connection to ground.
- the circuit 1713 provides one or more inputs into a differential AC amplifier 1714 .
- the embodiment of the differential AC amplifier 1714 shown in FIG. 17 accepts inputs from ground and from the input device 1712 .
- the differential AC amplifier 1714 outputs an amplified signal waveform having a substantial frequency component f s substantially equal to the frequency of the infrared source of the analyzer 910 .
- the signal waveform may also have other frequency components that correspond to, for example, noise.
- the output of the differential AC amplifier 1714 is supplied to a function device 1716 .
- the circuit 1713 also provides an input into a phase-lock loop (PLL) 1718 .
- the PLL 1718 conditions the signal and sends a reference waveform having a frequency component f r to the function device 1716 .
- the frequency f r of the predominant component of the reference waveform is substantially equal to the frequency of the infrared source of the analyzer 910 .
- the function device 1716 outputs the product of the inputs of the device.
- the function device 1716 multiplies the reference waveform and the signal waveform.
- the product of the reference waveform and the signal waveform has substantial components at frequencies f s ⁇ f r and f s +f r that are proportional to the amplitude of the signal at the detector output. Because f s and f r are substantially equal, the resulting waveform will have a 0 Hz component (e.g., a DC component or a demodulated component) and a component at about 2f r that are proportional to the amplitude of the signal at the detector output.
- a 0 Hz component e.g., a DC component or a demodulated component
- the resulting waveform is fed through a low pass filter 1720 , which can eliminate substantially all components of the signal above a threshold frequency. For example, if the pass-band of the filter 1720 is small enough (e.g., if the threshold frequency is low enough), the filter 1720 can substantially remove the portions of the waveform that are above the demodulated component at 0 Hz. This resulting waveform can then pass to a DC amplifier 1722 , which provides an output from the lock-in amplifier system 1710 .
- the center frequency of a lock-in amplifier system 1710 is 10 Hz and the integration time is 2 seconds, yielding an approximately 0.5 Hz wide pass band.
- the use of such a device excludes a significant amount of vibration induced noise outside of the 9.75 to 10.25 Hz frequency band.
- the energy present in the pass band generates the signal of the input device 1712 .
- the optical source can be a radiation source that is modulated at 10 Hz.
- the optical source can have all of its energy centered in the pass band, and it can proceed through the demodulator and filter unattenuated.
- a large part of the vibration-induced noise is eliminated because most of it falls outside of the pass band.
- FIG. 18 graphically depicts example signal levels, in volts, of the 10 Hz signal of an input device 1712 (e.g., from the analyzer 910 of FIG. 9 ).
- the illustrated signal can originate from a modulated infrared radiation source that passes through a medium and is then detected by an optical detector.
- the illustrated signal can be a detector signal before it reaches the differential AC amplifier 1714 of FIG. 17 .
- Signals such as these can be produced, for example, as an analyzer 910 cycles through various filters 915 in a filter wheel 918 .
- the example shown here depicts 13 different signal regions corresponding to 13 different filters. Shown are relatively large AC signals centered approximately around zero volts.
- FIG. 19 graphically depicts a signal that can be produced by the same analyzer 910 of FIG. 9 when the original modulated infrared radiation source is blocked.
- FIG. 19 also illustrates the signal output from a detector before that signal reaches the differential AC amplifier 1714 of FIG. 17 .
- the illustrated signal (output from a detector) is due only to vibrations of the detector.
- the top portion of FIG. 19 shows how, when a system (e.g., the apparatus 100 illustrated in FIG. 1 ) is not moving (stationary), the output signal is steady.
- the bottom portion of FIG. 19 shows how, when the system is placed on a rolling platform and is rolled across the floor, there is some signal due to the microphonics of the detector.
- the magnitude of the microphonic signal is generally much less than the magnitude of the optical signal depicted in FIG. 18 . Nevertheless, even this small amount of microphonic signal may produce an output that has an unacceptably low signal-to-noise ratio.
- FIG. 20 graphically depicts signals after they have passed through a lock-in amplifier such as the one depicted in FIG. 17 (e.g., a demodulator and filter), resulting in a DC output.
- the signals shown in FIG. 20 can originate from the analyzer 910 of FIG. 9 when the original modulated infrared radiation source is blocked. Thus, the optical signal is removed, leaving only noise or interference signals.
- the top signal shows, the noise level is very low when the apparatus 100 is stationary.
- the bottom signal shows a higher noise level resulting from microphonic interference when the apparatus 100 is non-stationary (e.g., rolling across a hospital floor).
- the two signals shown in FIG. 20 are derived from the two signals of FIG.
- a lock-in amplifier such as the one depicted in FIG. 17 (e.g., a demodulator and filter).
- the pre-filtering “rolling” microphonic signals of FIG. 19 have a peak-to-peak amplitude on the order of 5+ mV; as shown in FIG. 20 , however, after demodulation and filtering, the demodulated voltage is typically 0.5 mV.
- a lock-in amplifier system (or, for example, a sync demodulator and a digital filter) can result in a 10:1 reduction in microphonics.
- the optical signal level (S) can be approximately 3 V peak-to-peak.
- the typical microphonics free (no vibration) noise level (N) is approximately 0.03 mV after the signal passes through the lock-in amplifier system 1710 .
- the signal-to-noise ratio (S/N) is relatively low. For example, it meets the 100 dB S/N test, which can be a requirement in some embodiments.
- the apparatus 100 is not stationary (e.g., when it rolls across a floor), the lower portion of FIG.
- lock-in amplifier 1710 can condition the signal and remove some noise, as described above with respect to FIG. 17 , some microphonic effects can still persist in the signal even after the signal has passed through the lock-in amplifier 1710 . These effects can increase the error in the measurements provided by the apparatus 100 .
- microphonic effects can create such a noisy signal that the signal passes beyond a threshold of acceptability.
- rolling the apparatus 100 of FIG. 1 over a rough floor can result in an unacceptable error.
- an unacceptable standard error 1 ⁇
- An example of a threshold level, beyond which errors are unacceptable in some embodiments, is 5 mg/dL of glucose.
- Other threshold error levels can also be selected.
- additional steps can be advantageously taken to prevent the apparatus 100 from reporting an incorrect value.
- a “hold off” can be included with or in the hardware or software of an apparatus 100 .
- the “hold off” can prevent the system from taking and/or recording a measurement (e.g., a spectroscopic measurement) when a noise threshold is exceeded.
- a “hold off” can be activated when a non-noise threshold is exceeded.
- the threshold can be a value having the same units as the final output of the system (e.g., mg/dL).
- the system can be designed to track when a certain guaranteed accuracy level is or is not being met, and when it is not being met, the system can automatically stop providing data regarding glucose levels or concentrations, for example.
- the “hold off” can comprise sensing accelerations of the apparatus 100 and preventing the apparatus 100 from measuring during periods of excessive vibrations.
- vibrations are measured using one or more accelerometers that can be located, for example, near pyroelectric detectors (e.g., the detectors 930 and 936 ).
- the accelerometers sense a vibration that exceeds a pre-determined value
- the analyzer 910 is instructed to stop making spectroscopic measurements until the vibration has subsided to below the threshold level.
- the pre-determined threshold value of acceleration can be selected so that the error in analyte measurement is at an acceptable level for accelerations below the pre-determined threshold value.
- FIG. 21 shows a block diagram of a system 2110 for dealing with noise in an analyte detection environment.
- the illustrated system 2110 can be used with the apparatus 100 of FIG. 1 , for example.
- a debugging module 2114 can comprise an RJ-45 chip.
- the debugging module 2114 can be connected to an algorithm engine 2118 , which can comprise a computer chip, a digital signal processor, and/or a field-programmable gate array, for example.
- the algorithm engine 2118 can function as the algorithm processor 416 described above.
- the algorithm engine can communicate (via a custom interface 2120 , for example) with a spectrometer control device 2122 that can preferably, among other things, convert analog signals to digital format.
- FIG. 21 also shows how several other components of the system 2110 can connect to the spectrometer control device 2122 .
- a temperature sensor 2124 can all connect to (and, in some embodiments, be controlled by) the spectrometer control device 2122 .
- vibration sensor(s) 2140 (which can be accelerometers, for example) and detector(s) 2150 (which can be the input device 1712 , for example) can both feed signals to the spectrometer control device 2122 .
- these signals are fed in with little or no delay (e.g., in “real time”).
- the spectrometer control device 2122 can monitor and quantify the output from the vibration sensors 2140 , which can allow a user or system to detect the vibration. This can, in turn, indicate that a significant portion of the detector signal is due to vibration, which can trigger a “hold off” procedure such as that described above. After a “hold off” has occurred, the spectrometer control device 2122 preferably causes a re-measurement or a resumption of measurement.
- FIG. 22 when monitoring accelerometers, it can be useful to gather information related not only to a vibration above a predetermined threshold value but also to a signal within the pass band of the detector signal, the 10 Hz signal produced by the detector. This can be done as vibration sensors (e.g., the vibration sensor(s) 2140 and/or the accelerometer 2212 ) are monitored and quantified as described above. To gather such information, an accelerometer 2212 can feed a signal through a connection 2214 into a “lock in amplifier” style demodulator 2216 similar to (or exactly like) the lock-in amplifier 1710 used for the detector signal (see FIG. 17 ). FIG. 22 shows schematically how such a system can be set up.
- the input device 1712 can be a detector.
- FIG. 22 illustrates how the PLL output 2218 can be taken from the detector channel and fed directly into the accelerometer demodulation multiplier 2220 . This same result can be accomplished by software, rather than with the hardware configuration illustrated.
- a system is designed to detect when the vibration-induced error is so large that the instrument (e.g., the apparatus 100 ) will not be able to achieve a given accuracy specification. Including one or multiple accelerometers can help achieve this goal.
- an accelerometer measures forces in three axes. To convert that measurement into glucose error a linear regression calibration can be used to estimate an error (e.g., Optical Density error, or ODe) based on the three accelerometer signals for the three axes.
- an error e.g., Optical Density error, or ODe
- Some embodiments of a noise removal method include performing a regression analysis to determine when vibration-induced errors in parameter values exceed acceptable limits. For example, calibration computations based on vibration-induced glucose concentration errors can allow a monitoring device to predict when vibrations measured by the one or more accelerometers will produce errors in glucose concentration that exceed acceptable limits. In some embodiments, measurements taken during an initial period, such as the first 30 seconds of apparatus operation, can be used for calibration and the remaining measurements can be used to predict when calculated parameters have an unacceptable amount of error.
- Vibration errors can be more harmful during some measurements than others.
- vibration-induced noise can be especially difficult to filter and/or detect and correct.
- maximum error in the optical density (OD) measurement can coincide with the measurement of the 9.22 micron filter. This can occur, for example, when that filter has the highest coefficient in terms of mg/dL per ODe when compared to the other filters.
- spectroscopic measurements of a sample be completed within a certain period of time (the “maximum measurement period”).
- a maximum measurement period can be selected based on the stability of the sample.
- one system can produce accurate glucose measurements if a spectroscopic measurement can be completed within 90 seconds.
- the “maximum measurement period” in this case is 90 seconds.
- the cumulative hold-off time for making a complete measurement is specifically designed to be less than the maximum measurement period.
- vibrations prevent accurate measurement from being made for some short period of time
- measurement can be allowed to continue if that short period of time still permits scanning through all the filters (or a pre-determined sub-set of filters) within the maximum measurement period.
- vibrations While measurement may be periodically interrupted during times of vibration (e.g., when the accelerometers indicate acceleration beyond a threshold value), there may still be enough time to complete a full measurement cycle within the maximum measurement period of 90 seconds.
- the vibrations may be aborted and restarted later.
- a spectroscopic scan of analyzer 910 takes approximately 2-3 seconds per filter (see, e.g., filters 915 of FIG. 9 ), and there can be 25-30 filters in a system, in some embodiments. Some embodiments require that all filters be scanned without significant change in the sample to meet accuracy guidelines for analyte (e.g., glucose) computation. Typically, without any vibration interference, all filters are scanned in 50 to 90 seconds. In some embodiments, the sample meets certain stability requirements during this time to achieve a glucose measurement with error below a predetermined accuracy, which may be, for example, 5 mg/dL. If the acceleration is unacceptably high for a short period of time, for example less than 5 seconds, then the measurement using a particular filter can be temporarily stopped. If the acceleration is unacceptably high and a complete filter scan cannot be made within the period for a stable sample, for example 90 seconds, then the entire scan can be restarted and completed within the stability period of the sample.
- analyte e.g., glucose
- any single delay of up to approximately 5 seconds can be tolerated for a single filter before having to restart a measurement cycle.
- the measurement can “hold off” on that single filter measurement, subsequently completing the measurement and moving on to the next when the vibration stops. Holding off on only one or two filters in this way can add only a few seconds to the measurement time and is likely not unacceptable to a user.
- the system can wait (in some embodiments, the system is required to wait) until the vibration has stopped and then restart the entire 25-30 filter scans. Thus, another 50-90 seconds may be required before the analyte measurement can be completed.
- biological factors can also be pertinent to the timing settings of the system. For example, clotting, aggregation, or other biological processes can cause blood to plug the flow cell if allowed to remain stagnant in the flow cell for too long. In some embodiments, 10 minutes is too long. Therefore if a vibration period strong enough to trigger measurement “hold off” lasts so long that the 10 min. blood holding period will be exceeded, the entire measurement process (including drawing new blood for measurement, in some embodiments) is preferably re-started. In some embodiments, the minimum processing time of the instrument is approximately 10 minutes, so in such a case the glucose reading will be refreshed 10 minutes after the vibration ceases.
- each of the methods described herein may include a computer program accessible to and/or executable by a processing system, e.g., a one or more processors and memories that are part of an embedded system.
- a processing system e.g., a one or more processors and memories that are part of an embedded system.
- embodiments of the disclosed inventions may be embodied as a method, an apparatus such as a special purpose apparatus, an apparatus such as a data processing system, or a carrier medium, e.g., a computer program product.
- the carrier medium carries one or more computer readable code segments for controlling a processing system to implement a method.
- various ones of the disclosed inventions may take the form of a method, an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects.
- any one or more of the disclosed methods may be stored as one or more computer readable code segments or data compilations on a carrier medium.
- a carrier medium Any suitable computer readable carrier medium may be used including a magnetic storage device such as a diskette or a hard disk; a memory cartridge, module, card or chip (either alone or installed within a larger device); or an optical storage device such as a CD or DVD.
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Heart & Thoracic Surgery (AREA)
- Optics & Photonics (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Chemical & Material Sciences (AREA)
- Immunology (AREA)
- General Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Emergency Medicine (AREA)
- Investigating Or Analysing Biological Materials (AREA)
Abstract
Description
- This application claims priority to U.S. Provisional Patent Application No. 60/791,931, filed Apr. 11, 2006, the entirety of which is hereby incorporated by reference and made part of this specification.
- 1. Field
- Some embodiments disclosed herein relate to methods and apparatus for determining the concentration of an analyte in a sample, such as an analyte in a sample of bodily fluid, as well as methods and apparatus which can be used to support the making of such determinations.
- 2. Description of the Related Art
- It is advantageous to measure the levels of certain analytes, such as glucose, in a bodily fluid, such as blood. This can be done in a hospital or clinical setting when there is a risk that the levels of certain analytes may move outside a desired range, which in turn can jeopardize the health of a patient. Currently known systems for analyte monitoring in a hospital or clinical setting suffer from various drawbacks.
- Embodiments described herein have several features, no single one of which is solely responsible for their desirable attributes. Without limiting the scope of the invention as expressed by the claims, some of the advantageous features will now be discussed briefly.
- Some embodiments use a synchronous demodulator and digital filter to reduce microphonic signal content. Some embodiments monitor the microphonic signal content and “holds off” on making a measurement until vibrations subside. In some embodiments, monitoring is performed using an accelerometer or other vibration sensor. In another embodiment, an algorithm is used to examine the detector output signal and detect excessive microphonic components thereby eliminating the need for the accelerometers.
- In some embodiments, a method for determining the concentration of an analyte in a sample comprises: providing an optical detector signal from an optical detector, the signal having information relating to the concentration of an analyte in a sample; passing the optical detector signal through a first demodulator; providing a vibration sensor signal, the vibration signal having information relating to the vibration of the optical detector; passing the vibration sensor signal through a second demodulator; providing a threshold value for the vibration sensor signal that is calculated to correspond to an accuracy parameter of a system concentration output; and controlling measurement by the optical detector based on whether or not the threshold value is exceeded.
- Some embodiments comprise a system for improving accuracy of a mobile analyte concentration measurement apparatus. The system can comprise: a sample detector configured to provide a detector output signal; a first signal conditioner configured to receive the detector output signal; a vibration sensor mounted to detect vibration of the sample detector and configured to provide a sensor output signal; a second signal conditioner configured to receive the sensor output signal; and a controller configured to communicate with the sample detector and the vibration sensor and prevent the sample detector from detecting when the sensor output signal exceeds a threshold value.
- Some embodiments comprise a method of mitigating the effects of vibration on an optical analyte detection system measurements. The method can comprise: monitoring microphonic effects to determine when a threshold is exceeded; and automatically delaying optical measurement until after the threshold is no longer exceeded. In some embodiments, monitoring microphonic effects comprises processing a detector output signal. In some embodiments, monitoring microphonic effects comprises monitoring output from a vibration sensor that is physically associated with a detector. In some embodiments, monitoring output from a vibration sensor further comprises monitoring output from an accelerometer.
- The following drawings and the associated descriptions are provided to illustrate embodiments of the present disclosure and do not limit the scope of the claims.
-
FIG. 1 shows an embodiment of an apparatus for withdrawing and analyzing fluid samples; -
FIG. 2 illustrates how various other devices can be supported on or near an embodiment of apparatus illustrated inFIG. 1 ; -
FIG. 3 illustrates an embodiment of the apparatus inFIG. 1 connected to a patient; -
FIG. 4 is a block diagram of an embodiment of a system for withdrawing and analyzing fluid samples; -
FIG. 5 schematically illustrates an embodiment of a fluid system within a system for withdrawing and analyzing fluid samples; -
FIG. 6 is an oblique schematic depiction of an embodiment of a modular monitoring device; -
FIG. 7 shows a cut-away side view of an embodiment of a monitoring device; -
FIG. 8 illustrates an embodiment of a disposable cartridge that can interface with a fluid system; -
FIG. 9 schematically illustrates an embodiment of an optical system that comprises a spectroscopic analyzer adapted to measure spectra of a fluid sample; -
FIG. 10 is a flowchart that schematically illustrates an embodiment of a spectroscopic method for determining the concentration of an analyte of interest in a fluid sample; -
FIG. 11 is a flowchart that schematically illustrates an embodiment of a method for estimating the concentration of an analyte in the presence of interferents; -
FIG. 12 is a flowchart that schematically illustrates an embodiment of a method for performing a statistical comparison of the absorption spectrum of a sample with the spectrum of a sample population and combinations of individual library interferent spectra; -
FIG. 13 is a flowchart that schematically illustrates an example embodiment of a method for estimating analyte concentration in the presence of the possible interferents; -
FIGS. 14A and 14B schematically illustrate the visual appearance of embodiments of a user interface for a system for withdrawing and analyzing fluid samples; -
FIG. 15 schematically depicts various components and/or aspects of a patient monitoring system and the relationships among the components and/or aspects; -
FIG. 16 is a block diagram of a system for reducing noise. -
FIG. 17 is a diagram of a lock-in amplifier system that can help reduce noise. -
FIG. 18 shows example signal levels of a signal from an input device. -
FIG. 19 shows a signal output from a detector before that signal reaches a lock-in amplifier. -
FIG. 20 shows signals after they have passed through a lock-in amplifier system such as the one depicted inFIG. 17 . -
FIG. 21 shows a block diagram of a system for dealing with noise in an analyte detection environment. -
FIG. 22 illustrates an output can be taken from a detector channel and fed directly into an accelerometer demodulation multiplier. - Although certain preferred embodiments and examples are disclosed below, the inventive subject matter extends beyond the specifically disclosed embodiments to other alternative embodiments and/or uses of the invention, and to modifications and equivalents thereof. Thus, the scope of the inventions herein disclosed is not limited by any of the particular embodiments described below. For example, in any method or process disclosed herein, the acts or operations of the method or process may be performed in any suitable sequence and are not necessarily limited to any particular disclosed sequence. For purposes of contrasting various embodiments with the prior art, certain aspects and advantages of these embodiments are described. Of course, it is to be understood that not necessarily all such aspects or advantages are achieved by any particular embodiment. Thus, for example, it should be recognized that the various embodiments may be carried out in a manner that achieves or optimizes one advantage or group of advantages as taught herein without necessarily achieving other aspects or advantages as may be taught or suggested herein. The systems and methods discussed herein can be used anywhere, including, for example, in laboratories, hospitals, healthcare facilities, intensive care units (ICUs), or residences. Moreover, the systems and methods discussed herein can be used for invasive techniques, as well as non-invasive techniques or techniques that do not involve a body or a patient.
-
FIG. 1 shows an embodiment of anapparatus 100 for withdrawing and analyzing fluid samples. Theapparatus 100 includes amonitoring device 102. In some embodiments, themonitoring device 102 can be an “OptiScanner®,” available from OptiScan Biomedical Corporation of Hayward, Calif. In some embodiments, theapparatus 100 can measure one or more physiological parameters, such as the concentration of one or more substance(s) in a sample fluid. The sample fluid can be, for example, whole blood from a patient 302 (see, e.g.,FIG. 3 ). In some embodiments, theapparatus 100 can also deliver an infusion fluid to thepatient 302. - In the illustrated embodiment, the
monitoring device 102 includes adisplay 104 such as, for example, a touch-sensitive liquid crystal display. Thedisplay 104 can provide an interface that includes alerts, indicators, charts, and/or soft buttons. Thedevice 102 also can include one or more inputs and/oroutputs 106 that provide connectivity. - In the embodiment shown in
FIG. 1 , thedevice 102 is mounted on astand 108. Thestand 108 can be easily moved and includes one ormore poles 110 and/or hooks 112. Thepoles 110 and hooks 112 can be configured to accommodate other medical implements, including, for example, infusion pumps, saline bags, arterial pressure sensors, other monitors and medical devices, and so forth. -
FIG. 2 illustrates how various other devices can be supported on or near theapparatus 100 illustrated inFIG. 1 . For example, thepoles 110 of thestand 108 can be configured (e.g., of sufficient size and strength) to accommodatemultiple devices monitoring device 102 such that multiple devices function in concert for the benefit of one or multiple patients without the devices interfering with each other. -
FIG. 3 illustrates theapparatus 100 ofFIG. 1 as it can be connected to apatient 302. Themonitoring device 102 can be used to determine the concentration of one or more substances in a sample fluid. The sample fluid can come from a fluid container in a laboratory setting, or it can come from apatient 302, as illustrated here. In some preferred embodiments, the sample fluid is whole blood. - In some embodiments, the
monitoring device 102 can also deliver an infusion fluid to thepatient 302. An infusion fluid container 304 (e.g., a saline bag), which can contain infusion fluid (e.g., saline and/or medication), can be supported by thehook 112. Themonitoring device 102 can be in fluid communication with both thecontainer 304 and the sample fluid source (e.g., the patient 302), throughtubes 306. The infusion fluid can comprise any combination of fluids and/or chemicals. Some advantageous examples include (but are not limited to): water, saline, dextrose, lactated Ringer's solution, drugs, and insulin. - The illustrated
monitoring device 102 allows the infusion fluid to pass to thepatient 302 and/or uses the infusion fluid itself (e.g., as a flushing fluid or a standard with known optical properties, as discussed further below). In some embodiments, themonitoring device 102 may not employ infusion fluid. Themonitoring device 102 may thus draw samples without delivering any additional fluid to thepatient 302. Themonitoring device 102 can include, but is not limited to, fluid handling and analysis apparatuses, connectors, passageways, catheters, tubing, fluid control elements, valves, pumps, fluid sensors, pressure sensors, temperature sensors, hematocrit sensors, hemoglobin sensors, colorimetric sensors, gas (e.g., “bubble”) sensors, fluid conditioning elements, gas injectors, gas filters, blood plasma separators, and/or communication devices (e.g., wireless devices) to permit the transfer of information within themonitoring device 102 or between themonitoring device 102 and a network. - In some embodiments, one or more components of the
apparatus 100 can be located at another facility, room, or other suitable remote location. One or more components of themonitoring device 102 can communicate with one or more other components of the monitoring device 102 (or with other devices) by communication interface(s) such as, but not limited to, optical interfaces, electrical interfaces, and/or wireless interfaces. These interfaces can be part of a local network, internet, wireless network, or other suitable networks. - System Overview
-
FIG. 4 is a block diagram of asystem 400 for withdrawing and analyzing fluid samples. Themonitoring device 102 can comprise such a system. Thesystem 400 includes afluid source 402 connected to afluid system 404. Thefluid system 404 prepares fluid samples that are analyzed by anoptical system 412. Thesystem 400 includes adisplay controller 414 and analgorithm processor 416 that assist in fluid sample analysis and presentation of data. In some embodiments, the sampling andanalysis system 400 is a mobile point of care apparatus that monitors physiological parameters such as, for example, blood glucose concentration. Tubes and connectors within thesystem 400 can be coated with an antibacterial coating to reduce the risk of infection. Connectors between at least some components of thesystem 400 can include a self-sealing valve, such as a spring valve, in order to reduce the risk of contact between port openings and fluids, and to guard against fluid escaping from the system. -
Fluid Source 402 - The sampling and
analysis system 400 includes afluid source 402 that contains fluid to be sampled. Thefluid system 404 of the sampling andanalysis system 400 is connected to afluid source 402 from which fluid samples can be drawn. Thefluid source 402 can be, for example, a patient's blood vessel such as a vein or an artery, a container such as a decanter or a tube, or any other corporeal or extracorporeal fluid source. The fluid to be sampled can be, for example, blood, plasma, or another bodily fluid. -
Fluid System 404 - In some embodiments, the
fluid system 404 withdraws a sample of fluid from thefluid source 402 for analysis, centrifuges at least a portion of the sample, and prepares at least a portion of the sample for analysis by an optical sensor such as a spectrophotometer. In some embodiments, at least a portion of the sample is returned to thefluid source 402. At least some of the sample, such as portions of the sample that are mixed with other materials or portions that are otherwise altered during the sampling and analysis process, can also be placed in a waste bladder. The waste bladder can be integrated within thefluid system 404 or supplied by a user ofsystem 400. Thefluid system 404 can also be connected to a saline source, a detergent source, and/or an anticoagulant source, each of which can be supplied by a user or integrated withinfluid system 404. - Components of the
fluid system 404 can be modularized into one or more non-disposable, disposable, and/or replaceable subsystems. In the embodiment shown inFIG. 4 , components of thefluid system 404 are separated into anon-disposable subsystem 406, a firstdisposable subsystem 408, and a seconddisposable subsystem 410. - The
non-disposable subsystem 406 can include components that do not generally require regular replacement during the useful lifetime of thesystem 400. In some embodiments, thenon-disposable subsystem 406 of thefluid system 404 includes one or more reusable valves and sensors. For example, thenon-disposable subsystem 406 can include one or more pinch valves (or non-disposable portions thereof), ultrasonic bubble sensors, non-contact pressure sensors, and optical blood dilution sensors. Thenon-disposable subsystem 406 can also include one or more pumps (or non-disposable portions thereof). In some embodiments, the components of thenon-disposable subsystem 406 are not directly exposed to fluids and/or are not readily susceptible to contamination. - First and second
disposable subsystems system 400. For example, the firstdisposable subsystem 408 can be replaced after a certain period of use, such as a few days, has elapsed. Replacement may be necessary, for example, when a bladder within the firstdisposable subsystem 408 is filled to capacity. Such replacement may mitigate fluid system performance degradation associated with and/or contamination wear on system components. - In some embodiments, the first
disposable subsystem 408 includes components that may contact fluids such as patient blood, saline, flushing solutions, anticoagulants, and/or detergent solutions. For example, the firstdisposable subsystem 408 can include one or more tubes, fittings, cleaner pouches and/or waste bladders. The components of the firstdisposable subsystem 408 can be sterilized in order to decrease the risk of infection and can be configured to be easily replaceable. - In some embodiments, the second
disposable subsystem 410 can be designed to be replaced under certain circumstances. For example, the seconddisposable subsystem 410 can be replaced when the patient being monitored by thesystem 400 is changed. The components of the seconddisposable subsystem 410 may not need replacement at the same intervals as the components of the firstdisposable subsystem 408. For example, the seconddisposable subsystem 410 can include a flow cell and/or at least some components of a centrifuge, components that may not become filled or quickly worn during operation of thesystem 400. Replacement of the seconddisposable subsystem 410 can decrease or eliminate the risk of transferring fluids from one patient to another during operation of thesystem 400, enhance the measurement performance ofsystem 400, and/or reduce the risk of contamination or infection. - In some embodiments, the flow cell of the second
disposable subsystem 410 receives the sample obtained from thefluid source 402 via the fluidics of the firstdisposable subsystem 408. The flow cell is a container that can hold fluid for the centrifuge and provide a window to the sample for analysis by a spectrometer. In some embodiments, the flow cell includes windows that are made of a material that is substantially transparent to electromagnetic radiation in the mid-infrared range of the spectrum. For example, the flow cell windows can be made of calcium fluoride. - An injector can provide a fluidic connection between the first
disposable subsystem 408 and the flow cell. In some embodiments, the injector can be removed from the flow cell to allow for free spinning of the flow cell during centrifugation. - In some embodiments, the components of the sample are separated by centrifuging at a high speed for a period of time before measurements are performed by the
optical system 412. For example, a blood sample can be centrifuged at 7200 RPM for 2 minutes in order to separate plasma from other blood components for analysis. Separation of a sample into the components can permit measurement of solute (e.g., glucose) concentration in plasma, for example, without interference from other blood components. This kind of post-separation measurement, (sometimes referred to as a “direct measurement”) has advantages over a solute measurement taken from whole blood because the proportions of plasma to other components need not be known or estimated in order to infer plasma glucose concentration. - An anticoagulant, such as, for example, heparin can be added to the sample before centrifugation to prevent clotting. The
fluid system 404 can be used with a variety of anticoagulants, including anticoagulants supplied by a hospital or other user of themonitoring system 400. A detergent solution formed by mixing detergent powder from a pouch connected to thefluid system 404 with saline can be used to periodically clean residual protein and other sample remnants from one or more components of thefluid system 404, such as the flow cell. Sample fluid to which anticoagulant has been added and used detergent solution can be transferred into the waste bladder. -
Optical System 412 - The
system 400 shown inFIG. 4 includes anoptical system 412 that can measure optical properties (e.g., transmission) of a fluid sample (or a portion thereof). In some embodiments, theoptical system 412 measures transmission in the mid-infrared range of the spectrum. In some embodiments, theoptical system 412 includes a spectrometer that measures the transmission of broadband infrared light through a portion of a flow cell filled with fluid. The spectrometer need not come in direct contact with the sample. As used herein, the term “flow cell” is a broad term that carries its ordinary meaning as an object that can provide a place for fluid. The fluid can enter the flow cell by flowing. - In some embodiments, the
optical system 412 includes a filter wheel that contains one or more filters. In some embodiments, twenty-five filters are mounted on the filter wheel. Theoptical system 412 includes a light source that passes light through a filter and the flow cell to a detector. In some embodiments, a stepper motor moves the filter wheel in order to position a selected filter in the path of the light. An optical encoder can also be used to finely position one or more filters. -
Display Controller 414 - The
system 400 shown inFIG. 4 includes adisplay controller 414 that provides for communication of information to a user of thesystem 400. Thedisplay controller 414 can include a display processor that controls or produces an interface to communicate information to the user. Thedisplay controller 414 can include a display screen. One or more parameters such as, for example, blood glucose concentration,system 400 operating parameters, and/or other operating parameters can be displayed on a monitor (not shown) associated with thesystem 400. An example of one way such information can be displayed is shown inFIGS. 14A and 14B . In some embodiments, thedisplay controller 414 can communicate measured physiological parameters and/or operating parameters to a computer system over a communications connection. -
Algorithm Processor 416 - The
system 400 shown inFIG. 4 includes analgorithm processor 416 that can receive optical density (OD) values (or other analog or digital optical data) from theoptical system 412. In some embodiments, thealgorithm processor 416 calculates one or more physiological parameters by adjusting the coefficients of a model, if necessary, and computing the physiological parameters using an equation having the adjusted coefficients. Thealgorithm processor 416, thedisplay controller 414, and any embedded controllers withinsystem 400 can be connected to one another with a communications bus. - Fluidics System
-
FIG. 5 schematically illustrates afluid system 510. In addition to the reference numerals used below, the various portions of the illustratedfluid system 510 are labeled with letters to suggest their role as follows: T# indicates a section of tubing. C# indicates a connector that joins multiple tubing sections. V# indicates a valve. BS # indicates a bubble sensor or ultrasonic air detector. N# indicates a needle (e.g., a needle that injects sample into a flow cell). PS# indicates a pressure sensor (e.g., a reusable pressure sensor). Pump# indicates a fluid pump (e.g., a syringe pump with a disposable body and reusable drive). “Hb 12” indicates a sensor for hemoglobin (e.g., a dilution sensor that can detect hemoglobin optically). - At the start of a measurement cycle, most lines, including the patient tube 512 (T1), can be filled with saline that can be introduced into the system through the
tubes infusion pump 518 and/or asaline bag 520. Theinfusion pump 518 and thesaline bag 520 can be provided separately from thesystem 510. For example, a hospital can use existing saline bags and infusion pumps to interface with the described system. Thevalve 521 can be open to allow saline to flow into the tube 512 (T1). - To draw a sample, a first pump 522 (pump #1) draws sample fluid to be analyzed (e.g. blood) from a fluid source (e.g., a laboratory sample container, a living patient, etc.) up into the patient tube 512 (T1), through the open valve F23 (V0), through the first connector 524 (C1), past the hemoglobin sensor 526 (Hb12), and into the looped tube 528 (T4). During this process, the valve 529 (V7 a) is open to fluid flow, but the valves 531 (V1 a) and 533 (V3 a) can be closed and therefore block (or substantially block) fluid flow.
- Initially the lines are filled with saline and the hemoglobin (Hb) level is zero. The tubes that are filled with saline are in fluid communication with the a sample source (not shown). The sample source can be the vessels of a living human or a pool of liquid in a laboratory sample container, for example. When the saline is drawn toward the
first pump 522, fluid to be analyzed is also drawn into the system because of the suction forces in the closed fluid system. Thus, thefirst pump 522 draws a relatively continuous column of fluid that first comprises generally nondiluted saline, then a mixture of saline and sample fluid (e.g., blood), and then eventually nondiluted sample fluid. In the example illustrated here, the sample fluid is blood. - The hemoglobin sensor 526 (Hb12) detects the level of Hemoglobin in the sample fluid. As blood starts to arrive at the hemoglobin sensor 526 (Hb12), the hemoglobin level rises. When the hemoglobin level reaches a preset value (e.g., which can occur after a draw of approximately 2 mL depending on the size of the catheter used) there is a nondiluted sample present at the first connector 524 (C1). A nondiluted sample can be, for example, a blood sample that is not diluted with saline solution, but instead has the characteristics of the rest of the blood flowing through a patient's body. A loop of tubing 530 (e.g., a 1-mL loop) can be advantageously positioned as illustrated to help insure that undiluted fluid (e.g., undiluted blood) is present at the first connector 524 (C1) when the
hemoglobin sensor 526 registers that the preset Hb threshold is crossed. The loop oftubing 530 provides additional length to the tube 528 (T4) to make it less likely that the portion of the fluid column in the tubing at the first connector 524 (C1) has advanced all the way past the mixture of saline and sample fluid, and the nondiluted blood portion of that fluid has reached the first connector 524 (C1). - When nondiluted blood is present at the first connector 524 (C1), a second pump 532 (pump #0) draws four “slugs” of blood into the tubing 534 (T3). As used herein, the term “slug” refers to a continuous column of fluid. Slugs can be separated from one another by injecting (or sucking in) small amounts of air to create bubbles at intervals in the tube. In the illustrated embodiment, blood slugs are alternated with air bubbles by maintaining the valve 523 (V0) closed, maintaining the valve 533 (V3 a) open, and alternately closing and opening the valves 529 (V7 a) and 531 (V1 a) such that one is closed while the other one is open. This periodically pulls either one or the other of 1) blood from the tube 528 (T4) through the valve 529 (V7 a) and 2) air from the tube 536 (T2) through the valve 531 (V1 a). In some embodiments, four blood slugs are created. The first three blood slugs are approximately 15 μL and the fourth is approximately 35 μL.
- As, or after, the slugs are formed, heparin can be inserted into each slug. A heparin vial 538 (e.g., an insertable vial provided independently by the user of the system 510) can be connected to a
tube 540. Ashuttle valve 541 can connect to both thetube 540 and the tube 534 (T3). The valve can open thetube 540 to a suction force (e.g., created by the pump 532), allowing heparin to be drawn from thevial 538 into thevalve 541. Then, theshuttle valve 541 can slide the heparin over into fluid communication with thetube 534. Theshuttle valve 541 can then return to its previous position. Thus, heparin can be shuttled from thetube 540 into the tube 534 (T3) such that each blood slug contains a precisely controlled amount of heparin. - Following the formation of four blood slugs, the majority of the sampled blood is returned to the patient. The first pump 522 (pump #1) pushes the blood out of the tube 528 (T4) and back to the patient by opening the valve 523 (V0), closing the valves 531 (V1 a) and 533 (V3 a), and keeping the valve 529 (V7 a) open. The tube 528 (T4) is preferably flushed with approximately 2 mL of saline. This can be accomplished by closing the valve 529 (V7 a), opening the valve 542 (PV1), drawing saline from the
saline source 520 into thetube 544, closing the valve 542 (PV1), opening the valve 529 (V7 a), and forcing the saline down the tube 528 (T4) with thepump 522. - In some embodiments, less than two minutes elapses between the time that blood is drawn from the patient and the time that the blood is returned to the patient after formation of the blood slugs.
- Following return of the unused patient blood sample, the four slugs are pushed up the tube 534 (T3), through the second connector 546 (C2), and into the
flow cell 548, which can be located on thecentrifuge wheel 550. The bubble sensor 552 (BS14) can identify the fourth slug by identifying and counting how many air bubbles (or inter-slug spaces) pass by the sensor. The fourth slug can be identified, and thepump 522 can stop forcing the fluid column through thetube 534 so that the fourth slug remains within theflow cell 548. Thus, the first three blood slugs can serve to flush any residual saline out theflow cell 548. The three leading slugs can be deposited in thewaste bladder 554 by passing through the tube F56 (T6) and through the valve 557 (V4 a). - In some embodiments, the fourth blood slug is centrifuged for two minutes at 7200 RPM. This separates the whole blood into its components, isolates the plasma, and positions the plasma in the
flow cell 548 for measurement. Thecentrifuge 550 is stopped with theflow cell 548 in a beam of radiation (not shown) for analysis. The radiation, a detector, and logic can be used to analyze the plasma spectroscopically (e.g., for glucose and/or lactate concentration). - Following analysis, the second pump 532 (pump #0) flushes the
flow cell 548 and sends the flushed contents to thewaste bladder 554. This flush can be done with a cleaning solution from theterg tank 558. In some embodiments, thesecond pump 532 is in fluid communication with the tube 560 (T9) and theterg tank 558 because the valve 559 (V7 b) is open. Thesecond pump 532 forces cleaning solution from theterg tank 558 through theopen valve 561 and the tube 562 (T7) when thevalve 559 is open. The cleaning flush can pass through theflow cell 548, through thesecond connector 546, through the tube 564 (T5) and the open valve 563 (V2 b), and into thewaste bladder 554. Following this flush, - Subsequently, the first pump 522 (pump #1) can flush the cleaning solution out of the
flow cell 548 using saline in drawn from thesaline bag 520. This flush pushes saline through the tube 528 (T4), the tube 534 (T3), theflow cell 548, and the tube 556 (T6). Thus, in some embodiments, the following valves are open for this flush: 529 (V7 a), 533 (V3 a), 557 (V4 a), and the following valves are closed: 542 (PV1), 523 (V0), 531 (V1 a), 566 (V3 b), 563 (V2 b), and 561 (V4 b). - When the fluid source is a living entity such as a patient, in between measurements, a low flow of saline (e.g., 1-5 mL/hr) is preferably moved through the patient tube 512 (T1) and into the patient to keep the patient's vessel open (e.g., to establish a keep vessel open, or “KVO” flow). The source of this KVO flow can be the
infusion pump 518, the third pump 568 (pump #3), or the first pump 522 (pump #1). In some embodiments, theinfusion pump 518 can run continuously throughout the measurement cycle described above. This continuous flow can advantageously avoid any alarms that may be triggered if theinfusion pump 518 senses that the flow has stopped or changed in some other way. In some embodiments, when thevalve 521 closes to allow pump 522 (pump #1) to withdraw fluid from a fluid source (e.g., a patient), the third pump 568 (pump #3) can withdraw fluid through theconnector 570, thus allowing theinfusion pump 518 to continue pumping normally as if the fluid path was not blocked by thevalve 521. If the measurement cycle is about two minutes long, this withdrawal by thethird pump 568 can continue for approximately two minutes. Once thevalve 521 is open again, the third pump 568 (pump #3) can reverse and insert the saline back into the system at a low flow rate. Preferably, the time between measurement cycles is longer than the measurement cycle itself (e.g., longer than two minutes). Accordingly, thethird pump 568 can insert fluid back into the system at a lower rate than it withdrew that fluid. This can help prevent an alarm by the infusion pump. - Mechanical Fluidics Interface
-
FIG. 6 is an oblique schematic depiction of amodular monitoring device 600. Themodular monitoring device 600 includes abody portion 602 havingreceptacles receptacles disposable cassettes first cassette 610 can be used to store at least a portion of thefluid system 510 described previously, including portions that contact sample fluids, saline, detergent solution, and/or anticoagulant. - In some embodiments, a
non-disposable fluidics subsystem 608 is disposed within thebody portion 602 of themonitoring device 600. Thefirst cassette 610 can include one or more openings that allow portions of thenon-disposable fluidics subsystem 608 to interface with thecassette 610. For example, thenon-disposable fluidics subsystem 608 can include one or more pinch valves that are designed to extend through such openings to engage one or more sections of tubing. When thefirst cassette 610 is inserted into a correspondingfirst receptacle 604, actuation of the pinch valves can selectively close sections of tubing within the cassette. Thenon-disposable fluidics subsystem 608 can also include one or more sensors that interface with connectors, tubing sections, or pumps located within thefirst cassette 610. - In the embodiment shown in
FIG. 6 , themonitoring device 600 includes anoptical system 614 disposed within thebody portion 602. Theoptical system 614 can include a light source and a detector that are adapted to perform measurements on fluids within a flow cell. In some embodiments, the flow cell is disposed within asecond cassette 612. Thesecond cassette 612 can include an optical window through which theoptical system 614 can shine radiation for measuring properties of a fluid in the flow cell when the cassette is inserted into a correspondingsecond receptacle 606. Theoptical system 614 can include other components (some of which may interface with the second cassette 612) such as, for example, a power supply, a centrifuge motor, a filter wheel, and/or a beam splitter. - In some embodiments, the
first cassette 610 and thesecond cassette 612 are adapted to be in fluid communication with each other. For example, thefirst cassette 610 can include a retractable injector that injects fluids into a flow cell disposed in thesecond cassette 612. In some embodiments, the injector can be retracted to allow the centrifuge to rotate the flow cell freely. In other embodiments, a fluid communication path can be provided by components disposed within thebody portion 602 of themonitoring device 600. - The
body portion 602 of themonitoring device 600 can also include one or more connectors for an external battery (not shown). The external battery can serve as a backup emergency power source in the event that a primary emergency power source such as, for example, an internal battery (not shown) is exhausted. -
FIG. 7 shows a cut-away side view of a monitoring device 700 (which can correspond, for example, to thedevice 102 shown inFIG. 1 ). Thedevice 700 includes acasing 702 that can include one or more receptacles. Depicted inFIG. 7 are examples of ways in which components of thedevice 700 mounted within thecasing 702 can interact with components of thedevice 700 disposed within cassettes inserted into the receptacles. Not all components of thedevice 700 are shown inFIG. 7 . - A
first cassette 704 having a variety of components is shown inserted into a receptacle formed in thecasing 702. Asecond cassette 706 is also inserted into a receptacle. Components mounted within the cassettes are indicated with dashed lines inFIG. 7 , while components mounted within thecasing 702 are depicted with solid lines. - In some embodiments, one or
more actuators 708 housed within thecasing 702 operatesyringe pumps 710 located within thefirst cassette 704. Thepumps 710 are connected to sections oftubing 716 that move fluid among various components of the system. The movement of fluid is at least partially controlled by the action of one ormore pinch valves 712 positioned within thecasing 702. Thepinch valves 712 have arms 714 that extend within thefirst cassette 704. Movement of the arms 714 can constrict a section oftubing 716 in order to create an effective seal. - In some embodiments, the
second cassette 706 includes aflow cell 720 that engages acentrifuge motor 718 mounted within thecasing 702 of thedevice 700 when the cassette is inserted into a receptacle. Afilter wheel motor 722 disposed within thehousing 702 rotates afilter wheel 724 in order to align a filter with a window of theflow cell 720. An optical light path including alight source 726 within thehousing 702 routes a beam of infrared light through the filter and theflow cell 720. Adetector 728 measures the optical density of the light transmitted through the filter and flowcell 720. -
FIG. 8 illustrates adisposable cartridge 800 that can interface with a fluid system such as thefluid system 510 ofFIG. 5 . Thedisposable cartridge 800 can be configured for insertion into a receptacle of thedevice 700 shown inFIG. 7 . In some embodiments, thecartridge 800 includes one or more features that ease insertion of thecartridge 800 into a corresponding receptacle. For example, thecartridge 800 can be shaped so as to promote insertion of thecartridge 800 in the correct orientation. Thecartridge 800 can also include labeling or coloring affixed to or integrated with the cartridge's exterior casing that help a handler insert thecartridge 800 into a receptacle properly. - The
cartridge 800 can include one or more ports for connecting to material sources. For example, oneport 802 can be configured to attach to ananticoagulant source 804. Other ports can be provided to connect to, for example, a saline source, an infusion pump, a sample source, and/or a source of nitrogen gas. The ports can be connected to sections of tubing within thecartridge 800. In some embodiments, the sections of tubing are opaque or covered so that fluids within the tubing cannot be seen. - The
cartridge 800 shown inFIG. 8 includes one or more injector needles 806. The injector needles 806 can be configured to inject at least a portion of a sample into a flow cell (not shown). The housing of thecartridge 800 can include atubing space 808 for one or more sections of tubing. In some embodiments, the body of thecartridge 800 includes one ormore apertures 809 through which various components, such as, for example, pinch valves and sensors, can interface with the fluidics contained in thecartridge 800. The sections of tubing found in thetubing space 808 can be aligned with theapertures 809 in order to implement at least some of the functionality shown in thefluid system 510 ofFIG. 5 . - The
cartridge 800 can include a pouch space 810 for storing one or more components of thefluid system 510. For example, one or more pouches and/or bladders can be disposed in the pouch space 810. In some embodiments, a cleaner pouch and a waste bladder are housed in the pouch space 810. The waste bladder can be placed under the cleaner pouch such that, as detergent is removed from the cleaner pouch, the waste bladder has more room to fill. The components placed in the pouch space 810 can also be placed side-by-side or in any other suitable configuration. The pouch space 810 can be isolated from the rest of thecartridge 800 by one ormore walls 811. One ormore connectors fluid system 510. - The
cartridge 800 can include one ormore pumps 816 that facilitate movement of fluid within thefluid system 510. Each of thepumps 816 can be, for example, a syringe pump having a plunger. The plunger can include aportion 818 configured to interface with an actuator housed outside thecartridge 800. For example, theportion 818 of the pump that interfaces with an actuator can be exposed to the exterior of thecartridge 800 housing by one or more apertures in the housing. - In some embodiments, the
disposable cartridge 800 is designed for single patient use. Thecartridge 800 may also be designed for replacement after a period of operation. For example, in some embodiments, if thecartridge 800 is installed in a continuously operating monitoring device that performs four measurements per hour, the waste bladder may become filled or the detergent in the cleaner pouch depleted after about three days. Thecartridge 800 can be replaced before the detergent and waste bladder are exhausted. - The
cartridge 800 can be configured for easy replacement. For example, in some embodiments, thecartridge 800 is designed to have an installation time of only several minutes. For example, the cartridge can be designed to be installed in less than about five minutes. During installation, various portions of the fluidics contained in thecartridge 800 can be primed by automatically filling the fluidics with saline. The saline can be mixed with detergent powder from the cleaner pouch in order to create a cleaning solution. - The
cartridge 800 can also be designed to have a relatively brief shut down time. For example, the shut down process can be configured to take less than about five minutes. The shut down process can include flushing the patient line; sealing off the insulin pump connection, the saline source connection, and the sample source connection; and taking other steps to decrease the risk that fluids within the usedcartridge 800 will leak after disconnection from the monitoring device. - In some embodiments, the
cartridge 800 is designed to fit within standard waste containers found in a hospital, such as a standard biohazard container. For example, thecartridge 800 can be less than one foot long, less than one foot wide, and less than two inches thick. In some embodiments, thecartridge 800 is designed to withstand a substantial impact, such as that caused by hitting the ground after a four foot drop, without damage to the housing or internal components. In some embodiments, thecartridge 800 is designed to withstand significant clamping force applied to its casing. For example, thecartridge 800 can be built to withstand five pounds per square inch of force without damage. In some embodiments, thecartridge 800 is non pyrogenic and/or latex free. - Spectroscopy
- As described above with reference to
FIG. 4 , thesystem 400 comprises theoptical system 412 for analysis of a fluid sample. In various embodiments, theoptical system 412 comprises one or more optical components including, for example, a spectrometer, a photometer, a reflectometer, or any other suitable device for measuring optical properties of the fluid sample. Theoptical system 412 may perform one or more optical measurements on the fluid sample including, for example, measurements of transmittance, absorbance, reflectance, scattering, and/or polarization. The optical measurements may be performed in one or more wavelength ranges including, for example, infrared (IR) and/or optical wavelengths. As described with reference toFIG. 4 (and further described below), the measurements from theoptical system 412 are communicated to thealgorithm processor 416 for analysis. For example, in one embodiment thealgorithm processor 416 computes concentration of analyte(s) (and/or interferent(s)) of interest in the fluid sample. Analytes of interest include, e.g., glucose and lactate in whole blood or blood plasma. -
FIG. 9 schematically illustrates an embodiment of theoptical system 412 that comprises aspectroscopic analyzer 910 adapted to measure spectra of a fluid sample such as, for example, blood or blood plasma. Theanalyzer 910 comprises anenergy source 912 disposed along an optical axis X of theanalyzer 910. When activated, theenergy source 912 generates an electromagnetic energy beam E, which advances from theenergy source 912 along the optical axis X. In certain embodiments, theenergy source 912 comprises an infrared energy source, and the energy beam E comprises an infrared beam. In some embodiments, the infrared energy beam E comprises a mid-infrared energy beam or a near-infrared energy beam. In certain embodiments, the energy beam E may include optical and/or radio frequency wavelengths. - The
energy source 912 may comprise a broad-band and/or a narrow-band source of electromagnetic energy. In some embodiments, theenergy source 912 comprises optical elements such as, e.g., filters, collimators, lenses, mirrors, etc., that are adapted to produce a desired energy beam E. For example, in some embodiments, the energy beam E is an infrared beam in a wavelength range between about 2 μm and 20 μm. In certain embodiments, the energy beam E comprises an infrared beam in a wavelength range between about 4 μm and 10 μm. In the infrared wavelength range, water generally is the main contributor to the total absorption together with features from absorption of other blood components, particularly in the 6 μm-10 μm range. The 4 μm to 10 μm wavelength band has been found to be advantageous for determining glucose concentration, because glucose has a strong absorption peak structure from about 8.5 μm to 10 μm, whereas most other blood components have a relatively low and flat absorption spectrum in the 8.5 μm to 10 μm range. Two exceptions are water and hemoglobin, which are interferents in this range. - The energy beam E may be temporally modulated to provide increased signal-to-noise ratio (S/N) of the measurements provided by the
analyzer 910 as further described below. For example, in some embodiments, the beam E is modulated at a frequency of about 10 Hz or in a range from about 1 Hz to about 30 Hz. Asuitable energy source 912 may be an electrically modulated thin-film thermoresistive element such as the HawkEye IR-50 available from Hawkeye Technologies of Milford, Conn. - As depicted in
FIG. 9 , the energy beam E propagates along the optical axis X and passes through anaperture 914 and afilter 915 thereby providing a filtered energy beam Ef. Theaperture 914 helps collimate the energy beam E and may include one or more filters adapted to reduce the filtering burden of thefilter 915. For example, theaperture 914 may comprise a broadband filter that substantially attenuates beam energy outside a wavelength band between about 4 μm to about 10 μm. Thefilter 915 may comprise a narrow-band filter that substantially attenuates beam energy having wavelengths outside of a filter passband (which may be tunable or user-selectable in some embodiments). The filter passband may be specified by a half-power bandwidth (“HPBW”). In some embodiments, thefilter 915 may have an HPBW in a range from about 0.01 μm to about 1 μm. In one embodiment, the bandwidths are in a range from about 0.1 μm to 0.35 μm. Other filter bandwidths may be used. Thefilter 915 may comprise a varying-passband filter, an electronically tunable filter, a liquid crystal filter, an interference filter, and/or a gradient filter. In some embodiments, thefilter 915 comprises one or a combination of a grating, a prism, a monochrometer, a Fabry-Perot etalon, and/or a polarizer. Other optical elements as known in the art may be utilized as well. - In the embodiment shown in
FIG. 9 , theanalyzer 910 comprises afilter wheel assembly 921 configured to dispose one ormore filters 915 along the optical axis X. Thefilter wheel assembly 921 comprises afilter wheel 918, afilter wheel motor 916, and aposition sensor 920. Thefilter wheel 918 may be substantially circular and have one ormore filters 915 or other optical elements (e.g., apertures, gratings, polarizers, etc.) disposed around the circumference of thewheel 918. In some embodiments, the number offilters 915 in thefilter wheel 916 may be, for example, 1, 2, 5, 10, 15, 20, 25, or more. Themotor 916 is configured to rotate thefilter wheel 918 to dispose a desired filter 915 (or other optical element) in the energy beam E so as to produce the filtered beam Ef. In some embodiments, themotor 916 comprises a stepper motor. Theposition sensor 920 determines the angular position of thefilter wheel 916, and communicates a corresponding filter wheel position signal to thealgorithm processor 416, thereby indicating which filter 915 is in position on the optical axis X. In various embodiments, theposition sensor 920 may be a mechanical, optical, and/or magnetic encoder. An alternative to thefilter wheel 918 is a linear filter translated by a motor. The linear filter may include an array of separate filters or a single filter with properties that change along a linear dimension. - The
filter wheel motor 916 rotates thefilter wheel 918 to position thefilters 915 in the energy beam E to sequentially vary the wavelengths or the wavelength bands used to analyze the fluid sample. In some embodiments, eachindividual filter 915 is disposed in the energy beam E for a dwell time during which optical properties in the passband of the filter are measured for the sample. Thefilter wheel motor 916 then rotates thefilter wheel 918 to position anotherfilter 915 in the beam E. In one embodiment, 25 narrow-band filters are used in thefilter wheel 918, and the dwell time is about 2 seconds for eachfilter 915. A set of optical measurements for all the filters can be taken in about 2 minutes, including sampling time and filter wheel movement. In some embodiments, the dwell time may be different fordifferent filters 915, for example, to provide a substantially similar S/N ratio for each filter measurement. Accordingly, thefilter wheel assembly 921 functions as a varying-passband filter that allows optical properties of the sample to be analyzed at a number of wavelengths or wavelength bands in a sequential manner. - In certain embodiments of the
analyzer 910, thefilter wheel 918 includes 25 finite-bandwidth infrared filters having a Gaussian transmission profile and full-width half-maximum (FWHM) bandwidth of 28 cm−1 corresponding to a bandwidth that varies from 0.14 μm at 7.08 μm to 0.28 μm at 10 μm. The central wavelength of the filters are, in microns: 7.082, 7.158, 7.241, 7.331, 7.424, 7.513, 7.605, 7.704, 7.800, 7.905, 8.019, 8.150, 8.271, 8.598, 8.718, 8.834, 8.969, 9.099, 9.217, 9.346, 9.461, 9.579, 9.718, 9.862, and 9.990. - With further reference to
FIG. 9 , the filtered energy beam Ef propagates to abeamsplitter 922 disposed along the optical axis X. Thebeamsplitter 922 separates the filtered energy beam Ef into a sample beam Es and a reference beam Er. The reference beam Er propagates along a minor optical axis Y, which in this embodiment is substantially orthogonal to the optical axis X. The energies in the sample beam Es and the reference beam Er may comprise any suitable fraction of the energy in the filtered beam Ef. For example, in some embodiments, the sample beam Es comprises about 80%, and the reference beam Er comprises about 20%, of the filtered beam energy Ef.A reference detector 936 is positioned along the minor optical axis Y. Anoptical element 934, such as a lens, may be used to focus or collimate the reference beam Er onto thereference detector 936. Thereference detector 936 provides a reference signal, which can be used to monitor fluctuations in the intensity of the energy beam E emitted by thesource 912. Such fluctuations may be due to drift effects, aging, wear, or other imperfections in thesource 912. Thealgorithm processor 416 may utilize the reference signal to identify changes in properties of the sample beam Es that are attributable to changes in the emission from thesource 912 and not to the properties of the fluid sample. By so doing, theanalyzer 910 may advantageously reduce possible sources of error in the calculated properties of the fluid sample (e.g., concentration). In other embodiments of theanalyzer 910, thebeamsplitter 922 is not used, and substantially all of the filtered energy beam Ef propagates to the fluid sample. - As illustrated in
FIG. 9 , the sample beam Es propagates along the optical axis X, and arelay lens 924 transmits the sample beam Es into asample cell 948 so that at least a fraction of the sample beam Es is transmitted through at least a portion of the fluid sample in thesample cell 948. Asample detector 930 is positioned along the optical axis X to measure the sample beam Es that has passed through the portion of the fluid sample. Anoptical element 928, such as a lens, may be used to focus or collimate the sample beam Es onto thesample detector 930. Thesample detector 930 provides a sample signal that can be used by thealgorithm processor 416 as part of the sample analysis. - In the embodiment of the
analyzer 910 shown inFIG. 9 , thesample cell 948 comprises the flow cell 648 located toward the circumference of thecentrifuge wheel 650. The flow cell 648 comprises windows that are substantially transmissive to energy in the sample beam Es. For example, in implementations using mid-infrared energy, the windows may comprise calcium fluoride. As described herein with reference toFIG. 5 , the flow cell 648 is in fluid communication with an injector system that permits filling the flow cell 648 with a fluid sample (e.g., whole blood) and flushing the flow cell 648 (e.g., with saline or a detergent). The injector system may disconnect after filling the flow cell 648 with the fluid sample to permit free spinning of thecentrifuge wheel 650 bycentrifuge motor 926. In certain embodiments of theanalyzer 910, the fluid sample (e.g., a whole blood sample) is spun at about 7200 rpm for about 2 minutes to separate blood plasma for spectral analysis. In some embodiments, an anti-clotting agent such as heparin may be added to the fluid sample before centrifuging to reduce clotting. - The embodiment of the
analyzer 910 illustrated inFIG. 9 advantageously permits direct measurement of the concentration of analytes in the plasma sample rather than by inference of the concentration from measurements of a whole blood sample. An additional advantage is that relatively small volumes of fluid may be spectroscopically analyzed. For example, in certain embodiments the fluid sample volume is between about 1 μL and 80 μL and is about 25 μL in some embodiments. In certain embodiments, the flow cell 648 is disposable and is intended for use with a single patient or for a single measurement. - In certain embodiments, the
reference detector 936 and thesample detector 930 comprise broadband pyroelectric detectors. As known in the art, some pyroelectric detectors are sensitive to vibrations. Thus, for example, the output of a pyroelectric infrared detector is the sum of the exposure to infrared radiation and to vibrations of the detector. The sensitivity to vibrations, also known as “microphonics,” can introduce a noise component to the measurement of the reference and sample energy beams Er, Es using some pyroelectric infrared detectors. Because it may be desirable for theanalyzer 910 to provide high signal-to-noise ratio measurements, such as, e.g., S/N in excess of 100 dB, some embodiments of theanalyzer 910 utilize one or more vibrational noise reduction apparatus or methods. For example, theanalyzer 910 may be mechanically isolated so that high S/N spectroscopic measurements can be obtained for vibrations below an acceleration of about 1.5 G. - In some embodiments of the
analyzer 910, vibrational noise can be reduced by using a temporally modulatedenergy source 912 combined with an output filter. In certain embodiments, theenergy source 912 is modulated at a known source frequency, and measurements made by thedetectors source 912 is sinusoidally modulated at 10 Hz, and outputs of thedetectors algorithm processor 416. Synchronous demodulators, coupled with low pass filters, are often referred to as “lock in amplifiers.” - The
analyzer 910 may also include a vibration sensor 932 (e.g., one or more accelerometers) disposed near one (or both) of thedetectors vibration sensor 932 is monitored, and suitable actions are taken if the measured vibration exceeds a vibration threshold. For example, in some embodiments, if thevibration sensor 932 detects above-threshold vibrations, the system discards any ongoing measurement and “holds off” on performing further measurements until the vibrations drop below the threshold. Discarded measurements may be repeated after the vibrations drop below the vibration threshold. In some embodiments, if the duration of the “hold off” is sufficiently long, the fluid in thesample cell 930 is flushed, and a new fluid sample is delivered to thecell 930 for measurement. The vibration threshold may be selected so that the error in analyte measurement is at an acceptable level for vibrations below the threshold. In some embodiments, the threshold corresponds to an error in glucose concentration of 5 mg/dL. The vibration threshold may be determined individually for eachfilter 915. - Certain embodiments of the
analyzer 910 include a temperature system (not shown inFIG. 9 ) for monitoring and/or regulating the temperature of system components (such as thedetectors 936, 930) and/or the fluid sample. Such a temperature system may include temperature sensors, thermoelectrical heat pumps (e.g., a Peltier device), and/or thermistors, as well as a control system for monitoring and/or regulating temperature. In some embodiments, the control system comprises a proportional-plus-integral-plus-derivative (PID) control. For example, in certain embodiments, the temperature system is used to regulate the temperature of thedetectors - The
analyzer 910 illustrated inFIG. 9 can be used to determine optical properties of a substance in thesample cell 948. The substance may include whole blood, plasma, saline, water, air or other substances. In some embodiments, the optical properties include measurements of an absorbance, transmittance, and/or optical density in the wavelength passbands of some or all of thefilters 915 disposed in thefilter wheel 918. As described above, a measurement cycle comprises disposing one ormore filters 915 in the energy beam E for a dwell time and measuring a reference signal with thereference detector 936 and a sample signal with thesample detector 930. The number offilters 915 used in the measurement cycle will be denoted by N, and eachfilter 915 passes energy in a passband around a center wavelength ξi, where i is an index ranging over the number of filters (e.g., from 1 to N). The set of optical measurements from thesample detector 936 in the passbands of the N filters 915 provide a wavelength-dependent spectrum of the substance in thesample cell 948. The spectrum will be denoted by Cs(ξi), where Cs may be a transmittance, absorbance, optical density, or some other measure of an optical property of the substance. In some embodiments, the spectrum is normalized with respect to one or more of the reference signals measured by thereference detector 930 and/or with respect to spectra of a reference substance (e.g., air or saline). The measured spectra are communicated to thealgorithm processor 416 for calculation of the concentration of the analyte(s) of interest in the fluid sample. - In certain embodiments, the
analyzer 910 performs spectroscopic measurements on the fluid sample (known as a “wet” reading) and on one or more reference samples. For example, an “air” reading occurs when thesample detector 936 measures the sample signal without thesample cell 948 in place along the optical axis X. A “water” or “saline” reading occurs when thesample cell 948 is filled with water or saline, respectively. Thealgorithm processor 416 may be programmed to calculate analyte concentration using a combination of these spectral measurements. - In some embodiments, a pathlength corrected spectrum is calculated using wet, air, and reference readings. For example, the transmittance at wavelength ξi, denoted by Ti, may be calculated according to Ti=(Si(wet)/Ri(wet))/(Si(air)/Ri(air)), where Si denotes the sample signal from the
sample detector 936 and Ri denotes the corresponding reference signal from thereference detector 930. In certain embodiments, thealgorithm processor 416 calculates the optical density, ODi, as a logarithm of the transmittance, e.g., according to ODi=−Log(Ti). In one implementation, theanalyzer 910 takes a set of wet readings in each of the N filter passbands and then takes a set of air readings in each of the N filter passbands. In other embodiments, theanalyzer 910 may take an air reading before (or after) the corresponding wet reading. - The optical density ODi is the product of the absorption coefficient at wavelength ξi,αi, times the pathlength L over which the sample energy beam Es interacts with the substance in the
sample chamber 948, e.g., ODi=αiL. The absorption coefficient αi of a substance may be written as the product of an absorptivity per mole times a molar concentration of the substance.FIG. 9 schematically illustrates the pathlength L of thesample cell 948. The pathlength L may be determined from spectral measurements made when thesample cell 948 is filled with a reference substance. For example, because the absorption coefficient for water (or saline) is known, one or more water (or saline) readings can be used to determine the pathlength L from measurements of the transmittance (or optical density) through thecell 948. In some embodiments, several readings are taken in different wavelength passbands, and a curve-fitting procedure is used to estimate a best-fit pathlength L. The pathlength L may be estimated using other methods including, for example, measuring interference fringes of light passing through anempty sample cell 948. - The pathlength L may be used to determine the absorption coefficients of the fluid sample at each wavelength. Molar concentration of an analyte of interest can be determined from the absorption coefficient and the known molar absorptivity of the analyte. In one embodiment, a sample measurement cycle comprises a saline reading (at one or more wavelengths), a set of N wet readings, followed by a set of N air readings. As discussed above, the sample measurement cycle can be performed in about 2 minutes when the filter dwell times are about 2 seconds. After the sample measurement cycle is completed, a detergent cleaner may be flushed through the flow cell 648 to reduce buildup of organic matter (e.g., proteins) on the windows of the flow cell 648. The detergent is then flushed to a waste bladder.
- In some embodiments, the system stores information related to the spectral measurements so that the information is readily available for recall by a user. The stored information may include wavelength-dependent spectral measurements (including fluid sample, air, and/or saline readings), computed analyte values, system temperatures and electrical properties (e.g., voltages and currents), and any other data related to use of the system (e.g., system alerts, vibration readings, S/N ratios, etc.). The stored information may be retained in the system for a time period such as, for example, 30 days. After this time period, the stored information may be communicated to an archival data storage system and then deleted from the system. In certain embodiments, the stored information is communicated to the archival data storage system via wired or wireless methods, e.g., over a hospital information system (HIS).
- Algorithm(s)
- The algorithm processor 416 (
FIG. 4 ) (or any other suitable processor) may be configured to receive from theanalyzer 910 the wavelength-dependent optical measurements of the fluid sample. In some embodiments, the optical densities ODi in each of the N filter passbands centered around wavelengths ξi are communicated to theprocessor 416, which analyzes the optical densities to measure and quantify one or more analytes in the presence of interferents. Interferents can comprise components of a material sample being analyzed for an analyte, where the presence of the interferent affects the quantification of the analyte. Thus, for example, in the spectroscopic analysis of a sample to determine an analyte concentration, an interferent could be a compound having spectroscopic features that overlap with those of the analyte. The presence of such an interferent can introduce errors in the quantification of the analyte. More specifically, the presence of interferents can affect the sensitivity of a measurement technique to the concentration of analytes of interest in a material sample, especially when the system is calibrated in the absence of, or with an unknown amount of, the interferent. - Independently of or in combination with the attributes of interferents described above, interferents can be classified as being endogenous (i.e., originating within the body) or exogenous (i.e., introduced from or produced outside the body). As an example of these classes of interferents, consider the analysis of a blood sample (or a blood component sample or a blood plasma sample) for the analyte glucose. Endogenous interferents include those blood components having origins within the body that affect the quantification of glucose, and may include water, hemoglobin, blood cells, and any other component that naturally occurs in blood. Exogenous interferents include those blood components having origins outside of the body that affect the quantification of glucose, and can include items administered to a person, such as medicaments, drugs, foods or herbs, whether administered orally, intravenously, topically, etc.
- Independently of or in combination with the attributes of interferents described above, interferents can comprise components which are possibly, but not necessarily, present in the sample type under analysis. In the example of analyzing samples of blood or blood plasma drawn from patients who are receiving medical treatment, a medicament such as acetaminophen is possibly, but not necessarily, present in this sample type. In contrast, water is necessarily present in such blood or plasma samples.
-
FIG. 10 is a flowchart that schematically illustrates an embodiment of aspectroscopic method 1010 for determining the concentration of an analyte of interest in a fluid sample in the presence of one or more possible interferents. Inblock 1012, spectral measurements of the fluid sample are obtained. For example, as described above with reference toFIG. 9 , theanalyzer 910 may be used to obtain optical measurements Cs(ξi) of the fluid sample in a number N of filter passbands centered around wavelengths ξi. Inblock 1014, quality of the spectral measurements is determined regardless of the concentration of the analyte of interest of the presence of possible interferents. In some embodiments, one or more of poor quality spectral measurements Cs(ξi) may be rejected (e.g., as having a S/N ratio that is too low), and themethod 1010 performed on the remaining, sufficiently high-quality measurements. In other embodiments, additional spectral measurements of the fluid sample are obtained to replace one or more of the poor quality measurements. - In
block 1016, the spectral measurements are tested to determine the possible presence of interferents. For example, the system may utilize spectroscopic signatures of possible interferents to test for their presence. Inblock 1017, if the test determines that no interferents are present or that any possible interferents, if present, are at concentrations below suitable thresholds, themethod 1010 proceeds to block 1022 in which analyte concentration is determined. In one embodiment, analyte concentration is determined using a hybrid linear algorithm (HLA) in which analyte concentration is estimated from measured spectra using one or more calibration coefficients and an offset. If inblock 1017 the test determines that one or more interferents are present at concentrations above threshold, then, inblock 1018, the above-threshold interferents are identified. Themethod 1010 proceeds to block 1020 in which the analyte concentration algorithm is adapted to account for the presence of one or more of the identified interferents. For example, in embodiments using HLA, the calibration coefficients may be adjusted to compensate for the presence of some or all of the identified interferents. Themethod 1010 proceeds to block 1022 in which analyte concentration is determined as further described below. - Certain disclosed analysis methods are particularly effective if each analyte and interferent has a characteristic signature in the measurement (e.g., a characteristic spectroscopic feature), and if the measurement is approximately affine (e.g., includes a linear term and an offset) with respect to the concentration of each analyte and interferent. In such methods, a calibration process is used to determine a set of one or more calibration coefficients and one or more optional offset values that permits the quantitative estimation of an analyte. For example, the calibration coefficients and the offsets may be used to calculate an analyte concentration from spectroscopic measurements of a material sample (e.g., the concentration of glucose in blood plasma). In some of these methods, the concentration of the analyte is estimated by multiplying the calibration coefficient by a measurement value (e.g., an optical density) to estimate the concentration of the analyte. Both the calibration coefficient and measurement can comprise arrays of numbers. For example, in some embodiments, the measurement comprises the spectra Cs(ξi) measured at the wavelengths ξi, and the calibration coefficient and optional offset comprise an array of values corresponding to each wavelength ξi. As described with reference to blocks 1017-1020 of
FIG. 10 , in some embodiments a hybrid linear algorithm (HLA) is used to estimate analyte concentration in the presence of a set of interferents, while retaining a high degree of sensitivity to the desired analyte. The data used to accommodate the random set of interferents may include (a) signatures of each of the members of the family of potential additional substances and (b) the typical quantitative level at which each additional substance, if present, is likely to appear. As described with reference to block 1020, in some embodiments, the calibration constant (and optional offset) are adjusted to minimize or reduce the sensitivity of the calibration to the presence of interferents that are identified as possibly being present in the fluid sample. - In one embodiment, the analyte analysis method uses a set of training spectra each having known analyte concentration(s) and produces a calibration that minimizes the variation in estimated analyte concentration with interferent concentration. The resulting calibration coefficient measures sensitivity of the measurement to analyte concentration(s) and, on average, is not sensitive to interferent concentrations. The training spectra need not include a spectrum from the individual whose analyte concentration is to be determined. That is, the term “training” when used in reference to the disclosed methods does not require training using measurements from the individual whose analyte concentration will be estimated (e.g., by analyzing a bodily fluid sample drawn from the individual).
- Several terms are used herein to describe the analyte analysis process. The term “Sample Population” is a broad term and includes, without limitation, a large number of samples having measurements that are used in the computation of a calibration—in other words, used to train the method of generating a calibration. For an embodiment involving the spectroscopic determination of glucose concentration, the Sample Population measurements can each include a spectrum (analysis measurement) and a glucose concentration (analyte measurement). In one embodiment, the Sample Population measurements are stored in a database, referred to herein as a “Population Database.”
- The Sample Population may or may not be derived from measurements of material samples that contain interferents to the measurement of the analyte(s) of interest. One distinction made herein between different interferents is based on whether the interferent is present in both the Sample Population and the sample being measured, or only in the sample. As used herein, the term “Type-A interferent” refers to an interferent that is present in both the Sample Population and in the material sample being measured to determine an analyte concentration. In certain methods it is assumed that the Sample Population includes only interferents that are endogenous, and does not include any exogenous interferents, and thus Type-A interferents are endogenous. The number of Type-A interferents depends on the measurement and analyte(s) of interest, and may number, in general, from zero to a very large number (e.g., greater than 300). The material sample being measured, for example a fluid sample in the
sample cell 948, may also include interferents that are not present in the Sample Population. - As used herein, the term “Type-B interferent” refers to an interferent that is either: 1) not found in the Sample Population but that is found in the material sample being measured (e.g., an exogenous interferent), or 2) is found naturally in the Sample Population, but is at abnormally high concentrations in the material sample (e.g., an endogenous interferent). Examples of a Type-B exogenous interferent may include medications, and examples of Type-B endogenous interferents may include urea in persons suffering from renal failure. For example, in mid-infrared spectroscopic absorption measurements of glucose in blood (or blood plasma), water is present in all fluid samples, and is thus a Type-A interferent. For a Sample Population made up of individuals who are not taking intravenous drugs, and a material sample taken from a hospital patient who is being administered a selected intravenous drug, the selected drug is a Type-B interferent. In addition to components naturally found in the blood, the ingestion or injection of some medicines or illicit drugs can result in very high and rapidly changing concentrations of exogenous interferents.
- In some embodiment, a list of one or more possible Type-B Interferents is referred to herein as forming a “Library of Interferents,” and each interferent in the library is referred to as a “Library Interferent.” The Library Interferents include exogenous interferents and endogenous interferents that may be present in a material sample due, for example, to a medical condition causing abnormally high concentrations of the endogenous interferent.
-
FIG. 11 is a flowchart that schematically illustrates an embodiment of amethod 1100 for estimating the concentration of an analyte in the presence of interferents. Inblock 1110, a measurement of a sample is obtained, and inblock 1120 data relating to the obtained measurement is analyzed to identify possible interferents to the analyte. Inblock 1130, a model is generated for predicting the analyte concentration in the presence of the identified possible interferents, and inblock 1140 the model is used to estimate the analyte concentration in the sample from the measurement. In certain embodiments of themethod 1100, the model generated inblock 1130 is selected to reduce or minimize the effect of identified interferents that are not present in a general population of which the sample is a member. - An example embodiment of the
method 1100 ofFIG. 11 for the determination of an analyte (e.g., glucose) in a blood sample will now be described. This example embodiment is intended to illustrate various aspects of themethod 1100 but is not intended as a limitation on the scope of themethod 1100 or on the range of possible analytes. In this example, the sample measurement inblock 1110 is an absorption spectrum, Cs(ξi), of a measurement sample S that has, in general, one analyte of interest, glucose, and one or more interferents. As described with reference toFIG. 9 , the absorption spectrum may comprise the set of optical densities ODi measured by theanalyzer 910. In general, the sample S includes Type-A interferents, at concentrations preferably within the range of those found in the Sample Population. - In
block 1120, a statistical comparison of the absorption spectrum of the sample S with a spectrum of the Sample Population and combinations of individual Library Interferent spectra is performed. The statistical comparison provides a list of Library Interferents that are possibly contained in sample S and may include either no Library Interferents or one or more Library Interferents. In this example, inblock 1130, a set of spectra are generated using the spectra of the Sample Population and their respective known analyte concentrations and known spectra of the Library Interferents identified inblock 1120. Inblock 1130, the generated spectra are used to calculate a calibration coefficient κ(ξi) that can be used with the sample measurements Cs(ξi) to provide an estimate of the analyte concentration, gest. Inblock 1140, the estimated analyte concentration is determined. For example, in some embodiments of HLA, the estimated analyte concentration is calculated according to a linear formula: gest=κ(ξi)·Cs(ξi). Because the absorption measurements and calibration coefficients may represent arrays of numbers, the multiplication operation indicated in the preceding formula may comprise an inner product or a matrix product. In some embodiments, the calibration coefficient is determined so as to have reduced or minimal sensitivity to the presence of the identified Library Interferents. - An example embodiment of
block 1120 of themethod 1100 will now be described with reference toFIG. 12 . In this example,block 1120 includes forming a statistical Sample Population model (block 1210), assembling a library of interferent data (block 1220), comparing the obtained measurement and statistical Sample Population model with data for each interferent from an interferent library (block 1230), performing a statistical test for the presence of each interferent from the interferent library (block 1240), and identifying possible interferents that pass the statistical test (block 1250). The acts ofblock 1220 can be performed once or can be updated as necessary. The acts ofblocks - In this example, in
block 1210, a Sample Population Database is formed that includes a statistically large Sample Population of individual spectra taken over the same wavelength range as the sample spectrum, Cs(ξi). The Database also includes an analyte concentration corresponding to each spectrum. For example, if there are P Sample Population spectra, then the spectra in the Database can be represented as C={C1, C2, . . . , CP}, and the analyte concentration corresponding to each spectrum can be represented as g={g1, g2, . . . , gP}. In some embodiments, the Sample Population does not have any of the Library Interferents present, and the material sample has interferents contained in the Sample Population and one or more of the Library Interferents. Stated in terms of Type-A and Type-B interferents, the Sample Population has Type-A interferents, and the material sample has Type-A and may have Type-B interferents. - In some embodiments of
block 1210, the statistical sample model comprises a mean spectrum and a covariance matrix calculated for the Sample Population. For example, if each spectrum measured at N wavelengths ξi is represented by an N×1 array, C, then the mean spectrum, μ, is an N×1 array having values at each wavelength averaged over the range of spectra in the Sample Population. The covariance matrix, V, is calculated as the expected value of the deviation between C and μ and can be written as V=E((C−μ)(C−μ)T), where E(·) represents the expected value and the superscript T denotes transpose. In other embodiments, additional statistical parameters may be included in the statistical model of the Sample Population spectra. - Additionally, a Library of Interferents may be assembled in
block 1220. A number of possible interferents can be identified, for example, as a list of possible medications or foods that might be ingested by the population of patients at issue. Spectra of these interferents can be obtained, and a range of expected interferent concentrations in the blood, or other expected sample material, can be estimated. In certain embodiments, the Library of Interferents includes, for each of “M” interferents, the absorption spectrum of each interferent, IF={IF1, IF2, . . . , IFM}, and a maximum concentration for each interferent, Tmax={Tmax1, Tmax2, . . . , TmaxM). Information in the Library may be assembled once and accessed as needed. For example, the Library and the statistical model of the Sample Population may be stored in a storage device associated with the algorithm processor 416 (FIG. 4 ). - Continuing in
block 1230, the obtained measurement data (e.g., the sample spectrum) and the statistical Sample Population model (e.g., the mean spectrum and the covariance matrix) are compared with data for each interferent from the Library of Interferents in order to determine the presence of possible interferents in the sample (block 1240). In some embodiments, the statistical test for the presence of an interferent comprises the following actions. The measured spectrum of the fluid sample, Cs, is modified for each interferent of the library by analytically subtracting, wavelength-by-wavelength, the spectrum of the interferent. For any of the M interferents in the Library, having an absorption spectrum per unit of interferent concentration, IF, the modified spectrum is given by C′s(T)=Cs−IFT, where T is the interferent concentration. In some embodiments, the interferent concentration is assumed to be in a range from a minimum value, Tmin, to a maximum value, Tmax. The value of Tmin may be zero or, alternatively, be a value between zero and Tmax, such as some fraction of Tmax. - In certain embodiments, the statistical test for determining the presence of possible interferents in
block 1240 further comprises determining a Mahalanobis distance (MD) between the modified spectrum C′s(T) and the statistical model (μ, V) of the Sample Population. The Mahalanobis distance can be calculated from
MD 2(C s −IFT,μ;ρ )=(C′ s(T)−μ)T V −l(C′ s(T)−μ). Eq. (1)
The value of MD2 found from Eq. (1) is referred to herein as the “squared Mahalanobis distance” or the “MD2 score.” The MD2 score is used in various embodiments of the statistical test for determining the presence of an interferent. - In
block 1250, a list of possible interferents may be identified as the particular Library Interferents that pass one or more statistical tests for being present in the sample. One or more tests may be used, alone or in combination, to identify the possible interferents. For example, if a statistical test indicates that the interferent is present in negative concentrations, then this non-physical result is used to exclude the possible interferent from the list of possible interferents. In some embodiments, only the single most probable interferent is included on the list. - In one test embodiment, for each interferent, the concentration T is varied from Tmin to Tmax (e.g., evaluate C′s (T) over a range of values of T in Eq. (1)). If the minimum value of MD (or MD2) in this interval is below a minimum threshold, then the test indicates the probable presence of the interferent in the sample. In some embodiments, the minimum threshold MD2 is chosen relative to quantiles of a χ2 random variable having N degrees of freedom, where N is the number of wavelengths in the spectrum Cs. In some embodiments, the 95% quantile is used as the minimum threshold.
- In another test embodiment, if the MD2 score is above a maximum threshold, then it is probable that the interferent is not actually present or is not present in a large enough concentration to modify the analyte concentration estimate. The maximum threshold generally is empirically determined. In one embodiment, it is found that a maximum threshold value is in a range from about 50 to about 200.
- Another test embodiment includes calculating a probability density that combines a range of probable interferent concentrations and the MD2 score for that interferent. For interferents that are not indicated as being present at negative concentrations and that do not have an MD2 score above the maximum threshold, the probability density ρ(T) is computed, which is given by the product:
ρ(T)=ρχ2N (MD 2(C s −IFT))ρT(T), Eq. (2)
The right-hand-side of Eq. (2) is the product of two probability densities: (1) the χ2 distribution with N degrees of freedom (where N is the number of wavelengths present in the spectral measurements), evaluated at the Mahalanobis score for the difference spectrum Cs−IFT, and (2) the distribution of concentrations T for the interferent. In some embodiments, interferent concentration is assumed to have a log-normal distribution with a value of 95% at the assumed maximum interferent concentration in the fluid and a standard deviation of one half the mean. - An integral of ρ(T) is then computed over a range of possible concentrations T, for example from 0 to infinity, or a smaller range, such as from TMIN= 1/2TOPT to TMAX=2TOPT, to give a “raw probably score” (RPS) for the interferent. The RPS is then compared to a minimum value (Pmin). Possible interferents are identified as interferents having an RPS greater than Pmin. Possible interferents are denoted herein with the variable ξ. In some embodiments, the value of Pmin is empirically determined from an analysis of the measurements. For example, a value of 0.70 may result in a single possible interferent (a “single interferent identification”) and a value of 0.3 may result in three possible interferents (a multiple interferent identification).
- Accordingly, in
block 1250, one or more of the above statistical tests (or other tests as known in the art) are used to determine a list of possible interferents ξ that may be present in the fluid sample. In some embodiments, the list of possible interferents includes only the single most probable interferent. In other embodiments, the list of possible interferents ξ may include each of the interferents in the Library of Interferents. - Returning to
FIG. 11 , themethod 1100 continues inblock 1130 where analyte concentration is estimated in the presence of the possible interferents ξ determined inblock 1250.FIG. 13 is a flowchart that schematically illustrates an example embodiment of the acts ofblock 1130. In block 1310, synthesized Sample Population measurements are generated to form an Interferent Enhanced Spectral Database (IESD). Inblock 1320, the spectra in the IESD are partitioned into a calibration set and a test. Inblock 1330, the calibration set is used to generate a calibration coefficient, and in block 1340, the calibration coefficient is used to estimate the analyte concentration of the test set. Inblock 1350, errors in the estimated analyte concentration of the test set are calculated, and inblock 1360 an average calibration coefficient is calculated based on errors in the test set(s). Inblock 1370, the average calibration coefficient is applied to the measured spectra to determine an estimated single-interferent analyte concentration. - In certain embodiments, the blocks 1310-1360 are performed for each possible interferent ξ to provide a corresponding “single-interferent” average calibration coefficient for each particular interferent. In other embodiments, the blocks 1310-1360 are performed only for the single most probable interferent in the list identified in
block 1250. - In one example embodiment for block 1310, synthesized Sample Population spectra are generated by adding a random concentration of one of the possible interferents ξ to each Sample Population spectrum. These spectra are referred to herein as an Interferent-Enhanced Spectral Database or IESD. In one method, the IESD is formed as follows. A plurality of Randomly-Scaled Single Interferent Spectra (RSIS) are formed by combinations of the interferent ξ having spectrum IFξ multiplied by the maximum concentration Tmaxm, which is scaled by a random factor between zero and one. In certain embodiments, the scaling places the maximum concentration at the 95th percentile of a log-normal distribution in order to generate a wide range of concentrations. In one embodiment, the log-normal distribution has a standard deviation equal to half of its mean value.
- Individual RSIS are then combined independently and in random combinations to form a large family of Combination Interferent Spectra (CIS), with each spectrum within the CIS comprising a random combination of RSIS, selected from the full set of identified Library Interferents. An advantage of this method of selecting the CIS is that it produces adequate variability with respect to each interferent, independently across separate interferents.
- The CIS and replicates of the Sample Population spectra are combined to form the IESD. Since the interferent spectra and the Sample Population spectra may have been obtained from measurements having different optical pathlengths, the CIS may be scaled to the same pathlength as the Sample Population spectra. The Sample Population Database is then replicated R times, where R depends on factors including the size of the Database and the number of interferents. The IESD includes R copies of each of the Sample Population spectra, where one copy is the original Sample Population Data, and the remaining R−1 copies each have one randomly chosen CIS spectra added. Accordingly, each of the IESD spectra has an associated analyte concentration from the Sample Population spectra used to form the particular IESD spectrum. In one embodiment, a 10-fold replication of the Sample Population Database is used for 130 Sample Population spectra obtained from 58 different individuals and 18 Library Interferents. A smaller replication factor may be used if there is greater spectral variety among the Library Interferent spectra, and a larger replication factor may be used if there is a greater number of Library Interferents.
- After forming the IESD in block 1310, the blocks 1320-1350 may be executed to repeatedly combine different spectra of the IESD to statistically average out effects of the interferentξ. For example, in
block 1320, the IESD may be partitioned into two subsets: a calibration set and a test set. Repeated partitioning of the IESD into different calibration and test sets improves the statistical significance of the calibration constant. In some embodiments, the calibration set includes a random selection of some of the IESD spectra, and the test set includes the remaining unselected IESD spectra. In one embodiment, the calibration set includes approximately two-thirds of the IESD spectra. - In
block 1330, the calibration set is used to generate a calibration coefficient for estimating the analyte concentration from a sample measurement. In an implementation in which glucose concentration is to be determined from absorption measurements, a glucose absorption spectrum is obtained and indicated as αG. The calibration coefficient is calculated in some embodiments as follows. Using the calibration set having calibration spectra C={c1, c2, . . . , cn} and corresponding glucose concentration values G={g1, g2, . . . , gn}, glucose-free spectra C′={c′1, c′2, . . . , c′n} are calculated as: c′j=cj−αGgj. The calibration coefficient, κ, is calculated from C′ and αG, according to the following 5 actions: -
- 1) C′ is decomposed into C′=AC′ΔC′BC′, for example, by a singular value decomposition, where the A-factor is an orthonormal basis of column space, or span, of C′;
- 2) AC′is truncated to avoid overfitting to a particular column rank r, based on the sizes of the diagonal entries of Δ(the singular values of C′). The selection of r involves a trade-off between the precision and stability of the calibration, with a larger r resulting in a more precise but less stable solution. In one embodiment, each spectrum c includes 25 wavelengths, and r ranges from 15 to 19;
- 3) The first r columns of AC′ are taken as an orthonormal basis of span(C′);
- 4) The projection from the background is found as the product PC′=AC′AC′ T, e.g., the orthogonal projection onto the span of C′. The complementary, or nulling projection PC′ ⊥=1−PC′, which forms the projection onto the complementary subspace C′⊥, is calculated; and
- 5) The calibration coefficient κ is found by applying the nulling projection to the absorption spectrum of the analyte of interest: κRAW=PC′ ⊥αG and normalizing the calibration coefficient κ=κRAW/κRAW,αG , where the angle brackets , denote the standard inner (or dot) product of vectors. The normalized calibration coefficient produces a unit response for a unit αG spectral input for one particular calibration set.
- In block 1340, the calibration coefficient is used to estimate the analyte concentration for the spectra in the test set. For example, each spectrum of the test set is multiplied by the calibration coefficient κ to calculate an estimated glucose concentration. Since each spectrum in the test set has a known glucose concentration, the error between the calculated and known glucose concentration may be calculated, in
block 1350. - Blocks 1320-1350 may be repeated for a number of different random combinations of calibration sets. The number of combinations may be in a range from hundreds to thousands. In
block 1360, an average calibration constant is calculated from the calibration coefficient and the error from the many calibration and test sets. For example, the average calibration coefficient may be computed as a weighted average of the individual calibration coefficients from the combinations. In one embodiment the weighting is in proportion to an inverse root-mean-square (rms), κave=Σ(κ*rms−2)/Σ(rms−2) for all tests. - In summary, one embodiment of a method of computing a calibration constant based on an identified interferent ξ can be summarized as follows:
-
- 1. Generate synthesized Sample Population spectra by adding the RSIS to raw (interferent-free) Sample Population spectra, thus forming an Interferent Enhanced Spectral Database (IESD). Each spectrum of the IESD is synthesized from one spectrum of the Sample Population, and thus each spectrum of the IESD has at least one associated known analyte concentration
- 2. Separate the spectra of the IESD into a calibration set of spectra and a test set of spectra
- 3. Generate a calibration coefficient based on the calibration set spectra and their associated known analyte concentrations.
- 4. Use the calibration coefficient generated in (3) to calculate the error in the corresponding test set as follows (repeat for each spectrum in the test set):
- a. Multiply (the selected test set spectrum)×(average calibration constant generated in (3)) to generate an estimated glucose concentration
- b. Evaluate the difference between this estimated glucose concentration and the known glucose concentration associated with the selected test spectrum to generate an error associated with the selected test spectrum
- 5. Average the errors calculated in (4) to arrive at a weighted or average error for the current calibration set—test set pair
- 6. Repeat (2) through (5) a number n times, resulting in n calibration coefficients and n average errors
- 7. Compute a “grand average” error from the n average errors and an average calibration coefficient from the n calibration coefficient (preferably weighted averages wherein the largest average errors and calibration coefficient are discounted), to arrive at a calibration coefficient that has reduced or minimal sensitivity to the effect of the identified interferents
- The average calibration coefficient determined in
block 1360 corresponds to a single interferent ξ from the list of possible interferents and is denoted herein as a single-interferent calibration coefficient κavg(κ). Inblock 1370 ofFIG. 13 , the single-interferent calibration coefficient is applied to the measured spectra Cs to determine an estimated, single-interferent analyte concentration, g(ξ)=κavg(ξ)·Cs for the interferent ξ. The blocks 1310-1370 can be repeated for each of the interferents on the list of possible interferents, thereby providing an array of estimated, single-interferent analyte concentrations. As noted above, in some embodiments the blocks 1310-1360 are performed only once for the single most probable interferent on the list (e.g., the array of single-interferent analyte concentrations has a single member). - Returning to block 1140 of
FIG. 11 , the array of single-interferent concentrations, g(ξ) are combined to determine an estimated analyte concentration, gest, for the fluid sample. In certain embodiments, a weighting function p(ξ) is determined for each of the interferents on the list of possible interferents. The weighting function may be normalized to unity, e.g., Σp(ξ)=1. For example, in some embodiments, the Raw Probability Score (RPS) (described above following Eq. (2)) is used in determining the weighting function. In one embodiment, the RPS's determined for the interferents on the list of possible interferents are rescaled to unit probability. The weighting function p(ξ) equals the rescaled RPS and may be calculated according to p(ξ)=RPS(ξ)/(ΣRPS(ξ)), where the sum in the denominator is over all interferents ξ in the list. In other embodiments, different weighting functions can be used. For example, in one embodiment, the weighting function is the same constant value for each interferent. - In certain embodiments, the estimated analyte concentration, gest, is determined by combining the single-interferent estimates, g(ξ), and the weighting functions, p(ξ), to generate a likelihood-weighted average analyte concentration. The likelihood-weighted average concentration may be computed according to gest=Σg(ξ)p(ξ), where the sum is over all possible interferents. By testing the above described likelihood-weighted average method on simulated data, it has been found that the likelihood-weighted average analyte concentration advantageously has reduced errors compared to other methods (e.g., using only a single most probable interferent). In embodiments using a constant value for the weighting functions, the estimated analyte concentration is the arithmetic average of the single interferent concentrations.
- In some embodiments, block 1370 of
FIG. 13 is not performed and instead the estimated analyte concentration is determined inblock 1140 ofFIG. 11 by combining the single interferent calibration coefficients κavg(ξ) (determined in block 1360) into a likelihood weighted average calibration coefficient according to κavg=Σκ(ξ)p(ξ). The estimated analyte concentration is determined from the average calibration coefficient and the spectral sample measurement according to gest=κavg·Cs. These embodiments determine the same estimated analyte concentration because of the linearity of the likelihood weighted average method. - The
algorithm processor 416 may be configured, additionally or alternatively, to implement other methods for determining analyte concentration. For example, in certain embodiments, a parameter-free interferent rejection algorithm is implemented. In certain such embodiments, a sample measurement is obtained, substantially as described above in reference to block 1110 ofFIG. 11 . Thealgorithm processor 416, inblock 1120, analyzes the obtained measurement to identify possible interferents. For example, thealgorithm processor 416 may form a statistical sample population model and calculate statistical sample population parameters including mean spectra and covariance matrix (e.g., as described above with reference to block 1210 ofFIG. 12 ). Theprocessor 416 may then assemble a library of interferent data (e.g., as described above with reference to block 1220 ofFIG. 12 ). The library may include interferent spectra, maximum plasma concentration, and a common random concentration distribution function for each interferent. In some embodiments, theprocessor 416 calculates a common variance (denoted by v) of the common random concentration distribution function. - The library may be divided into groups comprising some or all combinations of a number K of the library interferents. The number K may be an integer such as 1, 2, 3, 4, 5, 6, 7, 8, 15, 20, or more. A statistical test may then be performed to determine how well some or all of the groups of K library interferents fits the statistical population model. For example, the statistical test may provide a value for the Mahalanobis distance (of distance squared) for each group and/or an estimate of the concentration of some or all of the library interferents. In some embodiments, groups in which one or more estimated concentrations are negative are eliminated as being unphysical. In other embodiments, some or all groups having negative estimated concentrations may be retained, because they may indicate that the estimated concentration is lower than a standard or reference concentration (e.g., due to dilution of the sample measurement by saline or another fluid). A subset of the remaining groups may be selected, which provide the most likely interferents. For example, the subset may comprise the groups having a number N of the smallest values of the Mahalanobis distance (or distance squared). In various embodiments, the number N may be 1, 2, 5, 10, 20, 100, 200, or more. In certain embodiments, the subset is used to form a model group comprising some or all combinations of a number L of the groups in the subset. For example, the model group may comprise some or all combinations of pairs of subset groups (e.g., L=2). Because each group in the subset comprises K interferents and each model group comprises L subset groups, there are K*L interferents in each model group. For example, in an embodiment in which the each subset group comprises three interferents (K=3), and pairs of subset groups are combined (L=2), then each model group will have 3*2=6 interferents. Because interferents may be repeated when combinations of subset groups are formed, each model group will have between K+1 and K12 distinct interferents. For example, in the preceding example (K=3, L=2), there may be 4, 5, or 6 distinct interferents in any particular model group. The number of model groups may be determined from the well know formula for the number of combinations of the number N of subset groups taken L at a time: CN L=N!/(L!*(N−L)!). For example, if N=100 subset groups are taken two at a time (e.g., pairs), then there will be 4950 model groups.
- The
algorithm processor 416 may then, inblock 1130 ofFIG. 11 , generate a model for predicting the analyte concentration from the obtained sample measurement. For example, in some implementations, for some or all of the model groups, an average group interferent calibration coefficient is calculated, which accounts for the presence of the distinct interferents in any particular model group. The group interferent calibration coefficient may be calculated according to blocks 1310-1360 ofFIG. 13 in some embodiments. In these embodiments, the group interferents are used, rather than a single interferent, in block 1310 to generate synthesized sample population spectra by adding random concentrations of each interferent present in the group to form an Interferent Enhanced Spectral Database (IESD). Inblock 1320, the IESD is partitioned into a calibration set and a test set. Inblock 1330, the calibration set is used to generate a calibration coefficient for estimating the analyte concentration in the presence of the interferents in the group. In block 1340, the calibration coefficient is used to estimate the analyte concentration of the test set, assuming the presence of that interferent group's interferents. Inblock 1350, the error is calculated in the estimated analyte concentration for the test set. Blocks 1320-1350 may be repeated one or more times to obtain group interferent calibration coefficients and errors for different combinations of calibration and test sets. Inblock 1360, an average group interferent calibration coefficient for each group is calculated from the results determined from blocks 1320-1350. - Returning to block 1140 shown in
FIG. 11 , thealgorithm processor 416 may then use an average calibration coefficient to estimate analyte concentration from the obtained sample measurement. For example, in certain embodiments the average calibration coefficient is determined from an average of the group interferent calibration coefficients determined inblock 1360. The average may be a straight average or a weighted average in various embodiments. The analyte concentration is determined by multiplying this average calibration coefficient by the measured spectra. - In other embodiments, the
algorithm processor 416 uses different algorithms inblock 1130 ofFIG. 11 to determine an average calibration coefficient. For example, in some embodiments, every IESD is used as a calibration set, and there is no partition of the IESD into a calibration set and a test set and no error estimate is calculated. Accordingly, in some of these embodiments, thealgorithm processor 416 may not performblocks - In another embodiment, in
block 1130, the average group calibration coefficient may be determined from the following actions. -
- 1. From the group's NIF interferents, form an interferent spectra matrix, IF, having a mean
IF . - 2. Calculate the covariance of the group's IF spectral set:
- 3. Calculate the group's covariance according to K=K0+ρvΦ, where: K0 is the covariance of the original sample population (from block 1120), ρ is a weighting function that depends on the number of interferents in the group (e.g., ρ=NIF/(NIF+1)), and v is a variance of the (scalar) random concentration function.
- 4. Calculate all eigenvectors of K and their corresponding eigenvalues and sort them by decreasing magnitude. Typically, there is one eigenvector (eigenvalue) for each wavelength measured in the sample. The number of wavelengths is denoted by NW.
- 5. Calculate a QR-decomposition of the matrix of sorted eigenvectors, yielding a matrix Q having orthonormal columns and rows.
- 6. For index n ranging from 2 to NW−1, calculate the product P∥ n=Q(:,1:n)Q(:,1:n)T, where Q(:,1:n) refers to a submatrix comprising the first n columns of the full matrix Q. Subtract P∥ n from the Nw·Nw identity matrix I, thereby yielding the orthogonal projection P⊥ n away from the space spanned by Q(:,1:n). The nth calibration vector may be determined from κn=P⊥ nαG/αG TP⊥ nαG, where αG represents the analyte absorption spectrum. The nth error variance Vn may be determined as the projection of the full covariance K onto the subspace spanned by κn as follows: Vn=κn TKκn
- 7. The average group calibration coefficient κ may be selected to be the mth calibration vector κm for the value of m at which the minimum value for the error variance Vm is attained.
- 1. From the group's NIF interferents, form an interferent spectra matrix, IF, having a mean
- A possible advantage of the foregoing algorithms is more rapid execution time by the
algorithm processor 416, because the calibration coefficient is computed directly, without synthesizing spectra or breaking the data into calibration and test sets. In other embodiments, a skilled artisan will recognize that regression, partial least squares, and/or principal component resolution techniques may be used to determine the average group calibration coefficient. - User Interface
- The
system 400 may include adisplay controller 414, for example, as depicted inFIG. 4 . Thedisplay controller 414 may comprise an input device including, for example, a keypad or a keyboard, a mouse, a touchscreen display, and/or any other suitable device for inputting commands and/or information. Thedisplay controller 414 may also include an output device including, for example, an LCD monitor, a CRT monitor, a touchscreen display, a printer, and/or any other suitable device for outputting text, graphics, images, videos, etc. In some embodiments, a touchscreen display is advantageously used for both input and output. - The
display controller 414 may include auser interface 1400 by which users can conveniently and efficiently interact with thesystem 400. Theuser interface 1400 may be displayed on the output device of the system 400 (e.g., the touchscreen display). -
FIGS. 14A and 14B schematically illustrate the visual appearance of embodiments of theuser interface 1400. Theuser interface 1400 may showpatient identification information 1402, which may include patient name and/or a patient ID number. Theuser interface 1400 also may include the current date andtime 1404. An operating graphic 1406 shows the operating status of thesystem 400. For example, as shown inFIGS. 14A and 14B , the operating status is “Running,” which indicates that thesystem 400 is fluidly connected to the patient (“Jill Doe”) and performing normal system functions such as infusing fluid and/or drawing blood. Theuser interface 1400 can include one or moreanalyte concentration graphics FIG. 14A shows “Glucose” concentration of 150 mg/dl, while the graphic 1412 shows “Lactate” concentration of 0.5 mmol/L. The particular analytes displayed and their measurement units (e.g., mg/dl, mmol/L, or other suitable unit) may be selected by the user. The size of thegraphics user interface 1400 may also include a next-reading graphic 1410 that indicates the time until the next analyte measurement is to be taken. InFIG. 14A , the time until next reading is 3 minutes, whereas inFIG. 14B , the time is 6 minutes, 13 seconds. - The
user interface 1400 may include an analyte concentration status graphic 1414 that indicates status of the patient's current analyte concentration compared with a reference standard. For example, the analyte may be glucose, and the reference standard may be a hospital ICU's tight glycemic control (TGC). InFIG. 14A , the status graphic 1414 displays “High Glucose,” because the glucose concentration (150 mg/dl) exceeds the maximum value of the reference standard. InFIG. 14B , the status graphic 1414 displays “Low Glucose,” because the current glucose concentration (79 mg/dl) is below the minimum reference standard. If the analyte concentration is within bounds of the reference standard, the status graphic 1414 may indicate normal (e.g., “Normal Glucose”), or it may not be displayed at all. The status graphic 1414 may have a background color (e.g., red) when the analyte concentration exceeds the acceptable bounds of the reference standard. - The
user interface 1400 may include one ormore trend indicators 1416 that provide a graphic indicating the time history of the concentration of an analyte of interest. InFIGS. 14A and 14B , thetrend indicator 1416 comprises a graph of the glucose concentration (in mg/dl) versus elapsed time (in hours) since the measurements started. The graph includes atrend line 1418 indicating the time-dependent glucose concentration. In other embodiments, thetrend line 1418 may include measurement error bars and may be displayed as a series of individual data points. InFIG. 14B , theglucose trend indicator 1416 is shown as well as atrend indicator 1430 andtrend line 1432 for the lactate concentration. In some embodiments, a user may select whether none, one, or bothtrend indicators trend indicators - The
user interface 1400 may include one or more buttons 1420-1426 that can be actuated by a user to provide additional functionality or to bring up suitable context-sensitive menus and/or screens. For example, in the embodiments shown inFIGS. 14A and 14B , four buttons 1420-1426 are shown, although fewer or more buttons are used in other embodiments. The button 1420 (“End Monitoring”) may be pressed when one or both of thedisposable cassettes 610, 612 (seeFIG. 6 ) are to be removed. In many embodiments, because thecassettes button 1420 is pressed. If the user is certain that monitoring should stop, the user can confirm this by actuating an affirmative button in the confirmation window. If thebutton 1420 were pushed by mistake, the user can select a negative button in the confirmation window. If “End Monitoring” is confirmed, thesystem 400 performs appropriate actions to cease fluid infusion and blood draw and to permit ejection of one (or both)cassettes - The button 1422 (“Pause”) may be actuated by the user if patient monitoring is to be interrupted but is not intended to end. For example, the “Pause”
button 1422 may be actuated if the patient is to be temporarily disconnected from the system 400 (e.g., by disconnecting the tubes 306). After the patient is reconnected, thebutton 1422 may be pressed again to resume monitoring. In some embodiments, after the “Pause”button 1422 has been pressed, thebutton 1422 displays “Resume.” - The button 1424 (“
Delay 5 Minutes”) causes thesystem 400 to delay the next measurement by a delay time period (e.g., 5 minutes in the depicted embodiments). Actuating thedelay button 1424 may be advantageous if taking a reading would be temporarily inconvenient, for example, because a health care professional is attending to other needs of the patient. Thedelay button 1424 may be pressed repeatedly to provide longer delays. In some embodiments, pressing thedelay button 1424 is ineffective if the accumulated delay exceeds a maximum threshold. The next-reading graphic 1410 automatically increases the displayed time until the next reading for every actuation of the delay button 1424 (up to the maximum delay). - The button 1426 (“Dose History”) may be actuated to bring up a dosing history window that displays patient dosing history for an analyte or medicament of interest. For example, in some embodiments, the dosing history window displays insulin dosing history of the patient and/or appropriate hospital dosing protocols. A nurse attending the patient can actuate the
dosing history button 1426 to determine the time when the patient last received an insulin dose, the last dosage amount, and/or the time and amount of the next dosage. Thesystem 400 may receive the patient dosing history via wired or wireless communications from a hospital information system. - In other embodiments, the
user interface 1400 may include additional and/or different buttons, menus, screens, graphics, etc. that are used to implement additional and/or different functionalities. - Related Components
-
FIG. 15 schematically depicts various components and/or aspects of apatient monitoring system 15130 and how those components and/or aspects relate to each other. Some of the depicted components can be included in a kit containing a plurality of components. Some of the depicted components, including, for example, the components represented within the dashedrounded rectangle 15140 ofFIG. 15 , are optional and/or can be sold separately from other components. - The
patient monitoring system 15130 shown inFIG. 15 includes amonitoring device 15132. Themonitoring device 15132 can provide monitoring of physiological parameters of a patient. In some embodiments, themonitoring device 15132 measures glucose and/or lactate concentrations in the patient's blood. In some embodiments, the measurement of such physiological parameters is substantially continuous. Themonitoring device 15132 may also measure other physiological parameters of the patient. In some embodiments, themonitoring device 15132 is used in an intensive care unit (ICU) environment. In some embodiments, onemonitoring device 15132 is allocated to each patient room in an ICU. - The
patient monitoring system 15130 can include anoptional interface cable 15142. In some embodiments, theinterface cable 15142 connects themonitoring device 15132 to a patient monitor (not shown). Theinterface cable 15142 can be used to transfer data from themonitoring device 15132 to the patient monitor for display. In some embodiments, the patient monitor is a bedside cardiac monitor having a display that is located in the patient room. In some embodiments, theinterface cable 15142 transfers data from themonitoring device 15132 to a central station monitor and/or to a hospital information system (HIS). The ability to transfer data to a central station monitor and/or to a HIS may depend on the capabilities of the patient monitor system. - In the embodiment shown in
FIG. 15 , an optionalbar code scanner 15144 is connected to themonitoring device 15132. In some embodiments, thebar code scanner 15144 is used to enter patient identification codes, nurse identification codes, and/or other identifiers into themonitoring device 15132. In some embodiments, thebar code scanner 15144 contains no moving parts. Thebar code scanner 15144 can be operated by manually sweeping thescanner 15144 across a printed bar code or by any other suitable means. In some embodiments, thebar code scanner 15144 includes an elongated housing in the shape of a wand. - The
patient monitoring system 15130 includes afluidic system kit 15134 connected to themonitoring device 15132. In some embodiments, thefluidic system kit 15134 includes fluidic tubes that connect a fluid source to an analytic subsystem. For example, the fluidic tubes can facilitate fluid communication between a blood source or a saline source and an assembly including a flow cell and/or a centrifuge. In some embodiments, thefluidic system kit 15134 includes many of the components that enable operation of themonitoring device 15132. In some embodiments, thefluidic system kit 15134 can be used with anti-clotting agents (such as heparin), saline, a saline infusion set, a patient catheter, a port sharing IV infusion pump, and/or an infusion set for an IV infusion pump, any or all of which may be made by a variety of manufacturers. In some embodiments, thefluidic system kit 15134 includes a monolithic housing that is sterile and disposable. In some embodiments, at least a portion of thefluidic system kit 15134 is designed for single patient use. For example, thefluidic system kit 15134 can be constructed such that it can be economically discarded and replaced with a newfluidic system kit 15134 for every new patient to which thepatient monitoring system 15130 is connected. In addition, at least a portion of thefluidic system kit 15134 can be designed to be discarded after a certain period of use, such as a day, several days, several hours, three days, a combination of hours and days such as, for example, three days and two hours, or some other period of time. Limiting the period of use of thefluidic system kit 15134 may decrease the risk of malfunction, infection, or other conditions that can result from use of a medical apparatus for an extended period of time. - In some embodiments, the
fluidic system kit 15134 includes a connector with a luer fitting for connection to a saline source. The connector may be, for example, a three-inch pigtail connector. In some embodiments, thefluidic system kit 15134 can be used with a variety of spikes and/or IV sets used to connect to a saline bag. In some embodiments, thefluidic system kit 15134 also includes a three-inch pigtail connector with a luer fitting for connection to one or more IV pumps. In some embodiments, thefluidic system kit 15134 can be used with one or more IV sets made by a variety of manufacturers, including IV sets obtained by a user of thefluidic system kit 15134 for use with an infusion pump. In some embodiments, thefluidic system kit 15134 includes a tube with a low dead volume luer connector for attachment to a patient vascular access point. For example, the tube can be approximately seven feet in length and can be configured to connect to a proximal port of a cardiovascular catheter. In some embodiments, thefluidic system kit 15134 can be used with a variety of cardiovascular catheters, which can be supplied, for example, by a user of thefluidic system kit 15134. - As shown in
FIG. 15 , themonitoring device 15132 is connected to asupport apparatus 15136, such as an IV pole. Thesupport apparatus 15136 can be customized for use with themonitoring device 15132. A vendor of themonitoring device 15132 may choose to bundle themonitoring device 15132 with acustom support apparatus 15136. In one embodiment, thesupport apparatus 15136 includes a mounting platform for themonitoring device 15132. The mounting platform can include mounts that are adapted to engage threaded inserts in themonitoring device 15132. Thesupport apparatus 15136 can also include one or more cylindrical sections having a diameter of a standard IV pole, for example, so that other medical devices, such as IV pumps, can be mounted to the support apparatus. Thesupport apparatus 15136 can also include a clamp adapted to secure the apparatus to a hospital bed, an ICU bed, or another variety of patient conveyance device. - In the embodiment shown in
FIG. 15 , themonitoring device 15132 is electrically connected to anoptional computer system 15146. Thecomputer system 15146 can be used to communicate with one or more monitoring devices. In an ICU environment, thecomputer system 15146 can be connected to at least some of the monitoring devices in the ICU. Thecomputer system 15146 can be used to control configurations and settings for multiple monitoring devices (for example, the system can be used to keep configurations and settings of a group of monitoring devices common). Thecomputer system 15146 can also run optional software, such asdata analysis software 15148, HISinterface software 15150, andinsulin dosing software 15152. - In some embodiments, the
computer system 15146 runs optionaldata analysis software 15148 that organizes and presents information obtained from one or more monitoring devices. In some embodiments, thedata analysis software 15148 collects and analyzes data from the monitoring devices in an ICU. Thedata analysis software 15148 can also present charts, graphs, and statistics to a user of thecomputer system 15146. - In some embodiments, the
computer system 15146 runs optional hospital information system (HIS)interface software 15150 that provides an interface point between one or more monitoring devices and an HIS. The HISinterface software 15150 may also be capable of communicating data between one or more monitoring devices and a laboratory information system (LIS). - In some embodiments, the
computer system 15146 runs optionalinsulin dosing software 15152 that provides a platform for implementation of an insulin dosing regimen. In some embodiments, the hospital tight glycemic control protocol is included in the software. The protocol allows computation of proper insulin doses for a patient connected to amonitoring device 15146. Theinsulin dosing software 15152 can communicate with themonitoring device 15146 to ensure that proper insulin doses are calculated. - Noise Reduction
-
FIG. 16 is a block diagram of an embodiment of asystem 1600 for reducing noise and/or unwanted elements in a signal. Thesystem 1600 includes adetector 1602, such as, for example, asample detector 930 or areference detector 936. One type of detector that can serve as asample detector 930 or areference detector 936 is a pyroelectric infrared detector. It is typical for pyroelectric detectors to also be sensitive to vibrations. Thus, for example, the output of a pyroelectric infrared detector is the sum of the exposure to infrared radiation and to vibrations of the detector. The sensitivity to vibrations, also known as “microphonics,” can introduce a large noise component to the measurement of radiation using a pyroelectric infrared detector. It is desirable for a spectrometer, such as theanalyte detection system 910, to have a high signal-to-noise ratio, such as a S/N in excess of 100 dB. It can be difficult to achieve this low noise level in the presence of vibrations, even with good mechanical isolation. The presently disclosed system includes one or more of the following techniques for reducing the vibrational noise component of measurements using detectors (e.g., thesample detector 930 or the reference detector 936). - Some embodiments for reducing vibrational noise include the use of a modulated infrared source combined with an output filter. In some embodiments, the
analyzer 910 is the infrared source. The infrared source of theanalyzer 910 can be modulated at a known frequency, and the detector output can be filtered using a narrow band filter centered about the known source frequency. Thus, in some embodiments, theenergy source 912 has an energy output that is sine-wave modulated at 10 Hz, and the output of the detector(s) (e.g., thedetectors 930, 936) is filtered using a narrow-band pass filter. The narrow-band pass filter can have a frequency of less than 1 Hz and be centered about 10 Hz, for example. Microphonic (also referred to as vibration-induced) signals that are not exactly 10 Hz can be significantly attenuated with this arrangement. - In some embodiments, the detector output is filtered using
signal conditioning 1604, including a synchronous demodulator and digital filter. The demodulator can be a software component implemented in the signal processing computer. Synchronous demodulators, coupled with low pass filters, can be referred to as “lock-in amplifiers.” -
FIG. 17 illustrates an example ofsignal conditioning 1604 including a lock-inamplifier system 1710. Aninput device 1712 such as, for example, theanalyzer 910 ofFIG. 9 or thedetector 1602 ofFIG. 16 , is shown at the left. Theinput device 1712 is part of a circuit 1713 that provides one or more inputs into the lock-inamplifier system 1710. The circuit may contain one or more resistors R1, R2 and may also have a connection to ground. The circuit 1713 provides one or more inputs into adifferential AC amplifier 1714. The embodiment of thedifferential AC amplifier 1714 shown inFIG. 17 accepts inputs from ground and from theinput device 1712. Thedifferential AC amplifier 1714 outputs an amplified signal waveform having a substantial frequency component fs substantially equal to the frequency of the infrared source of theanalyzer 910. The signal waveform may also have other frequency components that correspond to, for example, noise. The output of thedifferential AC amplifier 1714 is supplied to afunction device 1716. The circuit 1713 also provides an input into a phase-lock loop (PLL) 1718. ThePLL 1718 conditions the signal and sends a reference waveform having a frequency component fr to thefunction device 1716. The frequency fr of the predominant component of the reference waveform is substantially equal to the frequency of the infrared source of theanalyzer 910. - In the embodiment shown in
FIG. 17 , thefunction device 1716 outputs the product of the inputs of the device. Thefunction device 1716 multiplies the reference waveform and the signal waveform. The product of the reference waveform and the signal waveform has substantial components at frequencies fs−fr and fs+fr that are proportional to the amplitude of the signal at the detector output. Because fs and fr are substantially equal, the resulting waveform will have a 0 Hz component (e.g., a DC component or a demodulated component) and a component at about 2fr that are proportional to the amplitude of the signal at the detector output. The resulting waveform is fed through alow pass filter 1720, which can eliminate substantially all components of the signal above a threshold frequency. For example, if the pass-band of thefilter 1720 is small enough (e.g., if the threshold frequency is low enough), thefilter 1720 can substantially remove the portions of the waveform that are above the demodulated component at 0 Hz. This resulting waveform can then pass to aDC amplifier 1722, which provides an output from the lock-inamplifier system 1710. - In some embodiments, the center frequency of a lock-in
amplifier system 1710 is 10 Hz and the integration time is 2 seconds, yielding an approximately 0.5 Hz wide pass band. The use of such a device excludes a significant amount of vibration induced noise outside of the 9.75 to 10.25 Hz frequency band. In some embodiments, the energy present in the pass band generates the signal of theinput device 1712. In some embodiments, the optical source can be a radiation source that is modulated at 10 Hz. Thus, the optical source can have all of its energy centered in the pass band, and it can proceed through the demodulator and filter unattenuated. Thus, a large part of the vibration-induced noise is eliminated because most of it falls outside of the pass band. -
FIG. 18 graphically depicts example signal levels, in volts, of the 10 Hz signal of an input device 1712 (e.g., from theanalyzer 910 ofFIG. 9 ). The illustrated signal can originate from a modulated infrared radiation source that passes through a medium and is then detected by an optical detector. Thus, the illustrated signal can be a detector signal before it reaches thedifferential AC amplifier 1714 ofFIG. 17 . Signals such as these can be produced, for example, as ananalyzer 910 cycles throughvarious filters 915 in afilter wheel 918. The example shown here depicts 13 different signal regions corresponding to 13 different filters. Shown are relatively large AC signals centered approximately around zero volts. -
FIG. 19 graphically depicts a signal that can be produced by thesame analyzer 910 ofFIG. 9 when the original modulated infrared radiation source is blocked. Thus, likeFIG. 18 ,FIG. 19 also illustrates the signal output from a detector before that signal reaches thedifferential AC amplifier 1714 ofFIG. 17 . However, this time the illustrated signal (output from a detector) is due only to vibrations of the detector. The top portion ofFIG. 19 shows how, when a system (e.g., theapparatus 100 illustrated inFIG. 1 ) is not moving (stationary), the output signal is steady. In contrast, the bottom portion ofFIG. 19 shows how, when the system is placed on a rolling platform and is rolled across the floor, there is some signal due to the microphonics of the detector. As illustrated, the magnitude of the microphonic signal is generally much less than the magnitude of the optical signal depicted inFIG. 18 . Nevertheless, even this small amount of microphonic signal may produce an output that has an unacceptably low signal-to-noise ratio. -
FIG. 20 graphically depicts signals after they have passed through a lock-in amplifier such as the one depicted inFIG. 17 (e.g., a demodulator and filter), resulting in a DC output. As with the signals depicted inFIG. 19 , the signals shown inFIG. 20 can originate from theanalyzer 910 ofFIG. 9 when the original modulated infrared radiation source is blocked. Thus, the optical signal is removed, leaving only noise or interference signals. As the top signal shows, the noise level is very low when theapparatus 100 is stationary. In contrast, the bottom signal shows a higher noise level resulting from microphonic interference when theapparatus 100 is non-stationary (e.g., rolling across a hospital floor). The two signals shown inFIG. 20 are derived from the two signals ofFIG. 19 after passing those signals through a lock-in amplifier such as the one depicted inFIG. 17 (e.g., a demodulator and filter). The pre-filtering “rolling” microphonic signals ofFIG. 19 have a peak-to-peak amplitude on the order of 5+ mV; as shown inFIG. 20 , however, after demodulation and filtering, the demodulated voltage is typically 0.5 mV. Thus, a lock-in amplifier system (or, for example, a sync demodulator and a digital filter) can result in a 10:1 reduction in microphonics. - As shown in
FIG. 18 , the optical signal level (S) can be approximately 3 V peak-to-peak. As shown in the upper portion ofFIG. 20 , the typical microphonics free (no vibration) noise level (N) is approximately 0.03 mV after the signal passes through the lock-inamplifier system 1710. Thus, when theapparatus 100 is stationary, the signal-to-noise ratio (S/N) is relatively low. For example, it meets the 100 dB S/N test, which can be a requirement in some embodiments. However, when theapparatus 100 is not stationary (e.g., when it rolls across a floor), the lower portion ofFIG. 20 shows that even with signal filtering, microphonics can degrade the signal, resulting in a noise level (N) of approximately 0.5 mV. Thus, when theapparatus 100 is not stationary, the signal-to-noise ratio (S/N) is higher. Accordingly, rolling theapparatus 100 can drop the S/N level by approximately 25 dB to 75 dB. A S/N level of 75 dB does not meet the 100 dB S/N test, and thus, in some embodiments, a further improvement in S/N level may be required. - Although the lock-in
amplifier 1710 can condition the signal and remove some noise, as described above with respect toFIG. 17 , some microphonic effects can still persist in the signal even after the signal has passed through the lock-inamplifier 1710. These effects can increase the error in the measurements provided by theapparatus 100. - In some cases, microphonic effects can create such a noisy signal that the signal passes beyond a threshold of acceptability. For example, in some embodiments of a glucose monitoring system, rolling the
apparatus 100 ofFIG. 1 over a rough floor can result in an unacceptable error. One example of an unacceptable standard error (1σ) is an error of over 20 mg/dL of glucose. An example of a threshold level, beyond which errors are unacceptable in some embodiments, is 5 mg/dL of glucose. Other threshold error levels can also be selected. In some embodiments, when the error exceeds a predetermined level, additional steps can be advantageously taken to prevent theapparatus 100 from reporting an incorrect value. - In some embodiments, a “hold off” can be included with or in the hardware or software of an
apparatus 100. The “hold off” can prevent the system from taking and/or recording a measurement (e.g., a spectroscopic measurement) when a noise threshold is exceeded. In some embodiments, a “hold off” can be activated when a non-noise threshold is exceeded. For example, the threshold can be a value having the same units as the final output of the system (e.g., mg/dL). Thus, the system can be designed to track when a certain guaranteed accuracy level is or is not being met, and when it is not being met, the system can automatically stop providing data regarding glucose levels or concentrations, for example. - The “hold off” can comprise sensing accelerations of the
apparatus 100 and preventing theapparatus 100 from measuring during periods of excessive vibrations. In some embodiments, vibrations are measured using one or more accelerometers that can be located, for example, near pyroelectric detectors (e.g., thedetectors 930 and 936). When the accelerometers sense a vibration that exceeds a pre-determined value, theanalyzer 910 is instructed to stop making spectroscopic measurements until the vibration has subsided to below the threshold level. The pre-determined threshold value of acceleration can be selected so that the error in analyte measurement is at an acceptable level for accelerations below the pre-determined threshold value. -
FIG. 21 shows a block diagram of asystem 2110 for dealing with noise in an analyte detection environment. The illustratedsystem 2110 can be used with theapparatus 100 ofFIG. 1 , for example. Adebugging module 2114 can comprise an RJ-45 chip. Thedebugging module 2114 can be connected to analgorithm engine 2118, which can comprise a computer chip, a digital signal processor, and/or a field-programmable gate array, for example. Thealgorithm engine 2118 can function as thealgorithm processor 416 described above. The algorithm engine can communicate (via acustom interface 2120, for example) with aspectrometer control device 2122 that can preferably, among other things, convert analog signals to digital format. -
FIG. 21 also shows how several other components of thesystem 2110 can connect to thespectrometer control device 2122. For example, atemperature sensor 2124, aheater 2128, and afilter wheel 2130 can all connect to (and, in some embodiments, be controlled by) thespectrometer control device 2122. Moreover, in some embodiments, vibration sensor(s) 2140 (which can be accelerometers, for example) and detector(s) 2150 (which can be theinput device 1712, for example) can both feed signals to thespectrometer control device 2122. Preferably, these signals are fed in with little or no delay (e.g., in “real time”). Thespectrometer control device 2122 can monitor and quantify the output from thevibration sensors 2140, which can allow a user or system to detect the vibration. This can, in turn, indicate that a significant portion of the detector signal is due to vibration, which can trigger a “hold off” procedure such as that described above. After a “hold off” has occurred, thespectrometer control device 2122 preferably causes a re-measurement or a resumption of measurement. - With reference to
FIG. 22 , when monitoring accelerometers, it can be useful to gather information related not only to a vibration above a predetermined threshold value but also to a signal within the pass band of the detector signal, the 10 Hz signal produced by the detector. This can be done as vibration sensors (e.g., the vibration sensor(s) 2140 and/or the accelerometer 2212) are monitored and quantified as described above. To gather such information, anaccelerometer 2212 can feed a signal through aconnection 2214 into a “lock in amplifier”style demodulator 2216 similar to (or exactly like) the lock-inamplifier 1710 used for the detector signal (seeFIG. 17 ).FIG. 22 shows schematically how such a system can be set up. Theinput device 1712 can be a detector. Because theaccelerometer 2212 has no frequency output of its own, (and thus is not modulated at 10 Hz), the 10 Hz demodulation signal from the detector channel's PLL can be used.FIG. 22 illustrates how the PLL output 2218 can be taken from the detector channel and fed directly into theaccelerometer demodulation multiplier 2220. This same result can be accomplished by software, rather than with the hardware configuration illustrated. - In some embodiments, a system is designed to detect when the vibration-induced error is so large that the instrument (e.g., the apparatus 100) will not be able to achieve a given accuracy specification. Including one or multiple accelerometers can help achieve this goal. In some embodiments, an accelerometer measures forces in three axes. To convert that measurement into glucose error a linear regression calibration can be used to estimate an error (e.g., Optical Density error, or ODe) based on the three accelerometer signals for the three axes.
- Some embodiments of a noise removal method include performing a regression analysis to determine when vibration-induced errors in parameter values exceed acceptable limits. For example, calibration computations based on vibration-induced glucose concentration errors can allow a monitoring device to predict when vibrations measured by the one or more accelerometers will produce errors in glucose concentration that exceed acceptable limits. In some embodiments, measurements taken during an initial period, such as the first 30 seconds of apparatus operation, can be used for calibration and the remaining measurements can be used to predict when calculated parameters have an unacceptable amount of error.
- Vibration errors can be more harmful during some measurements than others. For example, when a filter wheel 2230 is allowing some wavelengths of radiation to propagate, vibration-induced noise can be especially difficult to filter and/or detect and correct. In some embodiments, maximum error in the optical density (OD) measurement can coincide with the measurement of the 9.22 micron filter. This can occur, for example, when that filter has the highest coefficient in terms of mg/dL per ODe when compared to the other filters.
- In some embodiments, it is desired that spectroscopic measurements of a sample be completed within a certain period of time (the “maximum measurement period”). In some embodiments, a maximum measurement period can be selected based on the stability of the sample. Thus, for example, one system can produce accurate glucose measurements if a spectroscopic measurement can be completed within 90 seconds. The “maximum measurement period” in this case is 90 seconds. In some embodiments, the cumulative hold-off time for making a complete measurement (for example by scanning all of the filters) is specifically designed to be less than the maximum measurement period. Thus, if vibrations prevent accurate measurement from being made for some short period of time, measurement can be allowed to continue if that short period of time still permits scanning through all the filters (or a pre-determined sub-set of filters) within the maximum measurement period. In this case, while measurement may be periodically interrupted during times of vibration (e.g., when the accelerometers indicate acceleration beyond a threshold value), there may still be enough time to complete a full measurement cycle within the maximum measurement period of 90 seconds. On the other hand, if the vibrations occur over a long time (either individually or in the aggregate), there may not be enough time remaining within the maximum measurement period to complete the measurement. In this case, the measurement may be aborted and restarted later.
- In some embodiments, a spectroscopic scan of
analyzer 910 takes approximately 2-3 seconds per filter (see, e.g., filters 915 ofFIG. 9 ), and there can be 25-30 filters in a system, in some embodiments. Some embodiments require that all filters be scanned without significant change in the sample to meet accuracy guidelines for analyte (e.g., glucose) computation. Typically, without any vibration interference, all filters are scanned in 50 to 90 seconds. In some embodiments, the sample meets certain stability requirements during this time to achieve a glucose measurement with error below a predetermined accuracy, which may be, for example, 5 mg/dL. If the acceleration is unacceptably high for a short period of time, for example less than 5 seconds, then the measurement using a particular filter can be temporarily stopped. If the acceleration is unacceptably high and a complete filter scan cannot be made within the period for a stable sample, for example 90 seconds, then the entire scan can be restarted and completed within the stability period of the sample. - In some embodiments, any single delay of up to approximately 5 seconds can be tolerated for a single filter before having to restart a measurement cycle. In the case of a brief vibration incident while scanning a filter (1-2 seconds) the measurement can “hold off” on that single filter measurement, subsequently completing the measurement and moving on to the next when the vibration stops. Holding off on only one or two filters in this way can add only a few seconds to the measurement time and is likely not unacceptable to a user.
- If the vibration lasts so long that holding off and re-scanning a single filter after the vibration stops would cause the entire filter scan interval to exceed approximately 90 seconds, the system can wait (in some embodiments, the system is required to wait) until the vibration has stopped and then restart the entire 25-30 filter scans. Thus, another 50-90 seconds may be required before the analyte measurement can be completed.
- In addition to convenience to a user, biological factors can also be pertinent to the timing settings of the system. For example, clotting, aggregation, or other biological processes can cause blood to plug the flow cell if allowed to remain stagnant in the flow cell for too long. In some embodiments, 10 minutes is too long. Therefore if a vibration period strong enough to trigger measurement “hold off” lasts so long that the 10 min. blood holding period will be exceeded, the entire measurement process (including drawing new blood for measurement, in some embodiments) is preferably re-started. In some embodiments, the minimum processing time of the instrument is approximately 10 minutes, so in such a case the glucose reading will be refreshed 10 minutes after the vibration ceases.
- Some embodiments of each of the methods described herein may include a computer program accessible to and/or executable by a processing system, e.g., a one or more processors and memories that are part of an embedded system. Thus, as will be appreciated by those skilled in the art, embodiments of the disclosed inventions may be embodied as a method, an apparatus such as a special purpose apparatus, an apparatus such as a data processing system, or a carrier medium, e.g., a computer program product. The carrier medium carries one or more computer readable code segments for controlling a processing system to implement a method. Accordingly, various ones of the disclosed inventions may take the form of a method, an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, any one or more of the disclosed methods (including but not limited to the disclosed methods of measurement analysis, interferent determination, and/or calibration constant generation) may be stored as one or more computer readable code segments or data compilations on a carrier medium. Any suitable computer readable carrier medium may be used including a magnetic storage device such as a diskette or a hard disk; a memory cartridge, module, card or chip (either alone or installed within a larger device); or an optical storage device such as a CD or DVD.
- Reference throughout this specification to “some embodiments” or “an embodiment” means that a particular feature, structure or characteristic described in connection with the embodiment is included in at least some embodiments. Thus, appearances of the phrases “in some embodiments” or “in an embodiment” in various places throughout this specification are not necessarily all referring to the same embodiment. Furthermore, the particular features, structures or characteristics may be combined in any suitable manner, as would be apparent to one of ordinary skill in the art from this disclosure, in one or more embodiments.
- Similarly, it should be appreciated that in the above description of embodiments, various features of the inventions are sometimes grouped together in a single embodiment, figure, or description thereof for the purpose of streamlining the disclosure and aiding in the understanding of one or more of the various inventive aspects. This method of disclosure, however, is not to be interpreted as reflecting an intention that any claim require more features than are expressly recited in that claim. Rather, inventive aspects lie in a combination of fewer than all features of any single foregoing disclosed embodiment.
- Further information on analyte detection systems, sample elements, algorithms and methods for computing analyte concentrations, and other related apparatus and methods can be found in U.S. Patent Application Publication No. 2003/0090649, published May 15, 2003, titled REAGENT-LESS WHOLE BLOOD GLUCOSE METER; U.S. Patent Application Publication No. 2003/0178569, published Sep. 25, 2003, titled PATHLENGTH-INDEPENDENT METHODS FOR OPTICALLY DETERMINING MATERIAL COMPOSITION; U.S. Patent Application Publication No. 2004/0019431, published Jan. 29, 2004, titled METHOD OF DETERMINING AN ANALYTE CONCENTRATION IN A SAMPLE FROM AN ABSORPTION SPECTRUM; U.S. Patent Application Publication No. 2005/0036147, published Feb. 17, 2005, titled METHOD OF DETERMINING ANALYTE CONCENTRATION IN A SAMPLE USING INFRARED TRANSMISSION DATA; and U.S. Patent Application Publication No. 2005/0038357, published on Feb. 17, 2005, titled SAMPLE ELEMENT WITH BARRIER MATERIAL. The entire contents of each of the above-mentioned publications are hereby incorporated by reference herein and are made a part of this specification.
- A number of applications, publications and external documents are incorporated by reference herein. Any conflict or contradiction between a statement in the bodily text of this specification and a statement in any of the incorporated documents is to be resolved in favor of the statement in the bodily text.
- Although the invention(s) presented herein have been disclosed in the context of certain preferred embodiments and examples, it will be understood by those skilled in the art that the invention(s) extend beyond the specifically disclosed embodiments to other alternative embodiments and/or uses of the invention(s) and obvious modifications and equivalents thereof. Thus, it is intended that the scope of the invention(s) herein disclosed should not be limited by the particular embodiments described above.
Claims (28)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US11/734,261 US20070258083A1 (en) | 2006-04-11 | 2007-04-11 | Noise reduction for analyte detection systems |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US79193106P | 2006-04-11 | 2006-04-11 | |
US11/734,261 US20070258083A1 (en) | 2006-04-11 | 2007-04-11 | Noise reduction for analyte detection systems |
Publications (1)
Publication Number | Publication Date |
---|---|
US20070258083A1 true US20070258083A1 (en) | 2007-11-08 |
Family
ID=38610168
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/734,261 Abandoned US20070258083A1 (en) | 2006-04-11 | 2007-04-11 | Noise reduction for analyte detection systems |
Country Status (2)
Country | Link |
---|---|
US (1) | US20070258083A1 (en) |
WO (1) | WO2007120745A2 (en) |
Cited By (44)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060200071A1 (en) * | 2005-02-14 | 2006-09-07 | Sterling Bernhard B | Method and apparatus for detection of multiple analytes |
US20070273881A1 (en) * | 2006-05-23 | 2007-11-29 | Str Systems, Inc. D/B/A Safety Technology And Rese | Biofiedlic displacement measuring system for an anthropomorphic testing device |
US20090043171A1 (en) * | 2007-07-16 | 2009-02-12 | Peter Rule | Systems And Methods For Determining Physiological Parameters Using Measured Analyte Values |
US20100064257A1 (en) * | 2008-09-05 | 2010-03-11 | Roche Diagnostics Operations, Inc. | Insulin pump configuration programming invalid settings notification and correction |
WO2011011462A1 (en) * | 2009-07-20 | 2011-01-27 | Optiscan Biomedical Corporation | Adjustable connector and dead space reduction |
US20110257593A1 (en) * | 2007-03-30 | 2011-10-20 | Medtronic, Inc. | Devices and methods for detecting catheter complications |
WO2013010142A1 (en) * | 2011-07-13 | 2013-01-17 | Bio-Rad Laboratories, Inc., A Delawae Corporation | Computation of real-world error using meta-analysis of replicates |
US20130053099A1 (en) * | 2011-08-25 | 2013-02-28 | Kyocera Corporation | Human body detecting system |
US8425444B2 (en) | 2006-04-11 | 2013-04-23 | Optiscan Biomedical Corporation | Anti-clotting apparatus and methods for fluid handling system |
EP2604998A1 (en) * | 2011-12-12 | 2013-06-19 | ABB Research Ltd. | Gas analysis device |
US8633015B2 (en) | 2008-09-23 | 2014-01-21 | Bio-Rad Laboratories, Inc. | Flow-based thermocycling system with thermoelectric cooler |
US8663920B2 (en) | 2011-07-29 | 2014-03-04 | Bio-Rad Laboratories, Inc. | Library characterization by digital assay |
US8709762B2 (en) | 2010-03-02 | 2014-04-29 | Bio-Rad Laboratories, Inc. | System for hot-start amplification via a multiple emulsion |
US8731639B2 (en) | 2009-07-20 | 2014-05-20 | Optiscan Biomedical Corporation | Adjustable connector, improved fluid flow and reduced clotting risk |
US8730479B2 (en) | 2010-03-25 | 2014-05-20 | Bio-Rad Laboratories, Inc. | Detection system for droplet-based assays |
US8951939B2 (en) | 2011-07-12 | 2015-02-10 | Bio-Rad Laboratories, Inc. | Digital assays with multiplexed detection of two or more targets in the same optical channel |
US9033920B2 (en) | 2003-10-02 | 2015-05-19 | Medtronic, Inc. | Determining catheter status |
US9089844B2 (en) | 2010-11-01 | 2015-07-28 | Bio-Rad Laboratories, Inc. | System for forming emulsions |
US9126160B2 (en) | 2008-09-23 | 2015-09-08 | Bio-Rad Laboratories, Inc. | System for forming an array of emulsions |
US9132394B2 (en) | 2008-09-23 | 2015-09-15 | Bio-Rad Laboratories, Inc. | System for detection of spaced droplets |
US9138537B2 (en) | 2003-10-02 | 2015-09-22 | Medtronic, Inc. | Determining catheter status |
US9194861B2 (en) | 2009-09-02 | 2015-11-24 | Bio-Rad Laboratories, Inc. | Method of mixing fluids by coalescence of multiple emulsions |
US9222128B2 (en) | 2011-03-18 | 2015-12-29 | Bio-Rad Laboratories, Inc. | Multiplexed digital assays with combinatorial use of signals |
US20150378511A1 (en) * | 2014-06-26 | 2015-12-31 | Sitronix Technology Corp. | Capacitive Voltage Information Sensing Circuit and Related Anti-Noise Touch Circuit |
US9347059B2 (en) | 2011-04-25 | 2016-05-24 | Bio-Rad Laboratories, Inc. | Methods and compositions for nucleic acid analysis |
US9393560B2 (en) | 2010-03-25 | 2016-07-19 | Bio-Rad Laboratories, Inc. | Droplet transport system for detection |
US9399215B2 (en) | 2012-04-13 | 2016-07-26 | Bio-Rad Laboratories, Inc. | Sample holder with a well having a wicking promoter |
US9417190B2 (en) | 2008-09-23 | 2016-08-16 | Bio-Rad Laboratories, Inc. | Calibrations and controls for droplet-based assays |
US9492797B2 (en) | 2008-09-23 | 2016-11-15 | Bio-Rad Laboratories, Inc. | System for detection of spaced droplets |
US9500664B2 (en) | 2010-03-25 | 2016-11-22 | Bio-Rad Laboratories, Inc. | Droplet generation for droplet-based assays |
US9598725B2 (en) | 2010-03-02 | 2017-03-21 | Bio-Rad Laboratories, Inc. | Emulsion chemistry for encapsulated droplets |
US20170188943A1 (en) * | 2007-10-10 | 2017-07-06 | Optiscan Biomedical Corporation | Fluid component analysis systems and methods for glucose monitoring and control |
US9764322B2 (en) | 2008-09-23 | 2017-09-19 | Bio-Rad Laboratories, Inc. | System for generating droplets with pressure monitoring |
US9872956B2 (en) | 2011-04-29 | 2018-01-23 | Medtronic, Inc. | Limiting pressure in an implanted catheter |
US20180217055A1 (en) * | 2011-12-05 | 2018-08-02 | Gassecure As | Method and system for gas detection |
US20180238794A1 (en) * | 2015-09-01 | 2018-08-23 | Apple Inc. | Reference switch architectures for noncontact sensing of substances |
CN109682958A (en) * | 2018-09-21 | 2019-04-26 | 深圳沃德生命科技有限公司 | A kind of acceleration transducer signals compensation method for thrombelastogram instrument |
US10512910B2 (en) | 2008-09-23 | 2019-12-24 | Bio-Rad Laboratories, Inc. | Droplet-based analysis method |
US11130128B2 (en) | 2008-09-23 | 2021-09-28 | Bio-Rad Laboratories, Inc. | Detection method for a target nucleic acid |
US11579080B2 (en) | 2017-09-29 | 2023-02-14 | Apple Inc. | Resolve path optical sampling architectures |
US11726036B2 (en) | 2014-12-23 | 2023-08-15 | Apple Inc. | Optical inspection system and method including accounting for variations of optical path length within a sample |
US11852318B2 (en) | 2020-09-09 | 2023-12-26 | Apple Inc. | Optical system for noise mitigation |
US11960131B2 (en) | 2018-02-13 | 2024-04-16 | Apple Inc. | Integrated photonics device having integrated edge outcouplers |
US12007275B2 (en) | 2016-04-21 | 2024-06-11 | Apple Inc. | Optical system for reference switching |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080112853A1 (en) * | 2006-08-15 | 2008-05-15 | Hall W Dale | Method and apparatus for analyte measurements in the presence of interferents |
US10475529B2 (en) | 2011-07-19 | 2019-11-12 | Optiscan Biomedical Corporation | Method and apparatus for analyte measurements using calibration sets |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6091504A (en) * | 1998-05-21 | 2000-07-18 | Square One Technology, Inc. | Method and apparatus for measuring gas concentration using a semiconductor laser |
US20030040683A1 (en) * | 2001-07-06 | 2003-02-27 | Peter Rule | Site selection for determining analyte concentration in living tissue |
US20030086073A1 (en) * | 2001-11-08 | 2003-05-08 | Braig James R. | Reagent-less whole-blood glucose meter |
US20030108976A1 (en) * | 2001-10-09 | 2003-06-12 | Braig James R. | Method and apparatus for improving clinical accuracy of analyte measurements |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS61178662A (en) * | 1985-02-04 | 1986-08-11 | Kanegafuchi Chem Ind Co Ltd | Method and sensor for measuring concentration of low molecular weight sugars |
JPS63275944A (en) * | 1987-05-07 | 1988-11-14 | Nec Corp | Blood sugar meter |
US5118404A (en) * | 1989-04-28 | 1992-06-02 | Nec Corporation | Enzyme electrode and a method of determining concentration of an analyte in a sample solution |
US5398681A (en) * | 1992-12-10 | 1995-03-21 | Sunshine Medical Instruments, Inc. | Pocket-type instrument for non-invasive measurement of blood glucose concentration |
US5615672A (en) * | 1993-01-28 | 1997-04-01 | Optiscan, Inc. | Self-emission noninvasive infrared spectrophotometer with body temperature compensation |
US6947781B2 (en) * | 2002-12-13 | 2005-09-20 | Massachusetts Institute Of Technology | Vibratory venous and arterial oximetry sensor |
-
2007
- 2007-04-11 WO PCT/US2007/008997 patent/WO2007120745A2/en active Application Filing
- 2007-04-11 US US11/734,261 patent/US20070258083A1/en not_active Abandoned
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6091504A (en) * | 1998-05-21 | 2000-07-18 | Square One Technology, Inc. | Method and apparatus for measuring gas concentration using a semiconductor laser |
US20030040683A1 (en) * | 2001-07-06 | 2003-02-27 | Peter Rule | Site selection for determining analyte concentration in living tissue |
US20030108976A1 (en) * | 2001-10-09 | 2003-06-12 | Braig James R. | Method and apparatus for improving clinical accuracy of analyte measurements |
US20030086073A1 (en) * | 2001-11-08 | 2003-05-08 | Braig James R. | Reagent-less whole-blood glucose meter |
Cited By (85)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10357620B2 (en) | 2003-10-02 | 2019-07-23 | Medtronic, Inc. | Determining catheter status |
US9033920B2 (en) | 2003-10-02 | 2015-05-19 | Medtronic, Inc. | Determining catheter status |
US9138537B2 (en) | 2003-10-02 | 2015-09-22 | Medtronic, Inc. | Determining catheter status |
US7722537B2 (en) | 2005-02-14 | 2010-05-25 | Optiscan Biomedical Corp. | Method and apparatus for detection of multiple analytes |
US20060200071A1 (en) * | 2005-02-14 | 2006-09-07 | Sterling Bernhard B | Method and apparatus for detection of multiple analytes |
US8425444B2 (en) | 2006-04-11 | 2013-04-23 | Optiscan Biomedical Corporation | Anti-clotting apparatus and methods for fluid handling system |
US20070273881A1 (en) * | 2006-05-23 | 2007-11-29 | Str Systems, Inc. D/B/A Safety Technology And Rese | Biofiedlic displacement measuring system for an anthropomorphic testing device |
US9044537B2 (en) * | 2007-03-30 | 2015-06-02 | Medtronic, Inc. | Devices and methods for detecting catheter complications |
US20110257593A1 (en) * | 2007-03-30 | 2011-10-20 | Medtronic, Inc. | Devices and methods for detecting catheter complications |
US8412293B2 (en) | 2007-07-16 | 2013-04-02 | Optiscan Biomedical Corporation | Systems and methods for determining physiological parameters using measured analyte values |
US20090043171A1 (en) * | 2007-07-16 | 2009-02-12 | Peter Rule | Systems And Methods For Determining Physiological Parameters Using Measured Analyte Values |
US20170188943A1 (en) * | 2007-10-10 | 2017-07-06 | Optiscan Biomedical Corporation | Fluid component analysis systems and methods for glucose monitoring and control |
US8231562B2 (en) * | 2008-09-05 | 2012-07-31 | Roche Diagnostics Operations, Inc. | Insulin pump configuration programming invalid settings notification and correction |
US20100064257A1 (en) * | 2008-09-05 | 2010-03-11 | Roche Diagnostics Operations, Inc. | Insulin pump configuration programming invalid settings notification and correction |
US10279350B2 (en) | 2008-09-23 | 2019-05-07 | Bio-Rad Laboratories, Inc. | Method of generating droplets |
US9156010B2 (en) | 2008-09-23 | 2015-10-13 | Bio-Rad Laboratories, Inc. | Droplet-based assay system |
US11633739B2 (en) | 2008-09-23 | 2023-04-25 | Bio-Rad Laboratories, Inc. | Droplet-based assay system |
US9636682B2 (en) | 2008-09-23 | 2017-05-02 | Bio-Rad Laboratories, Inc. | System for generating droplets—instruments and cassette |
US9649635B2 (en) | 2008-09-23 | 2017-05-16 | Bio-Rad Laboratories, Inc. | System for generating droplets with push-back to remove oil |
US11612892B2 (en) | 2008-09-23 | 2023-03-28 | Bio-Rad Laboratories, Inc. | Method of performing droplet-based assays |
US11130134B2 (en) | 2008-09-23 | 2021-09-28 | Bio-Rad Laboratories, Inc. | Method of performing droplet-based assays |
US9492797B2 (en) | 2008-09-23 | 2016-11-15 | Bio-Rad Laboratories, Inc. | System for detection of spaced droplets |
US11130128B2 (en) | 2008-09-23 | 2021-09-28 | Bio-Rad Laboratories, Inc. | Detection method for a target nucleic acid |
US8633015B2 (en) | 2008-09-23 | 2014-01-21 | Bio-Rad Laboratories, Inc. | Flow-based thermocycling system with thermoelectric cooler |
US9417190B2 (en) | 2008-09-23 | 2016-08-16 | Bio-Rad Laboratories, Inc. | Calibrations and controls for droplet-based assays |
US10512910B2 (en) | 2008-09-23 | 2019-12-24 | Bio-Rad Laboratories, Inc. | Droplet-based analysis method |
US9126160B2 (en) | 2008-09-23 | 2015-09-08 | Bio-Rad Laboratories, Inc. | System for forming an array of emulsions |
US9132394B2 (en) | 2008-09-23 | 2015-09-15 | Bio-Rad Laboratories, Inc. | System for detection of spaced droplets |
US9764322B2 (en) | 2008-09-23 | 2017-09-19 | Bio-Rad Laboratories, Inc. | System for generating droplets with pressure monitoring |
US9623384B2 (en) | 2008-09-23 | 2017-04-18 | Bio-Rad Laboratories, Inc. | System for transporting emulsions from an array to a detector |
US10258988B2 (en) | 2008-09-23 | 2019-04-16 | Bio-Rad Laboratories, Inc. | Device for generating droplets |
US9216392B2 (en) | 2008-09-23 | 2015-12-22 | Bio-Rad Laboratories, Inc. | System for forming an array of emulsions |
US9248417B2 (en) | 2008-09-23 | 2016-02-02 | Bio-Rad Laboratories, Inc. | System for droplet-based assays using an array of emulsions |
US10258989B2 (en) | 2008-09-23 | 2019-04-16 | Bio-Rad Laboratories, Inc. | Method of making a device for generating droplets |
US9243288B2 (en) | 2008-09-23 | 2016-01-26 | Bio-Rad Laboratories, Inc. | Cartridge with lysis chamber and droplet generator |
US20110092784A1 (en) * | 2009-07-20 | 2011-04-21 | Optiscan Biomedical Corporation | Adjustable connector and dead space reduction |
US9326717B2 (en) * | 2009-07-20 | 2016-05-03 | Optiscan Biomedical Corporation | Adjustable connector and dead space reduction |
US10028692B2 (en) | 2009-07-20 | 2018-07-24 | Optiscan Biomedical Corporation | Adjustable connector, improved fluid flow and reduced clotting risk |
WO2011011462A1 (en) * | 2009-07-20 | 2011-01-27 | Optiscan Biomedical Corporation | Adjustable connector and dead space reduction |
US20140228710A1 (en) * | 2009-07-20 | 2014-08-14 | Optiscan Biomedical Corporation | Adjustable connector and dead space reduction |
US8731639B2 (en) | 2009-07-20 | 2014-05-20 | Optiscan Biomedical Corporation | Adjustable connector, improved fluid flow and reduced clotting risk |
US8731638B2 (en) | 2009-07-20 | 2014-05-20 | Optiscan Biomedical Corporation | Adjustable connector and dead space reduction |
US9194861B2 (en) | 2009-09-02 | 2015-11-24 | Bio-Rad Laboratories, Inc. | Method of mixing fluids by coalescence of multiple emulsions |
US10166522B2 (en) | 2009-09-02 | 2019-01-01 | Bio-Rad Laboratories, Inc. | System for mixing fluids by coalescence of multiple emulsions |
US10677693B2 (en) | 2009-09-02 | 2020-06-09 | Bio-Rad Laboratories, Inc. | System for mixing fluids by coalescence of multiple emulsions |
US10378048B2 (en) | 2010-03-02 | 2019-08-13 | Bio-Rad Laboratories, Inc. | Emulsion chemistry for encapsulated droplets |
US11866771B2 (en) | 2010-03-02 | 2024-01-09 | Bio-Rad Laboratories, Inc. | Emulsion chemistry for encapsulated droplets |
US9598725B2 (en) | 2010-03-02 | 2017-03-21 | Bio-Rad Laboratories, Inc. | Emulsion chemistry for encapsulated droplets |
US11060136B2 (en) | 2010-03-02 | 2021-07-13 | Bio-Rad Laboratories, Inc. | Emulsion chemistry for encapsulated droplets |
US8709762B2 (en) | 2010-03-02 | 2014-04-29 | Bio-Rad Laboratories, Inc. | System for hot-start amplification via a multiple emulsion |
US10099219B2 (en) | 2010-03-25 | 2018-10-16 | Bio-Rad Laboratories, Inc. | Device for generating droplets |
US9500664B2 (en) | 2010-03-25 | 2016-11-22 | Bio-Rad Laboratories, Inc. | Droplet generation for droplet-based assays |
US9393560B2 (en) | 2010-03-25 | 2016-07-19 | Bio-Rad Laboratories, Inc. | Droplet transport system for detection |
US10272432B2 (en) | 2010-03-25 | 2019-04-30 | Bio-Rad Laboratories, Inc. | Device for generating droplets |
US8730479B2 (en) | 2010-03-25 | 2014-05-20 | Bio-Rad Laboratories, Inc. | Detection system for droplet-based assays |
US10744506B2 (en) | 2010-03-25 | 2020-08-18 | Bio-Rad Laboratories, Inc. | Device for generating droplets |
US9089844B2 (en) | 2010-11-01 | 2015-07-28 | Bio-Rad Laboratories, Inc. | System for forming emulsions |
US9222128B2 (en) | 2011-03-18 | 2015-12-29 | Bio-Rad Laboratories, Inc. | Multiplexed digital assays with combinatorial use of signals |
US10190115B2 (en) | 2011-04-25 | 2019-01-29 | Bio-Rad Laboratories, Inc. | Methods and compositions for nucleic acid analysis |
US9347059B2 (en) | 2011-04-25 | 2016-05-24 | Bio-Rad Laboratories, Inc. | Methods and compositions for nucleic acid analysis |
US11939573B2 (en) | 2011-04-25 | 2024-03-26 | Bio-Rad Laboratories, Inc. | Methods and compositions for nucleic acid analysis |
US10760073B2 (en) | 2011-04-25 | 2020-09-01 | Bio-Rad Laboratories, Inc. | Methods and compositions for nucleic acid analysis |
US9885034B2 (en) | 2011-04-25 | 2018-02-06 | Bio-Rad Laboratories, Inc. | Methods and compositions for nucleic acid analysis |
US9872956B2 (en) | 2011-04-29 | 2018-01-23 | Medtronic, Inc. | Limiting pressure in an implanted catheter |
US8951939B2 (en) | 2011-07-12 | 2015-02-10 | Bio-Rad Laboratories, Inc. | Digital assays with multiplexed detection of two or more targets in the same optical channel |
CN103930886A (en) * | 2011-07-13 | 2014-07-16 | 伯乐生命医学产品有限公司 | Computation of real-world error using META-analysis of replicates |
WO2013010142A1 (en) * | 2011-07-13 | 2013-01-17 | Bio-Rad Laboratories, Inc., A Delawae Corporation | Computation of real-world error using meta-analysis of replicates |
US8663920B2 (en) | 2011-07-29 | 2014-03-04 | Bio-Rad Laboratories, Inc. | Library characterization by digital assay |
US20130053099A1 (en) * | 2011-08-25 | 2013-02-28 | Kyocera Corporation | Human body detecting system |
US9377361B2 (en) * | 2011-08-25 | 2016-06-28 | Kyocera Corporation | Human body detecting system |
US10508988B2 (en) * | 2011-12-05 | 2019-12-17 | Gassecure As | Method and system for gas detection |
US20180217055A1 (en) * | 2011-12-05 | 2018-08-02 | Gassecure As | Method and system for gas detection |
EP2604998A1 (en) * | 2011-12-12 | 2013-06-19 | ABB Research Ltd. | Gas analysis device |
US9399215B2 (en) | 2012-04-13 | 2016-07-26 | Bio-Rad Laboratories, Inc. | Sample holder with a well having a wicking promoter |
US20150378511A1 (en) * | 2014-06-26 | 2015-12-31 | Sitronix Technology Corp. | Capacitive Voltage Information Sensing Circuit and Related Anti-Noise Touch Circuit |
US9524056B2 (en) * | 2014-06-26 | 2016-12-20 | Sitronix Technology Corp. | Capacitive voltage information sensing circuit and related anti-noise touch circuit |
US11726036B2 (en) | 2014-12-23 | 2023-08-15 | Apple Inc. | Optical inspection system and method including accounting for variations of optical path length within a sample |
US11585749B2 (en) | 2015-09-01 | 2023-02-21 | Apple Inc. | Reference switch architectures for noncontact sensing of substances |
US20180238794A1 (en) * | 2015-09-01 | 2018-08-23 | Apple Inc. | Reference switch architectures for noncontact sensing of substances |
US10801950B2 (en) * | 2015-09-01 | 2020-10-13 | Apple Inc. | Reference switch architectures for noncontact sensing of substances |
US12007275B2 (en) | 2016-04-21 | 2024-06-11 | Apple Inc. | Optical system for reference switching |
US11579080B2 (en) | 2017-09-29 | 2023-02-14 | Apple Inc. | Resolve path optical sampling architectures |
US11960131B2 (en) | 2018-02-13 | 2024-04-16 | Apple Inc. | Integrated photonics device having integrated edge outcouplers |
CN109682958A (en) * | 2018-09-21 | 2019-04-26 | 深圳沃德生命科技有限公司 | A kind of acceleration transducer signals compensation method for thrombelastogram instrument |
US11852318B2 (en) | 2020-09-09 | 2023-12-26 | Apple Inc. | Optical system for noise mitigation |
Also Published As
Publication number | Publication date |
---|---|
WO2007120745A2 (en) | 2007-10-25 |
WO2007120745A3 (en) | 2008-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20070258083A1 (en) | Noise reduction for analyte detection systems | |
US8425444B2 (en) | Anti-clotting apparatus and methods for fluid handling system | |
US20080072663A1 (en) | Accurate and timely body fluid analysis | |
US10677688B2 (en) | Fluid injection and safety system | |
EP2580589B1 (en) | Measuring analytes in a fluid sample drawn from a patient | |
EP2695573B1 (en) | Fluid component analysis system and method for glucose monitoring and control | |
US9863837B2 (en) | Systems and methods for detecting leaks | |
US20180263555A1 (en) | Systems and methods to reduce fluid contamination | |
US20160371451A1 (en) | Analyte monitoring systems and methods | |
US20080161723A1 (en) | Infusion flow interruption method and apparatus | |
WO2011011462A1 (en) | Adjustable connector and dead space reduction | |
WO2009049245A1 (en) | Synchronization and configuration of patient monitoring devices | |
US20090032111A1 (en) | Fluid mixing systems and methods | |
US20100145175A1 (en) | Systems and methods for verification of sample integrity | |
EP2729784A1 (en) | Sample cell for fluid analysis system | |
US20090156911A1 (en) | Low draw volume analyte detection systems | |
US20200227144A1 (en) | Method and apparatus for analyte measurements using calibration sets | |
US20090160656A1 (en) | Analyte monitoring system alarms |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: OPTISCAN BIOMEDICAL INCORPORATED, CALIFORNIA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:HEPPELL, KEVIN;BRAIG, JAMES R.;REEL/FRAME:019492/0753 Effective date: 20070626 |
|
AS | Assignment |
Owner name: HERCULES TECHNOLOGY II, L.P., CALIFORNIA Free format text: SECURITY AGREEMENT;ASSIGNOR:OPTISCAN BIOMEDICAL CORPORATION;REEL/FRAME:020995/0202 Effective date: 20080522 Owner name: HERCULES TECHNOLOGY II, L.P.,CALIFORNIA Free format text: SECURITY AGREEMENT;ASSIGNOR:OPTISCAN BIOMEDICAL CORPORATION;REEL/FRAME:020995/0202 Effective date: 20080522 |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |
|
AS | Assignment |
Owner name: OPTISCAN BIOMEDICAL CORPORATION, CALIFORNIA Free format text: RELEASE BY SECURED PARTY;ASSIGNOR:HERCULES TECHNOLOGY II, L.P.;REEL/FRAME:041344/0534 Effective date: 20170111 |