使用腕戴式可穿戴设备的加速度进行在线心率预测

摘要。

在本文中，我们研究了使用手腕佩戴的可穿戴设备根据加速度来预测心率。 尽管现有的光电体积描记法 (PPG) 心率传感器可提供可靠的测量，但它们使用的能量比加速度计多得多，并且对可穿戴设备的电池寿命有重大影响。 通过使用节能加速计来预测心率，可以显着节省能源。 此外，我们有兴趣了解心率干预后患者的恢复情况，我们预计心率会随着时间的推移而变化。 因此，我们提出了一种在线方法来解决随着时间的推移而产生的概念漂移。 我们根据三名患者几个月内大约 4 周的自由生活数据来评估这些方法。 我们表明，我们的方法可以在不频繁使用 PPG 心率传感器（例如 20.25% 的样本）的情况下实现良好的预测性能（例如，2.89 平均绝对误差）。

1. 介绍

医疗保健技术在消费者和临床应用中的重要性持续快速增长（等人，2015）。 消费者越来越希望能够通过使用技术来监控和了解他们的健康和福祉。 智能手表和健康手表都是常用的设备。 临床医生越来越多地使用这些设备来了解患者在较长时间内的行为，而不是在临床环境中进行单一测量。

健康手表通常以活动检测、定位、呼吸或心率的形式提供健康测量。 由于加速度计的低成本和低功耗，使用腕戴式设备进行活动检测（例如通过步数计数）的方式被广泛使用。

虽然加速度计的能耗较低，但常见的心率测量传感器消耗的电量要多得多。 最典型的方法之一是光电体积描记法 (PPG)，它通过反射光从微血管组织（Challoner 和 Ramsay，1974） 中观察到的血量脉冲来测量心率。 消费级 PPG 传感器的功耗通常是可穿戴设备中使用的加速度计（Elsts 等人，2017 年） 的 5000 倍，这阻碍了可穿戴技术所需的长电池寿命。 由于加速度计很普遍并且存在于任何可能包含心率传感器的设备中，因此我们有兴趣考虑加速度作为预测心率的手段的可行性。 我们注意到这其中固有的困难；也就是说，有许多已知因素会影响心率，但不会被加速度捕获。 此外，我们预计相同的加速度模式与心率的关系将随着时间的推移而改变，例如，由于衰老或健康水平的变化。 因此，我们认为，不适应这些变化的方法将无法根据加速度很好地预测心率。

2. 设想

EurValve项目正在构建主动脉瓣和二尖瓣疾病的决策支持系统（DSS）。 其中之一涉及广泛收集环境数据，用于患者在家中的行为分析（Pope 等人，2017；McConville 等人，2018b），在心脏瓣膜介入治疗之前和之后两次。 其中 42 名患者已被招募，目前正在其家中分阶段部署定制智能家居盒 (SHiB) （Pope 等人，2017 年）。 SHiB 包含四个网关，用于从手腕佩戴的可穿戴设备接收加速度计数据。 这些网关随后用于室内定位（McConville 等人，2018a），但也可用于本地数据处理和机器学习。 这款可穿戴设备仅包含一个加速度计，没有心率传感器，但电池寿命约为三周。 招募的患者还佩戴了 Phillips Health 手表（Hendrikx 等人，2017），该手表能够使用 PPG 测量心率，但电池寿命长达 4 天。 在参与研究期间，患者将佩戴这两种设备，因此我们将收集每位患者最多 6 周的免费生活数据。 这提供了一个独特的第一次机会来研究利用加速度来预测现实世界自由生活条件下不同恢复阶段的相当长一段时间内的心率的潜力。

2.1. 愿望

我们将首先产生一种基线方法，在该方法中，我们根据从加速度中提取的特征训练回归模型，与心率对齐，然后尝试仅根据加速度来预测未来的心率。 由于上述挑战，我们认为这种方法无法捕捉加速度和心率之间的关系。 为了改进这一点，我们将提出一种可以在 SHiB 网关上运行的主动学习方法，该网关与可穿戴设备进行双向通信。 这种方法将以在线方式根据流式加速计数据预测心率，并能够在需要时通过 PPG 请求测量真实心率。

3. 实验

实验评估将涉及 3 名患者，他们在心脏瓣膜介入治疗后不久佩戴该设备约 2 周，然后在大约 2 个月后再佩戴 2 周。 我们将单独的事件称为阶段。 加速度特征是在一秒窗口内提取的，如下所示：三个轴中每个轴的最小值、最大值、标准差、中值、平均值、第 25 个和第 75 个百分位数、四分位数、偏度、峰度和过零数。 还提取了每个轴的谱能量和谱熵。 由于飞利浦健康手表以一分钟的间隔记录心率，因此加速度计功能在一分钟的间隔内进行平均，以与心率进行时间对齐。

3.1. 离线基线

