本发明设计一种新兴的血糖浓度光学检测数据回归算法,特别是采用逐步回归分析,在有限数据中得到了最优回归。
背景技术:
糖尿病是一组以高血糖为特征的代谢性疾病,高血糖则是由于胰岛素分泌缺陷或其他生物作用损伤引起的。糖尿病时长期存在的高血糖,导致各种组织,特别是眼、肾、心脏、血管、神经的慢性损害、功能障碍。目前,随着我国生活质量的提高,这一疾病在大城市中越来越常见,因此这一疾病也得到了人们的极大关注。目前,越来越多的血糖仪出现在市场,为病人提供实时、便利的血糖浓度监测,为医护人员的治疗方案提供思路。通过新型的传感器、光谱分析、电渗透技术、光电检测技术以及信号处理等手段,血糖监测可以以无创的方式进行。
无创光电检测利用的是测量信号与人体组织之间的相互作用,通过反馈回啊立的幅度信息,我们可以分析患者体内的血糖含量。这一方法对内部结构非常敏感,因此我们需要通过复杂的算法对大量的数据进行分析,建立最有模型。
在实际光学检测数据与血糖浓度分析的过程中,传统的多元线性回归建立的线性关系可能存在相关性低、置信区间太大、包含了无统计意义的变量,最后导致建立的关系很难被应用在无创检测的设备中,降低了血糖浓度的检测结果。
技术实现要素:
为了克服以上不足,本发明利用了多元逐步回归算法,以在有限的数据量中建立最优的回归分析。
具体实施方式
以下通过具体的实例分析本算法的实施方式,本领域技术人员可以根据本说明书所揭示的内容轻易的了解到这一算法的优势。本算法亦可通过其他不同的计算机语言实现,本说明书中的各项细节亦可基于不同观点与应用,在不背离基本统计原理的情况下可以进行各种修饰与变更。
本实施中所淑的y变量包括血糖估计值,优化血糖数量,用来作为自变量的为无创光学测量的数据,包括特征光峰值,并配合心率,血流速度等参数,建立回归方程。
首先利用检测到的数据,将两个峰值设为x1和x2,心率设为x3,血流速度设为x4……而优活血糖数量设为y1,血糖浓度为y2,通过临床实验我们得到n组数据,接下来以y1(优活血糖数量)为例,进行回归分析我们将回归方程写为:
首先建立系数矩阵
y=[y1y2y3...yn]t
b=[blb2b3...bn]t
首先通过最小二乘法得出系数矩阵
b=(xtx)-1xty
接着利用计算好的系数矩阵参数估计残差平方和,再多个变量的情况下,应变量的总离均差平方和ss总可分解为回归平方和ss回归和残差平方和ss残差两个部分
上式中yi是观察值而
其中bj为回归系数,ljy为第j个自变量与应变量y的离均差积和。残差平方和为ss残差=∑(y-yj)2=ss总-ss回归
根据残差的公式我们要对引入的自变量进行判断,看其与应变量的关系如何
比如以引入的血流速x4为例,我们去掉x4回到第一步重新用最小二乘法进行判断,得到新的ss回归-4
uj=ss回归-ss回归-4
接着利用f统计进行判断
我们做出假设认为x4与y有线性关系,利用f统计的书指标,我们设定α=0.05,对比我们所计算的f与f0.05(1,n-m-1)值,如果大于则认为该变量对y有线性关系,否则认为没有,接下来去掉该变量重新回到第一步进行回归分析。
本算法具有如下有点:
(1)适合于自变量有很多时的血糖浓度回归分析,可以得到“最优”回归方程;
(2)本算法可以简单的通过各种计算机语言实现,适合于大数据分析;
(3)本算法可以与专用集成电路结合,便于硬件实现。
附图说明
图1为多元逐步回归算法流程图。