专利名称：：改善连续式葡萄糖传感器的准确度的制作方法
技术领域：
：本发明涉及葡萄糖监测(监控)的领域，并且更具体地涉及用于连续葡萄糖监测(监控)的方法和系统。
背景技术：
：i见有的i正才居，如才艮寺居美国泮唐尿病十办会(AmericanDiabetesAssociation)的"全国泮唐尿病调查(NationalDiabetesFactSheet)，，，表明目前在美国大约1820万人患有糖尿病；并且在美国，糖尿病是死亡的第六主要原因。在2000年出生的每三个美国人中有一人将发展成H型4唐尿病。由于大量的糖尿病患者，以及由于预期会增加的糖尿病发病率，因此对于用于监测葡萄糖水平的准确葡萄糖监测系统存在持续增长的需要。连续式葡萄糖传感器(continuousglucosesensor)一皮i殳计成不^f又在单个时间点4是供实时葡萄^唐水平，而且基于每一确定的时段(具有最少的指刺(finger-stick))进行的分析来提供人葡萄糖水平的趋势，这导致改善的血糖(glycemic)/if唐尿病控制。然而，由于难以在动力永或血管中直4妄测量血4唐(下文称作BG)，最现代的连续式葡萄糖传感器(下文称作CGS)是通过对在间质液中的间质葡萄糖(下文称作IG)进行取样来产生或获得血糖估计(值)、而不是对BG进行取样来产生血糖估计(值)。来自IG的典型葡萄糖(BG)估计产生自至少两个连续近似步骤1)血液至间质葡萄糖(BG-至-IG)输送；以及2)从由传感器记录的IG相关的电流导出BG值。因此，虽然CGS技术已取得显著的进展，但就校准、灵敏度、稳定性、以及血液与间质葡萄糖浓度之间的生理时滞而言，准确和可靠的CGS装置的发展继续面临许多挑战。BG和CGS读ft之间的差异起因于以下主要因素生理过程、传感器校准、噪声(干扰)、以及操作(设计，engineering)。生理时滞和梯度动态地随时间、随BG水平、以及交叉受治疗者(acrosssubject)而变化；并且IG的直接频繁的体内取样是极端困难的。因此，CGS的操作性能的评估留下一个中心问题分离起因于校准、传感器噪声、以及BG/IG梯度的BG/CGS误差部分。因此，期望一种方法和设备来改善CGS的准确度(精确性)和可靠性。
发明内容基于本4皮露内容，将明了本发明的优选具体实施方式的各种目的和优点。才艮据优选的具体实施方式，通过改善CGS传感器的才交准或纠正由生物时滞引起的误差或其组合，本发明可改善CGS的准确度和可靠性。为本发明的一种示例性具体实施方式，本文纟皮露了一种用于改善连续式葡萄糖传感器(CGS)的准确度的方法。该方法包括在第一时间才交准CGS;以及在第二时间变4匕CGS才交准，该第二时间是基于动态监测的CGS<i、CGS变4匕的速率和予贞定准则力o以确定的。作为本发明的另一种示例性具体实施方式，本文4皮露了一种用于改善连续式葡萄糖传感器(CGS)的准确度的方法。该方法包括利用第一血糖数据和不同于第一血糖凄t据的第二血糖数据对CGS进行校准。作为本发明的又一示例性具体实施方式，本文4皮露了一种连续式葡萄糖传感(CGS)装置(或连续式葡萄糖检测装置)。该装置包括第一装置，用于测量间质葡萄糖水平以便获得CGS输出；以及可访问CGS输出的校准模块，用于改善CGS的准确度，进一步包括可访问CGS输出的监测才莫块，用于动态地监测CGS和CGS的时间导数；以及基于动态CGS值、CGS值的时间导数和预定准则来指导(指示)另一个校准事件。作为本发明的又一种示例性具体实施方式，披露了一种计算机可读介质，该介质具有计算才几可执4亍指令以实施一种用于改善连续式葡萄糖传感器的准确度的方法，其中所述方法包括取得(retrieving)最初血糖值和在最初血糖值的测量中获得的CGS值；随时间监测CGS值和CGS值的时间导lt;基于监测的CGS值和CGS值的时间导凄t确定是否开始另一个才交准；以及如果确定开始所述的另一个4交准，则4交准CGS。作为本发明的又一种示例性具体实施方式，#皮露了一种计算枳^可读介质，该介质具有计算机可执行指令以实施一种用于改善连续式葡萄糖传感器的准确度的方法，其中所述方法包括在第一时间取得血糖值和在最初血糖值的测量中获得的CGS值；以及在第二时间才交准该CGS，该第二时间是通过基本上在该第二时间的一个CGS值、该CGS值的时间导数和预定准则来确定的。作为本发明的又一种示例性具体实施方式，本文4皮露了一种用于治疗与血糖有关的疾病的系统。该系统包括连续式葡萄糖装置；用于将CGS值递送到疾病治疗中心的装置，其中疾病治疗中心能够发出相应的治疗指令或采耳又相应的治疗4亍为(treatingaction)。尤点可以(在某些优选的实施例中)借助于本文所附的独立权利要求的特征来实现。在从属权利要求中进一步阐述了另外的优选具体实施方式。根据35U.S.C.§112,第六款，在权利要求中，仅由词组"用于...的装置"表示的元件用来解释为装置加功能性的权利要求。才艮据以下连同附图的详细描述，可以最好地理解本发明的优选具体实施方式，其中图1示意地图解说明了一种利用才艮据本发明一个实施例的连续式葡萄糖传感器来测量血糖水平的方法；图2是示出了包括图1的校准仪的连续式葡萄糖传感器的分解图3是图解说明了根据本发明一个实施例的校准仪和时滞校正器的功能模块的示意图4a是根据本发明一个实施例的校准模块的方块图；图4b是根据本发明一个实施例的校准过程的流程图5图解说明了其中在校准期间的血糖差异以及通过本发明的示例性方法改善的CGS的准确度；图6是示出了利用本发明的时滞才交正方法为改善CGS的准确度所执行的步骤的流程图7是示出了根据本发明的又一个实施例为校准CGS所执行的步骤的流禾呈图8是示出了利用本发明的一种示例性方法的改善的CGS测量结果的曲线图9是示出了利用本发明的一种示例性方法的改善的CGS测量结果的另一曲线图IO是示出了利用本发明的一种示例性方法的改善的CGS测量结果的又一曲线图11是示出了^4居本发明的又一个实施例为改善CGS的准确度所执行的步骤的流程图12是图解i兌明了一种其中可以实施本发明的实施例的系统的示意图；以及图13是示出了具有计算机可读指令的示例性计算装置的示意图，在该计算装置中可以实施本发明的实施例。具体实施例方式通过改善CGS的校准或通过校正起因于BG与IG之间的生理时滞的误差或它们的组合，本发明提供了一种用于改善连续式葡萄糖传感器的准确度的方法和装置。考虑到在本发明的精神范围内的许多可能的变型，在下文将参照特定的实施例来讨论本发明。然而，本领i或4支术人员应当明了，以下讨i仑是用于i兌明目的，并且不应当解释为限制本发明。在不偏离本发明的精神的情况下，其它变型也是可用的。大多数目前的CGS装置的不准确性(不准确度)主要归因于不良的CGS校准、生理时滞、以及随机误差。为了降低CGS的不准确性，提出了一种改善的校准程序。也在本发明中提出了，可以可替换地通过校正与生理时滞有关的误差来实现降低不准确性。事实上，改善的校准程序和时滞校正程序可以可替换地结合在一起以便获得更好的CGS性能。参照图1，其中示意性地图解i兌明了才艮据本发明的一个实施例测量体内葡萄冲唐水平的方法。由于在血管或动月永(血液)102中直接测量葡萄糖水平非常困难，所以CGSIOO通过测量与血管或动脉相互作用的间质液104中的葡萄^f唐7jc平来4企测血斗唐水平；并JM寻CGS的输出值与血糖水平联系起来(该过程经常称作校准)。这种间接测量基于已证明的基础，即血糖水平与间质液中的葡萄糖水平共同变化。为了改善CGS的准确度，如图2中所示，才是供了准确度改善^f义(才交准4义，accuracyimprover)106。通过经由最^f圭"尤4匕)才交准建议来改善CGS4交准(其将在下文讨论)、或通过校正与生理时滞有关的误差(其将在下文讨论)、或通过其组合，准确度改善4义可以改善CGS的准确度。虽然图2示出了准确度改善4义包括在CGS(100)中作为功能构件，4旦并不需要如此。在其它实施例中，准确度改善仪可以是与CGS分开的独立模块或方法。具体地说，为了改善CGS的准确度，可以在与CGS输出相连接的装置中实施(实现或设置，implement)准确度改善仪。此夕卜，可以以存储在计算机可读介质中的计算机可执行代码的序列的形式来实施准确度模块；或可以在装置的石更件中加以实施，其将随后详述。作为本发明的一个实施例，图3示意性地图解i兌明了图2中的准确度改善仪(106)的分解图。在该特定的实施例中，准确度改善仪(106)包括CGS校准模块108和时滞校正模块110。通过执行最佳校准建议循环(optimalcalibrationrecommendationcycle)，CGS校准模块可指定用于改善CGS的校准，从而改善CGS的总准确度。时滞校正才莫块(110)被指定用来校正由血糖水平与间质葡萄糖水平之间的生理时滞引起的误差。取决于如上所讨论的准确度改善4义的具体功能，准确度改善仪可以仅具有CGS校准模块和时滞校正模块其中之一或两者均具有。改善的CGS校准在本领域已知的，CGS校准的准确度取决于血糖(下文称作BG)变化的速率和在校准时刻(f)的BG值(SGW)。BG变化的速率可以数学上表示为BG(^)的时间导数《BG(U)/力。鉴于下述事实带有变量输入(variantinput)的4交准好于那些借助于单一或非变化的输入的校准，所以本发明的CGS校准使用变量输入。作为一个实施例，图4a是示出了图3的示例性CGS4交准的功能模块(118)的示意图。参照图4a，CGS校准仪118包括初始校准模块lll、存储模块113、再赋值模块(reassigningmodule)115、以及决策环才莫块117,其中决策环冲模块117进一步包4舌监测模块119、内部请求模块121、自检启动模块(自检起始模块)123、决定模块125、以及传感器4交准才莫块129。初始纟交准才莫块进行最初4交准以使>人SG(^)的最初测量结果获得初始校准数据对SG(^)和SGW,存储模块113耦合于初始41准模块111的输出并存储初始校准数据对SG(^)和SG(^)。监测模块119连接于存储模块的输出并动态地监测CGS输出SG^和SGr"变化的速率5"G'"戶^SG(U)/力。内部请求模块121连接于监测模块的输出并管理用于校准的内部请求。自检启动模块123连接于内部请求模块的输出并被指定用于启动自检程序以用于最佳校准定时。与自检模块123的输出连接，决定模块125决定是否在特定时间通过传感器校准模块129进行另一次校准。在决定以后，来自校准的校准数据通过再赋值模块115被再赋值(再分配)到CGS。应该注意，实际上一个或多个上述功能才莫块可以合并到其它功能模块中。尤其是，传感器校准模块129可以合并到初始校准模块111中，以进行CGS校准。可以将监测模块119、内部请求模块121、自检模块123、以及决定模块125中的任何一个合并到组合功能模块中。在图4b中示意性地图解说明了图4a中的功能模块的示例性操作，以1更完成本发明的所期望的最佳才交准建议过程。参照图4b，通过在最初时间f测量血糖水平BG(^)来进行初始校准，以^f更获得最初CGS输出SG(^(步骤112)。可以通过图4a中的初始校准才莫块111来执行该步骤。在步骤114记录最初数据对SG(^)和5"GYW。可以通过图4a中的存储模块113来执行此步骤。可以在由CGS制造商推荐的时间或在由用户确定的时间冲丸4亍上述初始4交准。然后才交准程序进入校准决定开^成环(calibrationdecisionmakingloop)118，如图4a所示该纟交准决定形成环是通过决策环才莫块117来执行的。才交准决定形成环开始于步骤120，该步骤120为随时间监测传感器值SC7(^和变化速率SG丫"，其中5U丫^定义为5"C7(^的一阶时间导数，即SG，f>」=^SGW/&。可以通过图4a中的监测模块119进行SG，W的动态监测和推导。在特定时间^接受另一个校准请求(步骤122)以后，则在步骤124启动用于最佳校准定时的自检程序。此步骤触发(启动)一系列确定在此时是否应进行校准的才喿作，例如图4b中的126和128。接受校准请求并传送这种请求以启动自枱r禾呈序可以通过图4a中的内部i青求才莫块121来完成；并且可以通过图4a中的自冲企启动才莫J夹123来完成自冲企启动。应当注意到，可以在由制造商规定的时间、或可替换地在由用户(如医生以至适当的患者)少见定的时间，启动对于另一次才交准的内部请求。在步骤124启动自检程序以后，则在步骤126确定是否ISC丫^)1〈1。可以通过在图4a中的决定才莫块125来才丸4亍it匕决定步骤。如果I丫。』21，则程序回到步骤120以继续监测5UW和，^>值。否则，程序进行另一决定是否ISG(^)-5"G(^1大于《mg/d/入其中《mg/d/)是初始SG(^)与在时间^时的COS"输出SG^"之间的预定差异阈4i。例3口，Wwg/(^可以是10wg/d/或更高，如15wg/d/或更高，并且更优选30wg/^/或更高。在步骤126和128的一种或两种决定可以通过如图4a中所示的决定模块125来执行。在步骤128,如果|SG(^)-SG(^)1等于或小于^mg々/仏则程序回到步骤120。否则，在步骤130,通过例如图4a中的传感器校准模块129进行另一校准。基于在步骤130的再校准，在步骤116再赋予CGS校准值，例如，通过分别用再校准值5"G(^)和BC7(^)代替校准值^GY4)和5Gf^)。在步骤114存储再赋予的校准值。可以通过图4a中的再赋《直才莫块115来完成上述再!^[直。在再U武^直和i己录以后，4交准过禾呈再次进入决定形成环118,并且重复上述步骤114、120、122、124、126、128、130、以及116。可以通过由CGS制造商建议的默认校准数目、或可替换地由用户来确定校准循环的数目。例如，在CGS运行(寿命，life)的最初24小时期间，可以进4亍多次4交准循环，例如2至10次，或更通常3至4次4交准循环。利用如上所讨i仑的最佳才交准方法的改善的CGS的准确度可以通过以下实-睑lt据和计算机才莫拟来证实，如图5所示。实验数据为了测试体现如上所讨论的精确度改善方法的该CGS的准确度，对39位患有I型糖尿病(T1DM)的受治疗者进行了测量。39位参与者具有以下统计资料平均年龄为42.5岁，其中标准差(SD)为12(SD42)，T1DM的平均持续时间为21.6年(SD=94)，平均HbAlc=7.4%(SD=0.8)，16位为男性。该研究由弗吉尼亚大学IRB受治疗者(UniversityofVirginiaIRBSubjects)认可。受治疗者在研究之前的晚上被允许进入一般临床研究中心(综合临床研究中心，generalclinicresearchcenter,GCRC)。整夜将参与者的BG水平控制在100-150wg/c//(5.5-8.3wwo///)的j&4唐正常范固内。将MinimedCGMSTM连4妻于每位受治疗者并在研究期间按照制造商的说明进行校准。所有CGMStm被插入腹中。在早上实施高月夷岛素4计夹(hyperinsulmemicclamp)。每个钳夹使用了1mt//&g/m/"的恒定胰岛素输注速率以及可变的葡萄糖车命注速率以达到和维持BG7jc平在约110wg/d/(约6wwo///)。其后，降低葡萄糖输注速率以允许BG的受控下降(约1直到BG达到50wgA//(约2.8mwo///)。然后重新开始葡萄并唐l命注以恢复到正常葡萄糖水平。钳夹研究的血糖正常部分以70至210分钟的时间长短变化；以及BG降低操作的持续时间在30至60分钟范围内变动。恢复时间在30至60分钟范围内变动。通过将手加热至50。C来获得动月永化血液(arterializedblood),然后每5分钟耳又样，用于参照BG水平。考虑到胰岛素要达到其稳定状态效应，在开始输注以后的最初15分钟的数据被忽略。使CGMSTM读数与参照BG同步。传感器最佳校准的计算机模拟利用在上述4甘夹研究期间获得的两个参照BG值对传感器的再校准进行了计算机模拟，如图5所示，其描述了作为两个BG值之间的差异的函数的传感器误差。模拟的再才交准使用了CGMSTM的标准线性校准函数。对来自再校准的结果与在实验期间显示的传感器本身的精确度以及与利用所有可获得的参照BG值的"完美"；艮准进行了比较。参照图5，X轴表示以BG单位(wgA//)的在两个模拟校准点之间的距离；而Y轴表示利用两点校准的传感器输出的平均绝对误差(MAE)。在该图中可以看到，如果两个校准BG的值接近，则MAE会较高。当差异接近20wg/d/时，MAE快速下降，并且其后緩十曼下降。图中的上水平线表示传感器本身4交准的MAE;而下水平线表示来自利用所有可获得的参照点的"完美"校准的MAE。在图中还可以看到，数值大于30wg/d/但小于40wgA//的BG校准差异d获得极好的结果；而数值大于40wg/W的差异d获得"几乎完美的"结果。值得指出的是，上述在实验期间的传感器校准总是在稳定的BG(保持在血糖正常)期间完成的，因此BG变化的速率的影响为最小。生理时滞的4交正除CGS的校准之外，BG与IG之间的生理时滞也引起CGS输出的不准确性。这来自以下事实大多^t目前的CGS装置并不直接测量血糖水平，而是以在间质液中的IG水平来代替。然后CGS装置将IG读数转换成BG的估计值。因此，从IG到BG的改善的转换方法将导致改善的CGS性能。本发明的一个目的是通过包括IG与BG水平之间的生理时滞来改善从IG到BG的转换。这样的改善是通过分析和结合IG与BG之间的时间相关性(timedependence)来完成的。具体i也i兌，建立lt学才莫型以用来描述BG与IG或CGS输出之间的时间相关性。基于已建立的模型，导出了数学方程以定量地表示CGS输出对BG的时间相关性，即CGS是BG的函数。然后转换此方程，以Y更表示作为CGS的函lt的BG。因此，该反向方禾呈(inverted叫uation)可以用来预测对于纟合出CGS丰釙出值的BG水平。在应用中，将反向方程应用于原始CGS凄t据以产生准确的BG估计。数学模型鉴于葡萄糖是相对小分子的事实，所以认为，葡萄糖可以自由地扩散通过毛细管壁，如血管和脂肪组织。脂肪组织是高度血管化的；并且间质液占据细胞之间相对较薄的层。该事实意味着，不存在远离细月包表面的体积、元素(volumeelement)，也不存在远离毛细管壁的体积元素。因此，葡萄糖在间质液中的摄取和扩散可以假设为在拓朴上是相对均匀的。才艮据本发明的IG和BG的运输行为描述在图1中。再次参照图1，葡萄糖在间质液与血管(或脂肪组织)之间的运输行为可以作为扩散来建立沖莫型。在间质液中的IG还经受消摔毛，其导致量和/或浓度降低。为了导出数学扩散方程，假定所考虑的特定局部间隙环境并不显著地对BG/时间曲线的展开有贡献，因此，BG水平SG(^)的时间相关性独立地逐渐形成，并且可以处理为系统中的外源变量。这种假设是特别安全的，尤其在高胰岛素钳夹的情况下，其中BG水平受到葡萄糖IV输注的最大控制。进一步假设，葡萄糖的摄取按照IG独立途径、或4姿照由米-门氏(Michaelis-Menten)动力学描述的途径，如分别表示在方程la和lb中/邵)方程la卿一d7邻)一磨一,卿爐方程lb在上述方程中，oc是每单位体积每单位时间内葡萄糖的摄取。应当注意到，方禾呈la和lb描述了明确的时间相关性。其它变量，如胰岛素水平、锻炼等(其可以直接影响葡萄糖摄取SG和/GO并没有从方程中排除。《w是方程lb中的常数；并且它并不引入另夕卜的可拟合(fitable)参fc实际上，ATw可以采用那些针对GLUI-4的；舌寸生戶斤发表的<直，3口在ZierlerK的"P^7w/e6cx^vg/wcosemd<36o//sm(全身葡萄糖代谢)"，J尸/yw'o/.276:E409-E426,1999中所陈述的，将其主题的全部内容以引用方式结合于本文。通过参照菲克定卿方程2其中P是毛细管壁对葡萄糖的渗透性。因为在间质液中没有葡萄糖的其它明确的来源或渗透(sink),所以葡萄糖的净变化可以通过加入方程la和lb来导出，其可以表示为以下方程3a和3b，其中方程3a对应于均匀摄取扩散模型，而方程3b对应于米-门氏动力学模型。<formula>formulaseeoriginaldocumentpage24</formula>方程3b方程3a和3b的数学解方程3a是具有解析解的常孩吏分方程；而方程3b是需要数值才莫拟的第二型阿贝耳(Abel)方程的非线性微分方程。方程3a的解析解表示在以下方禾呈4中<formula>formulaseeoriginaldocumentpage24</formula>方程4通过，支定a和(3是随时间恒定的，利用如表示在以下方程中的A-t表示法，方禾呈4可以^皮简^f匕成以下方禾呈6:<formula>formulaseeoriginaldocumentpage24</formula>方程<formula>formulaseeoriginaldocumentpage24</formula><formula>formulaseeoriginaldocumentpage25</formula>方程6通过在方程6中消除5G(^)的高阶导数项，于是方程6可以被简化成一种形式，而基于卡尔曼滤波/平滑技术的卡尔曼递归分析(Kalmanrecursionanalysis)贝'J可以应用于it开j式。这才羊的才支术的一个实例在"0/"w"/Cow/ro/(最佳4空制和估计)，，NewYork:DoverPublications,199437中陈述，将其主题的全部内容以引用方式结合于本文，其中^f吏用了"状态空间模型"和递归方法。该递归方法附在本文的附录B中。利用"状态空间;溪型"并借助于CGS读数被小时间《均匀隔开的假设，在以下方程9中可以描述用于CGS输出和IG的相关性的状态空间模型。<formula>formulaseeoriginaldocumentpage25</formula>用于包括BG展开的CGS输出的状态空间模型可以表示为以下方禾呈10。<formula>formulaseeoriginaldocumentpage26</formula>方程10用于包括BG展开和线'性才殳影(linearprojection)的CGS输出的状态空间才莫型可以表示为以下方程11。<formula>formulaseeoriginaldocumentpage26</formula>方程11给出消耗、渗透性、/G、以及的变化速率的估计值，则可以类似地进行方程3a和3b的反演。3a和3b的反向方程分别表示为以下方禾呈12a和12b:方程12a<formula>formulaseeoriginaldocumentpage26</formula>方程12b<formula>formulaseeoriginaldocumentpage26</formula>方程12a和12b表明，使用CGS对提供/G的变化速率的精确估计变得很重要。反向方程12a和12b的表达也是可能的，其表示为以下方禾呈13和14。<formula>formulaseeoriginaldocumentpage27</formula>方程13<formula>formulaseeoriginaldocumentpage27</formula>方程14应当注意到，在上述才莫型中的观测值是在方程7中定义的5G^」的函数。如果接受多项式平滑/插值公式来描述BG(t)的进程，那么它也可以净皮线性反演，如所附的附录C中所示。算法实现因而，上述数学才莫型和方程可以应用于CGS读数，以便利用CGS^T出通过预测BG水平来才交正/GW与SGW之间的生理时滞。本发明的一个示例性过程呈现在图6的流程图中。参照图6,在步骤132展开如上所讨论的数学模型。该模型描述了5G^、/G^、以及CGS输出5"GY^之间的动力学。基于在步骤132展开的模型，第一数学方程被简化，以便求出作为5G(^的函数的SGf^)。该第一方程优选利用菲克扩散定律和米-门氏假设加以展开。然后反演该第一方程以便在步骤136估计作为CGS输出的函数的5GW的动力学。然后反向方程可以应用于CGS原始数据以在步骤138改善CGS读数的精确度。在步骤138反向方程对于原始CGS凄t据的示例性应用图解"i兌明于图7中。参照图7，在步骤140获得原始数据集。数据集包括CGS输出SGfV、取样时间7，mfV、以及参数Par(a、/、戎凑)，其中参翁:戎雀是在方程9中表示为我凉的参数。在步骤142优选对原始CGS凄丈据进4于预处理(虽然不需要)。具体i也i兌，通过初始回归和/或滤波对CGS原始lt据进行处理以平滑CGS原始lt据。该步骤是重要的，因为研究表明，原始CGS数据目前受到噪声和随机尖峰(信号)的影响，其中噪声和随机尖峰(信号)需要滤掉以便改善反向才莫型方程的精确度。存在许多方式来回归和/或过滤CGS数据。例如，可以基于临床7见测的BG变化速率结合贝叶斯估计一支术(Bayesianestimationtechnique)来对原始CGS数据进行过滤。在另一实施例中，可以用方程10和11的基于卡尔曼滤波器的方法来处理原始CGS输出，其借助于BG的展开满足其限制的假设来产生最佳估计。在给出预处理的CGS数据(或直接的原始CGS数据而没有上述预处理)的情况下，在步骤144计算CGS输出的变化速率(时间导数)。CGS输出的独特特点是能够估计它们的输出的导数。然而，由于观测值和系统噪声以及漂移的(wandering)传感器灵敏度，观测到，这样的即时原始估计很少是最佳的。相反，利用指数权(exponentialweight)的短区间多项式光滑技术可以产生更好的结果，如实验所指出的。然后，在步骤146将IG比例系数应用于CGS,接着在步骤148应用如上所述的反向模型来估计BG水平。软件实现如参照图7所描述的过程可以以许多方式加以实施，其中之一是以软件形式。例如，上述过程可以用实时版本加以实施，其尤其可直接应用于CGS输出以便将原始CGS输出转换成BG估计值，并且产生实时结果。本文附带了实时实施程序代码的集作为附录D。可^,换;也，可以4乍为追溯片反本(retrospectiveversion)来实施该过禾呈，其中追溯片反本尤其可用于进行CGS4青确度的追溯改善，并且可应用于并不显示实时邀:据的CGS。本文附带了用于追溯实施的程序4戈码集作为附录E。应当注意到，为了便于参考，消除了每个行号的起始点。采用了参数的线性化版本，以有利于计算速度。本领域技术人员应当明了，在附录D和E中呈现的程序代码是基于说明的目的，并且不应当解释为限制本发明。在不偏离本发明的精神的情况下，还可以应用许多其它片反本。^寸方'法进^亍i式马全一项研究期间获得的数据加以评估，其中所述研究是在弗吉尼亚大学一^殳临床研究中心(UniversityofVirginiageneralClinicalResearchCenter(GCRC))进行的，其是正在进行的NIH研究基金(ROlDK51562，PrincipalInvestigatorBorisKovatchev)的土曾力口的工页目。i亥i曽力口的研究由AbbottDiabetesCare(P.I.WilliamClarke)资助以在两种CGS(AbbottNavigatorTM和MinimedCGMS)之间进4亍直4妻比專交。才莫型的展开和试验属于增加的研究的目的。用于研究的受治疗者16位患有TIDM的受治疗者参与了研究，其中ll位为男性，5位为女性，年龄为42岁(标准差(SD)为3年)，糖尿病的持续时间为20年(SD为3年)。/人每一位受治疗者获得知情同意。在体检以后，在研究之前的晚上受治疗者被许可进入一般临床研究中心。按照制造商的说明以及如建议的那样进行校准，在开始数据记录以前，爿寻一种CGS系统，即FreestyleNavigator丁M应用于每4立受治疗者约12小时。所有系统被插入腹中。直到第二天早上没有进行BG参考与CGS的比较。研究方案定义为，连续两天进行相同的高月夷岛素4计夹。每天高月夷岛素4甘夹^f吏用了恒定的40mt/Ag/附/w的月夷岛素l俞注速率以及可变的葡萄纟唐，lr注速率以达到和维持BG水平在约110wg/d。其后，降低葡萄糖输注速率以允许BG水平的受控下降(约1wg/W/m/")直到BG水平达到40mg/d/。该研究的血^唐正常4甘夹部分以70至210分钟的时间长短变化，而BG降低操作的持续时间在30至60分钟范围内变动。每5分钟采集动脉化血液并利用BeckmanGlucoseAnalyzer(BeckmanInstruments,Inc,Fullerton,CA)确定参考BG水平。每分4中i己录FreestyleNavigator葡萄泮唐读凄t并与参照BG同步，其中4青确度为30秒。以5分钟的间隔计算参照和NavigatorTM速率以及BG变化的方向。这种程序为研究中的16位参与者产生29个钳夹数据集。用于分析的软件利用R2.1.1进行数值分析，其是开放原始码的免费程序设计^吾言并且适合于统计分冲斤(http〃www.r-project.org)。除了基本禾呈序包以外，还使用了"odesolve"、"fields"、和"dsel"程序包以及它们与CRAN储存库的相关性，另外MicrosoftExcel用来产生图表。结果将方程12a应用于未过滤的NavigatorTM原始数据，其中参数是通过非线性最小二乘法而获得的。每个凄t据行程(lt据运行，datarun)在下降到^f氐血4唐的起点处开始。表1示出了所有29个4甘夹事件结果的总结。可以看到，Navigator的平均RMS误差减少3倍(fold)以上，并且o/。RMS误差减少5倍以上。此夕卜，通过该才莫型改善了参照与传感器BG之间的相关性表1<table>tableseeoriginaldocumentpage31</column></row><table>图8至图10描绘了这些结果。具体地r说，图8示出了29个事件(以5分钟数据间隔在最低点比较)的平均值。实心菱形符号是由Beckman分才斤4义记录的参照BG。实心正方形是Navigator的凄t据。空心正方形是通过本发明的方法校正的Navigator的数据。在图中可以看到，模型校正数据比原始的Navigator输出更接近参照BG。图9和图10示出了具有l分钟数据间隔的两个个体患者。在两种情况下，通过本发明的方法校正Navigator数据(空心正方形)导致参照BG的改善的跟踪(tracing)。组合的准确度改善方法和装置如上所讨i仑的，通过本发明的一个实施例经由改善的4交准方法可以改善CGS输出的准确度。可替换地，还可以通过校正BG与IG之间的生理时滞来改善CGS的准确度。在本发明的另一实施例中，可以结合上述两种4交正方法以《更进一步改善CGS准确度。图11图解说明了才艮据本发明的另一个实施例用于进行组合准确度改善过程的流程图。参照图11，组合过程开始于4交准CGS，其中利用如上文参照图4所描述的方法。可以消除或减小系统误差。在步-骤152,通过滤波和/或平滑校准的CGS输出，对来自校准的CGS的CGS输出进行预处理。可以利用如在图7中的步骤142所描述的相同或不同的方法执行此步骤。应当注意到，虽然优选^f旦并不需要此步驶《。由于该预处理，在CGS输出中的传感器噪声和/或较小的随机波动以及一阶时间导数可以消除或简化。在步骤144计算CGS输出SG(^)的一阶时间导数。然后，将IG比例系数应用于CGS输出，接着应用反向才莫型方程以估计BG水平，如参照图6所描述的，在这里将不重复。中通过比较由本发明的方法所获得的CGS输出和本领域中典型CGS(没有采用本发明的方法)的CGS输出。获自本领域典型CGS的CGS车命出可参见"Evaluatingtheaccuracyofcontinuousglucosemonitoringsensors:continuousglucoseerrorgridanalysisillustratedbytherasensefreestylenavigatordata,"来自BKovatchev，LGonderFrederick,DCox,和WClarke,飾6etoC脏，vol.27,pp1922-1928，2004。<table>tableseeoriginaldocumentpage32</column></row><table>表2的小组A提供了在上述钳夹研究期间MinimedCGMStm的准确度的连续式葡萄4唐误差删才各分冲斤(continuousglucoseerror-gridanalysis)(CG-EGA),分层为低血糖和血糖正常。临床上准确的传感器读数在低血糖期间为50.0%而在血糖正常期间为96.4%。这些百分比之间的较大差异主要起因于要求更高的临床准确度标准。对于低血糖事件虽然对于稳定的血糖正常状态对于传感器误差存在车交大的临床允i午i吴差(clinicaltolerance),^f旦是在临床上危险的和快速发展的状态期间，如低血糖，则希望传感器满足更高的标准以侵_提供用于适当的和及时的治疗决定的准确反馈。CG-EGA反映了这种区别。另夕卜，MAE和平均绝只十百分i吴差(MAPE)包4舌在表1中并且也由BG范围分层。在表2中的小组B提供了由两个参照BG再校准的传感器的CG-EGA、MAE以及MAPE，其中两个参照BG相隔是30wg/J/(例如，微分d是30mg/dl)，其是在研究的BG范围内的临床上合理的樣爻分。可以看到，CG-EGA准确读数的百分比从50%增加到86,7%,而在4氐血3f唐期间MAE乂人27.9wg/d/减小到10.9mgA/Z。在血4唐正常期间也7见察到MAE和MAPE的改善。表2的小组C提供了在传感器再校准以后估计的SIG与BG的CG-EGA、MAE以及MAPE。可以看到，在BG波动以后SIG的"准确度"较高，几乎100%，其意味着用于潜在的传感器准确度的极好的理i仑限度。可以以i午多方式来实施本发明的实施例。例30，它可以实施为连续式葡萄糖传感器的功能构件、或可以实施为与连续式葡萄糖传感器有关的独立的才莫块。在<壬<可一种情况下，本发明的实施例可以实施为安装在计算装置中的软件的一组程序代码、或耦合于连续式葡萄糖传感器的硬件装置的一组计算机可执行代码。与实施介质无关，本发明的实施例可以与单个连续式葡萄糖传感器有关，用于改善单个葡萄糖传感器的准确度。可替换地，可以以这样的方式来实施本发明的实施例，即实时CGS数据、连同相关的误差和误差校正参tt和凄t据可以^皮传送到误差处理中心。该传输可以或可以不带有葡萄糖数据。这样，可以建立通用化的连续式葡萄糖传感系统。图12扭克略;也图解i兌明一种示例性系统，在该系统中可以实施本发明的实施例。参照图12，临床机构(clinicsetup)158提供了用于医生(例如164)i貪断患有与葡萄并唐有关的疾病的患者(例如160)的地方。连续式葡萄糖传感器(或结合有葡萄糖测试功能的传感装置)162可以用来监测和/或测试患者的葡萄糖水平。这样的监测和/或测试可以是短期的(例如门_&就-〖貪)或长期的(例如病房停留或家庭)。该连续式葡萄糖传感器在其中结合有如上所讨论的准确度改善方法的实施例。具有改善的准确度的CGSllT出可以由医生用于适当的操作(action),如对患者进行胰岛素注射或食物供给、或其它适当的操作。可替换地，具有改善的准确度的CGS输出可以-波传送至计算才几纟冬端168,用于即时或未来分析。该传送可以通过电缆或无线的或4壬何其它适宜的方法。来自患者的具有改善的准确度的CGS输出还可以被传送至便携式装置，如PDA166。具有改善的准确度的CGS输出可以被传送至葡萄糖监测中心172，用于处理和/或分片斤。这才羊的传送可以通过"i午多方式来完成，如网全各连*接170，其可以是有线的或无线的。除具有改善的准确度的CGS输出以外，还可以将误差、用于准确度改善的参数、以及任何准确度相关的信息传送至如计算机168、和/或凝:据处理中心172,用于进行误差分析。这可以为葡萄糖监测中心4是供集中化准确度监测和/或准确度增强，这是起因于葡萄糖传感器的重要性。如早先所述，还可以在与目标连续式葡萄糖传感器有关的独立的计算装置中实施本发明的实施例。在图13中示意性地图解说明了一种示例性计算装置，在该计算装置中可以实施本发明的实施例。虽然这样的装置对于本领域技术人员来说是众所周知的，但为了方便其他读者，本文将提供简要说明。参照图13,在其最基本的构造中，计算装置174通常包括至少一个处理单元180和存储器176。取决于计算装置的精确构造和类型，存4诸器176可以是易失的(如RAM)、非易失的(》口ROM、闪存等)或两种的某种组合。另外，装置174还可以具有其它特征和/或功能性。例如，该装置还可以包括另外的可移动式存储器和/或不可移动式存储器，其包括但不限于磁盘或光盘或磁带、以及可写电存储介质。这样的另外的存储器表示为可移动式存储器182和不可移动式存储器178。计算机存储介质包括易失性介质和非易失性介质、可移动式介质和不可移动式介质，这些介质是以任何方法或技术来实施的，用于存储信息如计算机可读指令、数据结构、程序模块或其它数据。存储器、可移动式存储器以及不可移动式存储器都是计算机存储介质的实例。计算才几存4诸介质包括4旦不限于RAM、ROM、EEPROM、闪存或其它存卞者才支术、CDROM、凄t字多用途光石乘(DVD)或其它光学存储器、盒式磁带、磁带、磁盘存储器或其它磁存储装置、或可以用来存储期望的信息并且可以由装置访问的任何其它介质。任何这样的计算机存储介质可以是装置的一部分、或连同装置一起使用。该装置还可以包括一个或多个通信连接184，该通信连接使该装置可以与其它装置(例如其它计算装置)通信。通信连接携带在通信介质中的信息。通信介质通常体现为计算才几可读指令、数据结构、程序模块或在调制的数据信号(如载波或其它运输机构)中的其它数据，并且包括任何信息传送介质。术语"调制的数据信号"是指这样的信号，该信号具有一个或多个其特征集合或以这样的方式加以变化^f吏得编码在信号中的信息。作为实例H旦并不作为限制)，通信介质包4舌有线介质如有线网会备或直4妾有线连4妄，以及无线介质如声音、RF、红外线以及其它无线介质。如上所述，如在本文中所使用的术语计算机可读介质包括存储介质和通信介质两者。本领域4支术人员应当明了，本文已讨i仑了用于改善连续式葡萄糖传感装置和利用该装置的系统的准确度的新的和有用的方法。然而，考虑到可以应用本发明的原理的许多可能的具体实施方式，应当认识到，本文参照附图所描述的具体实施方式仅是示例性的并且不应当视为限制本发明的范围。本领域的技术人员应当明了，在不偏离本发明的精神的情况下，可以在安排和细节上对示出的具体实施方式加以改进。因此，如本文描述的本发明i殳想所有这样的具体实施方式可以属于所附权利要求以及其等同替换的范围。附录A将以下出版物和美国临时申请的主题的全部内容以引用方式结合于本文。1，SantiagoJV:LessonsfromtheDiabetesControlandComplicationsTrial.IHct&e&s1993,^21549-1554.2.B.Feldman，RBazgS，Schwartz,andR.Weinstdn，"Acontinuousglucosesensorbasedonwiredenzymetechnologyandresultsfroma3daytrialinpatientswithtype1diabetes,"Dia&fesr"/mo/77(er，vol5，pp769-7，78,2003，3.DKlonoff》"Continuousglucosemonitoring:Roadraapfor21stcenturydiabetestherapy,"DZa6e^sCwe,vo128,1231-1239,2005.4.ECheyne，D.Cavan，andD.Kerr，言，Perfo迈anceofcontinuousglucosemonitoringsystemduringcontrolledhypoglycemiainhealthyvohinteers，"Z)t'a6efesTecAwo/77w/'，vol4，pp60'7-613,20025，MBoyne，DSilver,JKaplan，andC,Saudek,"Timingofchangesininterstitialandvenousbloodglucosemeasuredwithacontinuoussubcutaneousglucosesensor;"JD/fl办e魄vol52,pp2790-2794,2003,6.P7"Stout，JR,Raccliini，andM.EHilgers,"Anovelapproachtomitigatingthephysiologicallagbetweenbloodandinterstitialfluidglucosemeasurements,"7k力wc^3T^r,vol'(5,pp,635-644,20047*E，Kulcu,J,Tamada，G.Reach,RPo汰s，andM,Lesho，"Physiologicaldifferencesbetweeninterstitialglucoseandbloodglucosemeasuredinhumansubjects,"Z)纟"Z^toCe,vol26，pp2405-2409，2003,8.ASchoonenandKWientjes，"Amodelfortransportofglucoseinadiposetissuetoamicrodialysisprobe,"jD/"^/es1JfecA"o/7Te,'vol.5，pp,589-598，20039.G,MSteil,KRebrin，FHariri，SJinagonda，STadros,CDarwin,andMFSaa4"Interstitialfluidglucosedynamicsduringinsulininducedhypoglycaemia，"D/"h咖/ogte,vol48，pp.1833-1840,2005,10.KRebrin肌dGSteil,"Caninterstitialglucoseassessmentreplacebloodglucosemeasurements"Z)/aZ^es7fec/moiT7眠vol,2,pp.461472,2000"11.Kovatchev,LGonderFrederick,DCox，andWClarke,"Evaluatingtheaccuracyofcontinuousglucosemonitoringsensors:contiauousghiooseerlorgridanalysisillustratedbytherasensefteestylenavigatordata，"Z)faZ^ayvol.27,pp1922-1928,200412.WClarke,D.Cox，LGonderFrederick,WCarter,andSPoM,"Evaluatingclinicalaccuracyofsystemsforself-monitoringofbloodglucose,"Z)f"6&esC"/'e，vol.10,pp622、628,1987，13.HanefbldM:Postprandialhyperglycemia:noxiouseffectsonthevessdwall/"紐〃础'<￡7//0"/7^/^/'07/"^//^'"^^2002,Supplement129:45-5014.EspositoKjGiuglianoD，NappoF，MartellaK>fortheCampanianPostprandialHyperglycemiaStudyGroup.Ocw/加'o",2004,110:214-21915.QuagliiuL,PiconiL，AssaloneR,MartinelliL，MotzE，CerieloA:IntermittentI-IighGlucoseEnhancesApoptosisRelatedtoOxidativeStressinHumanUmbilicalVeinEndothelialCells:TheRoleo每roteinKinaseCandNAD(P)H-OxidaseActivationD油魄2003，a:2795，4.16.VanderDoesFE.DeNeelingJN，SnoekFJ，KostensePI,GrootenhuisPA》BouterLM，andRJHeine:Symptomsandwell-beinginrelationtoglycemiccontrolintypeIIdiabetesZ)Zate&yCa/^1996，u:204-21017.DeSonnaviUeJJ.SnoekFJCoUyLP，DevilleW.WijkelDHeineRJ:Weil-beingandsymptomsinrelationtoinsulintherapyintype2diabetes.ZZaZWejCare,1998，a:919-24，18.Co;cDJ，Gonder曙R:ederickLA，McCallA，etalTheeffectsofglucosefluctuationoncognitivefimctionandQOL:thefunctionalcostsofhypoglycaemiaan3hypeglycaemiaamongadultswithtype1ortype2diabetes/"em加'cma//ow,'"a/q/C/imWJV"crice，2002，Supplement129:20-2619.HirshIB,BrownleeM:Shouldminimalbloodglucosevariabilitybecomethegoldstandardofgiycemiccontrol/o/Z)/"Z^eawe/TteCo/n/>//cariow，2005,D:178-18120.ReichardP，PhilMMortalityandtreatmentsideeffectsduringlong-termintensifiedconventionalinsulintreatmentintheStockholmDiabetesInterventionstudyDiabetes43:313-317,199421.TheDiabetesControlandComplicationsTrialResearchGroupTheeffectofintensivetreatmentofdiabetesonthedevelopmentandprogressionoflong-termcomplicationsofinsulin-dependentdiabetesindlitusNEnglJMed329:978-986，199322.GrossTM，BodeBW，EinhornD，KayneDM，ReedIH,WhiteNH,MasttototaroJJ:PerformanceevaluationoftheMiMedcontimiousglucosemonitoringsystemduringpatienthomeuseDiabetesTechnolThel2000;2:49-5623.TaviDsR，ShoaibiA:ThepublichealthimpactoftheMiniMedContinuousGlucoseMonitoringSystem(CGMSO):anassessmentoftheliteratureDiabetesTechnolTher2004;6:518-52224.ChicoA，Vidal-RiosP，SubiiaM，NovialsA:Thecontinuousglucosemonitoringsystemisusehlfordetectingunrecognizedhypoglycemiasinpatientswithtype1andtype2diabetesb\itisnotbetterthanfiequentcapillaryglucosemeasurementsforimprovingmetaboliccontrolDiabetesCare26:1153-，1157,200325,MasbototaroJ:TheMiniMedContinuousGlucoseMonitoringSystem(CGMS),JPediattEndocrinolMetab12:751-758，199926,FddmanB，BrazgR，SchwartzS,WeinsteinR:Aoontinuousglucosesensorbasedonwiredenzymetechnology"resultsfroma3-daytrialinpatientswithtype1diabetesDiabetesTechnolThel2003;5(5):769-7927,PickupJ:SensitiveglucosesensingindiabetesLancet355:426-427，200028,GerzitsenM，JansenJA，LuttermanJA:PerformanceofsubcutaneouslyimplantedglucosesensorsforcontinuousmonitoringNethJMed54:167-179,199929,Boyne，Silver,Kaplan,&Saudek:TimingofChangesinInterstitialandVenousBloodGlucoseMeasuredWithaContinuousSubcutaneousGlucoseSensorDiabetes52:2790-2794，200330,PhilipJStout，JoelRRacchini，MichaeEHilgers:ANovelApp^oachtoMitigatingthePhysiologicalLagBetweenBloodandInterstitialFluidGlucoseMeasurementsDiabetesTechnolTher20046:5,635-64431,Schoonen&Wien勺'es:AModelforTransportofGlucoseinAdiposeTissuetoaMicrodialysisProbeDiabetesTechnolTher20035:4,589-59%32,ASPopelAnalysisofcapillary-tissuediffusioninmulticapillarysystemsMathematicalBiosciences,39:187-211，1978,33,G.MSteil，K.Rebrin，FHariri,SJinagonda，STadros，CDarwinandMFSaad,，Interstitialfluidglucosedynamicsduringinsulin-inducedhypoglycaemiaDiabetologia48:7,publishedonline,DOI:10.1007/s00125-,005-1852-x34，ZierlerK:WholebodyglucosemetabolismAmJPhysiol276:E409-E426，199935,PalfreymanRW，ClarkAE,DeatonRM，HolmanGD,KozkaIS:KineticresolutionoftheseparateGLUT1andGLUT4glucosetransportactivitiesin313-LIcellsBiochemJ.284:275-,282，1992.36,StengelRF:OptimalControlandEstimationNewYork:DoverPublications,199437,KnobbeEJ,BuckinghamB:TheextendedKalmanfilterfocontinuousglucosemonitoringDiabetesTechnolTher2005;7:15-27"38,PalermCesxC.，WillisJohnP,，DesemoneJames，BequetteBWayne:HypoglycemiaPredictionandDetectionUsingOptimalEstimationDiabetesTechnolTher200亏；7:3、14，39,"GM.Stcil，KRebrin，FHariri，SJinagonda,STadros，CDarwinandMF*Saad,Interstitialfluidglueosedynamicsduinginsulin-inducedhypoglycaemiaDiabetologia48:7，publishedonline.DOI:101007/s00125-005-1852-x"40.InternationalPatentApplicationSerialNo.PCTiUS20051013792，entitled"Method:SystemandComputer.ProgtamProductforEvaluatingtheAccuracyofBloodGlucoseMonitoringSensorsiDevices,"filedApril21,2005，,41,USprovisionalapplicationserialnumber60/815,191toKovatcheva/,filedJune20:2006;附录B-卡尔曼递归W寸录B:卡尔曼递归(KalmanRecursion)如果我们假设我们的数据遵循形式为Al.l的系统，其中x是隐系统状态向量(矢量)、y是观测向量、Ws是具有协方差矩阵Q的系统噪声、w。是具有协方差矩阵R的观测噪声，以及F、G、H矩阵定义系统转变，那么该系统被称为状态空间模型。<formula>formulaseeoriginaldocumentpage41</formula>卡尔曼滤波器是两步过程。首先，进行一步向前预测。接着，基于新的测量结果对预测进行校正。让我们引入记号XW，其是指"在测量结果达到和包括j的情况下在i处的X的估计值"。*1/没我们具有X的目前估计值以及该估计值的误差协方差P,那么我们可以将第一步骤写成A1.2。为了对估计值进行一交正，然后我们在A1.3中计算了卡尔曼增益，并且在A1.4中利用新的测量结果y对估计值进行更新。<formula>formulaseeoriginaldocumentpage41</formula>附录c-递归方法<formula>formulaseeoriginaldocumentpage42</formula>可以如在A2.9中写对角算子，其将A2.8简化成A2.10。此外，通过再次并列(concatenate)第一q+l列向量，我们可以形成在A2.ll中的所期望的矩阵A2.1。<formula>formulaseeoriginaldocumentpage42</formula>最后，通过S夸A2.ll入A2.3,我们到达A2.12，一种可以通过诸如QR分解的方法来求解的线性系统。<formula>formulaseeoriginaldocumentpage42</formula>附录D-实时实施的禾呈序代码版本l:实时应用1GetNewBGPreds<-function(nav，times，par){#需要几分钟的时间使参数工作正常#消除漏测值2.times[!is，.na(Nav)]->times;Nav[!is.na(Nav)l画>Nav;#第二，应用校准参数3.Navipar[[31]-〉Nav;#第三，产生阵列，以保持IG，和点(dot)IG的估计值4.rep(0.,times=length(Nav))-〉Idot;Nav->IG;#第四，在大量数据已被收集之前的时间内产生点IG的估计值5.(IG[2]画IG[11)/(times[Zl-timesill)画〉Idot121;6.#线性内插法(线性插值)7dataframe(y=IG[l:3]，x=times[l:31)國〉td;8lm(y-x,data=td)-〉tx;txScoefficients[2]->tx;attr(tx，"names")<-NULL;tx->Idot[3];9dataframe(y=IG[l:4]，x=times[1:4])->td;101m(y-x，data=td)->tx;tx$coefficients[2]-〉tx;attr(tx，"namesw)<-NULL;tx-〉Idot[4];#第五，在CGS凄t据上4吏用样条来产生估计值11for(iin5:length(Nav》(12times[(i-4):(i)]->tt;Nav[(i画4):(i)1->td;13smoothspline(x=:tt，y=td，df=l)-〉temp;14predict(temp,times[i]，deriv=l)$y-2Idot[i];15predict(temp,times[I])Sy->IG[il;16}#第六，应用主方禾呈17IG+Idotkpar[[2]]+par[[1]]->NewBG;#第七，返回结果18NewBG)附录E-追溯实施的程序代码版本2:追溯4交正除了实时实施外，可以通过延长将每一样条内插入12行(从"(i-4):(i)""(i匿4):(i+4)，，)的时间来引入人工延迟。在完全的追溯分冲斤中，可以z使用下面乡会出的更简单的形式。19GetOldBGPreds<-function(Nav，times,par){弁消除漏测^直20times[！is.na(Nav)]->times;21Nav[！is.na(Nav)1->Nav;#应用才交准参凄史22par[[3]1*Nav-〉Nav;#在所有数据中产生平滑样条回归23sreg(times,Nav)->Navfit;#4吏用它来预测IG24predict(Navfit，tirnes)->IG;弁^f吏用它来预测点IG25predict(Navfit,times，deriv=l)-〉Idot;#应用主方禾呈26(Idot)*par[[2]]+par[[I]1+IG->BGpreds;弁返回予贞测27BGpreds)最后，如果期望方程3的全部米-门氏形式，则在17和27^f亍中可以简单地用"par[[l]]*IG/(126.+IG)，，代替"par[[l]]"。R包含许多其他平滑和内插程序，并且使用sreg、平滑样条、以及/w可以用它们中的多种来代替，虽然输入/输出格式经常不同。权利要求1.一种用于改善连续式葡萄糖传感器(CGS)的准确度的方法，包括在第一时间校准所述CGS；以及在第二时间校准所述CGS，所述第二时间是基于动态监测的CGS值和CGS变化的速率中的至少一个来确定的。2.根据权利要求1所述的方法，其中，所述动态监测的CGS值在进4亍所述第二次4交准以前相差至少15mg/dl。3.根据权利要求1所述的方法，其中，所述动态监测的CGS值在进行所述第二次校准以前相差至少20mg/dl。4.根据权利要求1所述的方法，其中，所述动态监测的CGS值在进4亍所述第二次4交准以前相差至少30mg/dl。5.根据权利要求1所述的方法，进一步包括随时间动态地监测所述CGS值和所述CGS变化的速率。6.才艮据权利要求5所述的方法，进一步包括计算所述监测的CGS值的一阶时间导数；确定所述CGS值的所述一阶时间导数的绝对值是否基本上小于1;以及当所述CGS值的所述一阶时间导凄t的所述绝对值基本上小于1时，保持在所述第一时间的对所述CGS进行的校准。7.根据权利要求6所述的方法，进一步包括当所述CGS值的所述一阶时间导凄t的所述绝对值基本上小于1时，确定在所述第二时间的所述CGS与在所述第一时当在所述第二时间的所述CGS与在所述第一时间4交准的所述CGS之间的所述绝对变化大于预定值时，^f呆持在所述第一时间的对所述CGS进行的校准。8.根据权利要求7所述的方法，其中，所述预定值是30wg/c//。9.根据权利要求7所述的方法，其中，所述预定值是40wg/W。10.根据权利要求7所述的方法，进一步包括当确定了在所述第二时间的所述CGS与在第一时间校准的所述CGS之间的所述绝对变化大于预定值时，对所述CGS进行再校准。11.根据权利要求1所述的方法，进一步包括通过4交正所述血4唐水平和在与所述血^唐相互作用的间质液中的葡萄糖水平之间的生理时滞来改善所述CGS输出的所述准确度。12.根据权利要求11所述的方法，其中，改善所述准确度的所述步骤进一步包括导出数学方程，所述数学方程描述CGS输出对所述血糖水平的时间相关性；以及将所述导出的时间相关性血糖水平函数应用于一组CGS原始输出以便在随后预测所述血糖水平或校正所述CGS输出。13.根据权利要求12所述的方法，其中，所述数学方程是基于一种模型而导出的，所述模型包括对所述血糖与在所述间质液中的葡萄糖之间的扩散相互作用的描述。14.根据权利要求12所述的方法，其中，所述数学方程是基于一种模型而导出的，所述模型包括对在所述间质液中的所述葡萄糖的假设的描述。15.根据权利要求1所述的方法，其中，所述第二时间是同时基于所述动态监测的CGS值和所述CGS变化的速率来确定的。16.—种用于改善连续式葡萄糖传感器的准确度的方法，包括利用第一血糖值和第二血糖值校准所述连续式葡萄糖传感器；当所述连续式葡萄糖传感器测量相应的第一间质葡萄糖<直时，测量所述第一io4唐^直；以及当所述连续式葡萄糖传感器测量相应的第二间质葡萄糖值时，测量所述第二血糖值，其中所述第二间质葡萄糖值与所述第一间质葡萄糖值相差预定量。17.才艮据4又利要求16所述的方法，其中，所述预定量为至少15mg/dl。18.才艮才居4又利要求16所述的方法，其中，所述预定量为至少20mg/dl。19.根据权利要求17所述的方法，其中，所述第一和第二葡萄糖数据是在不同时间获得的。20.根据权利要求17所述的方法，其中，所述第二葡萄糖数据是通过决定形成环获得的，所述决定形成环包括基于预定准则来确定CGS变4匕的速率。21.根据权利要求20所述的方法，其中，所述预定准则》见定如果所述CGS变化的速率等于或大于预定值，则通过所述第一葡萄糖数据保持所述CGS校准。22.根据权利要求21所述的方法，其中，所述预定值是l。23.根据权利要求22所述的方法，其中，所述预定准则进一步夫见定如果所述CGS变化的速率小于1，则基于另一个预定准则进4亍另一个决定。24.根据权利要求23所述的方法，其中，所述另一个预定准则规定如果所述CGS与通过所述第一血糖数据4交准的CGS之间的绝对差异等于或小于预定差值，则通过所述第一葡萄糖数据保持所述CGS校准。25.根据权利要求24所述的方法，其中，所述另一个预定准则进一步^L定如果所述CGS与通过所述第一血净唐lt据冲交准的CGS之间的绝对差异大于所述预定差值，则通过所述第二血糖凄丈据进4亍另一个冲交准。26.才艮据4又利要求17所述的方法，进一步包括通过4交正所述血4唐水平和在与所述血斗唐相互作用的间质液中的葡萄糖水平之间的生理时滞来改善所述CGS输出的所述准确度。27.根据权利要求26所述的方法，其中，改善所述准确度的步骤进一步包括导出数学方程，所述方程描述CGS输出对所述血糖水平的时间相关性；导出作为所述CGS输出的函数的时间相关性血糖水平；以及将所述导出的时间相关性血糖水平函数应用于一组CGS原始输出以便对所述CGS输出作固有时滞的校正。28.根据权利要求27所述的方法，其中，所述数学方程是基于一种才莫型而导出的，所述才莫型包括对所述血糖与在所述间质液中的葡萄糖之间的扩散相互作用的描述。29.根据权利要求27所述的方法，其中，所述数学方程是基于一种才莫型而导出的，所述才莫型包括对在所述间质液中的所述葡萄糖的假设的描述。30.根据权利要求27所述的方法，其中，作为所述CGS输出的函数的所述时间相关性血糖水平是描述所述CGS输出对所述血糖水平的时间相关性的数学方程的反函数。31.—种用于改善连续式葡萄糖传感器的准确度的系统，包括可访问所述连续式葡萄糖传感器的输出的准确度改善模块，进一步包括决定形成机构，所述决定形成机构能够确定在某一时间的校准，其中所述时间是基于所述连续式葡萄糖传感器的动态监测^r出中的至少一个。32.根据权利要求31所述的系统，其中，所述决定形成机构能够确定在所述时间的所述校准，其中所述时间是基于所述连续式葡萄糖传感器的所述输出和所述连续式葡萄糖传感器的所述输出的变化速率。33.根据权利要求31所述的系统，其中，所述准确度改善模块进一步包4舌能够以预定时间序列启动校准的初始请求才莫块。34.根据权利要求33所述的系统，其中，所述决定形成机构进一步包括能够对来自所述连续式葡萄糖传感器的输出进行动态监测的监测模块。35.根据权利要求34所述的系统，其中，所述决定形成机构进一步包括确定指令，其规定在进行所述第二次校准以前来自所述连续式葡萄糖传感器的所述动态监测的输出相差至少15mg/dl。36.根据权利要求34所述的系统，其中，所述决定形成机构进一步包括确定指令，其规定在进行所述第二次才交准以前来自所述连续式葡萄糖传感器的所述动态监测的输出相差至少20mg/dl。37.才艮据4又利要求34所述的系统，其中，所述决定形成4几构进一步包括确定指令，其^见定在进^f于所述第二次才交准以前绝乂于变4b速率基本上等于或大于1。38.根据权利要求34所述的系统，进一步包括时滞4交正，莫块，其通过4交正所述血糖水平和在与所述血外唐相互作用的间质液中的葡萄^f唐水平之间的生理时滞能够改善所述CGS输出的准确度。39.—种试剂盒，包括连续式葡萄糖传感器；与所述连续式葡萄糖传感器数据通信的校准单元，用于校准所述连续式葡萄糖传感器；以及指令集，其包括用于指示监测来自所述连续式葡萄糖传感器的受监测的输出的指令、以及当第一和第二连续式葡萄糖传感器输出值相差至少15mg/dl时用于进4亍第一次和第二次校准的指令。40.根据权利要求39所述的试剂盒，其中，所述校准单元可连接于所述连续式葡萄糖传感器。41.根据权利要求39所述的试剂盒，其中，所述校准单元是所述连续式葡萄糖传感器的一体构件。42.根据权利要求39所述的试剂盒，其中，所述指令进一步规定，在进行所述第二个校准以前，所述连续式葡萄糖传感器的所述动态监测的llT出相差至少20mg/dl。43.根据权利要求39所述的试剂盒，其中，所述指令进一步规定，在进行所述第二个校准以前，所述连续式葡萄糖传感器的输出的绝对变化速率小于1。44.才艮据权利要求39所述的试剂盒，其中，所述准确度改善才几构进一步包4舌时滞校正模块，其可访问所述连续式葡萄糖传感器的输出并且能够校正生理时滞。45.—种用于改善连续式葡萄糖传感器的准确度的方法，包括利用第一血糖水平和与第一葡萄糖数据相差预定量的第二血糖水平对所述CGS进行校准。46.根据权利要求45所述的方法，其中，所述预定量为至少15mg/dl。47.根据权利要求46所述的方法，其中，所述预定量为至少20mg/dl。48.根据权利要求47所述的方法，其中，第二葡萄糖数据是通过决定形成环获得的，所述决定形成环包括基于预定准则来确定CGS变4匕的速率。49.根据权利要求48所述的方法，其中，所述预定准则规定如果所述CGS变化的速率等于或大于1，则通过所述第一葡萄糖数据保持所述CGS校准。50.—种计算才几可读介质，所述介质具有用于实施在4又利要求1中提出的所述方法的计算机可执行指令。51.—种计算才几可读介质，所述介质具有用于实施在4又利要求16中提出的所述方法的计算机可执行指令。52.—种计算冲几可读介质，所述介质具有用于实施在片又利要求45中提出的所述方法的计算机可执行指令。53.—种试剂盒，包括用于校准连续葡萄糖的校准单元，所述4交准单元包括血糖传感器；以及指令集，所述指令集包括用于监测来自连续式葡萄糖传感器的输出的指令、以及当第一和第二连续式葡萄糖传感器输出^:相差至少15mg/dl时用于进4亍第一次和第二次才交准的指令，其中所述第一次和第二次4交准分别包括通过所述血冲唐传感器对血糖进行第一次和第二次确定。全文摘要一种利用连续式葡萄糖传感器的动态输出的能够改善CGS装置的准确度的方法、装置、以及试剂盒。文档编号A61B5/05GK101365374SQ200680040174公开日2009年2月11日申请日期2006年8月29日优先权日2005年8月31日发明者克里斯托弗·瑞安·金,鲍里斯·P·科瓦特舍维申请人:弗吉尼亚大学专利基金委员会;雅培糖尿病护理公司