基线方法将假设在两个单独的大约时间内收集所有加速度计数据和心率数据。 14 天的监测阶段。 目标是，当给定一分钟的加速度数据时，预测该时间点相应的心率。 在第一个实验中，一个阶段中一名患者的 60% 数据将用于构建模型，其余 40% 用于测试。 在这种情况下，我们将数据视为独立同分布。

图1显示了该实验的结果。 图 1(a) 和 1(c) 表明，平均绝对误差 (MAE) 介于 4 和 6 之间，均方误差 (MSE) 介于 4 和 6 之间，可实现合理的性能。 23 和 51 取决于评估的患者和阶段。 图 1(b) 显示了“虚拟”预测器的结果，该预测器仅输出每分钟的平均心率。 这可能比基于加速的模型差多达 2 倍。

如果我们希望建立一个可以预测更长时间的未来的模型，当前的评估还不够，因为它是在连续的大约时间内进行的。 14天期限。 因此，我们提出了另一个实验，其中模型建立在一个相（大约）上。 连续 14 天的数据）并在另一个阶段进行评估（即 2 个月后的另外 14 个连续天）。 在这种情况下，如图 2 所示，性能显着下降至 6 至 8 MAE 和 65 至 90 MSE。 因此，通过这种方法，任何根据加速度预测心率的模型似乎都会随着时间的推移而变得不太准确。

在图 3 中，我们可以看到单个患者的实际心率与预测心率的关系。 由此可见，预测并不充分，而且很少偏离均值太远，这在 MSE (94.8) 中得到了反映。

（一） 随机森林回归器 (MAE))

(b) 虚拟均值预测器 (MAE)

(c) 随机森林回归器 (MSE)

(d) 虚拟均值预测器 (MSE)

图 1。 在同一连续约 14 天的时间段内进行测试时的 MAE (a) MAE 介于 4 到 6 之间，比 (b) 中的虚拟预测器表现更好。 MSE (c) 的趋势相同，但值得注意的是，性能显然很差。

（一） MAE

(b) 均方误差

图 2。 MAE，在不同的连续大约 14 天的时间段上进行训练和测试时 (a) MAE 介于 6 和 8 之间，如图 1 所示，性能有所下降。 MSE方面的表现也进一步下降。

图 3。 使用患者 0130 的离线基线方法的真实心率与预测心率。 图表的左侧用于训练，而预测的心率值显示在右侧。

通过研究这些阶段的心率分布，就可以清楚为什么表现会如此下降。 在图4中我们可以看到心率的分布随时间明显变化。 此外，我们预计由于非基于加速度的特征，加速度和心率之间的概念漂移会在短时间内频繁变化。 考虑到这两个观察结果，我们可以将其建模为由于后验分布的变化而导致的真实概念漂移（Gama等人，2014）的形式$P(y|X)$ 通过开发一种更具体地捕获心率 ($y$) 和加速度 ($X$) 之间数据的时间变化的方法。

（一） 0130号病人

(b) 1231号病人

(c) 549号病人

图 4。 三名患者在两个不同的 14 天阶段（间隔大约 2 个月）的心率分布。

3.2. 主动学习

如果我们假设特定的加速模式（例如步行）表明心率会增加，我们可以将其重新表述为主动学习问题，在该问题中我们可以持续获得加速度，并且可以查询真实心率需要时标记当前加速度样本。 如果我们可以通过使用（能源廉价）加速度预测心率来减少（能源昂贵）心率测量的数量，我们就可以延长可穿戴设备的电池寿命，这对于可穿戴技术、数字健康、智能家居和物联网 (IoT) 应用。

通常，在主动学习问题中，可以查询样本的标签（即心率），特别是对于标签特别有用的数据样本。 在最常见的主动学习方法中，假设存在一个数据池，可以选择要标记的数据。 然而，这在我们的设置中是不可行的，因为加速度数据不断到达，我们希望一致地产生心率估计，并且不可能回顾性地测量心率。 因此，我们必须在每个时间点以在线方式决定是否要承担与获取真实心率相关的能量成本，或者是否可以可靠地使用加速度来预测该值。 $P(y|X)$ 可能会发生变化，因此必须考虑到自然的时间方面，这一事实削弱了这一决定。

3.3. 算法

我们将这些不同的属性建模到算法 1 中描述的过程中。 该算法需要回归器$L$的集合，它提供了预测不确定性的度量。 这种不确定性通过参数 $O$ 进行阈值化，该参数确定了可容忍的不确定性量（就整体方差而言）。 该算法维护大小为 $N$ 的历史记录以捕获新近度。 存储先前的 $N$ 预测方差，并使用它们的 $O$ 倍标准差设置阈值，以确定当前预测是否确实是异常值。 此参数 $N$ 还控制训练学习者集合 $L$ 的初始数据量，以及为在集合中重新训练学习者而保留的最新标记数据的大小。 当存在不确定性时 根据预测，测量真实心率，并使用最新标记的 $N$ 样本更新集合中的每个模型。 我们还借此机会检查集成的实际准确性，如果任何学习者的错误预测大于 $T$，那么他们将在最新标记的 $N$ 样本上进行重新训练。 整体可能对预测有确定性，但底层的 $P(y|X)$ 在整体不知情的情况下已经发生了变化。 因此，我们包含一个 $TTL$ 参数，它代表每个学习者的生存时间。 每次预测后，它们的生存时间都会增加，一旦学习器达到 $TTL$，就会使用最新的 $N$ 样本进行重新训练。

