CN101079082A - 用于处理自我监测血糖(smbg)数据从而提高糖尿病患者自我管理的方法、***和计算机程序产品 - Google Patents
用于处理自我监测血糖(smbg)数据从而提高糖尿病患者自我管理的方法、***和计算机程序产品 Download PDFInfo
- Publication number
- CN101079082A CN101079082A CNA2007101032423A CN200710103242A CN101079082A CN 101079082 A CN101079082 A CN 101079082A CN A2007101032423 A CNA2007101032423 A CN A2007101032423A CN 200710103242 A CN200710103242 A CN 200710103242A CN 101079082 A CN101079082 A CN 101079082A
- Authority
- CN
- China
- Prior art keywords
- lbgi
- data
- hba
- rlo
- months
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 103
- 206010012601 diabetes mellitus Diseases 0.000 title abstract description 77
- 239000008280 blood Substances 0.000 title abstract description 54
- 210000004369 blood Anatomy 0.000 title abstract description 53
- 239000008103 glucose Substances 0.000 title abstract description 39
- 238000012544 monitoring process Methods 0.000 title abstract description 36
- 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 title abstract description 34
- 238000004590 computer program Methods 0.000 title abstract description 28
- 238000012545 processing Methods 0.000 title description 5
- 208000013016 Hypoglycemia Diseases 0.000 claims abstract description 183
- 230000002218 hypoglycaemic effect Effects 0.000 claims abstract description 182
- 230000007774 longterm Effects 0.000 claims abstract description 31
- 230000002045 lasting effect Effects 0.000 claims description 51
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims description 48
- 230000007246 mechanism Effects 0.000 claims description 8
- 201000001421 hyperglycemia Diseases 0.000 abstract description 12
- 238000012423 maintenance Methods 0.000 abstract description 2
- 238000012806 monitoring device Methods 0.000 abstract 2
- 230000002411 adverse Effects 0.000 abstract 1
- 238000004422 calculation algorithm Methods 0.000 description 204
- 238000012360 testing method Methods 0.000 description 105
- 208000001072 type 2 diabetes mellitus Diseases 0.000 description 70
- 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 46
- 238000012549 training Methods 0.000 description 39
- 238000011160 research Methods 0.000 description 37
- 101100121625 Caenorhabditis elegans rhi-1 gene Proteins 0.000 description 35
- 230000003203 everyday effect Effects 0.000 description 32
- 230000006870 function Effects 0.000 description 26
- 102000004877 Insulin Human genes 0.000 description 23
- 108090001061 Insulin Proteins 0.000 description 23
- 230000002641 glycemic effect Effects 0.000 description 23
- 229940125396 insulin Drugs 0.000 description 23
- 238000009826 distribution Methods 0.000 description 22
- 238000004458 analytical method Methods 0.000 description 21
- 101000655352 Homo sapiens Telomerase reverse transcriptase Proteins 0.000 description 18
- 102100032938 Telomerase reverse transcriptase Human genes 0.000 description 18
- 238000004891 communication Methods 0.000 description 18
- 239000000126 substance Substances 0.000 description 18
- 238000005259 measurement Methods 0.000 description 17
- 206010067584 Type 1 diabetes mellitus Diseases 0.000 description 15
- 230000000694 effects Effects 0.000 description 15
- 238000005457 optimization Methods 0.000 description 15
- 230000008859 change Effects 0.000 description 14
- 230000002776 aggregation Effects 0.000 description 11
- 238000004220 aggregation Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 11
- 238000007726 management method Methods 0.000 description 11
- 238000003860 storage Methods 0.000 description 11
- 108010014663 Glycated Hemoglobin A Proteins 0.000 description 10
- 102000017011 Glycated Hemoglobin A Human genes 0.000 description 10
- 238000000540 analysis of variance Methods 0.000 description 10
- 230000008878 coupling Effects 0.000 description 10
- 238000010168 coupling process Methods 0.000 description 10
- 238000005859 coupling reaction Methods 0.000 description 10
- 238000011835 investigation Methods 0.000 description 10
- 230000008901 benefit Effects 0.000 description 9
- 230000002354 daily effect Effects 0.000 description 9
- 239000000243 solution Substances 0.000 description 9
- 238000011156 evaluation Methods 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 7
- 230000000875 corresponding effect Effects 0.000 description 7
- 238000012417 linear regression Methods 0.000 description 7
- 208000010513 Stupor Diseases 0.000 description 6
- 238000012937 correction Methods 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 6
- 108010041872 Islet Amyloid Polypeptide Proteins 0.000 description 5
- 102000036770 Islet Amyloid Polypeptide Human genes 0.000 description 5
- 238000011161 development Methods 0.000 description 5
- 235000012054 meals Nutrition 0.000 description 5
- 210000002381 plasma Anatomy 0.000 description 5
- 230000009466 transformation Effects 0.000 description 5
- 241001269238 Data Species 0.000 description 4
- 208000003443 Unconsciousness Diseases 0.000 description 4
- 235000021152 breakfast Nutrition 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000011545 laboratory measurement Methods 0.000 description 4
- 238000012886 linear function Methods 0.000 description 4
- 230000036961 partial effect Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 208000035408 type 1 diabetes mellitus 1 Diseases 0.000 description 4
- 230000003442 weekly effect Effects 0.000 description 4
- 238000012935 Averaging Methods 0.000 description 3
- 239000012491 analyte Substances 0.000 description 3
- 238000013459 approach Methods 0.000 description 3
- 238000003556 assay Methods 0.000 description 3
- 230000033228 biological regulation Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 230000000295 complement effect Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 239000003814 drug Substances 0.000 description 3
- 230000036541 health Effects 0.000 description 3
- 230000004060 metabolic process Effects 0.000 description 3
- 229940126701 oral medication Drugs 0.000 description 3
- 230000002441 reversible effect Effects 0.000 description 3
- 230000000630 rising effect Effects 0.000 description 3
- 210000004243 sweat Anatomy 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- GLDQAMYCGOIJDV-UHFFFAOYSA-N 2,3-dihydroxybenzoic acid Chemical compound OC(=O)C1=CC=CC(O)=C1O GLDQAMYCGOIJDV-UHFFFAOYSA-N 0.000 description 2
- 235000017166 Bambusa arundinacea Nutrition 0.000 description 2
- 235000017491 Bambusa tulda Nutrition 0.000 description 2
- 241001330002 Bambuseae Species 0.000 description 2
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 2
- 206010054998 Neuroglycopenia Diseases 0.000 description 2
- 235000015334 Phyllostachys viridis Nutrition 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 230000001154 acute effect Effects 0.000 description 2
- 239000011425 bamboo Substances 0.000 description 2
- 238000004140 cleaning Methods 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 230000008034 disappearance Effects 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 230000003345 hyperglycaemic effect Effects 0.000 description 2
- 230000008676 import Effects 0.000 description 2
- 230000000670 limiting effect Effects 0.000 description 2
- 230000014759 maintenance of location Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- XEBWQGVWTUSTLN-UHFFFAOYSA-M phenylmercury acetate Chemical compound CC(=O)O[Hg]C1=CC=CC=C1 XEBWQGVWTUSTLN-UHFFFAOYSA-M 0.000 description 2
- 229960003611 pramlintide Drugs 0.000 description 2
- 108010029667 pramlintide Proteins 0.000 description 2
- NRKVKVQDUCJPIZ-MKAGXXMWSA-N pramlintide acetate Chemical compound C([C@@H](C(=O)NCC(=O)N1CCC[C@H]1C(=O)N[C@@H]([C@@H](C)CC)C(=O)N[C@@H](CC(C)C)C(=O)N1[C@@H](CCC1)C(=O)N1[C@@H](CCC1)C(=O)N[C@@H]([C@@H](C)O)C(=O)N[C@@H](CC(N)=O)C(=O)N[C@@H](C(C)C)C(=O)NCC(=O)N[C@@H](CO)C(=O)N[C@@H](CC(N)=O)C(=O)N[C@@H]([C@@H](C)O)C(=O)N[C@@H](CC=1C=CC(O)=CC=1)C(N)=O)NC(=O)[C@H](CC(N)=O)NC(=O)[C@H](CC(N)=O)NC(=O)[C@H](CO)NC(=O)[C@H](CO)NC(=O)[C@H](CC=1NC=NC=1)NC(=O)[C@@H](NC(=O)[C@H](CC(C)C)NC(=O)[C@H](CC=1C=CC=CC=1)NC(=O)[C@H](CC(N)=O)NC(=O)[C@H](C)NC(=O)[C@H](CC(C)C)NC(=O)[C@H](CCCNC(N)=N)NC(=O)[C@H](CCC(N)=O)NC(=O)[C@@H](NC(=O)[C@H](C)NC(=O)[C@H](CS)NC(=O)[C@@H](NC(=O)[C@H](C)NC(=O)[C@@H](NC(=O)[C@H](CC(N)=O)NC(=O)[C@H](CS)NC(=O)[C@@H](N)CCCCN)[C@@H](C)O)[C@@H](C)O)[C@@H](C)O)C(C)C)C1=CC=CC=C1 NRKVKVQDUCJPIZ-MKAGXXMWSA-N 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 238000005728 strengthening Methods 0.000 description 2
- 208000024891 symptom Diseases 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 208000028698 Cognitive impairment Diseases 0.000 description 1
- 101001051478 Equus caballus Beta-lactoglobulin-1 Proteins 0.000 description 1
- 240000004859 Gamochaeta purpurea Species 0.000 description 1
- 206010018473 Glycosuria Diseases 0.000 description 1
- 102100030373 HSPB1-associated protein 1 Human genes 0.000 description 1
- 102100021467 Histone acetyltransferase type B catalytic subunit Human genes 0.000 description 1
- 101000843045 Homo sapiens HSPB1-associated protein 1 Proteins 0.000 description 1
- 101000898976 Homo sapiens Histone acetyltransferase type B catalytic subunit Proteins 0.000 description 1
- 238000004566 IR spectroscopy Methods 0.000 description 1
- 206010021703 Indifference Diseases 0.000 description 1
- 206010022489 Insulin Resistance Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000009098 adjuvant therapy Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 239000003833 bile salt Substances 0.000 description 1
- 238000010241 blood sampling Methods 0.000 description 1
- 208000027503 bloody stool Diseases 0.000 description 1
- 150000001720 carbohydrates Chemical class 0.000 description 1
- 235000014633 carbohydrates Nutrition 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000000718 cholinopositive effect Effects 0.000 description 1
- 208000010877 cognitive disease Diseases 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 238000005183 dynamical system Methods 0.000 description 1
- 230000002124 endocrine Effects 0.000 description 1
- 230000007717 exclusion Effects 0.000 description 1
- 210000003722 extracellular fluid Anatomy 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 235000013305 food Nutrition 0.000 description 1
- 150000002303 glucose derivatives Chemical class 0.000 description 1
- 208000035861 hematochezia Diseases 0.000 description 1
- 108010013766 hemoglobin A(0) Proteins 0.000 description 1
- 229940088597 hormone Drugs 0.000 description 1
- 239000005556 hormone Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000003834 intracellular effect Effects 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- XZWYZXLIPXDOLR-UHFFFAOYSA-N metformin Chemical compound CN(C)C(=N)NC(N)=N XZWYZXLIPXDOLR-UHFFFAOYSA-N 0.000 description 1
- 229960003105 metformin Drugs 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 239000006187 pill Substances 0.000 description 1
- 238000002203 pretreatment Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000012797 qualification Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000000306 recurrent effect Effects 0.000 description 1
- 230000002829 reductive effect Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000028327 secretion Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 210000002784 stomach Anatomy 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000009469 supplementation Effects 0.000 description 1
- 210000000106 sweat gland Anatomy 0.000 description 1
- 239000013076 target substance Substances 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 239000002023 wood Substances 0.000 description 1
Images
Landscapes
- Investigating Or Analysing Biological Materials (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种方法、***和计算机程序产品涉及根据由自我监测血糖设备收集的血糖数据保持对糖尿病的最佳控制,并致力于预测对高血糖症的长期暴露,以及糖尿病患者罹患严重或者中度低血糖症的长期和短期危险。本方法、***和计算机程序产品适合于通过引入一种能够同时预测HbA1c和低血糖症高危期的智能数据判读组件直接提高现有家用血糖监测设备的性能,并通过相同的特征提高将来生产的连续监测设备的性能。利用这些预测,糖尿病患者能够采取措施防止与高血糖症和低血糖症有关的负面后果。
Description
相关专利申请
本分案申请是基于申请号为03824009.2,申请日为2003年08月08日,发明名称为“用于处理自我监测血糖(SMBG)数据从而提高糖尿病患者自我管理的方法、***和计算机程序产品”的中国专利申请的分案申请。
本专利申请是于2003年8月提出申请的国际专利申请No.PCT/US2003/025053的国内申请阶段,其根据35U.S.C.条款119(e)要求获得如下专利申请的优先权,即美国临时专利申请系列No.60/402,976,其于2002年8月13日提出申请,标题为“用于处理自我监测血糖(SMBG)数据从而提高糖尿病患者自我管理的方法、***和计算机程序产品”,和No.60/478,377,其于2003年6月13日提出申请,标题为“用于处理自我监测血糖(SMBG)数据从而提高糖尿病患者自我管理的方法、***和计算机程序产品”,本文引用上述三个公开的全部公开内容作为参考。
专利申请涉及国际申请No.PCT/US01/09884,其于2001年3月29日提出申请(专利申请Nos.WO 01/72208 A2,WO 01/72208 A3),标题为“由自我监测数据评估糖尿病血糖控制的方法、***和计算机程序产品”,和美国专利申请系列No.:10/240,228,其于2002年9月26日提出申请,标题为“由自我监测数据评估糖尿病血糖控制的方法、***和计算机程序产品”,本文引用其全部的公开内容作为参考。
技术领域
本***总体上涉及糖尿病个体的血糖控制,更特别地,涉及一种基于计算机的***和方法,用于评估预测糖基化血红蛋白(HbA1c和HbA1)和发生低血糖症的危险。
背景技术
大量的研究反复证实,防止糖尿病长期并发症的最有效方法是采取胰岛素强化治疗将血糖(BG)水平严格控制在正常范围之内,这些研究包括糖尿病控制和并发症实验(DCCT)(见DCCT研究会:糖尿病强化治疗对胰岛素依赖型糖尿病长期并发症的发生和发展的影响(The Effect Of Intensive Treatment Of Diabetes On The DevelopmentAnd Progresion Of Long-Term Complications Of Insulin-DependentDiabetes Mellitus).New England Journal ofMedicine,329:978-986,1993)),斯德哥尔摩糖尿病干预研究(见Reichard P,Phil M:斯德哥尔摩糖尿病干预研究中在传统胰岛素长期强化治疗期间的死亡率和治疗副作用(Mortality and TreatmentSide Effects Druing Long-term Intensified Conventional InsulinTreatment in the Stockholm Diabetes Intervention Study).Diabetes,3:313-317,1994)),和英国前瞻性糖尿病研究(见英国前瞻性糖尿病研究会:用二甲双胍进行强化血糖控制对于2型糖尿病患者并发症的影响(Effect of Intensive Blood Glucose Control With Metformin OnComplications In Patients With Type 2 Diabetes)(UKPDS 34).
Lancet,352:837-853,1998)。
然而,相同的研究也证实了胰岛素强化治疗的一些副作用,最严重的是会增加频繁严重低血糖症(SH)的危险,这是一种无法进行自我治疗和需要外界帮助才能恢复的神经低血糖症事件(见DCCT研究会:糖尿病控制和并发症实验中严重低血糖症的流行病学(Epidemiology of Severe Hypoglycemia In The Diabetes Control andComplications Trial).American Journal of Medicine,90:450-459,1991,和DCCT研究会:糖尿病控制和并发症实验中的低血糖症(Hypoglycemia in the Diabetes Control and ComplicationsTrial).Diabetes,46:271-286,1997)。因为SH会导致意外事故、昏迷甚至死亡,所以患者和康护人员对继续进行强化治疗感到灰心。结果,低血糖症被认定为是提高血糖控制的一个主要障碍(Cryer PE:低血糖症是糖尿病管理的限制因素(Hypoglycemia is the Limiting Factor inthe Management Of Diabetes).Diabetes Metab ResRev,15:42-46,1999)。
因此,糖尿病患者面临着一个终生的优化问题,即在保持对血糖严格控制的同时又不增加低血糖症的危险。与这个问题有关的主要挑战是产生一种能够同时评估患者的血糖控制及其低血糖症危险,并且能够在日常环境中使用的简单而可靠的方法。
二十年来人们已经熟知,糖基化血红蛋白是糖尿病(1型或2型)个体血糖控制的标志。大量的研究人员对这一关系进行了研究并且发现,糖基化血红蛋白基本上反映了患者过去两个月内的平均BG水平。因为在大多数糖尿病患者体内,BG水平在一段时间内会有相当大的波动,所以建议整体血糖控制与HbA1c之间的实际关联只能在已知患者在一个较长时期内处于稳定的血糖控制的情况下进行观察。
对该类患者的早期研究建立了前5周内平均BG水平与HbA1c之间的几乎确定的关系,并且该曲线关系产生了大小为0.98的相关系数(见Aaby Svendsen P,Lauritzen T,Soegard U,Nerup J(1982).1型(胰岛素依赖型)糖尿病中的糖基化血红蛋白与稳态平均血糖浓度(Glycosylated Hemoglobin and Steady-State Mean Blood GlucoseConcentration in Type 1(Insulin-Dependent)Diabetes).
Diabetologia,23,403-405)。1993年,DCCT得出结论,HbA1c是黄金标准糖基化血红蛋白化验(gold-standard glycosylated hemoglobin assay)的“合理推荐”,并且DCCT确定了先前平均BG与HbA1c之间的线性关系(见Santiago JV(1993).来自糖尿病控制与并发症实验的教训(Lessonsfrom the Diabetes Control and Complications Trial),
Diabetes,42,1549-1554)。
已提出的指导方针表明,7%的HbA1c相应于8.3mM(150mg/dl)的平均BG,9%的HbA1c相应于11.7mM(210mg/dl)的平均BG,并且HbA1c增加1%相应于平均BG增加1.7mM(30mg/dl,2)。DCCT还建议,因为直接测量平均BG并不现实,所以可以用单一的简单测试,即HbA1c,评估患者的血糖控制。然而,研究清晰地表明,HbA1c对低血糖症并不敏感。
确实,由任何数据都不能得到患者SH直接危险的可靠预测值。DCCT得出结论,只有大约8%的将来SH能够由已知的参数,例如SH历史、低HbA1c和低血糖症昏迷(unawareness),预报出来。一篇近期的综述详细介绍了该问题的当前临床状况,并且为预防SH提供了患者及其康护人员能够获得的选择(见Bolli,GB:如何在1型糖尿病强化及非强化治疗中改善低血糖症问题(How To Ameliorate ThePreblem of Hypoglycemia In Intensive As Well As NonintensiveTreatment Of Type I Diabetes).Diabetes Care,22,Supplement2:B43-B52,1999)。
现代家用BG监测器通过自我监测BG(SMBG)提供了进行频繁BG测量的装置。然而,SMBG的问题在于,通过BG监测器收集的数据与HbA1c和低血糖症之间缺乏联系。换言之,目前还没有可靠的方法能够根据SMBG读数估计HbA1c和识别即将发生的低血糖症(见Bremer T和Gough DA:血糖能够由先前值预测出来吗?数据的引发(Is blood glucose predictable from previous values?A solicitationfor data).Diabetes 48:445-451,1999)。
因此,本发明的一个目的是通过提出三种不同但相互兼容的算法给出这一缺失的联系,三种算法用于由SMBG数据估计HbA1c和低血糖症危险,以预测低血糖症的长期和短期危险以及高血糖症的长期危险。
发明人在先前曾经报道过,SMBG数据的日常获得与估计HbA1c和低血糖症危险之间缺乏联系的原因在于,在糖尿病研究中使用的数据收集与临床评估的精密方法很少获得糖尿病特异的和数学上精密的统计处理的支持。
出于对能够顾及BG特殊分布的统计分析的需要,发明人提出了一种对血糖测量值范围的对称变换(见Kovatchev BP,Cox DJ,Gonder-Frederick LA和WL Clarke(1997).血糖测量值范围的对称化及其应用(Symmetization of the Blood Glucose Measurement Scaleand Its Applications),
Diabetes Care,
20,1655-1658),其操作如下。BG水平在美国以mg/dl为单位进行测量,在其他大多数国家中以mmol/L(或mM)为单位。两种标度的直接关系为18mg/dl=1mM。在大多数参考文献中给出的整体BG值范围是1.1-33.3mM,并且认为它实际地覆盖了全部的观察值。根据DCCT的推荐(见DCCT研究会(1993),糖尿病强化治疗对胰岛素依赖型糖尿病并发症的发生和发展的影响(TheEffect Of Intensive Treatment Of Diabetes On The Development AndProgresion Of Long-Term Complications Of Insulin-DependentDiabetes Mellitus).New England Journal ofMedicine,329:978-986,1993)),糖尿病人的目标BG值范围——也称作正常血糖范围——是3.9-10mM,当BG降到低于3.9mM时则出现低血糖症,当BG升高到10mM以上时则出现高血糖症。不幸的是,该范围在数值上并不对称——高血糖范围(10-33.3mM)宽于低血糖范围(1.1-3.9Mm),且正常血糖范围(3.9-10mM)不在该范围的中心。发明人通过引入一种变换f(BG)修正了该对称性,其中f(BG)是定义域为BG范围[1.1,33.3]上的连续函数,具有双参数解析形式:
f(BG,α,β)=[(ln(BG))α-β],α,β>0
它满足假设:
A1:f(33.3,α,β)=-f(1.1,α,β)和
A2:f(10.0,α,β)=-f(3.9,α,β)。
接着,f(.)乘以第三比例参数从而将经过变换的BG范围的最小和最大值分别固定在
和
。这些数值很方便,因为一个具有标准正态分布的随意变量在区间
内具有99.8%的数值。如果BG以mmol/L为单位测量,则根据假设A1和A2进行数值解时,函数f(BG,α,β)的参数为α=1.026,β=1.861,比例参数γ=1.794。如果BG以mg/dl为单位测量,则计算出的参数为α=1.084,β=5.381,γ=1.509。
因此,当以mmol/L为单位测量BG时,对称变换是f(BG)=1.794[(ln(BG))1.026-1.861],当以mg/dl为单位测量BG时,对称变换为f(BG)=1.509[(ln(BG))1.084-5.381]。
根据对称变换f(.),发明人引入了低BG指数——一个用于由SMBG读数估计低血糖症危险的新测量值(见Cox DJ,Kovatchev BP,Julian DM,Gonder-Frederick LA,Polonsky WH,Schlundt DG,Clarke WL:IDDM中的严重低血糖症频率能够由自我监测血糖数据进行预测(Frequency of Severe Hypoglycemia In IDDM Can BePredicted From Self-Monitoring Blood Glucose Data).Journal ofClinical Endocrinology and Metabolism,79:1659-1662,1994,和Kovatchev BP,Cox DJ,Gonder-Frederick LA,Young-Hyman D,Schlundt D,Clarke WL.IDDM成年人中严重低血糖症危险估计:低血糖指数的有效性(Assessment of Risk for Severe HypoglycemiaAmong Adults With IDDM:Validation of the Low Blood GlucoseIndex),Diabetes Care 21:1870-1875,1998)。给定一系列SMBG数据,当f(BG)<0时,低BG指数计算为10.f(BG)2的平均值,否则为0。还提出了高BG指数,按照和低BG指数对称的方式计算而得,但是该指数还没有发现其实际用途。
在回归模型中使用低BG指数,发明人能够根据SH历史和SMBG数据解释随后6个月中SH事件40%的变异,后来将该预测提高到了46%(见Kovatchev BP,Straume M,Farhi LS,Cox DJ:估计血糖转变速度及其与严重低血糖的关系(Estimating the Speed ofBlood Glucose Transitions and its Relationship With SevereHypoglycemia),Diabetes,48:Supplement 1,A363,1999)。
此外,发明人还提出了一些关于HbA1c和SMBG的数据(见Kovatchev BP,Cox DJ,Straume M,Farhi LS.自我监测血糖曲线与糖基化血红蛋白的关系(Association of Self-monitoring Blood GlucoseProfiles with Glycosylated Hemoglobin),
In:Methods in Enzymology, vol.321:Numerical Computer Methods,Part C,Mechael Johnson和Ludvig Brand,Eds.,Academic Press,NY;2000)。
这些研究成果成为本发明理论背景的一部分。为了将该理论付诸于实践,添加了几个关键的理论参量,在下面的部分中将进行说明。特别地,提出了三种方法用于估计HbA1c、低血糖症的长期和短期危险。这些方法的提出是基于,但不仅限于,对867个糖尿病个体超过300,000个SMBG读数、严重低血糖症记录和HbA1c结果的详细分析。
因此,发明人试图改进与传统方法有关的前述限制,借此提供简单而可靠的方法,从而能够用于同时评估患者的血糖控制及其低血糖症的危险,并能够在其日常环境中使用。
发明内容
本发明包括一种数据分析方法和基于计算机的***,用于由日常收集的SMBG数据同时估计糖尿病血糖控制中两种最重要的分量:HbA1c和低血糖症危险。出于本文献的目的,BG自我监测(SMBG)被定义为用于在糖尿病患者的自然条件下确定血糖的方法,并且包括当前通常存储200-250个BG读数的SMBG设备使用的方法,以及将来生产的连续检测技术使用的方法。通过给出SMBG的这一广泛定义,本发明直接致力于通过引入一种能够同时预测HbA1c和低血糖症高危期的智能数据判读组件提高(但不仅限于)现有家用血糖监测设备的性能,以及通过相同的部件提高将来生产的连续监测设备的性能。
本发明的一个方面包括一种方法、***和计算机程序产品,用于由在预定时期内,例如大约4-6周,收集的SMBG数据估计HbA1c。在一个实施例中,本发明提供了一种计算机化的方法和***,用于根据在预定时期内收集的BG数据估计患者的HbA1c。本方法(或***或计算机可用介质)包括根据在第一预定时期内收集的BG数据估计患者的HbA1c。本方法包括:使用如下定义的预定序列数学公式准备用于估计HbA1c的数据:数据预处理;使用四个预定公式中的至少一个估计HbA1c;和通过样本选择标准验证估计的有效性。
本发明的另一个方面包括一种方法、***和计算机程序产品,用于估计严重低血糖症(SH)的长期概率。本方法使用预定时期内,例如4-6周,的SMBG读数,并预测随后大约6个月内的SH危险。在一个实施例中,本发明提供了一种计算机化的方法和***,用于根据在预定时期内收集的BG数据估计严重低血糖症(SH)的长期概率。本方法(或***或计算机可用介质)包括根据在预定持续时间内收集的BG数据估计患者严重低血糖症(SH)或者中度低血糖症(MH)的长期概率。本方法包括:根据所收集的BG数据计算LBGI;和根据计算出的LBGI利用预定的数学公式估计将来SH事件的数目。
本发明的再一个方面包括一种方法、***和计算机程序产品,用于识别24小时内(或者其他选择时期)的低血糖症高危期。其是通过使用在先前24小时收集的SMBG读数计算低血糖症的短期危险而实现的。在一个实施例中,本发明提供了一种计算机化的方法和***,用于根据在预定时期内收集的BG数据估计严重低血糖症(SH)的短期危险。本方法(或***或计算机可用介质)包括根据在预定持续时间内收集的BG数据估计患者严重低血糖症(SH)的短期概率。本方法包括:根据所收集的BG数据计算scale值;并为每个BG数据计算低BG危险值(RLO)。
本发明实施例的一个方面包括根据在第一预定持续时间内收集的BG数据估计患者HbA1c的方法(或者可替换的计算机程序)。该方法包括使用预定序列的数学公式准备用于估计HbA1c的数据。该数学公式定义为:数据预处理;通过样本选择标准验证BG数据样本;和如果样本有效则估计HbA1c。
本发明实施例的一个方面包括一种根据在第一预定持续时间内收集的BG数据估计患者HbA1c的***。该***包括可操作用于保持识别所述BG数据的数据库组件,和一个处理器,其中该处理器被编程用于使用预定序列的数学公式准备用于估计HbA1c的数据。该数学公式定义为:数据预处理;通过样本选择标准验证BG数据样本;和如果样本有效则估计HbA1c。
本发明实施例的一个方面包括一种根据在第一预定持续时间内收集的BG数据估计患者HbA1c的***。该***包括:BG获取机构,其用于从患者获取BG数据;可操作用于保持识别所述BG数据的数据库组件;和一个处理器。该处理器被编程用于使用预定序列的数学公式准备用于估计HbA1c的数据。该数学公式定义为:数据预处理;通过样本选择标准验证BG数据样本;和如果样本有效则估计HbA1c。
本发明实施例的一个方面包括一种不需要先前HbA1c的信息根据在第一预定持续时间内收集的BG数据估计患者HbA1c的方法(或者可替代的计算机程序)。该方法包括用预定序列的数学公式准备用于估计HbA1c的数据。该数学公式定义为:数据预处理;通过样本选择标准验证BG数据样本;和如果样本有效则估计HbA1c。
本发明实施例的一个方面包括一种不需要先前HbA1c的信息根据在第一预定持续时间内收集的BG数据估计患者HbA1c的***。该***包括可操作用于保持识别所述BG数据的数据库组件,和一个处理器。该处理器被编程用于使用预定序列的数学公式准备用于估计HbA1c的数据。该数学公式定义为:数据预处理;通过样本选择标准验证BG数据样本;和如果样本有效则估计HbA1c。
本发明实施例的一个方面包括一种不需要先前HbA1c的信息根据在第一预定持续时间内收集的BG数据估计患者HbA1c的***。该***包括:BG获取机构,其用于从患者获取BG数据;可操作用于保持识别所述BG数据的数据库组件;和一个处理器。该处理器被编程用于使用预定序列的数学公式准备用于估计HbA1c的数据。该数学公式定义为:数据预处理;通过样本选择标准验证BG数据样本;和如果样本有效则估计HbA1c。
本发明的这些方面,以及本文献全文讨论的其他方面,能够合并在一起从而提供有关糖尿病患者血糖控制的连续信息,并提高对低血糖症危险的监视。
本文中,本发明的这些及其他目标以及其优点和特点从下面的说明、附图和权利要求中将变得更加显而易见。
附图说明
通过联系附图一起阅读如下的优选实施例说明将对本发明的前述及其他目的、特点和优点以及本发明自身具有更全面的理解,其中:
图1图示地提供了在对15个由实例No.1的低BG指数定义的危险水平范围的每一个进行SMBG估计之后1个月内发生中度(虚线)和严重(实线)低血糖症的经验和理论概率。
图2图示地提供了在对15个由实例No.1的低BG指数定义的危险水平范围的每一个进行SMBG估计之后3个月内发生中度(虚线)和严重(实线)低血糖症的经验与理论概率。
图3图示地提供了在对15个由实例No.1的低BG指数定义的危险水平范围的每一个进行SMBG估计之后6个月内发生中度(虚线)和严重(实线)低血糖症的经验与理论概率。
图4图示地提供了在对15个由实例No.1的低BG指数定义的危险水平范围的每一个进行SMBG估计之后3个月内发生2次或更多次中度(虚线)和严重(实线)低血糖症的经验与理论概率。
图5图示地提供了在对15个由实例No.1的低BG指数定义的危险水平范围的每一个进行SMBG估计之后6个月内发生2次或更多次中度(虚线)和严重(实线)低血糖症的经验与理论概率。
图6是用于实现本发明的计算机***的功能框图。
图7-9是本发明相关处理器、通信连接和***的可替代变型的示意性框图。
图10图示地提供了在对15个由实例No.1的低BG指数定义的危险水平范围的每一个进行SMBG估计之后6个月内发生3次或更多次中度(虚线)和严重(实线)低血糖症的经验与理论概率。
图11图示地显示了本模型的残差分析,显示出与实例No.1训练数据组1残差的正态分布相近。
图12图示地显示了本模型的残差分析,显示出与实例No.1残差的正态分布相近。
图13图示地显示了由实例No.1的正态概率图1给出的统计证据。
图14图示地提供了实例No.1中用百分率表示的命中率和比值Rud之间的平滑依赖性。
图15图示地提供了实例No.1中预测时期与相应命中率之间的依赖性。
图16(A)-(B)图示地提供了由LBGI预测的T1DM显著低血糖症1个月内的危险,在实例No.2中各危险组严重低血糖症事件数目的ANOVA(F=7.2,p<0.001)和各危险组中度低血糖症事件数目的ANOVA(F=13.9,p<0.001)。
图17(A)-(B)图示地提供了由LBGI预测的T1DM显著低血糖症3个月内的危险,在实例No.2中各危险组严重低血糖症事件数目的ANOVA(F=9.2,p<0.001)和各危险组中度低血糖症事件数目的ANOVA(F=14.7,p<0.001)。
图18(A)-(B)图示地提供了由LBGI预测的T2DM显著低血糖症1个月内的危险,在实例No.2中各危险组严重低血糖症事件数目的ANOVA(F=6.0,p<0.001)和各危险组中度低血糖症事件数目的ANOVA(F=25.1,p<0.001)。
图19(A)-(B)图示地提供了由LBGI预测的T2DM显著低血糖症3个月内的危险,在实例No.2中各危险组严重低血糖症事件数目的ANOVA(F=5.3,p<0.01)和各危险组中度低血糖症事件数目的ANOVA(F=20.1,p<0.001)。
具体实施方式
本发明使得,但不仅限于,有可能产生用于估计糖尿病患者血糖控制的精确方法,并且包括在计算该方法的关键分量时使用的固件和软件编码。用于估计HbA1c、SH长期概率和低血糖短期危险的发明方法也能够根据所收集的大量数据加以验证,并且在本文的后面将进行讨论。最后,这些方法的方案能够被组合成结构化显示或者矩阵。
I.估计HbA1c
本发明的一个方面包括一种方法、***和计算机程序软件,用于由在预定时期内,例如4-6周,收集的SMBG数据估计HbA1c。在一个实施例中,本发明提供了一种计算机化(或者其他类型)的方法和***,用于根据在预定持续时间内收集的BG数据估计患者的HbA1c。本方法包括根据在第一预定持续时间内收集的BG数据估计患者的HbA1c,本方法包括:利用预定序列的数学公式准备用于估计HbA1c的数据。该数学公式被定义为:数据的预处理;使用四个预定公式中的至少一个估计HbA1c;和通过样本选择标准验证估计的有效性。笫一预定持续时间能够为大约60天,或者选择地,第一预定持续时间范围是大约45天到大约75天,或者大约45天到大约90天,或者根据期望。每个患者的数据预处理包括:将血浆BG转换成全血BG mg/dl;将以mg/dl为单位测量的BG转换成mmol/l的单位;和计算低血糖指数(RLO1)和高血糖指数(RHI1)。对每个患者的数据预处理使用如下定义的预定数学公式:通过BG=PLASBG(mg/d1)/1.12将血浆BG转换成全血BG mg/dl;通过BGMM=BG/18将以mg/dl测量的BG转换成mmol/l单位;和计算低血糖指数(RLO1)和高血糖指数(RHIl)。数据的预处理进一步使用如下定义的预定数学公式:Scale=[ln(BG)]1.0845-5.381,其中BG是以mg/dl为单位测量的;Risk1=22.765(Scale)2,其中RiskLO=Risk1,如果(BG小于大约112.5),因此存在LBGI的危险,否则RiskLO=0;RiskHI=Risk1,如果(BG大于大约112.5),因此存在HGBI的危险,否则RiskHI=0;BGMM1=每个患者的平均BGMM;RLO1=每个患者的平均RsikLO;RHI1=每个患者的平均RiskHI;L06=只对夜间读数计算的平均RiskLO,如果没有夜间读数则缺省;N06,N12,N24是各时间间隔中SMBG读数的百分率;NC1=第一预定持续时间内SMBG读数的总数;NDAYS=第一预定持续时间内具有SMBG读数的天数。N06、N12、N24分别如下时间间隔中SMBG读书的百分率,即大约0-6:59、大约7-12:59和大约18-23:59,或者其他的期望百分率和间隔数目。
本方法进一步包括根据用预定数学公式计算的患者高BG指数为组赋值。该公式可以定义为:如果(RHI1≤大约5.25或者如果RHI1≥大约16),则赋值group=0;如果(RHI1>大约5.25并且如果RHI1<大约7.0),则赋值group=1;如果(RHI1≥大约7.0并且如果RHI1<大约8.5),则赋值group=2;和如果(RHI1≥大约8.5并且如果RHI1<大约16),则赋值group=3。
接着,本方法可以进一步包括使用如下定义的预定数学公式给出估计:
E0=0.55555*BGMM1+2.95;E1=0.50567*BGMM1+0.074*L06+2.69;
E2=0.5555*BGMM1-0.074*L06+2.96;
E3=0.44000*BGMM1+0.035*L06+3.65;并且如果(group=1),则EST2=E1,或者如果(group=2)则EST2=E2,或者如果(group=3)则EST2=E3,否则EST2=E0。
本方法包括使用如下定义的预定数学公式对估计进行进一步的修正:如果(缺省(L06)),EST2=E0,如果(RLO1≤大约0.5并且RHI1≤大约2.0),则EST2=E0-0.25;如果(RLO1≤大约2.5并且RHI1>大约26),则EST2=E0-1.5*RLO1;并且如果((RLO1/RHI1)≤大约0.25并且L06>大约1.3)则EST2=EST2-0.08。
根据在第一预定持续时间内收集的BG数据估计患者的HbA1c能够通过使用四个预定数学公式中的至少一个估计HbA1c而实现,四个公式的定义如下:
a)HbA1c=如上面定义的或者如上面修正的EST2;
b)HbA1c=0.809098*BGMM1+0.064540*RLO1-0.151673*RHI1+1.873325,其中BGMM1是BG的平均值(mmol/l),RLO1是低BG指数,RHI1是高BG指数;
c)HbA1c=0.682742*HBA0+0.054377*RHI1+1.553277,其中HBA0是在估计之前的大约第二预定时期内采用的先前参考HbA1c读数,其中RHI1是高BG指数;或者
d)HbA1c=0.41046*BGMM+4.0775,其中BGMM1是BG的平均值(mmol/l)。第二预定持续时间能够为大约3个月;大约2.5个月到大约3.5个月;或者大约2.5个月到6个月;或者根据期望。
只有当第一预定持续时间样本满足如下四个标准中的至少一个时,才用HbA1c估计的样本选择标准验证估计的有效性:
a)测试频率标准,其中第一预定持续时间样本平均每天至少测试大
约1.5到大约2.5次;b)可选择测试频率标准,其中预定持续时间样本在第三预定采样时期内读数的平均频率为大约1.8个读数/天(或者其他期望的平均频率);
c)数据随机化标准-1,其中只有当比值(RLO1/RHI1)>=大约0.005时
才验证和显示HbA1c估计,其中RLO1是低BG指数,RHI1是高BG指数;或者
d)数据随机化标准,其中只有当比值(NO6>=大约3%)时才验证和显示HbA1c估计,且其中NO6是夜间读数的平均值。第三预定持续时间能够为至少35天,范围从大约35天到大约40天,或者从大约35天到大约和第一预定持续时间一样长,或者根据期望。
II.严重低血糖症(SH)的长期概率
本发明的另一个方面包括一种方法、***和计算机程序产品,用于估计严重低血糖症(SH)的长期概率。该方法使用预定时期,例如大约4-6周,的SMBG读数并预测在随后大约6个月内SH的危险。在一个实施例中,本发明提供了一种计算机化的(或者其他类型)方法和***,用于根据在预定持续时间内收集的BG数据估计患者严重低血糖症(SH)的长期概率。根据在预定持续时间内收集的BG数据估计患者严重低血糖症(SH)或者中度低血糖症(MH)的长期概率的方法包括:根据收集的BG数据计算LBGI;和根据计算的LBGI利用预定的数学公式估计将来SH事件的数目。LBGI的计算由在时间点t1,t2,...,tn采集一系列BG读数x1,x2,...,xn数学地定义:
我们定义了预定危险范围(risk category)(RCAT),借此每个危险范围(RCAT)表示LBGI的一个数值范围;并且将LBGI赋值给所述危险范围(RCAT)中的至少一个。危险范围(RCAT)的定义如下:
范围1,其中所述LBGI小于大约0.25;
范围2,其中所述LBGI介于大约0.25-大约0.50;
范围3,其中所述LBGI介于大约0.50-大约0.75;
范围4,其中所述LBGI介于大约0.75-大约1.0;
范围5,其中所述LBGI介于大约1.0-大约1.25;
范围6,其中所述LBGI介于大约1.25-大约1.50;
范围7,其中所述LBGI介于大约1.50-大约1.75;
范围8,其中所述LBGI介于大约1.75-大约2.0;
范围9,其中所述LBGI介于大约2.0-大约2.5;
范围10,其中所述LBGI介于大约2.5-大约3.0;
范围11,其中所述LBGI介于大约3.0-大约3.5;
范围12,其中所述LBGI介于大约3.5-大约4.25;
范围13,其中所述LBGI介于大约4.25-大约5.0;
范围14,其中所述LBGI介于大约5.0-大约6.5;和
范围15,其中所述LBGI大于大约6.5。
接着,分别为每个所述指定危险范围(RCAT)限定发生所选数目SH事件的概率。利用如下的公式分别为每个所述指定危险范围(RCAT)限定在下一个第一预定持续时间内发生所选数目SH事件的概率:F(x)=1-exp(-a.xb),x>0,否则为0,其中:a≈-4.19,b≈1.75(a和/或b可以是其他的期望值)。第一预定持续时间能够是大约1个月,范围从0.5个月到大约1.5个月,或者范围从大约0.5个月到大约3个月,或者根据期望。
此外,利用如下的公式分别为每个所述指定危险范围(RCAT)限定在下一个第二预定持续时间内发生所选数目SH事件的概率:F(x)=1-exp(-a.xb),x>0,否则为0,其中:a≈-3.28,b≈1.50(a和/或b可以是其他的期望值)。第二预定持续时间能够为大约3个月,范围从大约2个月到大约4个月,或者从大约3个月到大约6个月,或者根据期望。
进一步,利用如下的公式分别为每个所述指定危险范围(RCAT)限定在下一个第三预定持续时间内发生所选数目SH事件的概率:F(x)=1-exp(-a.xb),x>0,否则为0,其中:a≈-3.06,b≈1.45(a和/或b可以是其他的期值)。第三预定持续时间能够为大约6个月,范围从大约5个月到大约7个月,或者从大约3个月到大约9个月,或者根据期望。
选择地,利用如下的公式分别为每个所述指定危险范围(RCAT)限定在下一个第一预定时期内(范围是大约1个月,大约0.5-1.5个月,大约0.5-3个月,或者根据期望)发生所选数目MH事件的概率:F(x)=1-exp(-a.xb),x>0,否则为0,其中:a≈-1.58,b≈1.05(a和/或b可以是其他的期望值)。
选择地,利用如下的公式分别为每个所述指定危险范围(RCAT)限定在下一个第二预定持续时间内(范围是大约3个月,大约2-4个月,大约3-6个月,或者根据期望)发生所选数目MH事件的概率:F(x)=1-exp(-a.xb),x>0,否则为0,其中:a≈-1.37,b≈1.14(a和/或b可以是其他的期望值)。
选择地,利用如下的公式分别为每个所述指定危险范围(RCAT)限定在下一个第三预定持续时间内(范围是大约6个月,大约5-7个月,大约3-9个月,或者根据期望)发生所选数目MH事件的概率:F(x)=1-exp(-a.xb),x>0,否则为0,其中:a≈-1.37,b≈1.35(a和/或b可以是其他的期数值)。
而且,指定了患者在将来发生显著低血糖症的危险分类。该分类的定义如下:最小危险,其中所述LBGI小于大约1.25;低危险,其中所述LBGI为大约1.25-大约2.50;中度危险,其中所述LBGI为大约2.50-大约5之间;和高危险,其中所述LBGI大于大约5(也能够根据期望实现其他的分类范围)。
III.严重低血糖症(SH)短期概率
本发明的再一个方面包括一种用于识别24小时内(或者其他选择时期)低血糖症高危期的方法、***和计算机程序产品。这通过使用在先前24小时收集的SMBG读数计算低血糖症的短期危险而实现。在一个实施例中,本发明提供了一种计算机化方法和***,用于根据在预定持续时间内收集的BG数据估计患者严重低血糖症(SH)的短期危险。根据在预定持续时间内收集的BG数据估计患者严重低血糖症(SH)短期危险的方法包括:根据所述收集的BG数据计算scale值;和为每个BG数据计算低BG危险值(RLO)。RLO(BG)的计算被数学定义为:Scale=[ln(BG)]1.0845-5.381,其中BG是以mg/dl为单位测量的;Risk=22.765(Scale)2,如果(BG小于大约112.5),则:RLO(BG)=Risk,否则RLO(BG)=0。选择地,RLO(BG)的计算被数学地定义为:Scale=[ln(BG)]1.026-1.861,其中BG是以mmol/l为单位测量的;Risk=32.184(Scale)2,如果(BG小于大约112.5),则:RLO(BG)=Rsik,否则RLO(BG)=0。
根据所收集的BG数据能够计算LBGI。LBGI的计算由在时间点t1,t2,...,tn采集一系列BG读数x1,x2,...,xn数学地定义:
根据已收集的BG数据能够计算临时LBGI。临时LBGI的计算被数学地定义为:
LBGI(1)=RLO(x1);RLO2(1)=0;LBGI(j)=((j-1)/j)*LBGI(j-1)+(1/j)*RLO(xj);和RLO2(j)=((j-1)/j)*RLO2(j-1)+(1/j)*RLO(xj)-LBGI(j))2。
SBGI能够用如下定义的数学公式加以计算:
接着,本发明提供了对即将发生的短期SH的认证(qualification)和报警。在下述情况下将进行认证和报警,如果:(LBGI(150)≥2.5且(LBGI(50)≥(1.5*LBGI(150)且SBGI(50)≥SBGI(150)),则确认或者发出所述报警,或者RLO≥(LBGI(150)+1.5*(SBGI(150)),则确认或者发出所述报警,否则不需要认证或者提供报警。
接着选择地,本发明提供了对即将发生的短期SH的认证或报警。在下述情况下将进行认证和报警,如果:(LBGI(n)≥α且LBGI(n)ge(β)),则确认或者发出所述报警,和/或(RLO(n)≥(LBGI(n)+γ*SBGI(n))),则确认或者发出所述报警,否则不需要认证或者提供报警,其中α、β和γ是阈值参数。
阈值参数α、β和γ定义为α≈5,β≈7.5,γ≈1.5。在下面的表格中给出了其他可能的参数组合。该数值可以与下面给出的值相近,也可以是下面表格中数值的任意中间组合。
α | β | γ | α | β | γ |
6.4 | 8.2 | 1.5 | 5.0 | 7.5 | 1.3 |
6.0 | 7.5 | 1.5 | 4.9 | 7.0 | 1.2 |
5.5 | 7.5 | 1.5 | 4.9 | 7.0 | 1.2 |
IV.示例***
本发明的方法可以用硬件、软件或者其组合加以实现,并且能够在一个或多个计算机***或者其他处理***中实现,例如个人数字助理(PDA),或者直接在具有足够存储和处理能力的血糖自我监测设备(SMBG存储仪表)中实现。在一个示例实施例中,本发明是在如图6所示的通用计算机900上运行的软件。计算机***600包括一个或多个处理器,例如处理器604。处理器604连接通信基础设备606(例如,通信总线,交叉棒(cross-over bar)或者网络)。计算机***600可以包括显示接口602,其从通信基础设备606(或者从未显示的帧缓冲器)传送图表、文本和其他数据用于在显示单元630上进行显示。
计算机***600还包括主存储器608,优选地随机存取存储器(RAM),并还可以包括一个二级存储器610。二级存储器610可以包括,例如,硬盘驱动612和/或可移动存储驱动器614,其代表软盘驱动、磁带驱动、光盘驱动、闪存等。可移动存储驱动614以众所周知的方式从和/或向可移动存储单元618读和/或写。可移动存储单元618代表软盘、磁带、光盘等,其由可移动存储驱动614读出和写入。可以意识到,可移动存储单元618包括计算机可用的存储介质,其中存储有计算机软件和/或数据。
在可选择实施例中,二级存储器610可以包括其他的装置,用于允许计算机程序或者其他指令装载到计算机***600中。这类装置可以包括,例如,可移动存储单元622和接口620。这种可移动存储单元/接口的实例包括一个程序匣(cartridge)和匣接口(例如在视频游戏设备中见到的),可移动存储芯片(例如ROM、PROM、EPROM或者EEPROM)及相关插槽,和其他可移动存储单元622和接口,其允许插槽和数据从可移动存储单元62转移到计算机***600。
计算机***600还可以包括通信接口624。通信接口624允许软件和数据在计算机***600和外部设备之间进行传递。通信接口624的实例包括一个调制解调器,一个网络接口(例如以太网卡),一个通信端口(例如串联或者并联等),PCMCIA插槽和卡、一个调制解调器等。通过通信接口624传递的软件和数据呈信号628的形式,其可以是电子、电磁、光或者其他能够被通信接口624接收的信号。信号628通过通信路径(也就是通道)626提供给通信接口624。通道626携带信号628,并且可以用导线或者电缆、光纤、电话线、便携式电话连接、RF连接、红外连接和其他通信通道加以实现。
在本文献中,术语“计算机程序介质”和“计算机可用介质”用于一般地指如下的介质,例如可移动存储驱动614,安装在硬盘驱动612中的硬盘,和信号628。这些计算机程序产品是用于提供软件到计算机***600的装置。本发明包括这种计算机程序产品。
计算机程序(也称作计算机控制逻辑)被存储在主存储器608和/或二级存储器610中。计算机程序还可以通过通信接口624被接收。这种计算机程序在执行时能够使计算机***600执行如下文中讨论的本发明的特征。特别地,计算机程序在被执行时使处理器604能够执行本发明的功能。因此,这种计算机程序代表了计算机***600的控制器。
在用软件实现本发明的一个实施例中,软件可以存储在计算机程序产品中,并用可移动存储驱动器614、硬盘驱动612或者通信接口624装载到计算机***600中。控制逻辑(软件)在被处理器604执行时使处理器604执行如下文所述的本发明的功能。
在另一个实施例中,本发明主要以硬件形式实现,其使用例如硬件组件如设备专用集成电路(application specific integrated circuit)(ASIC)。执行硬件状态机器从而执行这里说明的功能,对于相关领域的技术人员而言是显而易见的。
在另一个实施例中,本发明使用硬件与软件的组合加以实现。
在本发明的一个示例软件实施例中,上述的方法用SPSS控制语言实现,但也能够用其他的程序实现,例如,但不仅限于,C++程序语言或者其他本领域技术人员能够获得的程序。
图7-9显示了代表本发明可选择实施例的框图。参考图7,显示了代表***710的框图,其主要包括由患者712使用的葡萄糖计728用于记录,尤其是,胰岛素剂量读数和测量的血糖(“BG”)水平。由葡萄糖计728获得的数据优选地通过合适的通信连接714或者数据调制解调器732转移到处理站或者芯片,例如个人计算机740、PDA,或者通过便携式电话或通过合适的互连网端口。例如,所存储的数据可以存储在葡萄糖计728中并通过合适的接口电缆直接下载到个人计算机,然后通过互连网发送到处理站。一个实例是由LifeScan有限公司生产的ONE TOUCH监测***或测量计,其与IN TOUCH软件兼容,并包括接口电缆以便下载到个人计算机。
血糖计是行业中通用的,并且主要包括任何能够起BG获取机构功能的设备。BG计或者获取机构、设备、工具或者***包括各种用于为每次测试提取血样(例如通过刺手指)的传统方法,并包括使用通过电机或者claorimetric方法读出葡萄糖浓度的装置确定葡萄糖的水平。最近,开发出了多种无需取血便可确定血液分析物浓度的方法。例如授予Yang等的美国专利No.5,267,152(本文引用作为参考),说明了一种利用近IR辐射漫反射激光光谱测量血糖浓度的非侵入技术。类似的近IR光谱设备在授予Rosenthal等的美国专利No.5,086,229,和授予Robinson等的美国专利No.4,975,581中有所说明(本文引用作为参考)。
授予Stanley的美国专利No.5,139,023(本文引用作为参考)说明了一种经皮血糖监测装置,其依靠渗透性增强器(例如胆汁盐)方便葡萄糖沿着在间隙液和接收介质之间建立的浓度梯度经皮运动。授予Sembrowich的美国专利No.5,036,861(本文引用作为参考)说明了一种被动葡萄糖监测器,其通过皮肤贴收集汗液,其中使用类胆碱能的制剂刺激汗腺分泌汗液。类似的汗液收集设备在授予Schoendorfer的美国专利No.5,076,273和授予Schroeder的美国专利No.5,140,985中有说明(本文引用作为参考)。
此外,授予Glikfeld的美国专利No.5,279,543(本文引用作为参考)说明了使用离子电渗疗法非侵入地将物质通过皮肤采样到皮肤表面的容器中。Glikfeld讲解了,该采样过程能够与葡萄糖特异的生物传感器或者葡萄糖特异的电极耦合从而监测血糖。而且,授予Tamada的国际公开No.WO 96/00110(本文引用作为参考)说明了一种用于经皮监测目标物质的离子电渗疗装置,其中离子电渗电极用于将分析物移动进入收集器中,并使用生物传感器探测容器中的目标分析物。最后,授予Berner的美国专利No.6,144,869(本文引用作为参考)说明了一种用于测量所存在分析物的浓度的采样***。
进一步,BG计或者获取结构还可以包括内在导管和皮下组织液采样。
计算机或者PDA 740包括软件和硬件,它们是根据预定的流程顺序(如上面详细说明的)处理、分析和判读自我记录的糖尿病患者数据并产生合适的数据判读输出所必需的。优选地,根据由计算机740存储的患者数据执行的数据分析和判读的结果通过与个人计算机740相连的打印机生成报表加以显示。选择地,数据判读程序的结果可以直接显示在与计算机740相连的视频显示单元上。
图8显示的框图代表了具有糖尿病管理***的可选择实施例,其是由患者操作的装置810,其外壳(housing)优选地足够紧凑从而使装置810能够用手持并由患者携带。在外壳816的表面上装载有一个用于接收血糖测试条(未显示)的条导向装置。测试条用于接收来自患者的血液样本。该装置包括一个微处理器822和与微处理器相连的存储器824。微处理器822被设计成执行存储在存储器824中的计算机程序从而执行各种计算和控制功能,如上面详细说明的。袖珍键盘816通过标准袖珍键盘解码器826与微处理器822相连。显示器814通过显示驱动器830与微处理器822相连。微处理器822通过接口与显示驱动器830通信,并且显示驱动器830在微处理器822的控制下修改和刷新显示器814。扬声器854和时钟856也连接微处理器822。扬声器854在微处理器822的控制下工作,从而发出可听的语音,使患者对将来可能的低血糖症警觉。时钟856向微处理器822提供当前的日期和时间。
存储器824也存储患者812的血糖值、胰岛素剂量值、胰岛素类型和微处理器822使用的参数值,用以计算将来血糖值、补充胰岛素剂量和碳水化合物补充。每个血糖值和胰岛素剂量值都与相应的日期和时间一起存储在存储器824中,存储器824优选地是非易失性存储器,例如电擦除只读存储器(EEPROM)。
装置810还包括一个与微处理器822相连的血糖计828。葡萄糖计828被设计用于测量被接收在血糖测试条上的血液样品,并产生血液样品测量的血糖值。如前所述,这种葡萄糖计在本领域是熟知的。葡萄糖计828的类型优选地是产生直接输出到微处理器822的数字值。选择地,血糖计828的类型可以是产生模拟值。在该选择实施例中,血糖计828通过模拟-数字转换器(未显示)与微处理器822相连。
装置810进一步包括一个输入/输出端口834,优选地,包括一系列连接微处理器822的端口。端口834通过接口,优选地,通过标准RS232接口连接调制解调器832。调制解调器832用于通过通信网络836在装置810与个人计算机840,或者康护人员的计算机838,之间建立通信。通过连接电缆连接电子设备的专用技术在本领域是熟知的。另一种可选择实例是“蓝牙”技术通信。
选择地,图9显示的框图代表了具有糖尿病管理***的可选择实施例,其是由患者操作的装置910,与图8所示类似,其外壳(housing)优选地足够紧凑从而使装置910能够用手持并由患者携带。例如,分离式或者可差卸式葡萄糖计或者BG获取机构/模块928。早已经有了自我监测设备,其能够直接计算算法1、2、3并向患者显示结果而无需将数据发送给其他事物。这种设备的实例是Lifescan有限公司生产的ULTRA SMART,Therasense,Alameda,CA生产的Milpitas,CA和FREESTYLE TRACKER。
因此,这里说明的实施例能够在数据通信网络例如互连网上实现,使得任何遥远位置的任何处理器和计算机都能够获得该评价、估计和信息,如图6-9和/或授予Wood的美国专利No.5,851,186中描述的,本文引用其内容作为参考。选择地,遥远位置处的患者可以将BG数据发送到中心康护人员或者医疗所或者不同的遥远位置。
总之,本发明提出了一种计算机化的(或者非计算机化的)数据分析方法和***,用于同时估计糖尿病个体血糖控制中的两个最重要的分量:HbA1c和低血糖症危险。本方法尽管只使用日常的SMBG数据,但是除了其他的之外,提供了三组输出。
本发明方法、***和计算机程序产品的潜力在于提供了如下的优点,但不仅限于此。首先,本发明通过执行和显示如下的内容提高了现有家用BG监测设备的性能:1)估计HbA1c的范围,2)估计在随后6个月内SH的概率,和3)估计低血糖症的短期危险(也就是,今后24小时)。后者可以包括警告,例如报警,表明即将发生低血糖症事件。这三个部分还能够组合在一起,从而提供有关糖尿病个体血糖控制的连续信息,进而提高对低血糖症危险的监测。
作为附加的优点,本发明提高了现有SMBG数据检索软件或者硬件的性能。几乎每个家用BG监测设备的制造商都生产这种软件或者硬件,并且患者和康护人员通常使用它判读SMBG数据。本发明的方法和***能够直接合并到现有家用血糖监测器中,或者通过引入能够同时预测HbA1c和低血糖症高危期的数据判读组件提高SMBG数据检索软件的性能。
另外一个优点是,本发明能够同时在低和高BG范围内以及在BG的整个数值范围内对家用BG监测设备的精确性进行评估。
而且,另一个优点,本发明能够评估各种糖尿病疗法的有效性。
再进一步,因为糖尿病患者终生面临着在保持严格血糖控制的同时而不使其低血糖症危险增加的优化问题,所以本发明通过使用其简单而可靠的方法缓和了这一相关问题,也就是,本发明能够同时评估患者的血糖控制及其低血糖症的危险,同时能够在患者的日常条件下加以应用。
另外,本发明通过提出三种截然不同但相互兼容的算法给出了缺失的联系,该三种算法均用于由SMBG数据估计HbA1c和低血糖症的危险,从而用于预测低血糖症的短期和长期危险以及高血糖症的长期危险。
另一个优点,本发明能够评估新型胰岛素或者胰岛素投放(delivery)设备的有效性。任何胰岛素或胰岛素投放设备的制造商或者研究人员都能够利用本发明的实施例测试其提出或者检测的胰岛素类型或者设备投放设计的相对成功性。
最后,另一个优点,本发明能够评估胰岛素辅助治疗药物的有效性。
发明实例
I.实例No.1
实例No.1包括三种算法,用于由日常SMBG数据同时估计糖尿病血糖控制中的两个最重要的分量,HbA1c和低血糖症危险。该方法直接适合于通过引入能够同时预测HbA1c和低血糖症高危期的智能数据判读组件提高现有家用BG监测设备的性能。该数据分析方法具有三个部分(算法):
●算法1:估计HbA1c
●算法2:估计严重低血糖症(SH)的长期危险,和
●算法3:估计低血糖症的短期(24-48小时内)危险。
算法1和2提供了有关1型或者2型糖尿病(T1DM,T2DM)个体整体血糖控制的不间断监测和信息,同时覆盖BG数值范围的上限和下限。当算法2指示低血糖症的长期危险增加时,算法3被激活。一旦激活,算法3要求更加频繁的监测(每天4次),并提供对中度/严重低血糖症的24-48小时预报。
实例1的另一个重要目标是用现有数据检验大量的假说和观点,这些假说和观点有可能产生其他其他的算法,它们以概念上不同于本发明公开所提供的方式估计HbA1c和计算低血糖症的危险。目标是找到潜在的更好的解决方法,或者简单地证实某些观点不能够产生更好的结果,其主要用于优化和改进对目前在本研究实例No.2中收集的数据的分析。
数据组(data sets)
为了保证我们的最优化结果能够推广到更广泛的水平,算法1和2首先用训练数据组加以优化,然后用不相关的测试数据组检验其精度。对于算法3,我们目前只有一组含有平行SMBG和SH记录的数据。患者人口调查(population)的详细说明如下:
(1)训练数据组1:96名患者至少在本研究2年之前被诊断患有TIDM。其中有43名患者报告在过去一年里有过至少2次严重低血糖症事件,53名患者报告在相同的时间里没有这种事件。其中有38名男性和58名妇女。平均年龄为35±8岁,糖尿病的平均持续时间为16±10年,每天的平均胰岛素为0.58±0.19单位/kg,平均HbA1c为8.6±1.8%。这些对象在40-45天的时间内收集了大约13,000个SMBG读数。SMBG的频率为大约3个读数/天。这些数据收集继续进行了6个月,按月记录中度和严重低血糖症事件。该数据组用作算法1(无先前HbA1c)和算法2的训练数据组。
(2)训练数据组:85名患者至少在本研究2年之前被诊断患有TIDM,所有患者均报告在过去一年里有SH事件。其中有44名男性和41名妇女。平均年龄为44±10岁,糖尿病的平均持续时间为26±11年,每天的平均胰岛素为0.6±0.2单位/kg,平均基准HbA1c为7.7±1.1%,且平均6个月HbA1c为7.4±1%(60名对象具有6个月HbA1c)。这些对象在两次HbA1c化验之间的6个月内收集了大约75,500个SMBG读数。数据组2中的SMBG频率更高,为4-5个读数/天。此外,在SMBG的6个月期间,这些对象保持记录中度和严重低血糖症事件及其发生的日期和时间,结果有399次SH事件。该数据组用作算法1(有先前HbA1c)和所有与算法3有关的分析的训练数据组。
(3)测试数据组:我们使用N=600个对象的数据,其中277名具有T1DM,323名具有T2DM,他们全部使用胰岛素治疗糖尿病。这些数据由美国中部圣地亚哥Amylin公司的药剂师收集,包括6-8个月的SMBG数据(大约300,000个读数),同时具有基准和6个月HbA1c结果以及一些人口统计数据。这些对象参加了pramlintide(剂量为60-120微克)对代谢控制的临床实验调查。T1DM和T2DM组中的对象使用pramlintide是随机的(表1)。
表1.测试数据组中对象的人口统计特征
变量 | T1DM平均值(SD) | T1DM平均值(SD) | P水平 |
年龄(岁) | 38.0(13.4) | 58.1(9.4) | <0.001 |
性别:男/女 | 136/41 | 157/166 | Ns |
基准HbA1c | 9.74(1.3) | 9.85(1.3) | Ns |
6个月HbA1c | 8.77(1.1) | 8.98(1.3) | 0.04 |
糖尿病持续时间 | 14.6(9.8) | 13.5(7.6) | Ns |
开始年龄(岁) | 23.4(12.8) | 44.6(10.4) | <0.001 |
#SMBG读数/对象/天 | 3.2(1.1) | 2.9(0.9) | <0.005 |
表1给出了人口统计特征以及T1DM与T2DM对象的比较。在研究最初的6个月,T1DM和T2DM组的平均HbA1c都显著下降,这可能是由于药物的使用,这超出了本报告书(presentation)(表1)的范围。HbA1c相对快速的变化允许更好地评估算法1的预测能力。在所有数据组中,SMBG用Lifescan公司的ONE TOUCH II或ONETOUCH PROFILE测量计执行。
算法1:估计HbA1c
实例No.1给出了,但不仅限于,HbA1c预测(算法1)的优化,其通过:(1)SMBG越接近中心,权重越高;(2)高BG事件时间越长,权重越高;(3)用较早前的HbA1c校正高BG指数;和(4)合并其他的患者变量,例如年龄、性别和疾病持续时间。
算法1包括SMBG数据的优化函数,其估计随后的HbA1c,并推荐数据收集周期的最佳持续时间和该周期内自我监测的最佳频率。然而要非常注意,算法1的更广泛目的是评估患者血糖控制的状态。尽管HbA1c被接受作为评估血糖控制的“黄金标准”,但是目前并不清楚其他的测量量,例如平均SMBG或者高BG指数是不是比HbA1c更好的糖尿病长期并发症的预报因子。在该问题被澄清之前,算法1的目的都是估计HbA1c。为了尽可能地接近将来算法1的实际应用,我们如下进行:
(1)首先,由我们在先前对T1DM患者的研究中收集的两个训练数据组1和2推导出多个使用不同独立变量的优化函数,最佳持续时间和最佳SMBG频率;
(2)然后,固定所有的系数,将算法1应用于大得多的测试数据组,该测试数据组同时包括在由Amylin药剂师指导的临床实验中在非常不同的条件下收集的T1DM和T2DM对象的数据;
(3)仅用测试数据组详细评估算法1对各种优化函数的精度。
训练和测试数据组分离允许我们声明,算法1的估计精度能够一般化到T1DM或者T2DM对象任何其他的数据。而且,因为Amylin数据(测试数据组)是从正接受治疗以降低其HbA1c的对象中收集的,因此他们的HbA1c在6个月的观察周期内显示出了不正常的大变异,我们能够声明,算法1不仅能够预测相对恒定的HbA1c,而且能够预测大而不正常快速变化的HbA1c。在相同的行中,算法1对于希望优化其HbA1c的患者是最有用的,可以推测,该患者群体很可能对于寻求具有先进特征的测量计,例如连续估计HbA1c,感兴趣。
结果摘要
●最佳SMBG数据收集周期为45天;
●最佳SMBG频率是3次读数/天;
●提出了两个最佳HbA1c估计函数:F1-只使用SMBG数据,和F2-使用SMBG数据加上在所预测HbA1c之前大约6个月获得的HbA1c读数;
●通过在下页(表2)中详细说明的多个标准对测试数据组(N=573个对象)中HbA1c预测的精度进行评估。我们在这里说明,在T1DM中,F1的总精度(HbA1c测量值的20%之内)为96.5%,而F2的总精度为95.7%。对于T2DM,F1的总精度为95.9%,F2的总精度为98.4%。因此,F1和F2的精度均与HbA1c的直接测量结果相当;
●最重要的,对于HbA1c从其基准读数改变2个或者更多个单位的患者(N=68),F1无论对于T1DM还是T2DM,预测该变化的精度均为100%,而F2对于T1DM和T2DM的精度分别为71%和85%;
●F1和F2对于6个月HbA1c的预测都比0月时的最初HbA1c估计要精确得多。但使用平均BG作为HbA1c的直接估计并不准确;
●试验了大量可选择的方法,例如选择每天的特殊时间点(餐后读数)估计HbA1c,根据每次SMBG读数和HbA1c确定之间间隔的时间给SMBG读数不同的权重,对具有不同平均血糖/HbA1c比值的对象分别进行评估等。尽管这些可选择方法中有一些获得了比上面提出的两个函数好一些的结果,但在整体上没有一个更好。我们能够得出结论:在算法1的将来应用中将使用优化函数F1和F2。
详细结果-测试数据组
算法1评估最重要的部分是评估它对于如下数据的性能,即与用于提出并优化它的数据不相关的数据。由测试数据组,573个对象的数据,其中T1DM的N=254,T2DM的N=319,足以用于评估算法1。
最佳算法1:对于每个对象,选择了他/她SMBG读数45天的子集。该子集的开始日期大约在对象6个月HbA1c化验之前的75天,结束日期大约为该化验前的30天。因为在该数据组中,HbA1c化验的时间只是大概地知道,所以在上一次用于分析的SMBG读数与HbA1c之间间隔的时间并不确切。该时间周期是通过连续优化其持续时间和其结束点(HbA1c之前的时间)选择出来的。最佳持续时间为45天。最佳结束时间为HbA1c之前1个月。换言之,45天的SMBG能够提前大约1个月预测HbA1c的数值。然而,预测45-75天之间任何一个HbA1c数值都几乎一样好——差别只是数字而非临床的显著度。类似地,45天的监测周期与60天的监测周期之间的差异并不大。但是,小于45天的监测周期会使得预测能力快速下降。
最佳估计函数是线性的,并且由下给出:
估计1——不知道先前的HbA1c:
F1=0.809098*BGMM1+0.064540*LBGI1-0.151673*RHI1+1.873325
估计2——已知先前的HbA1c(大约6个月之前)
F2=0.682742*HBA0+0.054377*RHI1+1.553277
在这些公式中,BGMM1是由45天的SMBG读数计算出的平均血糖;LBGI1和RHI1是由相同读数计算出的低和高BG指数,HBA0是只用于估计2的基准HbA1c读数。系数值用训练数据组进行优化,并且在详细结果——训练数据组部分中给出了相关的统计和绘图。
函数F1和F2产生了HbA1c的点估计,也就是说,每个函数产生了一个HbA1c估计。利用详细结果——训练数据组部分中提供的回归误差估计能够获得区间估计。然而,对于测试数据组,这些区间估计并不是HbA1c真正的90%或者95%置信区间,因为它们最初得自于训练数据组并且只是应用于测试数据(见下一部分的统计注意)。
算法1精度评估:表2A和2B给出了分别用T1DM和T2DM对象的测试数据组数据评估最佳算法1的结果。使用了多个标准:
(1)从HbA1c测量值估计的绝对偏差(AERR);
(2)从HbA1c测量值估计的绝对百分比偏差(PERR);
(3)HbA1c测量值20%内的百分比估计(HIT20);
(4)HbA1c测量值10%内的百分比读数(HIT10);和
(5)HbA1c测量值周围25%区间之外的百分比读数(MISS 25)。
表2A:算法1在T1DM中的精度(N=254个对象)
F1 | F2 | 平均BG | 先前HbA1c | P值 | |
AERR | 0.77 | 0.61 | 1.68 | 1.1 | <0.001 |
PERR(%) | 8.3 | 7.1 | 19.4 | 12.8 | <0.001 |
HIT20(%) | 96.5 | 95.7 | 61.0 | 81.0 | <0.001 |
HIT10(10) | 65.4 | 75.5 | 29.9 | 48.2 | <0.001 |
MISS 25(%) | 2.4 | 1.6 | 28.4 | 9.9 |
表2B.算法1在T2DM中的精度(N=319个对象)
F1 | F2 | 平均BG | 先前HbA1c | P值 | |
AERR | 0.72 | 0.57 | 1.92 | 0.87 | <0.001 |
PERR(%) | 7.6 | 6.4 | 20.9 | 11.7 | <0.001 |
HIT20(%) | 95.9 | 98.4 | 56.4 | 82.8 | <0.001 |
HIT10(10) | 70.2 | 79.3 | 29.5 | 53.3 | <0.001 |
MISS 25(%) | 1.2 | 0.6 | 36.7 | 8.2 | <0.001 |
表2A和2B的前两列分别给出了优化函数F1和F2的结果。第三行给出了当采用平均BG(mmol/L)作为HbA1c估计时估计的精度。第四行给出了用0时的HbA1c化验作为6个月HbA1c估计计算的相同精度测量结果。显然,对于T1DM和T2DM,F2对于估计HbA1c在整体上都比F1略好。最重要的是,F1和F2对于HbA1c的估计都比用其较早数值或者平均BG进行估计要好得多。对于落在25%精度区间之外的%估计,这特别地真。执行F1和F2与由先前HbA1c化验进行估计的差异是高度显著的(第4列)。
统计注意:要重点注意到,使用传统的回归型标准,例如从ANOVA表中获得的R2或F和p值,评估算法1的精度是不合适的。这是因为,参数估计得自于另一个不相关的数据组(训练数据),而只是应用于该测试数据组。因此,基本(underlying)模型的统计假定被违反了(例如在测试数据组中,残差的总和不是零),因此R2或F和p值失去了它们的统计学意义。
通过检查SMBG基准读数与随后6个月的读数有显著变化的T1DM和T2DM对象,对算法1在测试数据组中的精度进行了进一步的评估。表3A和3B给出了HbA1c的绝对改变等于或者大于2个单位的T1DM和T2DM对象的列表。在每个对象组中,有34个对象的HbA1c具有这种变化。算法1,函数F1,在T1DM和T2DM中都100%地预测出了这种改变。由于公式中含有基准HbA1c(其部分地将估计值拉回到HbA1c的基准值),所以F2的预测力变小,在T1DM中为71%,在T2DM中为85%。除两个对象之外,基准HbA1c落在6个月HbA1c值20%区间的外部:
表3A:HbA1c变化>=2单位的T1DM对象
IDHBA0 | HBA0 | HBA6 | DHBA | F1 | F2 | HIT F1 | HIT F2 | HIT |
6504661340036204370947013614360260083723701062086202392482116012391367012307351658082201401062104904670966193921660374156515361137327409 | 12.010.512.411.013.012.811.911.511.313.012.711.510.69.911.09.311.011.210.611.89.611.812.411.911.310.39.510.911.010.69.810.313.210.0 | 7.06.88.97.59.79.58.78.38.310.19.88.77.87.28.36.78.48.68.19.37.29.510.19.69.18.17.38.88.98.57.88.311.28.0 | 5.003.703.503.503.303.303.203.203.002.902.902.802.802.702.702.602.602.602.502.502.402.302.302.302.202.202.202.102.102.102.002.002.002.00 | 6.828.028.457.298.999.508.247.939.308.808.098.427.917.718.767.827.888.757.957.767.618.908.578.338.637.837.647.208.187.947.138.369.307.99 | 9.909.3710.739.4511.5411.6110.309.9410.5311.4610.8910.099.378.7210.328.359.5410.079.2710.038.5210.7111.1510.4010.299.048.579.199.899.278.549.2311.999.04 | 100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00 | .00.00.00.00100.00.00100.00100.00.00100.00100.00100.00.00.00.00.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00 | .00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00100.00.00 |
表3B:HbA1c变化>=2单位的T1DM对象
IDHBA0 | HBA0 | HBA6 | DHBA | F1 | F2 | HIT F1 | HIT F2 | HIT |
675463616270626463553961655580525356396690865542353406463517551635839659141603170837613768326109150139644352785842564752655665628255 | 10.811.312.011.111.810.811.111.79.710.39.510.711.111.310.112.28.410.111.110.210.812.411.210.39.311.913.712.212.110.610.111.17.910.9 | 7.07.68.67.88.68.08.38.97.07.76.98.18.78.97.79.86.17.88.87.98.610.29.08.27.29.811.610.110.08.68.19.15.98.9 | 3.803.703.403.303.202.802.802.802.702.602.602.602.402.402.402.402.302.302.302.302.202.202.202.102.102.102.102.102.102.002.002.002.002.00 | 6.908.517.858.317.999.138.117.686.758.087.478.168.997.897.929.177.007.839.578.027.629.138.297.457.708.5210.089.519.538.768.518.727.078.90 | 9.0310.2010.039.709.909.739.559.808.209.078.239.429.909.888.6311.027.328.6410.338.889.2410.869.748.788.1810.1812.6511.1411.019.698.879.687.049.87 | 100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00 | .00.00100.00.00100.00.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00100.00 | .00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00100.00.00.00.00.00.00.00.00 |
在表3A和3B中:
ID——对象的ID数字;
HBA0——基准HbA1c;
HBA6——HbA1c的6个月测量值;
DHBA——HbA1c基准值与6个月值之间的绝对差值;
F1——由函数F1估计的HbA1c,只使用SMBG数据;
F2——由函数F2估计的HbA1c,使用先前的HbA1c化验结果;
HitF1=100,如果F1在6个月HbA1c值的20%之内,否则为0;
HitF2=100,如果F2在6个月HbA1c值的20%之内,否则为0;和
Hit HbA0=100,如果基准HbA1c在6个月HbA1c读数的20%之内,否则为0。
详细结果——训练数据组
这一部分说明了优化算法1的步骤。该优化包括两个部分:(1)假定不能获得先前的HbA1c读数,和(2)假定能够使用先前的HbA1c预测HbA1c。
我们考虑了多个不同的函数用于描述SMBG数据与HbA1c的关系。根据计算的精度和简单性,如果不使用先前的HbA1c读数,最佳函数似乎是SMBG读数的平均值、低和高BG指数的线性函数,另一个是先前HbA1c与高B指数的线性函数。非线性关系并没有提高该模型的匹配良好度,因此不考虑用于实际应用。
训练数据组1-无先前的HbA1c使用线性回归模型优化函数F1的系数。最佳系数在前面的部分中已经给出。这里,我们给出有关模型匹配良好度的数据:
倍数R(Multiple R).71461
R方 .51067
方差分析
DF 平方和 均方
回归 3 154.57097 51.52366
残差 90 148.10903 1.64566
F=31.30889 显著度F=.0000
该模型的残差分析显示出接近于残差的正态分布(见图11)。残差的SD为1.2(根据定义其平均值为0)。因此,我们能够接受,该模型良好地描述该数据。
训练数据组2-先前的HbA1c再一次,使用线性回归模型优化函数F3的系数。最佳系数在前面的部分中已经给出。这里,我们给出有关模型匹配良好度的数据:
倍数R .86907
R方 .75528
方差分析
DF 平方和 均方
回归 4 38.70237 9.67559
残差 54 12.54000 .23222
F=41.66522 显著度F=.0000
该模型的残差分析显示出接近残差的正态分布(见图12)。残差的SD为0.47(根据定义其平均值为0)。因此,我们能够接受,该模型良好地描述了该数据。
此外,比较没有和有先前HbA1c的模型,我们能够得出结论,如果在计算中包括先前的HbA1c,则最终的模型比根据R2和根据残差误差都要好得多。
然而,正如我们在先前部分中见到的,在不相关数据组中,先前的HbA1c对于预测的总精度没有贡献,在某些情况下,HbA1c变化显著甚至由于变化快速而妨碍了算法的能力。因此,我们能够得出结论,即使先前的HbA1c可能从统计的观点看更好,但是不可能有足够的实际效用,用于校正将来的测量计的读数输入。我们还不知道HbA1c化验与SMBG检测(profile)之间间隔的时间,但仍然使HbA1c的输入有用。或许,这取决于该时间周期中HbA1c的改变——正如我们在前一部分中看到的,2个HbA1c单位的变化就使先前的HbA1c读数完全无用。
SMBG/HbA1c比值
我们现在提供了一种可选择的方法,以改进模型匹配的统计精度并保持相当高的临床应用性。可见,45天SMBG读数的平均值与HbA1c的比值是一个具有最完全正态分布的测量量(由Kolmogorov-Simov实验可以证实),并且最重要地,区分了三组对象,其比值为<1.0、1.0-1.2和>1.2。前两组的每一个都能解释大约40%的对象,第三组能够解释大约20%的对象。这对于T1DM和T2DM都是有效的,并且在训练数据组以及测试数据组中都能观察到。此外,该比值对于时间似乎特别稳定,因此可能是反映患者SMBG***均值将低估HbA1c,从而相应的比值将小于1.0)。注意这只是一个假说,并不能够用可获得的数据加以证实,我们进行了一些分析,似乎证实能够在一定程度上知道每个人在某个时间点的比例。这似乎相当于知道先前的HbA1c,并且有可能等价于数据输入,但是该比值的应用与先前HbA1c的使用非常不同。除了直接用于该预测公式之外,该比值还用于将患者进行分类,对它们应用三个不同预测公式中的某一个。这些新的公式不直接包括HbA1c,因此不会受包含HbA1c导致的惯性(inertia)的影响。另外,由该比值限定的三个组之间平均HbA1c差异并不大,并且与该比值不相关,所以不同人中比例不同的原因一定与HbA1c不相关。
如果我们首先根据他们的比值将对象分成3组并在训练数据组中分别进行回归,则回归模型的配合良好度显著增加:(1)在组1(比值<1.0)中,我们得到倍数R=0.86,R2=0.73;(2)在组2(比值为1.0-1.2)中,匹配度最佳,R=0.97,R2=0.94,和(3)在组3(比值>1.2)中,匹配度最差,R=0.69,R2=0.47。因为所有三种回归模型都不包括先前的HbA1c,所以我们得出结论,对于大约80%的对象,匹配良好度显著增加,剩余20%的对象,匹配良好度保持相同,匹配度将变差的对象能够预先被识别出来。
进一步,根据对象的比值将测试数据组分成3组。我们获得的预测精度与先前获得的精度相似(表4A和4B):
表4A:算法1在T1DM(N=254个对象)中的精度
比值<1.0 | 1.0≤比值≤1.2 | 比值>1.2 | |
AERR | 0.70 | 0.63 | 0.74 |
PERR(%) | 7.8 | 7.4 | 7.9 |
HIT 20(%) | 93.8 | 93.0 | 95.5 |
HIT 10(%) | 68.8 | 73.4 | 72.7 |
MISS 25(%) | 3.1 | 2.6 | 0.0 |
表4B:算法1在T2DM(N=319个对象)中的精度
比值<1.0 | 1.0≤比值≤1.2 | 比值>1.2 | |
AERR | 0.63 | 0.68 | 0.89 |
PERR(%) | 7.6 | 7.8 | 8.8 |
HIT 20(%) | 97.4 | 95.0 | 95.3 |
HIT 10(%) | 67.2 | 65.3 | 57.7 |
MISS 25(%) | 0.0 | 1.7 | 0.0 |
简而言之,知道每个对象的SMBG/HBA1c比值并相应地分别进行估计,似乎能够提高该模型的统计性能而不损失临床精度。
被测试的其他假说和思路
我们测试了大量其他的假说和思路,它们被证实至少对于促进和更集中地分析由实例No.2收集的数据是有用的。简要的结果说明如下:
(1)HbA1c与下午获取的SMBG读数——中午12点至下午6点——最关联(相关),与空腹SMBG读数(4a.m.-8a.m.)最少关联。然而,并不是只采集餐后SMBG读数就会提高HbA1c的预测,相反,如果忽略全天所有小时(相对较小但重要)的贡献,则预测会变得更糟。如果给全天其它小时的读数进行不同的加权,则有可能提高HbA1c的预测,但是该提高并不足以弥补模型附加的复杂性;
(2)在T2DM中,HbA1c与平均SMBG的关系比在T1DM中更强,即使两组的HbA1c相互匹配。根据直接相关性,在整个研究中,T1DM中的系数为大约0.6,而T2DM中的系数大约为0.75;
(3)根据SMBG与HbA1c化验之间间隔的时间对SMBG读数进行不同的加权(例如结果越接近中间,加权越高)并不能获得对HbA1c更好的预测;
(4)包含人口统计变量,例如年龄、糖尿病持续时间、性别等,并不能提高HbA1c的预测;
(5)HbA1c与平均SMBG(以mmol/l为单位测量)之间可能最简单的线性关系由如下公式给出:HbA1c=0.41046*BGMM+4.0775。尽管与F1和F2相比在统计上较差,但是该公式给出的HbA1c估计在T1DM和T2DM中具有大约95%的精度(根据与HbA1c化验的偏差小于20%),并且如果通过计算低和高BG指数将其合并在一个测量计中,则可能是有用的(但是,不计算低BG指数是不能够实现低血糖症的预测的,因此该公式可能只对包括算法1但不包括算法2和3的测量计有用)。
算法2:SH长期危险的评估
实例No.1提供了,但不仅限于,算法2的扩展,包括估计个体生化显著低血糖症(biochemical significant hyoglycemia)(BSH,定义为BG读数<=39mg/dl)或者生化中度低血糖症(BMH,定义为39mg/dl<BG<=55mg/dl)的概率。另外,我们计划评估,与日间SH相比,算法2预测夜间(午夜至7:00am)SH的发生率是否更好。
算法2是一个分类算法。也就是说,根据对象的SMBG数据,它将对象的将来BSH或者MSH分类成特定的危险范围。为了尽可能接近算法2将来真实的应用,我们如下地进行:
(4)首先,从训练数据组1得出多个最佳分类变量和最佳分类范围、最佳持续时间和最佳SMBG频率;
(5)然后,将测试数据组分成两个部分:最初45天和其余的数据部分。将算法2的最佳参数应用于数据的最初45天部分,从而使用将来BSH或者MSH的概率估计值预测数据第二部分的BSH和MSH;
(6)只用测试数据对算法2精度进行详细评估。
训练与测试数据组分离允许我们声明,算法2的估计精度能够一般化到T1DM或者T2DM患者任何其他的数据。而且,因为Amylin数据是从经过强化治疗的对象中收集的,所以我们能够推测,算法2在低血糖症危险不断变化且不断增加的对象中进行试验并且证实是有效的。
结果摘要
●估计将来BSH或者BMH概率所需的最佳SMBG数据收集周期为40-45天。SMBG的最佳频率为3-4个读数/天。大量的读数并不能够导致算法2预测能力的显著增加。如果每天的读数小于3个,则预测能力下降。然而,该要求是参考45天观察周期内每天的平均读数,并非意味着,每天都需要执行3-4次读数;
●预测子变量与将来SH和MH之间的关系严格地非线性。因此,线性方法不能用于优化预测,尽管通过直接线性模型能够获得R2=50%(相比之下,DCCT的最佳结果为预测8%的将来SH);
●单独预测夜间SH一般比预测日间SH更弱;
●确定了15个将来BSH和BMH的危险范围。范围的最佳分离只是根据低BG指数获得的,尽管低BG指数与其它变量的组合能够发挥类似良好的作用;
●尽管BSH和BMH的频率在T1DM和T2DM中不同(见表5),但是当给定一个危险范围时,条件频率在T1DM和T2DM之间并无差异。这允许获得统一的预测SH和MH危险的方法;
●为15个危险范围计算了将来的各经验概率并进行了比较。所有的比较都具有高的显著度,p’s<0.0005。
●这些经验概率用双参数Weibull分布(two-parameter Weibulldistribution)加以近似,产生了每个危险范围中将来BSH和BMH的理论概率。
●这些近似的匹配良好度非常好——所有的确定系数都大于85%,一些高达98%(见图1-5和9-10)。
详细结果——测试数据组
确定SH/MH的个人危险范围:总共600个对象的数据用于该分析。由他/她前45天的SMBG数据收集为每个对象计算低BG指数(LBGI)。然后,将LBGI分类到15个最佳危险范围中的其中一个(变量RCAT的范围是0-14),如在训练数据组1中得出的。这些危险范围通过如下的不等式定义:
if(LBGI≤0.25),RCAT=0
if(0.25<LBGI≤0.5),RCAT=1
if(0.50<LBGI≤0.75),RCAT=2
if(0.75<LBGI≤1.00),RCAT=3
if(1.00<LBGI≤1.25),RCAT=4
if(1.25<LBGI≤1.50),RCAT=5
if(1.50<LBGI≤1.75),RCAT=6
if(1.75<LBGI≤2.00),RCAT=7
if(2.00<LBGI≤2.50),RCAT=8
if(3.00<LBGI≤3.50),RCAT=9
if(3.50<LBGI≤4.00),RCAT=10
if(4.00<LBGI≤4.50),RCAT=11
if(4.50<LBGI≤5.25),RCAT=12
if(5.25<LBGI≤6.50),RCAT=13
if(LBGI>6.50),RCAT=14
BSH和BMH的观察频率:对于每个对象,在最初45天数据收集之后的1个月、3个月和6个月内计数由SMBG指示的BSH和BMH的发生。表5A给出了在T1DM中观察到0、>=1、>=2、>=3次BSH和BMH的频率,表5B给出了在T2DM中观察到的相同数据:
表5A:在T1DM中观察到的BSH和BMH频率
BSH(BG≤39mg/dl) | BMH(39mg/dl<BG≤55mg/dl) | |||||
1个月 | 3个月 | 6个月 | 1个月 | 3个月 | 6个月 | |
平均#/对象 | 0.82 | 1.77 | 2.74 | 3.64 | 8.33 | 12.93 |
0次事件的%Ss | 62.8 | 50.8 | 46.6 | 25.2 | 18.0 | 17.7 |
≥2次事件的%Ss | 18.8 | 33.1 | 38 | 64.3 | 75.6 | 77.1 |
≥3次事件的%Ss | 9.8 | 23.3 | 28.2 | 50.8 | 68.0 | 71.1 |
表5B:在T2DM中观察到的BSH和BMH频率
BSH(BG≤39mg/dl) | BMH(39mg/dl<BG≤55mg/dl) | |||||
1个月 | 3个月 | 6个月 | 1个月 | 3个月 | 6个月 | |
平均#/对象 | 0.18 | 0.53 | 0.76 | 1.11 | 2.93 | 4.59 |
0次事件的%Ss | 91.4 | 84.9 | 81.3 | 73.0 | 61.1 | 55.8 |
≥2次事件的%Ss | 3.6 | 8.6 | 10.1 | 18.1 | 26.7 | 30.3 |
≥3次事件的%Ss | 1.5 | 5.9 | 7.4 | 13.9 | 21.4 | 25.8 |
夜间BSH和BMH占由SMBG指示的全部事件的大约15%。在训练数据组中,夜间事件与全部预测子变量之间的相关性较弱。我们得出结论,夜间事件的目标预测是无效的。
将来BSH和BMH的经验概率:我们计算了15个危险范围中每一个的将来BSH和BMH的特定经验概率。这些概率包括:(1)在之后1个月、3个月和6个月内发生至少一次BSH或者BMH的概率;(2)在之后3个月和6个月内发生至少二次BSH或者BMH的概率;和(3)在之后6个月内发生至少三次BSH或者BMH的概率。当然,有可能根据要求计算任何其他的组合概率。
从该分析得出的最重要的结论是,给定一个危险范围,将来BSH和BMH的概率在T1DM和T2DM之间没有显著差异。这允许获得统一的经验性和理论性预测T1DM和T2DM中这些概率的方法。结果,T1DM和T2DM患者的数据被组合用于如下的分析。
图1-5和9-10给出了按照15个危险范围绘制的6个计算经验概率的散点图。BSH的经验概率用黑三角表示,而BMH的经验概率用红方形表示。
用单变量ANOVA在15个危险范围内对所有组的经验概率进行比较,所有的p水平都小于0.0005。因此,我们观察到,在不同的危险范围内BSH和BMH事件之间的差异高度显著。
将来BSH和BMH的理论概率:为了能够使用直接的公式估计将来BSH和BMH的概率,我们用双参数Weibull概率分布对经验概率进行了近似。Weibull概率函数由如下公式给出:
F(x)=1-exp(-a,xb),x>0;否则F(x)=0
统计注意:参数a和b大于0,并且分别称作尺度参数和形状参数。在b=1的特殊实例中,Weibull分布变成指数型。该分布经常在工程问题中使用,因为随机发生技术失败的分布彼此并不完全无关(如果失败完全无关,则它们将形成Piosson处理,这将在指数分布中加以说明,例如b=1)。这里的情况很不相同——我们需要描述不完全独立并且趋向于群集发生的事件(失效)的分布,例如我们先前的研究所证实的。
每组经验概率都用上面给出的理论公式加以近似。用非线性最小平方法对参数进行估计(用线性双对数模型给出的初始参数估计)。每种模型的匹配良好度用其确定系数(coefficient of determination)(D2)进行评估。该统计的意义与线性回归中的R2相似,但是R2不能应用于非线性模型。
图1-6中给出了模型匹配,黑线用于BSH的概率,虚线用于BMH的概率。在每个图上,我们给出了对相应模型的参数估计,因此我们给出了直接的公式,用于计算初始SMBG之后1个月、3个月和6个月内发生0、>=1、>=2、>=3次BSH和BMH的频率。在监测设备或者软件中能够包括其中一些公式或其变型作为SH和MH危险的指示子。
每个图的下方给出了D2值,作为近似精度的指示子。所有的数值都大于85%,一些达到了98%,这证实该近似度非常好,并且证实在将来的研究/应用中能够使用理论概率,而不是经验概率。
发生一次或多次中度或严重低血糖症事件的理论概率由图1所示的公式给出:
P(MH>=1)=1-exp(-exp(-1.5839)*Risk**1.0483)
P(SH>=1)=1-exp(-exp(-4.1947)*Risk**1.7472)
图1给出了在15个由低BG指数定义的危险水平的每一个内在SMBG估计之后1个月内中度(虚线)和严重(黑线)低血糖症的经验和理论概率。因为模型是非线性的,所以用它们的确定系数D2估计匹配良好度,D2是线性模型中R2的类似物。确定系数及其平方根如下:
SH模型:D2=96%,D=98%
MH模型:D2=87%,D=93%
发生一次或多次中度或严重低血糖症事件的理论概率由图2所示的公式给出:
P(MH>=1)=1-exp(-exp(-1.3731)*Risk**1.1351)
P(SH>=1)=1-exp(-exp(-3.2802)*Risk**1.5050)
图2给出了在15个由低BG指数定义的危险水平的每一个内在SMBG估计之后3个月内中度(虚线)和严重(黑线)低血糖症的经验和理论概率。
确定系数及其平方根如下:
SH模型:D2=93%,D=97%
MH模型:D2=87%,D=93%
发生一次或多次中度或严重低血糖症事件的理论概率由图3所示的公式给出:
P(MH>=1)=1-exp(-exp(-1.3721)*Risk**1.3511)
P(SH>=1)=1-exp(-exp(-3.0591)*Risk**1.4549)
图3给出了在15个由低BG指数定义的危险水平的每一个内在SMBG估计之后6个月内中度(虚线)和严重(黑线)低血糖症的经验和理论概率。
确定系数及其平方根如下:
SH模型:D2=86%,D=93%
MH模型:D2=89%,D=95%
发生两次或多次中度或严重低血糖症事件的理论概率由图4所示的公式给出:
P(MH>=2)=1-exp(-exp(-1.6209)*Risk**1.0515)
P(SH>=2)=1-exp(-exp(-4.6862)*Risk**1.8580)
图4给出了在15个由低BG指数定义的危险水平的每一个内在SMBG估计之后3个月内发生2次或多次中度(虚线)和严重(黑线)低血糖症的经验和理论概率。
确定系数及其平方根如下:
SH模型:D2=98%,D=99%
MH模型:D2=90%,D=95%
发生两次或多次中度或严重低血糖症事件的理论概率由图5所示的公式给出:
P(MH>=2)=1-exp(-exp(-1.7081)*Risk**1.19555)
P(SH>=2)=1-exp(-exp(-4.5241)*Risk**1.9402)
图5给出了在15个由低BG指数定义的危险水平的每一个内在SMBG估计之后6个月内发生2次或多次中度(虚线)和严重(黑线)低血糖症的经验和理论概率。
确定系数及其平方根如下:
SH模型:D2=98%,D=99%
MH模型:D2=89%,D=95%
发生三次或多次中度或严重低血糖症事件的理论概率由图9所示的公式给出:
P(MH>=3)=1-exp(-exp(-2.0222)*Risk**1.2091)
P(SH>=3)=1-exp(-exp(-5.5777)*Risk**2.2467)
图10给出了在15个由低BG指数定义的危险水平的每一个内在SMBG估计之后6个月内发生3次或多次中度(虚线)和严重(黑线)低血糖症的经验和理论概率。
确定系数及其平方根如下:
SH模型:D2=97%,D=99%
MH模型:D2=90%,D=95%。
详细结果——训练数据组
训练数据组包括SMBG数据及之后按月记录的严重低血糖症记录。与用截止(cutoff)BG数值确定BSH和BMH的测试数据组相对,按月的记录所包括的严重症状的报告被定义为由低血糖症导致的无意识、昏迷、无法自我治疗或者显著认知损伤。在SMBG之后的6个月内,研究对象报告称每个人有平均2.24次该事件,其中67%的对象报告没有该事件。从统计的观点看,这只会使SH事件的分布显著偏颇,并且不适合于应用线性方法。尽管线性回归能够用于估计各种变量的相对分布从而预测SH,但是并不能用于构建最终的模型。我们进行了如下的三种分析:
(1)不知道SH历史:忽略SH的任何历史知识,我们由基准HbA1c和SMBG特征,例如平均BG、低BG指数和变化的估计BG危险率,通过回归预测将来的SH(所有的变量在本发明最初的公开中已经说明)。和以前发现的一样,HbA1c和平均BG对预测SH没有任何贡献。最终的回归模型包括低BG指数和BG危险率的变化,并具有如下的匹配良好度:
倍数R .61548
S平方 .37882
方差分析
F=27.74772 显著度F=.0000
——等式中的变量——
变量 | B | SE B | β | T | 显著度T |
LBGI | 4.173259 | .649189 | 2.104085 | 6.428 | .0000 |
比值 | -5.749637 | 1.091007 | -1.724931 | -5.270 | .0000 |
(常数) | -2.032859 | .790491 | -2.572 | .0117 |
(1)知道先前的SH:当我们包括了前一年SH事件的数目时,如筛选调查表中报告的,这一变量能够解释另外11%的将来SH变异:
倍数R .70328
S平方 .49461
方差分析
F=29.35999 显著度F=.0000
——等式中的变量——
变量 | B | SE B | β | T | 显著度T |
LBGI | .337323 | .704286 | .375299 | 4.541 | .0000 |
LDR | -4.350779 | 1.036380 | -1.305264 | -4.198 | .0001 |
RLO | 3.134519 | .631684 | 1.580371 | 4.962 | .0000 |
(常数) | -2.136619 | .717334 | -2.979 | .0037 |
(3)不知道先前SH的数目,只知道某人先前是有还是没有过SH,我们只用SMBG变量能够解释45%的将来SH变异;
(4)最后,两个单独的线性模型能够解释55%的日间SH变异和25%的夜间SH变异。所有预测子变量与夜间SH的直接相关性也微弱。夜间事件占全部SH的30%。
我们得出结论,线性预测模型能够直接解释大约40-50%的将来SH变异。然而,根据其残余误差,该模型并未良好平衡(这是由于在糖尿病患者人口统计中SH事件数目的高度偏颇分布)。图13的正态概率图给出了其统计证据,该图显示出了标准残差与其预测值的显著偏差。
因此,我们采用了另一种预测SH的方法,利用他们的SMBG数据将对象分成各种危险范围,并估计这些范围中随后SH的概率。我们尝试了各种分类模型,使危险范围之间的差异最大并试图获得最大的危险估计分辨率(根据范围的最大数目)。
通过只用低BG指数进行分类得到15个危险范围(在先一部分的开始部分给出)获得了最佳的结果。
除了范围之间最佳分离之外,该结果还具有其他的优点:(1)不需要先前的SH历史知识;(2)计算相对简单且不需要跟踪时间变量,例如BG的变化速度,和(3)分类似乎能够同等地应用于T1DM和T2DM患者(不需要任何先前SH的知识)。
算法3:低血糖短期危险的评估
实例No.1给出了,但不仅限于,算法3的优化:
(1)利用基准长期危险(来自算法2)和HbA1c(来自算法1);
(2)低血糖症报警的危险标准/阈值;
(3)SMBG频率;
(4)如果检测到低血糖症危险增加并且已经有一段时间没有SMBG,则是否应当发出低血糖症报警,和
(5)人口统计学变量的贡献,例如严重低血糖症历史。
引言
与具有长期发展历史的算法1和算法2不同,算法3处理的问题直到最近都被认为是不可能的。实际上,人们仍然普遍认为根据先前已知的数值不可能预测将来的BG值(特别是低血糖症)(Bremer T和Gough DA.能够从先前的数值预测血糖吗?数据引发.(Is bloodglucose predictable from previous values?A solicitation for data)Diabetes,1999,
48:445-451)。我们先前的工作,其在能够从Lifescan有限公司获得的一个原稿中有过报道并且在本发明公开中被详细地提出,质疑了这个普遍的看法。为了解释该质疑的基础以及阐明算法3背后的原理,我们包含了下面的段落。
我们量化糖尿病特征的“哲学”:
根据所研究的内分泌***,激素间的相互作用是由动态调控生化网络控制的,该网络具有由主要节点和管道组成的复杂或者简单的结构。糖尿病能够在各种水平上扰乱调控胰岛素-葡萄糖动力学的网络。例如,在T1DM中,胰岛素的自然产生被完全消除,而在T2DM中,胰岛素在细胞内的使用受到更强大的胰岛素抵抗的阻碍。在T1DM中(在T2DM中也常见),需要某种形式的外部胰岛素替代品,这使得调控***对于有缺陷的外部因素变得脆弱,这些外部因素包括丸药和胰岛素注射的时限和剂量,所吃的食物,体育活动等。这通常会导致BG极度偏差出现低血糖症和高血糖症。在许多但不是所有病例中,低血糖症激发内分泌反应,称之为反向调节。因此,在数学领域,BG在一段时间内的波动是受大量内部和外部因素影响的复杂动力学***活动的可测量结果。然而,根据众所周知的动力学***理论,当调控的复杂性增加时,一个纯粹的确定性***会发展成显示随机的宏观行为。因此,在短时期内(分钟),在人的水平上观察到的BG波动将是接近确定的,然而从长期看,波动将是接近随机的,这包括极端过渡,例如SH事件。因此,随机建模和统计推论最适于在较长的时期内分析***——算法1和算法2采用的范例,其在特定的观察周期之后利用我们最初提出的手段,例如LBGI和HBGI,预测数值范围和一个事件的概率。BG在短期内的波动能够用确定性网络加以建模和预测,这将在能够进行连续检测的将来智能胰岛素投递设备中实现。
算法3在数小时到数天的中间时间范围内工作,所以要求结合统计推论和确定性建模。前者用于估计个体SH的基准危险,而后者用于动态跟踪个体参数和在SH事件出现之前进行预报。当在一个装置中执行时,算法3将如下地工作:
(1)装置收集研究对象的某些基准信息,建立个人基准参数;
(2)然后,装置开始追踪SMBG数据的某组特征;
(3)装置具有判定规则,确定何时建立即将发生SH的标志,以及当数据表明危险减弱时何时降低该标志;
(4)当建立了标志时,我们认为对象在随后24小时内(预测时间)将接受到SH报警。
这种动态预测同时在模型参数优化水平上和最佳方法的精度评估水平上产生了理论上的问题。我们开始阐述第二个问题,因为它对于理解算法3的作用是最重要的。
算法3精度评估:虽然算法1和算法2使用静态预测,并且评估这些算法的标准在理论上是显然的——预测值更好,但是对于算法3,优化的标准不再是直接。这是因为在提高被预测SH事件的百分率的同时,我们不可避免地增加了“已建立标志”的数目,这反过来增加了潜在的“错误报警”的数目。因为“错误报警”并没有被清晰地定义,所以使问题进一步复杂化。在其纯粹意义上,错误报警是指建立了标志但随后并未发生SH事件。然而,如果人察觉到了症状并且采取了合适的行动,则SH是可以避免的。因此,即使可能存在SH的生化潜力,但事件可能并不出现。为了解决这个问题,我们采取了如下的优化标准:
(1)使24小时内将发生SH的预测最大化;
(2)“标记升高”和“标记降低”的持续时期的比值Rud最小化。
尽管这两点的第一点是清晰的,但是第二点可能需要额外的解释。从在测量计中执行算法3的前景来看,在每次确定SMBG时,测量计都确定是否建立或者不建立即将发生SH的标记。当标记被建立时,它可能会持续一段时间(伴随着随后的数个SMBG读数)直到做出降低标记的决定为止。因此,我们将具有交替的“标志升高”和“标志降低”过程,其在SMBG点发生变化。参考上面点(2)的比值Rud是一个人标志升高时的平均时间除以标志降低时的平均时间。
在发明公开中给出的我们先前的最好结果是预测了44%的24小时内SH事件,且Rud=1∶7,即1天高危险警报与7天没有警报相交替。因为当时我们假定报警周期至少为24小时,所以算法被优化为建立标志的频率不超过每周1次。假如该分析是用具有高比率SH事件的对象的数据实现的,则这一比值是可以接受的。
在本研究的实例No.1中,我们不得不使用相同的数据组完善算法3,因为没有其他可用的同时包括SMBG记录和SH记录的数据。我们也使用了相似的标准评估算法3的精度。然而,我们基本上改变了其他所有的事情。数据的跟踪、参数估计、全部阈值和判定规则都不再一样。这些变化是由于一个新的思路导致的,即在SH之前机体对反向调节的储备(reserve)会有一定的“损耗”,并且该损耗可以用SMBG数据进行跟踪。这一思路的确切执行在“判定规则”部分中有描述。因为判定规则包括一个连续的标准和一个在某种程度上的人为中止,所以存在多个解,我们选择了最佳的一个用于进一步的研究。然而,根据这些结果的表述,我们可以决定选择另一个解,从而在算法3的将来应用中加以执行。
结果摘要
首先,要重点注意到下面给出的全部结果恰好超过统计的显著性。正如我们将在下一部分的几个实例中见到的,所观察到的差异一直非常显著(p值低于任何可能的显著水平)。算法3的观点是在个人的基础上预测SH事件的发生。结果是:
(1)最小的基准观察周期是以每天3-4个读数的频率在大约两周的时间内获得50个SMBG读数。之后将每个对象分类到两个危险组中的一个,后者使用不同的判定规则;
(2)从我们具有的六个月数据中我们发现,在观察开始时足以进行该组的分配。因此,我们假定大约每六个月测量计将使用50个读数重新评价它自己的组分配;
(3)SMBG跟踪的最佳延迟是以每天3-4个读数的频率采集100-150个读数。换言之,最佳的判定标准是根据使用测量计存储器中全部150个读数进行的计算。这通过模拟ONE TOUCH ULTRA的存储容量而实现。一般地,使用在一周内获得的只有20个读数的延迟就能够获得良好的结果,但更长的延迟能够产生更好的预测;
(4)判定规则是根据新的计算程序,其使用“临时平均值”计算跟踪对象的低BG指数和其他相关参数。我们设计了专用软件用以执行该程序并处理我们能够得到的数据。从编程的观点看,执行该程序所需的编码只有大约20行,包括LBGI的计算;
(5)我们调查了多个判定规则(使用各种参数)。忽略SMBG频率,这些规则获得的24小时内SH预测从43.4%,Rud=1∶25到53.4%,Rud=1∶7。因此,与我们先前的结果相比,24小时内SH的预测提高了10%;
(6)作为进一步研究的最佳解,我们选择的判定规则可预测50%的24小时内SH,且Rud=1∶10。下面的结果参考了该最佳解在不同条件下的结果:
(7)SMBG的最佳频率是每天4个读数。如果达到了该频率,则24小时内SH的预测会增加到57.2%,并具有相同的Rud=1∶10。其他的SMBG频率也有调查和报道;
(8)如果我们将预测周期延长到36或48小时,则SH的预测分别增加到57%和63%,且具有同样的Rud=1∶10;
(9)利用基准信息能够显著增加SH的预测。实际上,比算法3先前的版本增加10%完全是由于使用了基准跟踪。但是,该基准跟踪现在被建模作为不使用患者任何额外输入的测量计在两周时期内的自我校正;
(10)个人从口统计信息,例如SH的历史或先前的HbA1c,对于更好的SH短期预测没有贡献;
(11)任何时候当长时间没有SMBG活动时建立标志是不恰当的。只有测量计发出即将发生SH警报的次数才是使用次数。这是因为SH预测的主要部分是根据非常低BGs的重现(群集)。对该重现的估计在为2002年6月ADA会议准备的摘要中给出(Kovatchev等.具有大量SH历史的T1DM患者中周期性发生的低血糖症和严重低血糖症(SH))(见附录)。
数据处理的详细说明
测量计将SMBG读数与每个读数的日期和准确时间(小时,分,秒)存储在一起。因此,在训练数据组中,我们具有每个对象一定时间顺序的SMBG记录。在研究期间,从参与者的存储测量计中总共下载了75,495个SMBG读数(每人每天平均4.0±1.5个)。由对象的月记中,我们获得了所发生SH事件的日期和时间。对象报告了399次SH事件(每人4.7±6.0次)。有68位参与者(占80%)曾经历一次或者多次SH事件。根据他们的人口统计学特征,这些对象与那些没有经历过SH的人(剩余20%的对象)没有不同。
数据预处理:
我们开发了专用的软件用于数据预处理。这包括:(1)将每个对象的存储测量计数据与BG读数的连续6-8个月序列组合在一起,和(2)按照日期和时间使每位对象的SH记录与该序列相匹配。后者的执行如下:为每个SMBG读数计算到最近一次SH事件的时间(小时/分)和从上一次SH事件经历的时间。因此有可能:(1)每次SH事件之前或之后的24小时或48小时等时间周期,(2)SMBG读数之间的时间周期。由于SH(昏迷,无意识)的本质,在SH时不会有SMBG,因此出于算法3的目的,SH事件不包括用于算法2的生化显著低血糖症。平均每个SH事件在SH与之前最近SMBG读数之间的最小时间间隔为5.2±4.1小时;有29次SH事件(占7%)在之前15分钟内具有SMBG读数。对每个SH事件,我们计数了在该事件之前24h、36h、48h和72h内执行SMBG读数的次数。
基准危险值的计算和自我校正:
根据他/她的第一SMBG读数为每个对象计算低BG指数。已经确定,计算基准LBGI所需的最少读数是在大约2周内采集50个。因此对于每个新的测量计,我们需要期望一段最初两周的自我校正时期,在此期间测量计将扫描其拥有者的SH总体危险。初始期之后,人们被分到两个危险组中的一个:低中度危险(LBGI≤3.5,LM组)或者中高度危险(LBGI>3.5,MH组)。我们的测试数据显示,更精确的分组是不必要的。这种分组允许在LM和MH组中使用不同的判定规则,从而与本发明公开中给出的其最初命中率相比,算法的命中率提高了大约10%。
使用测试数据,不需要重新校正基准危险。因此我们能够假定,如果人在治疗中没有经历变化,则大约每六个月进行一次重新校正。这与算法2的结果是一致的,其显示在初始观察周期之后的6个月,SH的长期预测是非常准确的。
然而,如果人在其血糖控制中经历了迅速的变化,则可能需要更加频繁地重新校正。重新校正的判定可能是自动的,并且根据所观察到的运行危险值(见下段)与基准LBGI之间的差异不断增加。然而,可获得的数据并不能允许我们阐明这个问题,因为我们所观察对象的低血糖症危险并没有显著变化。
计算SMBG参数:在预处理步骤之后,我们设计了另一个软件来计算用于预测即将发生的SH的SMBG参数。该软件包括:
(1)为每个BG读数计算低BG危险值(RLO),这通过如下的编码实现(这里BG以mg/dl为单位测量,如果单位是mmol/l则系数会不同):
scale=(ln(bg))**1.08405-5.381
risk=22.765*scale*scale
if(bg_1≤112.5)then
RLO=risk
else
RLO=0
endif
(2)为每个具有序列数n的SMBG读数,BG(n),计算LBGI(n)的运行值,和另一个统计量,SBGI(n),其是低BG危险值的标准偏差。这两个参数用每个SMBG读数之后的某个延迟(k)加以计算,例如包括该读数,BG(n),和BG(n)之前采集的(k-1)个读数。
(3)LBGI(n)和SBGI(n)的计算采用一个新临时方法程序(newprovisional means procedure),其基于如下的递归式编码:
初值在n-k(或者精确地在最大值(1,n-k),以便解释序数小于k的测量计读数):
LBGI(n-k)=rlo(n-k)
rlo2(n-k)=0
对于n-k与n之间任何连续迭代j的数值:
LBGI(j)=((j-1)/j)*LBGI(j-1)+(1/j)*RLO(j)
rlo2(j)=(j-1)/j)*rlo2(j-1)+(1/j)*(RLO(j)-LBGI(j))**2
在该循环完成之后,我们得到了LBGI(n)的值,并计算
SBGI(n)=sqrt(rlo2(n))
因为对于ONE TOUCH ULTRA测量计,最大值n是150,所以在k=10到k=150范围内寻找最佳的延迟k。虽然性能上的差异并不显著,但是最佳的延迟确定为k=150(见例如下一部分)。
判定规则:在每次SMBG读数时,程序判定是否建立警告即将发生SH的标志。如果已经建立了标志,则程序确定是否使其降低。这些判定是根据三个阈值参数:α、β、γ,其运行如下:
对于对低中度危险(LM组)的对象:
FLAG=0
if(LBGI(n)≥αand SBGI(n)≥β)FLAG=1
if(RLO(n)≥(LBGI(n)+γ*SBGI(n)))FLAG=1
对于中高度危险组的对象,只有第二if叙述起作用。换言之,如果LBGI(n)的运行值和它的标准偏差SBGI(n)都超过了某个阈值,则建立标志(即等于1),而且如果低BG危险RLO(n)的当前值超过LBGI(n)加γ标准偏差的值则标志也建立。
一个启发性的解释:LBGI(n)和SBGI(n)的数值反映了低血糖症危险的更缓慢变化——要经过数天的SMBG才能显著改变这些数值。因为LBGI(n)平均值越高,近期低血糖症越频繁和极端,所以我们可得出结论,LBGI(n)和SBGI(n)反映了反向调节储备在数天的时期内持续的损耗(或者缺乏补给)。此外,SBGI(n)是***稳定性的一个标志——较大的SBGI(n)表明对象的BG波动增加,因此控制***变得不稳定和容易受到极端异常的影响。因此,第一逻辑表达式反映的观点是,无论何时当反向调节防御被耗尽且控制(外部或内部)变得不稳定时就会发生SH。第二逻辑表达式说明了低BG危险的急剧改变,无论何时当目前低BG危险值突然大于其运行平均值时则会发出标志。实际上对于中高度危险组的对象,只有第二逻辑表达式与这些对象的最终“永久耗尽”和“永久不稳定”状态相对一致。因为这些对象连续运行低BG值,并且他们的BG是不稳定的,所以任何急剧的低血糖事件都可能激发SH。大体上,严重低血糖症标志或者低不稳定BG周期之后建立,或者在急性低血糖事件之后建立,急性低血糖事件严重偏离了(在危险空间内)最近的运行危险平均值(可能早经很高)。这些之前没有任何报警信号的SH事件仍然不能用本算法加以解释。在下面的表5C中,我们给出了解释算法3对多个对象的作用的样本输出:
表5C:解释算法3对多个对象的作用的样本输出:
ID135135135135135135135135…135135135135135 | BG70771245150664997130597641 | SH.00.00.00.00.00.00.001.00.00.00.00.00.00 | FLAG.00.00.001.001.00<24h1.00.00.00.00.001.00 | TIME53.7541.0935.0230.4414.7210.608.30140.0525.1720.205.23.62 | 这里α=5,β=7.5,γ=1.5对于135号对象,第一次标志在SH之前30小时建立下一次读数,其在16小时之后采集并且在SH之前14小时。后面的读数及之后的两个读数均在SH之前的24小时内,因此我们认为预测到了该事件。无论如何,该标志在SH之前大约8小时20分时被再次建立。该对象的第二个SH事件被提前36分钟标记。 |
135219219219219219219219219222222222222222222223223223223223223 | 200644322538437515617683664222814941110 | 1.00.00.00.00.00.00.00.001.00.00.00.00.00.001.00.00.00.001.00.001.00 | .00.001.00.001.00<24h<24h.00.00.00.001.00.00.001.00<24H | 40.7237.8828.7316.2211.1810.874.5219.0813.239.727.834.7518.8014.155.856.00 | 该对象在此次SH事件之前大约28.7小时和11.2小时获得两次警告。该对象在此次SH事件之前大约4小时45分获得警告。该对象在12小时内经历了两次重现的SH事件。标志在第一个事件之前大约6小时建立,因此我们认为两个事件均在24小时的高危期内被预测到。 |
该输出的每一行都给出了一个SMBG读数,或一个SH事件(没有读数)。ID是对象的ID号码,BG是以mg/dl计的BG水平,当发生SH事件时SH=1。如果算法3决定建立flag则FLAG=1;TIME是以小时计的最近一次SH事件的时间。
临时方法程序延迟的优化:在先前的出版物中,我们报道了在SH之前的48到24小时周期内,平均BG水平降低而BG方差增加。在SH之前的24小时周期内,平均值BG水平进一步降低,BG的方差继续增加,LBGI急剧增加。在SH之后的24小时周期内,平均BG水平正常化,而BG方差仍大幅增加。在SH之后的48小时内,BG的平均值和方差返回到基准水平(见Kovatchev et,al.1型糖尿病中严重低血糖症事件之前和之后48小时内可测得的血糖扰动(Episodes ofSevere Hypoglycemia in Type 1 Diabetes are Preceded,and Followed,within 48 Hours by Measurable Disturbances in Blood Glucose).
J of Clinical Endocrinology and Metabolism,85:4287-4292,2000)。我们现在使用这些观察值根据在SH之前24小时内观察到的LBGI(n)和SBGI(n)的平均值优化算法3所采用的临时方法程序的延迟,k。简言之,用于计算LBGI(n)和SBGI(n))的延迟被选择使SH之前24小时内的这些测量显示与其余的研究相比最大化,但不包括SH之后***仍处于不平衡的时期。发现最佳延迟是k=150。表6A和6B给出了LBGI(n)和SBGI(n)对于参数k的多个值和对于低中度危险组和中高度危险组两个对象组的平均值。显然,k的各个值之间的差异不大,因此实际使用任何k≥10的数值都是合适的。然而,根据目前的数据,我们推荐k=150,并且所有进一步的计算使用该延迟。该推荐还是根据LBGI(n)和SBGI(n)的方差岁延迟值的增大而减小,这能够由下面更大的t-值加以反映:
表6A:不同延迟下SH之前24小时内与其余时间的LBGI(n)
LBGI | 低中危险(LM组) | 中高危险(MH组) | ||||||
SH之前24h | 其余时间 | t | p | SH之前24h | 其余时间 | t | P | |
K=10 | 4.99 | 3.32 | 9.2 | <.000... | 6.94 | 5.28 | 11.0 | <.000... |
K=20 | 4.73 | 3.34 | 9.4 | <.000... | 6.56 | 5.32 | 11.0 | <.000... |
K=30 | 4.54 | 3.35 | 9.1 | <.000... | 6.50 | 5.34 | 11.3 | <.000... |
K=50 | 4.53 | 3.34 | 9.8 | <.000... | 6.54 | 5.36 | 11.7 | <.000... |
K=100 | 4.45 | 3.29 | 11.2 | <.000... | 6.53 | 5.40 | 12.3 | <.000... |
K=150* | 4.46 | 3.26 | 12.1 | <.000... | 6.56 | 5.41 | 12.9 | <.000... |
表6B:不同延迟下SH之前24小时内与其余时间的SBGI(n)
LBGI | 低中危险(LM组) | 中高危险(MH组) | ||||||
SH之前24h | 其余时间 | t | p | SH之前24h | 其余时间 | t | P | |
K=10 | 7.50 | 5.35 | 9.0 | <.000... | 9.64 | 7.78 | 11.4 | <.000... |
K=20 | 8.31 | 6.28 | 10.1 | <.000... | 10.36 | 8.91 | 10.7 | <.000... |
K=30 | 8.50 | 6.68 | 9.5 | <.000... | 10.82 | 9.38 | 11.5 | <.000... |
K=50 | 8.78 | 7.04 | 10.0 | <.000... | 11.30 | 9.80 | 12.3 | <.000... |
K=100 | 9.19 | 7.33 | 11.7 | <.000... | 11.69 | 10.18 | 13.5 | <.000... |
K=150* | 9.42 | 7.42 | 13.7 | <.000... | 11.88 | 10.33 | 15.1 | <.000... |
*最佳解
从表6A和6B可以看出,LBGI和SBGI在SH之前的24小时内都高度显著地增加。因此,人们试图运行一种直接的判别式或者对数模型来预测即将发生的SH。不幸的是,这种标准统计的效果并不很好,尽管两种模型在统计上都非常显著。判别式模型(其效果比对数回归要好)能够正确预测52.6%的即将发生的SH事件。然而,其标志升高(flag-up)与标志降低(flag-down)的比非常差——Rud=1∶4。因此,该模型倾向于更大数量的数据点,这是在任何统计程序中都会预见到的倾向。因此,我们不得不必须采用上面给出的判定规则。
严重低血糖症的预测精度
阈值参数α、β和γ的优化 下面我们详细说明了使用阈值参数α、β和γ的各种组合的算法3的预测能力。因为这些参数与所期望结果(高度预测SH且比值Rud最小)之间的关系非常复杂,因此我们所采用的优化程序不能够获得单一的解。此外,似乎也不需要获得单一的解。它似乎是一个商业的而非数学的判定,即在给定的“标志升高”与“标志降低”的比值之下,获得可以接受的SH预测百分率。因此,我们并不声明下面给出的任何解是最优解。但是,为了进一步探索这个主题,我们接受能预测50%的将来SH且Rud=1∶10作为基准,用于研究除24小时之外的预测周期以及为了获得更好的危险分布(profile)所需的每天SMBG读数个数的各种要求。
表7给出了算法3在α、β和γ的各种数值组合下的性能,这些组合代表了SH的预测百分率(命中率)与我们称之为“烦扰指数”(annoyance index)的比值Rud之间的关系。表7还包括在研究期间每个对象在被警告与无警告状态下经历的平均总时间(以天为单位),也就是说,用阐明了比值Rud意义的算法对研究对象经历的被警告-无警告周期交替过程的摘要结果。
表7:SH预测:命中率、烦扰指数和平均时间
总研究(天数) | ||||||
α | β | γ | %命中 | Rud | 标志升高 | 标志降低 |
6.4 | 8.2 | 1.5 | 43.4 | 1∶25 | 7.8 | 198.9 |
6.0 | 7.5 | 1.5 | 45.2 | 1∶20 | 9.6 | 197.3 |
5.5 | 7.5 | 1.5 | 47.2 | 1∶15 | 12.9 | 194.1 |
5.0 | 7.5 | 1.5 | 49.9 | 1∶10 | 19.0 | 190.1 |
5.0 | 7.5 | 1.3 | 51.3 | 1∶9.5 | 19.5 | 185.7 |
4.9 | 7.0 | 1.2 | 53.1 | 1∶8.4 | 21.6 | 182.0 |
4.8 | 7.0 | 1.2 | 53.4 | 1∶7 | 25.5 | 178.2 |
最好的解用于进一步的分析。假定本研究中的参与者平均经历4.7次SH事件,那么如果报警可以预防50%的SH,则19天的高报警期似乎是可以接受的。此外,高报警期趋向于群集地(in clusters)来临。因此,我们能够假设,在实践中,长而且相对平静的时期会与少数几天高危险报警期交替。表7中的最后一行给出了Rud=1∶7的解,其与本发明公开中给出的解相当。然而,当前的解具有大约高10%的命中率:53.4%对我们先前算法的44%。当命中率与我们先前的算法相当时,烦扰比小于1∶20,也就是说,好了3倍。
图14给出了以百分数表示的命中率与比值Rud之间的平滑依赖关系。显然,当算法3的命中率增加时,“标志升高”与“标志降低”之间的比值迅速增加。因此,给定这些数据,寻找能够产生高于50%命中率的参数组合是不恰当的:
可选择的预测周期:在开始说明算法3时,我们制定了基本假设,如果在SH事件之前的24小时内建立了标志,则认为SH事件被预测到了。该假设产生了在前一部分报告的命中率。现在我们根据范围从12小时到72小时的其他预测周期计算命中率。在整个实验中,参数α、β和γ被分别固定为5.0、7.5和1.5,也就是,固定在它们在表7中最佳解的数值。因此,标志升高比值与该解保持相同,Rud=1∶10,只是命中率发生了变化,因为我们改变了命中的定义。图15给出了预测周期与相应命中率之间的依赖关系。
显然,随着预测周期增加到大约24小时,命中率迅速增加,然后命中率的增加缓慢下降。因此,我们能够得出结论,提前24小时是最佳的且合理的预报周期。
每天SMBG读数的最佳数目最后,我们进行了实验,研究为了产生最佳的SH预报每天需要多少个读数。
正如我们在开始时说的,总共报告了399个SH事件。这些事件中,343个具有24小时之前的SMBG读数(还有3个事件具有之前48小时内的读数,另有4个事件具有之前72小时的读数)。另外有超过50个SH事件(14%)没有任何有助于预测的合理预先SMBG读数。在先前24小时内至少具有一个SMBG读数的343个事件用于计算前一部分中的命中率。其余事件自然被排斥在计算之外。
进一步的分析显示,随着SH事件之前采集的读数的增加,命中率快速增加。然而,如果我们强行严格要求一定数目的可获得读数以便考虑一个SH事件,则我们会发现,满足这一要求的SH事件数目会迅速下降(表8)。这是由于对象不遵守研究的要求,并且可能是在将来的测量计中出现某种警报的良好理由,即如果没有以合适的速度采集SMBG读数,则算法3将不再有用并且测量计会被关闭。
表8给出了具有一定数目先前SMBG读数的SH事件的数目以及算法3对这些事件的命中率。表中最佳的行包括了表7的最优解,其作为后面所有计算的基础。所有的命中率都是以24小时预测周期给出的,也就是,SH之前24小时内的标志。我们能够得出结论,随着对象顺从性的增加,算法3预测SH的精度也显著增加。每天进行5次SMBG读数,则精度从基准的50%命中率提高10%:
表8:当给出一定数目的先前SMBG读数时,算法3的性能
先前SMBG读数的数目 | 满足第1列要求的SH事件(占SH总数的%) | 命中率 |
24小时内至少有1个 | 343(86%) | 49.9% |
24小时内至少有3个 | 260(65%) | 54.2% |
24小时内至少有4个 | 180(45%) | 57.2% |
24小时内至少有5个 | 103(26%) | 64.1% |
36小时内至少有4个 | 268(67%) | 52.6% |
36小时内至少有5个 | 205(51%) | 54.6% |
36小时内至少有6个 | 146(37%) | 60.3% |
36小时内至少有7个 | 107(27%) | 60.7% |
48小时内至少有6个 | 227(57%) | 53.3% |
48小时内至少有7个 | 187(47%) | 54.0% |
48小时内至少有8个 | 143(36%) | 55.9% |
48小时内至少有9个 | 107(27%) | 59.8% |
被测试的其他可能的改善
通过包含外部参数,例如前一年SH的数目或者HbA1c的基准值,提高算法3预测能力的尝试并不成功。显然,SH的短期预测主要取决于当前或近期的事件。然而,该研究的限制条件是,所有的参与者在前一年经历≥2次SH。
最后,我们检验了当探测到低血糖症的危险增加时并且当时已经有一段时间没有SMBG,则是否应当发出SH的报警。这样做是想预测至少一部分之前没有SMBG读数的SH事件。但并不成功,所产生的主要是错误的报警。这一结果更加证实了遵守SMBG协议的重要性,该协议要有足够频繁的SMBG读数。
附录:摘要
实例No.1估计了低血糖症(BG<3.9mmol/l)事件之后再发生低血糖症和SH(定义为无法进行自我治疗的昏迷或者无意识)的频率。
85名在上一年经历>2次SH事件的T1DM患者(41名女性),每天进行3-5次SMBG达6-8个月,并根据日期和时间记录SH事件。对象的平均年龄为44±10岁,糖尿病的持续时间为26±11年,HbA1c为7.7±1.1%。
所有的SMBG读数(75,495)按照日期和时间与对象的SH事件结合(n=399;SH事件一般没有相应的SMBG读数)。为每个SMBG读数或者SH事件计算从最近的上一次低BG(<3.9mmol/l)之后经历的时间。下面的表9给出了3个低血糖范围:BG<1.9mmol/l、1.9-2.8mmol/l和2.8-3.9mmol/l中读数的百分比,以及在之前24小时、24-48小时、48-72小时以及超过72小时具有低BG读数(BG<3.9mmol/l)的SH事件的百分比。最后一列给出了运行检验,其拒绝了如下的假设,即包含低BG读数(或者SH事件)的日期在整个时间范围上随机分布。检验的负Z值显示出具有和没有低血糖读数或者SH事件的日期“群集出现”。
表9:之前具有低BG的低血糖症/SH的百分比:
BG | <24h | 24-48h | 48-72h | >72h | 运行检验 | |
Z | P水平 | |||||
<3.9mmol/l | 50% | 21% | 10% | 19% | -13.6 | <.0001 |
2.8-3.9mmol/l | 52% | 20% | 10% | 18% | -18.3 | <.0001 |
1.9-2.8mmol/l | 55% | 20% | 7% | 18% | -14.7 | <.0001 |
SH | 64% | 11% | 6% | 19% | -11.1 | <.0001 |
我们得出结论,有超过一半的低血糖SMBG读数和大约2/3的SH事件在其之前24小时内有至少一个低血糖读数。此外,低血糖事件趋向于群集地出现。因此,最初的低血糖事件可能是即将再次出现低血糖症的报警信号。
II.实例No.2
本发明采用日常的自我监测血糖(SMBG)数据,并且直接适合于通过引入能够预测HbA1c和显著低血糖症高危期的智能数据判读逻辑提高家用SMBG设备的性能。该方法包括两个部分:(1)算法1估计HbA1c,和(2)算法2&3分别预测长期和短期(24小时内)的显著低血糖症。在该报告中,我们描述了提出、优化和验证HbA1c估计算法1的步骤以及其在估计实验室获得的HbA1c的精度。
目标:
我们的主要目的是达到有95%的测量结果位于实验室参考值±1HbA1c单位内的精度,这是HbA1c化验精度的国家糖基化血红蛋白标准化大纲(NGSP)准则(National Glycohemoglobin StandardizationProgram Criterion)。
方法:
对象:SMBG数据是由100名1型糖尿病对象和100名2型糖尿病对象(T1DM,T2DM)分别经过6个月和4个月获得的,并且在T1DM中于第0、3和6个月进行HbA1c检测,而在T2DM中于第0、2和4个月进行检测。
算法1的提出和优化:训练数据组包括由T1DM收集了3个月的SMBG和HbA1c数据,和由T2DM收集了2个月的SMBG和HbA1c数据。这些训练数据用于优化算法1,和用于评估大量的用于保证更高精度的样本选择标准。样本选择标准是对通过测量计收集的SMBG样本的要求,如果能够满足则能够保证由样本精确地估计HbA1c。因此,测量计会扫描每一个SMBG样本,如果满足样本选择标准,则计算和显示HbA1c估计。在经过分析大量的切点(cut point)之后,我们选择了如下的标准:
1.测试频率:为了产生HbA1c的估计,测量计需要在过去的60天内平均每天进行2.5次或者更多次测试,也就是,在过去的两个月内进行总共要有150个SMBG读数。注意这是每天的平均数,并未要求每天进行测试,这一点是重要的。
2.数据随机化:60天样本中只具有餐后测试或者夜间测试不够(<3%样本)的样本将被排除。此外还包括避免高度集中在每天的一个最常用的时间内进行的测试。这些标准本报告中有详细说明。
结果:算法的前瞻性验证和精度:
然后,将算法,包括样本选择标准,应用于测试数据组1和独立测试数据组2,其中测试数据组1包括T1DM和T2DM对象在上一次HbA1c测试之前2个月内的SMBG和HbA1c数据,而独立测试数据组2由参加了先前NIH研究的60名T1DM对象的数据构成。出于验证的目的,将由算法1获得的估计与参考HbA1c水平进行比较。在测试数据组1中,算法达到了NGSP标准,位于实验室参考值±1HbA1c单位内的精度为95.1%。在测试数据组2中,算法达到了NGSP标准,同样位于实验室参考值±1HbA1c单位内的精度为95.1%。样本选择标准的研究显示,有72.5%的对象能够每天产生一个这样精确的估计,有94%的对象能够大约每5天产生1个这样精确的估计。
结论:日常SMBG数据允许精确地估计HbA1c,并且满足直接HbA1c化验精度的NGSP标准。
对象&取舍标准
我们得到了100名患有1型糖尿病(T1DM)的对象和100名患有2型糖尿病(T2DM)的对象的允许。179名对象完成了SMBG数据收集的主要部分,其中90名患有T1DM,89名患有T2DM。这179名对象的数据用于检验算法2和3。然而,检验算法1要求对象不仅具有SMBG数据,而且要有在SMBG之前60天内采集的HbA1c数据和SMBG记录。在本研究的第3个月(对于T2DM是死2个月),153名对象(78患有T1DM)完成了满足上述标准的HbA1c数据和SMBG记录。此外,我们使用N=60名患有T1DM的对象的数据检验算法1,这些对象参与了我们先前的NIH研究(NIH)。表10给出了所有对象的人口统计学特征。
表10:对象的人口统计学特征
变量 | T1DM | T2DM | NIH |
年龄(岁) | 41.5(11.6) | 50.9(8.1) | 44.3(10.0) |
性别:%女性 | 41% | 43% | 46% |
糖尿病的持续时间(年) | 20.1(10.1) | 11.7(8.2) | 26.4(10.7) |
体质量指数 | 25.4(4.7) | 34.2(8.1) | 24.3(3.4) |
基准HbA1c | 7.5(1.1) | 8.5(2.1) | 7.6(1.0) |
第二HbA1c | 7.3(1.2) | 7.9(1.6) | 7.4(0.8) |
第三HbA1c | 7.0(0.9) | 7.5(1.1) | - |
#SMBG读数/对象/天 | 5.4(2.3) | 3.5(0.8) | 4.1(1.9) |
在第二次HbA1c之前2个月内具有SMBG读数的天数 | 56.9(5.4) | 57.3(4.3) | 37.5(14.3) |
观察到的测量计误差
我们的调查显示,HbA1c化验或者其他化验之前60天内没有完成数据的主要原因不是对象不顺从,而是测量计失效。显然如果患者按下“M”按钮太长时间,则ONE TOUCH ULTRA测量计的时间和日期会“跳转”到随机的日期/时间(例如2017年11月)。一返回我们就检查每个测量计的日期/时间,我们发现在整个研究过程中有60个测量计发生了这种事件。时间/日期的偏移影响了15,280个读数,或者将近全部读数的10%。我们单独保存这些读数并让一名学生检查它们。在许多情况下,但并非在所有的情况下,他都能够恢复读数的日期/时间序列。该误差以及在邮寄过程中丢失了少数测量计,导致具有用于分析算法1的良好数据的对象数目从179减少到了140。恢复了12名对象的数据,最终采集到了153名对象的数据,其中78名患有T1DM,75名患有T2DM,他们在HbA1c之前的数据的时间序列没有被扰乱,因此适合于检验算法1。
程序
所有的对象都填写了IRB推荐的准许表格并参加了指导会(orientation meeting),在会上给他们介绍了ONE TOUCH ULTRA测量计并让他们完成了屏幕问卷。指导会之后,所有的对象立即参观UVA实验室并被采血获得了基准HbA1c。T1DM对象进行了6个月,于第3和第6月进行实验室HbA1c化验;T2DM对象进行了4个月,于第2和第4月进行实验室HbA1c化验。自我监测(SMBG)数据被规律地从测量计下载下来并保存在数据库中。每两周通过自动e-mail/电话跟踪***对显著低血糖和高血糖事件进行平行记录。
数据存储和清理
T1DM和T2DM对象的ONE TOUCH ULTRA粗数据被分别存储在InTouch数据库中。使用定制开发的软件清理这些粗数据的对象和测量计误差,在某些情况下进行手动数据清理(见上面的测量计误差)。当不可能进行修正时,则放弃该数据。
为了保证我们的优化结果能够一般化到全部人口的水平,首先用训练数据组优化算法,然后用测试数据组进行检验。
训练数据组包括T1DM对象在确定第3个月HbA1c之前采集的60天的SMBG数据。该数据组用于优化算法1的公式。T2DM对象在确定他们第2个月HbA1c之前收集的数据用于鉴定样本选择标准,该标准在T1DM数据中不明显。然而,T2DM对象的数据不用于优化算法1的公式。含有这些数据的文件是PASS01.DAT。
测试数据组1包括T1DM对象在确定第6个月HbA1c之前,T2DM对象在确定笫4个月HbA1c之前采集的60天SMBG数据。下面我们将这些数据称作数据组1。含有这些数据的文件是PASS02.DAT。
测试数据组2包括先前NIH研究中N=60名T1DM对象的数据。这些数据用ONE TOUCH PROFILE测量计收集。下面我们将这些数据称作数据组2。含有这些数据的文件是HAT0.XLS。
PASS01.DAT、PASS02.DAT和HAT0.XLS中的变量如下:
ID、月、日、时、年——自我解释的ID数字和读数时间。
PLASBG——由One Touch Ultra记录的BG(HAT0.DAT中的N/A,因为使用了One Touch Profile)。
RISKLO、RISKHI——表示数据转换结果的控制变量(见下文)。
BG和BGMM——将BG转换成全血BG,然后以mmol/l表达(见下文)。
(每个对象的)合计数据、HbA1c、其估计值以及估计误差被存储在Excel文件PASS01.XLS和PASS02.XLS中。
PASS01.XLS、PASS02.XLS和HAT1.XLS中的变量如下:
ID、(糖尿病)类型
HBA1——参考基准HbA1c值
HBA2——第3个月的参考HbA1c(T2DM为2个月)——这需要加以预测;
EST2和ERR2——HbA1c及其误差的估计;
控制变量(所有被算法1使用的变量):
BGMM1——以mmol/l为单位的平均BG(见下面的第2部分);
RLO1、RHI1——低和高BG指数(见下面的第2部分);
L06-夜间低BG指数——限据午夜-6:59a.m的读数加以计算(也就是,if(0.1e.HOUR.1e.6));
NC1=过去60天中SMBG读数的个数;
NDAYS=过去60天中具有SMBG读数的天数。
N06-时间间隔0-6:59;7-12:59中SMBG读数的百分数;
ECLUDE=0,1——如果ECLUDE=1,则建议算法排除该样本。
文件PASS01.DAT和PASS1.XLS能够通过对象的ID数字加以匹配。类似地,文件PASS02.DAT和PASS2.XLS,HAT0.XLS和AHAT1.XLS也能够通过对象的ID数字加以匹配。粗数据和全部第二代数据文件被发送给LifeScan有限公司。
算法1的提出
公式推导:
在本项目的实例No.1中给出了算法1大部分的解释和提出。实例No.1不包括数据收集。相反,我们使用由Amylin药剂师在临床试验中收集的数据组。实例No.1提出了三种可能的公式用于由SMBG数据估计HbA1c:(1)使用平均SMBG、低和高BG指数的公式;(2)使用平均SMBG和先前参考HbA1c读数的公式;和(3)只使用平均SMBG的简单线性公式(见实例No.1)。
还提出了另一个用于评估HbA1c精度的客观标准(在实例No.1中,我们使用最小平方估计、%误差和绝对误差估计每个公式的精度)。这一新的要求被转变成算法1的一个不同的优化标准,即公式的优化不再是产生最小的误差平方和(最小平方估计),而是使估计值固定在HbA1c参考值±1的均匀范围内。
为了这样做,我们只用T1DM对象的训练数据分析我们最初的线性模型(实例No.1的公式)相对于该均匀匹配的误差。我们发现,这些误差与对象的高BG指数正相关(r=0.3),并且我们使用这一关系修正我们最初的线性模型。我们发现,最好使用高BG指数作为分组变量,将对象样本分成高BG指数渐增的组,并且在每个组中将修正引入该线性模型。我们的想法是将使用低BG指数的修正引入到每个特殊的组内,而不是象实例No.1中建议的那样引入到所有的样本。这一变化用基于NBSP标准的不同优化方案加以指示。
因此,根据T1DM对象的训练数据,我们完成了如下的算法1:
部分1-数据预处理
BG=PLASBG/1.12(将血浆BG转换成全血BG,其是通用的)。
BGMM=BG/18(将BG转换成mmol/l)。
下面各行计算每个SMBG读数的低和高BG指数:
COM SCALE=(in(G))**1.08405-5.381.
COM RISK1=22.765*SCALE*SCALE.
COM RISKLO=0.
IF(BG≤112.5)RISKLO=RISK1.
COM RISKHI=0.
IF(BG>112.5)RISKHI=RISK1.
下面各行合计了每个对象的数据:
BGMM1=每个对象的平均(BGMM);
RLO1=每个对象的平均(RISKLO);
RHI1=每个对象的平均(RISKHI);
L06=只为夜间读数计算的平均(RISKLO),如果没有夜间读数则缺省.
N06,N12,N24——分别是时间间隔0-6:59;7-12:59和18-23:59中SMBG读数的百分数,例如if(0≤HOUR≤6);if(7≤HOUR≤2)和if(18≤HOUR≤24).
NC1=过去60天中SMBG读数的总数;
NDAYS=过去60天中具有SMBG读数的天数。
部分2-估计程序:
该估计程序是基于我们实例No.1的线性模型:
HbA1c=0.41046*BGMM+4.0775.
分析该公式的误差我们发现,误差取决于高BG指数。因此,我们根据对象的高BG指数对全部对象进行分类,然后在每个组中对线性模型进行修正,如下:
A.每个对象根据他/她的高BG指数指派一个组:
if(RHI1≤5.25 or RHI1≥16)GRP=0.
if(RHI1>5.25 and RHI1<7.0)GRP=1.
if(RHI1≥7.0 or RHI1<8.5)GRP=2.
if(RHI1≥8.5 or RHI1<16)GRP=3.
B.对于每个组,我们进行如下估计:
E0=0.55555*BGMM1+2.95.
E1=0.50567*BGMM1+0.074*L06+2.69.
E2=0.55555*BGMM1-0.074*L06+2.96
E3=0.44000*BGMM1+0.035*L06+3.65.
EST2=E0
if(GRP=1)EST2=E1.
if(GRP=2)EST2=E2.
if(GRP=3)EST2=E3.
C.对极少数局外值(outlier)进行修正:
if(missing(L06))EST2=E0.
if(RL01≤0.5 and RHI1≤2.0)EST2=E0-0.25.
if(RL01≤2.5 and RHI1>26)EST2=E0-1.5*RLO1.
if((RLO1/RHI1)≤0.25 and L06>1.3)EST2=EST2-0.08.
精度标准
为了评估算法1的精度,我们使用了多个标准判据:
1)NGSP精度标准要求全部估计中有至少95%位于HbA1c参考值±1HbA1c单位内。
2)HbA1c估计值与测量值的平均绝对偏差;
3)HbA1c估计值与测量值的平均百分比偏差。
重点注意:NGSP精度标准是为检验直接测量HbA1c的设备而设计的。这里,我们应用该标准到由SMBG数据对HbA1c的估计。但是,这种估计的目的不是取代HbA1c实验室测量,而是帮助患者和医生进行糖尿病的日常管理。与实验室测量相对,该估计采用能够以任何方法获得的并且能够在日常获得的数据,而不要求特殊的设备或者去医生诊所。
为了证实HbA1c的其他直接测量是与传统的实验室测量是否一致,我们检测了21名IDDM患者的血样并同时用DCA2000和临床化验分析HbA1c。这21个检测结果中,有一个大的误差,为2.5个HbA1c单位。表11给出了该FDA推荐的办公设备的精度结果:
表11:DCA2000在T1DM中的精度
DCA2000 | |
NGSP标准——位于±1HbA1c单位内的百分数 | 95.2% |
平均绝对误差(HbA1c单位) | 0.45 |
平均百分比误差 | 5.7% |
样本选择标准
公式推导:
HbA1c的估计使用连续60天的SMBG。我们以这连续60天的SMBG数据作为样本。在他/她的SMBG期间,每个人会产生大量的样本。实际上,每次新的测量都会产生一个与先前一个略微不同的新样本。因此,自然会假定测量计具有一些控制点,以保证用于估计HbA1c的SMBG样本数据的质量。
因此,一般化算法公式被优化之后,便被应用于整个训练数据组(T1DM和T2DM对象数据),以调查SMBG样本会产生不准确HbA1c估计的条件。
该调查集中于如下的在SMBG中出现的导致不准确估计的模式:
1)SMBG不频繁——在两个月中需要一定数目的读数以便估计HbA1c,如果该数目没有达到,则估计可能不准确;
2)当对象主要在餐后进行测试或者服用口服药物而主要关注于高BG时,会出现SMBG偏向于高血糖症的模式;
3)SMBG的时间偏颇模式,即主要在每天的固定时间进行测试,使得对象的BG波动没有良好的日内分布。
调查完这些模式之后,我们根据最精确的和最少排除的切点选择了最佳的样本选择标准。关于该程序逻辑的详细说明和编码目的的叙述请参考附录A。
最终样本选择标准:
标准1.测试频率:算法要求60天样本包括平均每天至少2.5次测试,即在过去的60天内至少有150个SMBG读数以产生HbA1c估计值(NC1>=150)。
标准2.数据随机化:
2a)口服治疗/餐后检测:(RLO1/RHI1>=0.005)。在一些SMBG样本中,SMBG的分布似乎非常偏向于高血糖症。这主要在T2DM对象中发生,他们似乎只在夜间测量BG。我们假定这些样本不包括低血糖范围内的测试。我们的调查显示,这些样本有大约1/3会对HbA1c产生过高的估计(2/3仍然会产生准确的估计)。据此我们推荐,如果遇到偏颇的样本,则测量计不显示结果,该计算被公式化为LBGI至少为HBGI的1/2%。
2b)夜间测试:(NO6>=3%)。该标准保证至少有一部分夜间低血糖症能够得到解释。该标准要求全部读数的3%在夜间进行(午夜-7:00am)。换言之,如果2个月内采集的150个读数中有至少5个在夜间,则该样本是可以接受的。注意,患者通常被建议在夜间进行测试,因此该标准能够促进良好的管理。
2c)防止高度异常的测试模式:如果超过3/4的读数在每天的任意6小时间隔内进行,则该样本不能够产生估计值。例如,如果样本中有80%的测试是恰在早饭之后进行的,则不进行估计。该标准是LifeScan有限公司所要求的,以防止人们试图“困惑算法”(beat thealgorithm),借此允许我们确保有效,特别是对临床医生。
顺序采用选择标准在训练数据组中的精度
下面的表说明了所选样本选择标准对在训练数据组中的精度和排除数目的影响。注意提出作为本研究(最终算法)一部分的算法1最终版本的精度和在实例No.1中提出的并包含在实例No.1中的最简线性函数的精度(见第一线性模型)。
我们给出了每个模型在没有任何样本选择标准时和顺序应用了样本选择标准1——检测频率,#读数NR≥150,和标准2——数据随机化,时的精度,如上所述。
如在所有表中见到的,算法1的精度随着顺序采用样本选择标准而提高,并且在应用了所有标准之后达到了NGSP要求的95%。后一结果在表中是突出的。
表12A:训练数据组中的最终样本选择标准——所有对象:
最终算法 | 最初线性模型 | |||||
无样本排除 | 标准1 | 标准1和2 | 无样本排除 | 标准1 | 标准1和2 | |
NGSP标准-位于±1HbA1c单位内的百分比 | 93% | 93% | 95.5% | 83% | 83% | 90% |
平均绝对误差(HbA1c单位) | 0.54 | 0.53 | 0.47 | 0.61 | 0.59 | 0.52 |
平均百分比误差 | 7.2% | 7.2% | 6.8% | 8.2% | 8.2% | 7.6% |
绝对误差>1的对象# | 11 | 9 | 5 | 26 | 22 | 11 |
表12B:训练数据组中的最终样本选择标准——T1DM:在该样本中优化了算法1的系数,其解释了即使没有样本选择时的高精度。
最终算法 | 最初线性模型 | |||||
无样本排除 | 标准1 | 标准1和2 | 无样本排除 | 标准1 | 标准1和2 | |
NGSP标准-位于±1HbA1c单位内的百分比 | 96% | 96% | 96% | 86% | 88% | 90% |
平均绝对误差(HbA1c单位) | 0.45 | 0.46 | 0.45 | 0.54 | 0.53 | 0.51 |
平均百分比误差 | 6.3% | 6.6% | 6.5% | 7.7% | 7.8% | 7.6% |
表12C:训练数据组中的最终样本选择标准——T2DM:样本选择标准2(数据随机化)主要用这一样本提出,其解释了当应用该标准时增加的5%精度。
最终算法 | 最初线性模型 | |||||
无样本排除 | 标准1 | 标准1和2 | 无样本排除 | 标准1 | 标准1和2 | |
NGSP标准-位于±1HbA1c单位内的百分比 | 89% | 90% | 95% | 80% | 79% | 91% |
平均绝对误差(HbA1c单位) | 0.63 | 0.62 | 0.52 | 0.68 | 0.66 | 0.53 |
平均百分比误差 | 8.2% | 8.0% | 7.3% | 8.7% | 8.5% | 7.5% |
训练数据中的样本排除频率:
在每次新读数时测量计都有机会估计HbA1c。如果样本没有满足选择标准,则测量计不会显示HbA1c的估计,并且会:
(a)等待直到收集到合适的样本,或者
(b)如果没有收集到合适的样本,例如某人具有永久的偏颇测量方式,则测量计会发出修正SMBG模式的提示。
我们的调查显示,大部分的对象(>95%)在60天内会获得至少10个HbA1c估计(只要他们的测量频率足够),而只有2%的对象会由于偏颇的测量方式而没有获得估计。需要提示这2%的对象修正他们的测量方式。该调查的最终结果在下文给出:
我们计算了(在60天中)有多少天测量计会由于样本不满足选择标准而不能向用户显示HbA1c结果:
1)对于全部对象的72.5%,测量计能够每天报告HbA1c;
2)对于全部对象的另外7.5%,测量计能够报告(60天中)45-59天的HbA1c;
3)对于全部对象的另外10%,测量计能够报告(60天中)12-44天的HbA1c;
4)对于9名对象(5.9%),测量计不能够报告HbA1c,除非他们改变SMBG方式。
重点注意:这些对象中的大多数都不会获得估计,因为他们没有满足测试频率标准1,即他们的样本总是小于150个读数。因此,全部对象中至少有94%能够大约每5天获得至少一个HbA1c,而无需改变他们的测量方式(这包括T1DM和T2DM)。
如果我们要求在60天内至少有150个读数,则只有3名对象不能获得HbA1c估计:
1)95.6%在60天内得到至少10个HbA1c估计;
2)2.2%不会获得任何估计。
因此,在60天内有大约98%的平均每天测量2.5次的对象会得到HbA1c估计,有>95%的将每周至少获得1个估计。我们得出结论,样本选择标准2——数据随机化在一段时期内对显示HbA1c估计的影响很微小。只有大约2%的对象需要被提示改进他们的SMBG方式。
应当注意,样本选择标准能够用于提高任何估计HbA1c的公式的精度。该选择标准独立于任何特殊的算法/公式,并且在估计开始之前应用。例如,当被使用时,样本选择标准将提高最新提出的作为本研究一部分的算法1的精度,以及在实例No.1中提出的我们最初的线性模型的精度。
此外,检查其他一些样本选择标准的效果显示,我们能够进一步提高精度,这是被期望的。例如,当将其中一个原始测试频率标准应用于该数据时,能够证实具有更多的效力。该标准在附录E中有进一步的说明。
算法1的前瞻性验证:
在测试数据组1中的精度:
随后将算法,包括最终的样本选择标准,应用于测试数据组1(T1DM 1和T2DM对象最后一次HbA1c之前2个月的SMBG)以产生HbA1c估计。然后将这些估计值与HbA1c参考值进行比较,从而前瞻性验证算法1。表13给出了该验证的概要。关于每个样本选择标准对算法的影响的更详细说明请见附录C。
表13:算法前瞻性应用的精度:
具有标准1和2的最终算法 | |||
全部对象 | T1DM | T2DM | |
NGSP标准——±1HbA1c单位内% | 95.1% | 97% | 93% |
平均绝对误差(HbA1c单位) | 0.45 | 0.38 | 0.54 |
平均百分比误差 | 6.2% | 5.4% | 7.4% |
测试数据组2中的精度:
另一个独立的NIH数据组(N=60名患有T1DM的对象)证实了结果具有类似的精度,即有95.5%处于实验室参考值±1HbA1c单位之内(表14):
表14:算法1在独立NIH数据组中的精度:
全部对象 | |
NGSP标准——±1HbA1c单位内的百分率 | 95.5% |
平均绝对误差(HbA1c单位) | 0.42 |
平均百分比误差 | 5.9% |
算法1的精度与FDA推荐的办公设备精度的比较:
如下面的表15中显示的,算法1的精度与医生诊所中使用的HbA1c化验精度相当。如在精度标准部分中说明的,DCA2000数据用于证实HbA1c的其他直接测量结果是否与实验室测量结果一致。我们同时用DCA2000和临床化验分析了21名T1DM患者的HbA1c血液样本。在这21次测试中有较大的误差,为2.5HbA1c单位:
表15:DCA2000在T1DM中的精度与算法1比较:
DCA2000 | 测试数据组1 | 测试数据组2 | |
NGSP标准——±1HbA1c单位内的百分率 | 95.2% | 95.1% | 95.5% |
平均绝对误差(HbA1c单位) | 0.45 | 0.45 | 0.42 |
平均百分比误差 | 5.7% | 6.2% | 5.9% |
测试数据组中样本排除的频率:
正如我们在提出算法1的部分中讨论的,在每次新读数的时候,测量计都有机会估计HbA1c。如果样本不满足选择标准,那么测量计就不会显示HbA1c。
我们使用测试数据组1和2前瞻性估计样本排除的频率。为此,我们的计算是根据测量计能够有多少天(60天中)向人显示HbA1c,也就是人们能够有多少天具有满足样本选择标准的样本。表16A和16B给出了测试数据组1和2中这些结果的摘要。我们包括了所有对象的数据,这些数据分成平均测量1.5次/天的对象(在60天中有90个SMBG读数)和2.5次/天的对象的数据:
表16A:测试数据组1中样本排除的频率
全部对象(N=148) | 平均测量≥1.5次/天的对象(N=146) | 平均测量≥2.5次/天的对象(N=130) | |
测量计能够每天报告HbA1c的对象的百分率 | 69.6% | 72.6% | 77.7% |
测量计能够每3天报告一次HbA1c的对象的百分率 | 87.8% | 91.1% | 93.1% |
测量计能够每周报告一次HbA1c的对象的百分率 | 91.9% | 95.5% | 96.9% |
表16B:测试数据组2中样本排除的频率
全部对象(N=60) | 平均测量≥1.5次/天的对象(N=55) | 平均测量≥2.5次/天的对象(N=30) | |
测量计能够每天报告HbA1c的对象的百分率 | 51.7% | 83.6% | 80.0% |
测量计能够每3天报告一次HbA1c的对象的百分率 | 95.0% | 100.0% | 100.0% |
测量计能够每周报告一次HbA1c的对象的百分率 | 96.7% | 100.0% | 100.0% |
结论:
表13-16证实,测量计能够平均一周产生一个满足95%NGSP精度标准的HbA1c精确估计,对于平均每天测量2.5次的人则为>96%。
附录-样本选择标准的软件逻辑
向对象发送样本选择标准——建议通过算法排除一些SMBG样本,或者消息。该样本选择标准被编程如下:
标准1.测试频率:本算法要求60天样本要包括平均每天至少测量2.5次,即在过去的60天中要有至少150个SMBG读数以产生HbA1c估计:
EXCLUDE=0
if(NC1>=150)EXCLUDE=1
标准2.数据随机化:
2a)口服治疗/餐后测试:在一些SMBG样本中,SMBG的分布似乎非常偏向于高血糖症。这主要在T2MD对象中发生,他们似乎只在高BG时进行测量。我们假定,这些样本不包括在低血糖范围测量的结果。我们的调查显示,这些对象中大约有1/3会过高地估计HbA1c(2/3的对象仍然产生精确的估计)。据此,我们建议如果样本出现了偏颇,则测量计不显示结果,该计算被公式化为LBGI至少为HBGI的1/2%。
if(RLO1-RHI1<0.005)EXCLUDE=1。
2b)夜间测试:(NO6>=3%)。该标准保证至少能够解释一部分的夜间低血糖症。该标准要求全部读数的3%在夜间进行(午夜-7:00am)。换言之,如果在两个月中获得的150个读数中至少有5个在夜间,则该样本可以接受。注意,通常建议患者在夜间进行测试,所以该标准能够促进良好的管理。
if(NO6≤3.0)EXCLUDE=1。
2c)确保防止高度异常的测试模式:如果有超过3/4的读数一天中的任意6小时间隔内进行,则样本不会产生估计值。例如,如果样本中有80%的测试恰好在早饭之后进行,则不会产生估计值。该标准是LifeScan公司要求的以便保证人们试图“困惑算法”,借此允许我们确保有效性,特别是对于医师。在我们的数据中,没有样本如此高度地异常以致于触发该标准(详细信息见附录B——标准2c)。依靠软件执行,需要从SMBG数据中计算如下的频率:
M12-6:00am-中午的SMBG读数%(早餐)
M18-中午-6:00pm的SMBG读数%(午餐)
M24-6:00pm-12:00的SMBG读数%(晚餐)
M06-12:00-6:00am的SMBG读数%(夜间)
M15-9:00am-3:00pm的SMBG读数%
M21-3:00pm-9:00pm的SMBG读数%
M03-9:00pm-3:00am的SMBG读数%
M09-3:00am-9:00am的SMBG读数%
然后对于上述的任何组合(i,j):
if(Mij>75.0)EXCLUDE=1。
附录B-样本选择标准2C
本标准是LifeScan公司要求的,以确保防止高度异常的测试模式。该标准的目的是防止人们“困惑算法”。
基本上,本标准规定:如果你的读数的3/4(或者其他期望的数目)在一天的任意6小时间隔或者其他期望的间隔内进行,则你将不能够得到估计值。
因此,例如,如果有超过3/4测试在晚饭后进行,则不会得到估计值。这将更加支持我们的一般声明:本计算中不包括测试不随机的人。我认为,对此的特殊计算和编码可能看上去复杂,但关键的是我们可以只输入“你在一天中必须随机进行测试”或者其他类似的话作为面板声明,就可以覆盖我们全部的排除标准(排除测试频率)。如果我们需要,则我们能够以精练的语句(fine print)简单地输入更加精确的定义“一天中任何6小时间隔内的读数都不可以超过全部读数的75%”。这保证,遵循我们的标准,可以提高算法的临床接受度。
更详细内容:
4个6小时间隔定义如下:
6:00am-中午(早饭)
中午-6:00pm(午饭)
6:00pm-12:00(晚饭)
12:00-6:00am(夜间)
该标准能够以不同的时间间隔运行两次,从而防止人们集中在6小时间隔的交界点附近进行测试但仍然不正确地满足该第一标准。例如,如果人们在11:50pm有40%,在12:10pm有40%,则依然是集中的测试,虽然满足间隔的第一关(pass),但不满足间隔的第二关。
第二组间隔:
9:00am-3:00pm
3:00pm-9:00pm
9:00pm-3:00am
3:00am-9:00pm
注意选择地,从编码的观点来看,人们可以如下地运行而得到相同的结果:
任意18小时周期的读数不得少于总读数的25%。
你必须在相互重叠3个小时的18小时周期内运行它:
9:00am-3:00am
12:00中午-6:00am
3:00pm-9:00am
6:00pm-12:00中午
9:00pm-3:00pm
12:00午夜-6:00pm
3:00am-9:00pm
6:00am-12:00午夜
附录C-测试数据组1中样本选择标准对算法精度的增加效应
如在算法提出部分中说明的,如下的表格涉及了被提出作为本研究一部分的算法1(最终算法)和在实例No.1中提出的并且包含在实例No.1中的最简线性函数的精度(见最初线性模型)。该表格给出了算法1在测试数据组1中无样本排除并依次应用两个样本选择标准的精度:
标准1——测试频率,#读数NR≥150,和
标准2——数据随机化,如在样本选择标准中说明的:
表17A:算法1的精度——全部对象:
最终算法 | 最初线性模型 | |||||
无样本排除 | 标准1 | 标准1和2 | 无样本排除 | 标准1 | 标准1和2 | |
NGSP标准——±1HbA1c单位内% | 86% | 91% | 95.1% | 83% | 85% | 89% |
平均绝对误差(单位HbA1c) | 0.56 | 0.48 | 0.45 | 0.59 | 0.54 | 0.49 |
平均百分率误差 | 7.4% | 6.7% | 6.2% | 7.8% | 7.5% | 7.0% |
#绝对误差>1的对象 | 21 | 10 | 5 | 25 | 16 | 11 |
表17B:算法1在T1DM中的精度:
最终算法 | 最初线性模型 | |||||
无样本排除 | 标准1 | 标准1和2 | 无样本排除 | 标准1 | 标准1和2 | |
NGSP标准——±1HbA1c单位内% | 90% | 95% | 97% | 90% | 93% | 95% |
平均绝对误差(单位HbA1c) | 0.45 | 0.40 | 0.38 | 0.49 | 0.45 | 0.43 |
平均百分率误差 | 6.1% | 5.6% | 5.4% | 6.9% | 6.7% | 6.4% |
表17C:算法1在T2DM中的精度:
最终算法 | 最初线性模型 | |||||
无样本排除 | 标准1 | 标准1和2 | 无样本排除 | 标准1 | 标准1和2 | |
NGSP标准-±1HbA1c单位内% | 81% | 86% | 93% | 76% | 76% | 81% |
平均绝对误差(单位HbA1c) | 0.68 | 0.59 | 0.54 | 0.69 | 0.64 | 0.57 |
平均百分率误差 | 8.6% | 8.0% | 7.4% | 8.7% | 8.6% | 7.8% |
附录D-可选择测试频率标准
高级的测试频率标准有利于更显著地提高算法1的精度。这是因为,采用测试频率标准1并不仅是基于数据分析,而且还基于其他的考虑。如果发现要求2个月内有150个读数的标准1太严格,则可以采用可替换的解决方法。就是原始的测试频率标准,其要求有35天(60天中)的SMBG读数频率为1.8个读数/天,也就是60天中有35天的总读数为63个。表18用该原始的宽松测试频率标准加上标准2(数据随机化)证实这一点,算法1的精度超过95%:
表18:使用选择测试频率标准(35天的读数为1.8个读数/天)和数据随机化标准的算法1的精度:
训练数据组 | 测试数据组 | 测试数据组2 | |||||
全部对象 | T1DM | T2DM | 全部对象 | T1DM | T2DM | ||
NGSP标准-±1HbA1c单位内% | 95% | 96% | 94% | 95% | 97% | 92% | 100% |
平均绝对误差(单位HbA1c) | 0.48 | 0.44 | 0.55 | 0.46 | 0.39 | 0.55 | 0.40 |
平均百分率误差 | 6.8% | 6.3% | 7.6% | 6.3% | 5.5% | 7.2% | 5.5% |
重点注意:除此之外,该可选择标准可以筛选掉具有大量缺失数据的样本,例如,如果SMBG间断达4周之后再恢复,则不会显示HbA1c估计。这种模式的一个明显实例出现在测试数据组2中——他/她的HbA1c估计值具有最大误差的对象只在60天的30天中收集了159个读数。因此,通过在少数几天内快速收集读数,对象仍可以满足收集150个读数的要求,但会导致不准确的HbA1c估计。
实例No.2的典型定义(但对本文没有限制)
1)严重低血糖症(SH)定义为导致无法自我治疗的昏迷、疾病发作或者无意识的低血糖(BG);
2)中度低血糖症(MH)定为扰乱对象的活动但不妨碍自我治疗的严重神经低血糖症;
3)生化严重低血糖症(BSH)定义为血浆BG读数<=39mg/dl;
4)生化中度低血糖症(BMH)定义为血浆BG读数为39-55mg/dl;
5)上面全部病症均称作显著低血糖症。
附加目的
本实例的数据用于前瞻性验证如下的算法:
算法2——使用某一对象30-45天SMBG数据的分类算法,将对象分成发生将来显著低血糖症的一定危险范围。该分类是临时的,例如,当对象的SMBG方式改变时,分类也改变。
算法3——数据跟踪/判定算法,其使用一定序列的SMBG数据判断是否为即将发生的(24小时)显著低血糖症建立标志。现在我们详细说明算法1&2及其测试结果。
对象
我们获得了100名患有1型糖尿病(T1DM)的对象和100名患有2形糖尿病(T2DM)的对象的许可。170名对象,其中90名患有T1DM,89名患有T2DM,完成了SMBG数据收集的主要部分。
程序
全部对象都分配了一个IRB推荐的许可表格并参加了指导会,在会上向他们介绍了ONE TOUCH ULTRA测量计并完成了屏幕问卷调查。指导会之后,全部对象立即参观了UVA实验室并采血测量了HbA1c基准值。T1DM对象在随后的6个月里于第3和6月进行实验室HbA1c化验;T2DM对象在随后的4个月里于第2和4月进行实验室HbA1c化验。自我监测(SMBG)数据被规则地从测量计下载下来并保存在数据库中。通过定制设计的自动e-mail/电话跟踪***平行地记录显著低血糖症和高血糖症事件,其每2周与全部参加者进行联系。表19给出了SMBG和严重低血糖症/中度低血糖症[SH/MH]数据收集的摘要。
表19:数据收集摘要
变量 | T1MD(N=90对象) | T2DM(N=89对象) |
#SH事件 | 88 | 24 |
#MH事件 | 1,660 | 190 |
#SMBG读数 | 92,737 | 35,306 |
#BSH读数 | 1,039 | 39 |
#BMH读数 | 5,179 | 283 |
算法2和3的公式没有显著改变。这些公式与2002年3月实例No.1的报告中给出的公式特别相似。只有两个变化:(a)修正了SH/MH(实例No.1)危险范围列表中的类型和(b)算法3中有一行改变。后者的原因在下面进行解释。
因为算法1和2保持未变,所以我们能够将整个实例No.2的数据收集看作这些算法的前瞻性测试。
算法2的公式
算法2的执行如下:
1)根据一个月的SMBG数据,根据每个对象的低血糖指数如下地将他/她分类成15个危险范围(RCAT)中的一个:
if(LBGI≤0.25),RCAT=0
if(0.25<LBGI≤0.5),RCAT=1
if(0.50<LBGI≤0.75),RCAT=2
if(0.75<LBGI≤1.00),RCAT=3
if(1.00<LBGI≤1.25),RCAT=4
if(1.25<LBGI≤1.50),RCAT=5
if(1.50<LBGI≤1.75),RCAT=6
if(1.75<LBGI≤2.00),RCAT=7
if(2.00<LBGI≤2.50),RCAT=8
if(2.50<LBGI≤3.00),RCAT=9
if(3.00<LBGI≤3.50),RCAT=10
if(3.50<LBGI≤4.25),RCAT=11
if(4.25<LBGI≤5.00),RCAT=12
if(5.00<LBGI≤6.50),RCAT=13
if(LBGI>6.50),RCAT=14
2)通过双参数Weibull概率分布并利用如下给出的分布函数计算将来显著低血糖症的理论概率:
F(x)=1-exp(-a,xb)对于任何x>0;否则0。该分布的参数取决于期望的预测持续时间,并且在实例No.1的报告中有说明。如果用测量计执行,该步骤将提供显著低血糖症危险的连续型估计,例如“在下个月中50%”。
3)将每个对象分类到将来显著低血糖症的最小、低、中或者高危险组中:这些范围的定义如下:最小危险(LBGI≤1.25);低危险(1.25<LBGI≤2.5);中度危险(2.5<LBGI≤5)和高危险(LBGI>5)。如果用测量计执行,则该步骤将提供显著低血糖症危险的离散型估计,例如“在下个月为高危险”。
算法3的公式
首先,为了避免计算在算法3实例No.1报告说明书中给出的基准危险值,我们修改了编码中的一行。现在,算法3转而使用算法2的结果。我们引入该变化用于介绍2002年10月28日两个对象的样本结果。此时,显然用简单的Excel电子数据表例证算法3的行动是方便的,并且如果避免了基准值的计算则是可能的。该步骤没有改变算法3的精度,所以保留作为方便算法3编程的永久改变。2002年10月28日之后不再给算法3引入改变。这里,我们给出了与实例No.1报告中相同的算法3公式,被改变的行已经标记出来。
1)通过如下的编码为每个BG读数计算低BG危险值(RLO)(这里BG以mg/dl进行测量,如果单位是mmol/l则系数会不同):
scale=(In(bg))**1.08405-5.381
risk=22.765*scale*scale
if(bg_1≤112.5)then
RLO=risk
else
RLO=0
endif
2)对于每个SMBG读数,我们计算了运行值LBGI(n),和另一个统计量SBGI(n),其是低BG危险值的标准差。这两个参数的计算是用特定的标记(n)从每个SMBG读数逆推算的,也就是包括该读数及该读数之前的(n-1)个读数。
3)LBGI(n)和SBGI(n)的计算使用临时平均程序(provisionalmeans procedure),其使根据如下的递归编码:
初值在n(或者精确地在最大值(1,n-k),以便解释序数小于k的测量计读数):
LBGI(n)=rlo(n)
Rlo2(n)=0
倒着计数n与1之间的任何连续迭代j值:
LBGI(j)=((j-1)/j)*BLGI(j-1)+(1/j)*RLO(j)
rlo2(j)=((j-1)/j)*rlo2(j-1)+(1/j)*(RLO(j)-LBGI(j))**2
完成该循环之后,我们获得了LBGI(n)的数值,接着计算
SBGI(n)=sqrt(rlo2(n))
由该计算,我们保存了两组数据:n=150和n=50(例如前150个和前50个观察值)。
4)判定规则:在每次SMBG读数,程序都判定是否建立即将发生SH的警告标志。如果标志已经建立,则程序判定是否降低该标志。这些判定取决于三个阈值参数α,β,γ,运行如下:
对于低中危险的对象(LM组):
FLAG=0.
if(LBGI(150)≥2.5 and LBGI(50)≥(1.5*LGI(150)and SBGI(50)≥SBGI(150))FLAG=1.
if(RLO≥(LBGI(150)+1.5*SBGI(150))FLAG=1.
换言之,在每次SMBG读数,如果满足两个条件中的一条则建立标志:
1)根据算法2由前150次试验进行的分类对象处于中度高危险SH,且在前50次试验中LBGI和LBGI的SD增加;
2)或者,通过第二个不等式确定出低BG指数突增。
这些叙述的启发性解释在实例No.1的报告中给出。如上面描述的,第一个“if”的陈述已经改变了其最初的形式,以便避免使用基准LBGI,从而使用算法2的输出。
如实例No.1的报告中说明的,一旦建立起标志,则它将保持24小时。为了评估算法3的精度,我们使用以前提出的技术——计算两个测量值:
1)24小时内即将发生SH/MH事件的预测%,和
2)“标志升高”与“标志降低”持续周期的比值Rud(烦扰指数)。
SH事件的预测%需要比较高,而比值Rud需要比较低。这是因为通过增加预测SH事件的百分率,我们会不可避免地增加“建立标志”的数目,其反过来增加潜在“错误报警”的数目。因为“错误报警”没有被清晰地定义(见实例No.1的报告),所以我们使用Rud作为算法3有效性的指标。
在实例No.1的报告中给出的我们先前最佳的结果是,预测出50%的24小时内SH/MH事件,Rud=1∶10,即在高危险报警一天之后,紧接着10天没有报警。这里我们将保持相同的标志升高/标志降低比值,并分别为T1DM和T2DM对象计算24小时内SH和MH事件的%预测。对于该预测,我们不使用BSH和BMH事件,因为这是通过测量计记录的,因此是预测函数的一部分。
估计1-3个月内显著低血糖症的危险:算法2的精度
我们如下地估计了算法2的预测能力:
1)首先,我们由一个月的SMBG数据计算了LBGI,并且如上所述地将每个对象分类到显著低血糖症的最小、低、中和高危险组。
2)然后,在随后的1-3个月中,我们计数每个对象前瞻性记录的SH、BSH、MH和BMH事件的数目。
下面的图16-19分别为T1DM和T2DM给出了在一个月的SMBG之后每个对象在未来1个月或者3个月内观察到的SH、BSH、MH和BMH事件的数目。同时还包括了统计比较。
此外,直接线性回归使用LBGI、在屏幕问卷中根据过去一年中SH事件数目报告的SH历史,和基准HbA1c,显著地预测了(R2=0.62,f=48,p<0.0001)在之后3个月内即将发生显著低血糖症事件的总数(SH+MH+BSH+BMH)。预测变量按照其重要性的顺序为:1)LBGI(t=8.2,p<0.0001),可单独解释将来显著低血糖症55%的变异(例如R2=0.55);2)SH历史(t=3.6,p<0.0005),可解释另外5%的变异,和HbA1c(t=2.2,p<0.03),可解释另外2%的变异。这证实了先前的结果,即LBGI是将来低血糖症最重要的预测指标,而HbA1c对该预测的贡献是中等的。
通过Weibull模型计算的将来显著低血糖症的理论概率与未来观察到的显著低血糖症事件非常一致——对于严重和中度事件,确定系数均超过90%。
即将发生的(24小时内)显著低血糖症的预测:算法3的精度
下面的表分别给出了T1DM和T2DM对象SH和MH事件短期预测(24小时内)的精度。如果在24小时周期内可以获得一定数目的SMBG读数用于预测,则表20和21的每一行给出了预测到的事件的百分率。例如,每个表的第一行给出了被预测到的事件的百分率,而不管在某个事件之前24小时内是否有SMBG读数。可见,预测的精度随着事件之前读数数目的增加而增加。因此,如果某人每天测量3次或者更多次,则测量计能够报警且有可能帮助避免超过一半的显著低血糖症事件。
重点注意:出于评估算法3精度的目的,我们只使用了通过e-mail/电话***报告的独立于SMBG的SH和MH事件,该***要求参加者每两周报告SH和MH的日期和时间。如我们的调查所显示的,有时参加者在他们的报告中使用某一次事件之前的上一个SMBG读数的时间和日期,而不是该事件的精确时间和日期,因为查询测量计有助于帮助他们进行回忆。结果,有大量的事件,其之前最近一次SMBG读数的时刻与该事件时刻之间的时间间隔接近于零。为了解释这种可疑的时间记录,每个表的第3列给出了只严格限制于如下事件的算法3的精度,即前导报警时间至少为15分钟。假定平均前导报警时间为11小时,我们得出结论,在大多数情况下,报警的出现足够早从而有利于充分地自我治疗。
在表20和21中,烦扰指数被设定为Rud>=10以便和实例No.1的报告相匹配。
表20:算法3在T1DM中的精度
Rud=10.2 | SH+MH事件预测% | SH事件预测% | 报警时间>15分钟时SH+MH事件的预测% | |
无限制 | 53% | 48% | 49% | |
事件之前24小时内SMBG读数的最小数目 | 3 | 55% | 58% | 52% |
4 | 59% | 60% | 55% | |
5 | 63% | 60% | 59% |
表20:算法3在T2DM中的精度
Rud=10.0 | SH+MH事件预测% | SH事件预测% | 报警时间>15分钟时SH+MH事件的预测% | |
无限制 | 52% | 38% | 48% | |
事件之前24小时内SMBG读数的最小数目 | 3 | 57% | 60% | 53% |
4 | 64% | 64% | 59% | |
5 | 73% | 73% | 68% |
本发明可以用其他的特殊形式加以实现,而不背离其精神或者基本特征。因此从任何方面讲,都应认为前述的实施例只是例证,而非对本文说明的发明具有限制。因此本发明的范围由附属权利要求指定,而不是由前述的说明,因此在与权利要求等价的意义和范围内的任何变化都应当包含在本文的范围之内。
Claims (46)
1.一种根据在预定持续时间内收集的BG数据估计患者严重低血糖症(SH)和/或中度低血糖症(MH)长期概率的***,所述***包括:
数据库组件,其用于保存识别所述BG数据的数据库;和
处理器,其被编程用于:
根据所述收集的BG数据计算LBGI;和
根据所述计算的LBGI利用预定的数学公式估计将来SH事件的数目。
2.根据权利要求1的***,其中:
所述LBGI的计算是由在时间点t1,t2,...,tn获取的一系列BG读数x1,x2,...,xn数学地加以定义:
其中:
lbgi(BG;a)=10.f(BG)a,如果f(BG)>0,否则为0,
a≈2,表示权重参数。
3.根据权利要求1的***,进一步包括:
定义预定危险范围(RCAT),每个所述危险范围(RCAT)代表LBGI的一个数值范围;和
将所述LBGI指定给至少一个所述危险范围(RCAT)。
4.根据权利要求3的***,其中所述危险范围(RCAT)定义如下:
范围1,其中所述LBGI小于大约0.25;
范围2,其中所述LBGI介于大约0.25到大约0.50;
范围3,其中所述LBGI介于大约0.50到大约0.75;
范围4,其中所述LBGI介于大约0.75到大约1.0;
范围5,其中所述LBGI介于大约1.0到大约1.25;
范围6,其中所述LBGI介于大约1.25到大约1.50;
范围7,其中所述LBGI介于大约1.50到大约1.75;
范围8,其中所述LBGI介于大约1.75到大约2.0;
范围9,其中所述LBGI介于大约2.0到大约2.5;
范围10,其中所述LBGI介于大约2.5到大约3.0;
范围11,其中所述LBGI介于大约3.0到大约3.5;
范围12,其中所述LBGI介于大约3.5到大约4.25;
范围13,其中所述LBGI介于大约4.25到大约5.0;
范围14,其中所述LBGI介于大约5.0到大约6.5;和
范围15,其中所述LBGI大于大约6.5。
5.根据权利要求3的***,进一步包括:
分别为每个所述指定的危险范围(RCAT)限定发生所选数目SH事件的概率。
6.根据权利要求3的***,进一步包括:
分别为每个所述指定的危险范围(RCAT)限定在下一个第一预定持续时间内发生所选数目SH事件的概率,使用公式:
F(x)=1-exp(-a.xb),对于任何x>0,否则为0,其中:
a≈-4.19,
b≈1.75。
7.根据权利要求6的***,其中所述第一预定持续时间为大约1个月。
8.根据权利要求6的***,其中所述第一预定持续时间范围从大约0.5个月到大约1.5个月。
9.根据权利要求6的***,其中所述第一预定持续时间范围从大约0.5个月到大约3个月。
10.根据权利要求3的***,进一步包括:
分别为每个所述指定的危险范围(RCAT)限定在下一个第二预定持续时间内发生所选数目SH事件的概率,使用公式:
F(x)=1-exp(-a.xb),对于任何x>0,否则为0,其中:
a≈-3.28,
b≈1.50。
11.根据权利要求10的***,其中所述第二预定持续时间为大约3个月。
12.根据权利要求10的***,其中所述第二预定持续时间范围从大约2个月到大约4个月。
13.根据权利要求10的***,其中所述第二预定持续时间范围从大约3个月到大约6个月。
14.根据权利要求3的***,进一步包括:
分别为每个所述指定的危险范围(RCAT)限定在下一个第三预定持续时间内发生所选数目SH事件的概率,使用公式:
F(x)=1-exp(-a.xb),对于任何x>0,否则为0,其中:
a≈-3.06,
b≈1.45。
15.根据权利要求14的***,其中所述第三预定持续时间为大约6个月。
16.根据权利要求14的***,其中所述第三预定持续时间范围从大约5个月到大约7个月。
17.根据权利要求14的***,其中所述第三预定持续时间范围从大约3个月到大约9个月。
18.根据权利要求3的***,进一步包括:
分别为每个所述指定的危险范围(RCAT)限定在下个月内发生所选数目MH事件的概率,使用公式:
F(x)=1-exp(-a.xb),对于任何x>0,否则为0,其中:
a≈-1.58,
b≈1.05。
19.根据权利要求3的***,进一步包括:
分别为每个所述指定的危险范围(RCAT)限定在下3个月内发生所选数目MH事件的概率,使用公式:
F(x)=1-exp(-a.xb),对于任何x>0,否则为0,其中:
a≈-1.37,
b≈1.14。
20.根据权利要求3的***,进一步包括:
分别为每个所述指定的危险范围(RCAT)限定在下6个月内发生所选数目MH事件的概率,使用公式:
F(x)-1-exp(-a.xb),对于任何x>0,否则为0,其中:
a≈-1.37,
b≈1.35。
21.根据权利要求1的***,进一步包括:
为患者的将来显著低血糖症指定危险类别。
22.根据权利要求21的方法,其中所述类别如下定义:
最小危险,其中所述LBGI小于大约1.25;
低危险,其中所述LBGI介于大约1.25到大约2.50;
中度危险,其中所述LBGI介于大约2.50到大约5;和
高危险,其中所述LBGI大约5.0。
23.一种根据在预定持续时间内收集的BG数据估计患者严重低血糖症(SH)和/或中度低血糖症(MH)长期概率的***,所述***包括:
BG获取机构,所述获取机构用于从患者获取BG数据;
数据库组件,其用于保存识别所述BG数据的数据库;和
处理器,其被编程用于:
根据所述收集的BG数据计算LBGI;和
根据所述计算的LBGI利用预定的数学公式估计将来SH事件的数目。
24.一种根据在预定持续时间内收集的BG数据估计患者严重低血糖症(SH)短期概率的方法,所述方法包括:
根据所述收集的BG数据计算scale数值;和
为每个BG数据计算低BG危险值(RLO)。
25.根据权利要求24的方法,其中:
所述RLO(BG)的计算被数学地定义为:
scale=[ln(BG)]1.0845-5.381,其中BG以mg/dl为单位测量;
Risk=22.765(Scale)2
如果(BG<大约112.5)则:
RLO(BG)=Risk,否则
RLO(BG)=0。
26.根据权利要求24的方法,其中
所述RLO(BG)的计算被数学地定义为:
Scale=[ln(BG)]1.026-1.861,其中BG以mmol/l为单位测量;
Risk=32.184(Scale)2
如果(BG小于大约112.5)则:
RLO(BG)=Risk,否则
RLO(BG)=0。
27.根据权利要求24的方法,其中
根据收集的BG数据计算LBGI,所述LBGI的计算由在时间点t1,t2,...,tn采集一系列BG读数x1,x2,...,xn数学地定义:
其中:
lbgi(BG;a)=BLO(BG)。
28.根据权利要求24的方法,其中
根据所述收集的BG数据计算临时LBGI,所述临时LBGI的计算由如下定义的数学公式数学地定义:
LBGI(1)=RLO(X1)
RLO2(1)=0
LBGI(j)=(j-1)/j)*LBGI(j-1)+(1/j)*RLO(xj)
RLO2(j)=((j-1)/j)*RLO2(j-1)+(1/j)*RLO(xj)-LBGI(j))2
29.根据权利要求28的方法,其中
计算SBGI,所述SBGI的计算能够如下地数学定义:
30.根据权利要求29的方法,包括确认或者发出将发生短期SH的报警,其中如果:
(LBGI(150)≥2.5且(LBGI(50)≥(1.5*LBGI(150)且SBGI(50)≥SBGI(150)),则确认或者发出所述报警,或者
RLO≥(LBGI(150)+1.5*(SBGI(150)),则确认或者发出所述报警;
否则:
不需要确认或者提供报警。
31.根据权利要求29的方法,包括确认或者发出将发生短期SH的报警,其中如果:
(LBGI(n)≥α且LBGI(n)≥(β)),则确认或者发出所述报警,和/或
(RLO(n)≥(LBGI(n)+γ*SBGI(n))),则确认或者发出所述报警;
否则:
不需要确认或者发出报警,其中α、β和γ是阈值参数。
32.根据权利要求31的方法,所述阈值参数α、β和γ被定义为α≈5,β≈7.5,γ≈1.5。
33.根据权利要求31的方法,其中所述阈值参数α、β和γ定义为a、b和/或c的任意组合,和/或在下面的a、b和/或c中所述参数值之间任意的中间组合:
a)α=6.4,β=8.2,γ=1.5,α=5.0,β=7.5,γ=1.3;
b)α=6.0,β=7.5,γ=1.5,α=4.9,β=7.0,γ=1.2;和/或
c)α=5.5,β=7.5,γ=1.5,α=4.8,β=7.0,γ=1.2
34.根据权利要求31的方法,其中所述阈值参数α、β和γ定义为a、b和/或c的任意组合,和/或在下面的a、b和/或c中所述参数值之间任意的中间组合:
a)α≈6.4,β≈8.2,γ≈1.5,α≈5.0,β≈7.5,γ≈1.3;
b)α≈6.0,β≈7.5,γ≈1.5,α≈4.9,β≈7.0,γ≈1.2;和/或
c)α≈5.5,β≈7.5,γ≈1.5,α≈48,β≈7.0,γ≈1.2。
35.一种根据在预定持续时间内收集的BG数据估计患者严重低血糖症(SH)短期概率的***,所述***包括:
数据库组件,其用于保存识别所述BG数据的数据库;和
处理器,其被编程用于:
根据所述收集的BG数据计算scale数值;和
为每个BG数据计算低BG危险值(RLO)。
36.根据权利要求35的***,其中:
所述RLO(BG)的计算被数学地定义为:
scale=[ln(BG)]1.0845-5.381,其中BG以mg/dl为单位测量;
Risk=22.765(Scale)2
如果(BG小于大约112.5)则:
RLO(BG)=Risk,否则
RLO(BG)=0。
37.根据权利要求35的***,其中
所述RLO(BG)的计算被数学地定义为:
Scale=[ln(BG)]1.026-1.861,其中BG以mmol/l为单位测量;
Risk=32.184(Scale)2
如果(BG小于大约112.5)则:
RLO(BG)=Risk,否则
RLO(BG)=0。
38.根据权利要求35的***,其中
根据收集的BG数据计算LBGI,所述LBGI的计算由在时间点t1,t2,...,tn采集一系列BG读数x1,x2,...,xn数学地定义:
其中:
lbgi(BG;a)=RLO(BG)。
39.根据权利要求35的***,其中
根据所述收集的BG数据计算临时LBGI,所述临时LBGI的计算由如下定义的数学公式数学地定义:
LBGI(1)=RLO(X1)
RLO2(1)=0
LBGI(j)=((j-1)/j)*LBGI(j-1)+(1/j)*RLO(xj)
RLO2(j)=((j-1)/j)*RLO2(j-1)+(1/j)*RLO(xj)-LBGI(j))2
40.根据权利要求39的***,其中
计算SBGI,所述SBGI的计算能够如下地数学定义:
41.根据权利要求40的***,包括确认或者发出将发生短期SH的报警,其中如果:
(LBGI(150)≥2.5且(LBGI(50)≥(1.5*LBGI(150)且SBGI(50)≥SBGI(150)),则确认或者发出所述报警,或者
RLO≥(LBGI(150)+1.5*(SBGI(150)),则确认或者发出所述报警;
否则:
不需要确认或者提供报警。
42.根据权利要求40的***,包括确认或者发出将发生短期SH的报警,其中如果:
(LBGI(n)≥α且LBGI(n)≥(β)),则确认或者发出所述报警,和/或
(RLO(n)≥(LBGI(n)+γ*SBGI(n))),则确认或者发出所述报警;
否则:
不需要确认或者发出报警,其中α、β和γ是阈值参数。
43.根据权利要求42的***,所述阈值参数α、β和γ被定义为α≈5,β≈7.5,γ≈1.5。
44.根据权利要求42的***,其中所述阈值参数α、β和γ定义为a、b和/或c的任意组合,和/或在下面的a、b和/或c中所述参数值之间任意的中间组合:
a)α=6.4,β=8.2,γ≈=1.5,α=5.0,β=7.5,γ=1.3;
b)α=6.0,β=7.5,γ=1.5,α=4.9,β=7.0,γ=1.2;和/或
c)α=5.5,β=7.5,γ=1.5,α=4.8,β=7.0,γ=1.2
45.根据权利要求42的***,其中所述阈值参数α、β和γ定义为a、b和/或c的任意组合,和/或在下面的a、b和/或c中所述参数值之间任意的中间组合:
a)α≈6.4,β≈8.2,γ≈1.5,α≈5.0,β≈7.5,γ≈1.3;
b)α≈6.0,β≈7.5,γ≈1.5,α≈4.9,β≈7.0,γ≈1.2;和/或
c)α≈5.5,β≈7.5,γ≈1.5,α≈48,β≈7.0,γ≈1.2。
46.一种根据在预定持续时间内收集的BG数据估计患者严重低血糖症(SH)短期概率的***,所述***包括:
BG获取机构,所述获取机构用于从患者获取BG数据;
数据库组件,其用于保存识别所述BG数据的数据库;和
处理器,其被编程用于:
根据所述收集的BG数据计算scale值;和
为每个BG数据计算低BG危险值(RLO)。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US40297602P | 2002-08-13 | 2002-08-13 | |
US60/402,976 | 2002-08-13 | ||
US60/478,377 | 2003-06-13 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB038240092A Division CN100466965C (zh) | 2002-08-13 | 2003-08-08 | 用于处理自我监测血糖(smbg)数据从而提高糖尿病患者自我管理的方法、***和计算机程序产品 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101079082A true CN101079082A (zh) | 2007-11-28 |
Family
ID=35306237
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNA2007101032423A Pending CN101079082A (zh) | 2002-08-13 | 2003-08-08 | 用于处理自我监测血糖(smbg)数据从而提高糖尿病患者自我管理的方法、***和计算机程序产品 |
CNB038240092A Expired - Lifetime CN100466965C (zh) | 2002-08-13 | 2003-08-08 | 用于处理自我监测血糖(smbg)数据从而提高糖尿病患者自我管理的方法、***和计算机程序产品 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB038240092A Expired - Lifetime CN100466965C (zh) | 2002-08-13 | 2003-08-08 | 用于处理自我监测血糖(smbg)数据从而提高糖尿病患者自我管理的方法、***和计算机程序产品 |
Country Status (1)
Country | Link |
---|---|
CN (2) | CN101079082A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102549587A (zh) * | 2009-08-10 | 2012-07-04 | 糖尿病工具瑞典股份公司 | 用于处理血糖数据的装置和方法 |
CN102802508A (zh) * | 2009-06-26 | 2012-11-28 | 霍夫曼-拉罗奇有限公司 | 用于根据血糖的结构化单点测量结果来提供估计真实平均血糖值和估计糖化血红蛋白(HbA1c)值两者的方法、***和计算机程序产品 |
CN102946804A (zh) * | 2010-06-18 | 2013-02-27 | 霍夫曼-拉罗奇有限公司 | 具有考虑与生物标志相关联的***噪声的经过调节的退出标准的胰岛素优化***和测试方法 |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101785701B (zh) * | 2009-01-23 | 2015-10-07 | 理康互联科技(北京)有限公司 | 信息处理的装置、方法、设备及试剂载体 |
CN106955391B (zh) | 2009-02-04 | 2020-11-20 | 赛诺菲-安万特德国有限公司 | 用于提供用于血糖控制的信息的医疗设备和方法 |
JP5997453B2 (ja) * | 2011-04-25 | 2016-09-28 | アークレイ株式会社 | 情報処理装置およびユーザ端末 |
EP2723229B1 (en) * | 2011-06-23 | 2018-09-26 | University Of Virginia Patent Foundation | Unified platform for monitoring and control of blood glucose levels in diabetic patients |
US20140012511A1 (en) * | 2012-07-09 | 2014-01-09 | Dexcom, Inc. | Systems and methods for leveraging smartphone features in continuous glucose monitoring |
WO2014130841A1 (en) * | 2013-02-21 | 2014-08-28 | University Of Virginia Patent Foundation | Tracking changes in average glycemia in diabetics |
CN104183104A (zh) * | 2014-08-13 | 2014-12-03 | 上海电机学院 | 一种血糖测试结果的通信方法及*** |
US20160113594A1 (en) | 2014-10-22 | 2016-04-28 | Dexcom, Inc. | User interfaces for continuous glucose monitoring |
CN107003315B (zh) * | 2014-12-25 | 2018-09-25 | 株式会社日立制作所 | 胰岛素分泌能力分析装置、具备该装置的胰岛素分泌能力分析***以及胰岛素分泌能力分析方法 |
US11004566B2 (en) * | 2015-07-06 | 2021-05-11 | Abbott Diabetes Care Inc. | Systems, devices, and methods for episode detection and evaluation |
CN113450906A (zh) * | 2020-03-26 | 2021-09-28 | 苏州佳世达光电有限公司 | 临床检测诊断***及临床检测诊断方法 |
CN113505925B (zh) * | 2021-07-09 | 2022-07-15 | 重庆邮电大学 | 一种基于anfis的实验室危险化学品异常信息预警方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5206144A (en) * | 1985-03-29 | 1993-04-27 | Novo Industri A/S | Determination of glycated (glycosylated) hemoglobin in blood |
JP3150857B2 (ja) * | 1994-10-19 | 2001-03-26 | 富士写真フイルム株式会社 | グリコヘモグロビン含有比測定用分析要素及び分析方法 |
US5695949A (en) * | 1995-04-07 | 1997-12-09 | Lxn Corp. | Combined assay for current glucose level and intermediate or long-term glycemic control |
JPH11296598A (ja) * | 1998-04-07 | 1999-10-29 | Seizaburo Arita | 血糖値の予測システム及び予測方法並びにこの方法を記録した記録媒体 |
-
2003
- 2003-08-08 CN CNA2007101032423A patent/CN101079082A/zh active Pending
- 2003-08-08 CN CNB038240092A patent/CN100466965C/zh not_active Expired - Lifetime
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102802508A (zh) * | 2009-06-26 | 2012-11-28 | 霍夫曼-拉罗奇有限公司 | 用于根据血糖的结构化单点测量结果来提供估计真实平均血糖值和估计糖化血红蛋白(HbA1c)值两者的方法、***和计算机程序产品 |
CN102802508B (zh) * | 2009-06-26 | 2015-05-13 | 霍夫曼-拉罗奇有限公司 | 用于根据血糖的结构化单点测量结果来提供估计真实平均血糖值和估计糖化血红蛋白(HbA1c)值两者的方法、***和计算机程序产品 |
CN102549587A (zh) * | 2009-08-10 | 2012-07-04 | 糖尿病工具瑞典股份公司 | 用于处理血糖数据的装置和方法 |
US9336355B2 (en) | 2009-08-10 | 2016-05-10 | Diabetes Tools Sweden Ab | Apparatus and method for generating a condition indication |
CN102549587B (zh) * | 2009-08-10 | 2016-08-03 | 糖尿病工具瑞典股份公司 | 用于处理葡萄糖值的装置和方法 |
US9533097B2 (en) | 2009-08-10 | 2017-01-03 | Diabetes Tools Sweden Ab | Apparatus and method for processing a set of data values |
CN102946804A (zh) * | 2010-06-18 | 2013-02-27 | 霍夫曼-拉罗奇有限公司 | 具有考虑与生物标志相关联的***噪声的经过调节的退出标准的胰岛素优化***和测试方法 |
Also Published As
Publication number | Publication date |
---|---|
CN1688244A (zh) | 2005-10-26 |
CN100466965C (zh) | 2009-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101079082A (zh) | 用于处理自我监测血糖(smbg)数据从而提高糖尿病患者自我管理的方法、***和计算机程序产品 | |
CN1423789A (zh) | 自动诊断***和方法 | |
Patrick et al. | Genetic signature to provide robust risk assessment of psoriatic arthritis development in psoriasis patients | |
Crawford et al. | Prevalence of obesity, type II diabetes mellitus, hyperlipidemia, and hypertension in the United States: findings from the GE Centricity Electronic Medical Record database | |
CN1679443A (zh) | 健康管理*** | |
Noble et al. | Risk models and scores for type 2 diabetes: systematic review | |
CN1862551A (zh) | 治疗效果预测***、治疗效果预测方法及其计算机程序产品 | |
Wu et al. | The association of psoriasis with autoimmune diseases | |
Leong et al. | Systematic review and meta-analysis of validation studies on a diabetes case definition from health administrative records | |
CN1902650A (zh) | 医疗信息电子化***、程序以及记录媒体 | |
Pantalone et al. | The probability of A1C goal attainment in patients with uncontrolled type 2 diabetes in a large integrated delivery system: a prediction model | |
CN113614850A (zh) | 预测血糖浓度 | |
CN1227645A (zh) | 迭代问题求解技术 | |
CN1102036A (zh) | 脉波分析装置和使用该装置的诊断装置 | |
Karpati et al. | Patient clusters based on HbA1c trajectories: a step toward individualized medicine in type 2 diabetes | |
CN1244104A (zh) | 病人生理参数的远程监控装置和*** | |
CN1472536A (zh) | 分析物检验设备 | |
CN101077307A (zh) | 生物器官功能模拟*** | |
CN1516562A (zh) | 投药量决定支援装置、注射器及健康管理支援*** | |
CN1660012A (zh) | 血糖值测定装置 | |
Shi et al. | Application of back propagation artificial neural network on genetic variants in adiponectin ADIPOQ, peroxisome proliferator-activated receptor-γ, and retinoid X receptor-α genes and type 2 diabetes risk in a Chinese Han population | |
CN1297343A (zh) | 耳式妇女体温计 | |
Liu et al. | Regulation and competition in the Taiwanese pharmaceutical market under national health insurance | |
Li et al. | Integration of genetic and clinical information to improve imputation of data missing from electronic health records | |
US11762950B1 (en) | Automatic and blind segmentation of diverse data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20071128 |