Input: $L$ learners, $N$ size of history, $O$ uncertainty threshold, $T$ error threshold, $TTL$ model time to live

Create $L$ learners in ensemble;

Train each learner on first $N$ instances;

while True do

$Acc$ = CollectAccDataFor60Seconds();

$HeartRate$, $Variance$ = PredictHeartRate(Acc);

if IsOutlier($O$, $Variance$) then

$TrueHeartRate$ = QueryTrueHeartRate();

TeachLearners($Acc$, $TrueHeartRate$);

for $l$ $\in$ $L$ do

if abs($HeartRate$-$TrueHeartRate$) ¿ $T$ then

retrain $l$ on most recent $N$ ;

end if

end for

end if

for $l$ $\in$ $L$ do

if older than $TTL$ then

retrain $l$ on most recent $N$ ;

end if

end for

Output $HeartRate$;

end while

Algorithm 1 PPAW (Predicting Pulse from Acceleration on Wearables)

3.4. 结果

我们将所提出的算法应用于三名患者的自由生活数据。 我们使用 MAE 来测量算法的性能，MAE 可以解释为每分钟平均心跳误差、MSE 以及测量真实心率以标记加速度数据样本的样本百分比。 虽然我们的算法有许多参数，但由于空间限制，我们将把调查限制在最有影响力的 $O$ 参数上，该参数处理当前预测的不确定性。 这里，$N$固定为5，$TTL$固定为10。

3.5. 变化的不确定性 ($O$)

在图 5 中，我们可以看到 MAE w.r.t。 不确定性参数$O$的三个明智选择。 虽然 $O$ 值较高，性能较低，但在所有情况下，MAE 都在 3 到 4 之间。 图6显示了相应的查询标签百分比。 对于较高的 $O$ 值，对心率传感器的查询数量会减少。 当$O$参数为1时，大约50%的时间请求真实心率，而当$O$参数值为1时，请求真实心率大约在15%到20之间所有患者的时间百分比。 这是心率传感器使用量的显着减少，而性能仅略有下降。 在所有情况下，与基线方法相比，我们的在线方法的性能都有显着提高。

（一） 0130号病人

(b) 1231号病人

(c) 549号病人

图 5。 $O$ 参数（不确定性）对 MAE 的影响。

（一） 0130号病人

(b) 1231号病人

(c) 549号病人

图 6。 $O$ 参数（不确定性）对查询标签百分比的影响。 在图 5 中，查询较高 $O$ 值的百分比显着下降，但相应 $O$ 的性能仅略有下降。

我们还在图 7 中绘制了所有阶段的真实心率（蓝色）和预测心率（绿色）。 我们可以看到我们的方法相对于图 3 中的基线有了显着的改进。 其重要性直观地反映在心率峰值和谷值之后的预测中，这在我们的在线方法中显然更好。 其他患者也发现了类似的结果，所有病例的 MAE 和 MSE 均显着改善。

图 7。 使用我们的方法($0$=3、$N$=5、$TTL$=10）对 549 号患者进行真实心率与预测心率的概述。 我们可以看到，预测心率与真实心率的趋势密切相关（MAE=2.89，MSE=21.07），而仅在 20.25% 的时间内请求测量心率。

4. 结论

在本文中，我们研究了根据加速度值预测心率的可行性。 一个关键动机是了解心率与加速度的关系，以及加速度是否可以成功预测长时间内的心率。 在这项工作中，我们特别感兴趣的是从心脏瓣膜介入治疗中恢复的患者心率随时间的变化。 另一个动机是通过利用更节能的加速计来预测心率，减少使用昂贵的心率传感器，从而从可穿戴设备中获得更大的节能效果。 我们首先提出了一种基线方法，该方法学习离线模型以根据加速度预测心率。 然而，我们发现，当尝试预测几个月后同一患者的心率时，这种方法表现不佳。 因此，我们提出了一种解决此问题的方法，并使用在线主动学习方法来适应基于非加速度的心率变化。 我们对三名患者评估了这种方法，每名患者都有大约 4 周的自由生活数据，这些数据在相隔几个月的两个阶段中收集，结果表明，我们可以实现每分钟 2.5 到 5 次心跳之间的 MAE，同时查询真实心率相对较少。 在一个示例中，MAE 为 2.89，而心率传感器仅有 20.25% 的时间被查询。