Урок 20 из 20 · Уровень 4: Аналитик

Как читать исследование ВСР или проверку устройства

Статья о ВСР ценна ровно настолько, насколько полно она описывает цепочку из урока 16 и настройки из уроков 17–19. Два опубликованных контрольных списка перечисляют, что статья обязана сообщить. Проверка устройства добавляет четыре статистических показателя, и каждый отвечает на свой вопрос. Урок показывает, как читать и то и другое, на примере исследования носимых устройств 2025 года, а в конце объясняет, что может и чего не может значить слово «валидировано» на коробке.

Содержание 6 разделов

Главное

  • Статью, которую нельзя воспроизвести по описанным в ней методам, нельзя и оценить. Что в ней должно быть, перечисляют 13 пунктов GRAPH и 30 пунктов Catai.
  • Корреляция не измеряет согласие. Bland-Altman даёт смещение и диапазон, в который попадают 95 % отдельных разностей.
  • ICC и CCC снижаются и от смещения, и от разброса; r Пирсона на смещение не реагирует. MAPE — это обычная средняя процентная ошибка.
  • В исследовании 2025 года, где эталоном была ЭКГ нагрудного датчика, ночной RMSSD лучших носимых устройств в 95 % ночей отличался от эталона примерно на 12 мс или меньше. Для многих людей этот диапазон шире изменения, на которое они реагируют.
  • «Валидировано» значит: относительно этого эталона, на этих людях, для этого показателя и окна, с этой прошивкой, по этому статистическому показателю. Спрашивайте обо всех шести пунктах.
После урока вы сможете
  • Сверить статью о ВСР с контрольными списками GRAPH и Catai и назвать, что в ней не указано.
  • Прочитать график Bland-Altman: смещение, границы согласия и что они значат для одного измерения.
  • Различать ICC, коэффициент корреляции конкордантности и MAPE и сказать, что каждый из них может и не может показать.
  • Объяснить, почему результат для группы не переносится на число одного человека, и как оценить согласие носимого устройства по сравнению с вашими собственными колебаниями от дня ко дню.
  • Расшифровать слово «валидировано» на странице продукта.

Статья о вариабельности сердечного ритма (ВСР) ценна ровно настолько, насколько полно она описывает цепочку из урока 16 и настройки из уроков 17–19. Два опубликованных контрольных списка перечисляют, что статья обязана сообщить. Проверка устройства добавляет четыре статистических показателя, и каждый отвечает на свой вопрос. Урок показывает, как читать и то и другое, на примере исследования носимых устройств 2025 года, а в конце объясняет, что может и чего не может значить слово «валидировано» на коробке.

Часть 1: можно ли воспроизвести исследование?

Начинайте с раздела о методах, а не с аннотации. Если непонятно, как удары записывали, очищали и анализировали, результат нельзя сравнить ни с каким другим, что бы он ни утверждал.

Контрольный список GRAPH (Quintana, Alvares и Heathers, 2016) написан для психиатрии, но подходит для любого исследования ВСР. В нём 13 пунктов в четырёх разделах. Отбор участников: как набирали и обследовали каждую группу, критерии включения, характеристики болезни и лекарства, демографические данные, включая физическую активность, алкоголь и никотин. Сбор межударных интервалов: оборудование и расположение электродов, исходная частота дискретизации Частота дискретизации Сколько раз в секунду датчик измеряет сигнал, в герцах (Гц). От неё зависит точность определения момента каждого удара: около ±1 мс при 1000 Гц, ±4 мс при 250 Гц и ±20 мс у камеры с частотой 25 Гц. Подробнее: Частота дискретизации , длительность записи, время суток, фильтрация, положение тела и инструкции. Подготовка данных: как рассчитывали интервалы и как их передискретизировали, как выявляли артефакты Артефакт Ошибка в записи сердечного ритма: пропущенный, лишний или неточно отмеченный удар, часто из-за движения или плохого контакта датчика. Артефакты могут сильно исказить показатели ВСР. Подробнее: Артефакт , почему данные были потеряны и как их очищали, с долей исправленных ударов. Расчёт ВСР: показатели, программы и любое логарифмическое преобразование, а также частотные диапазоны и их интерпретация. В статье отмечено, что порог потери данных 5 % распространён и не имеет эмпирического обоснования, а оценки HF Мощность HF Дисперсия интервалов между ударами в диапазоне от 0,15 до 0,40 Гц, в мс²: это диапазон обычного дыхания. Отражает в основном вагусную активность при частоте дыхания от 9 до 24 вдохов в минуту. Подробнее: Мощность HF , полученные быстрым преобразованием Фурье и по авторегрессионной модели, коррелируют с r = 0,96, поэтому выбор метода важен меньше, чем его подробное описание.

Контрольный список Catai (Catai и коллеги, 2020) содержит 30 пунктов, рассчитанных на кардиологию и кардиореабилитацию, и подробнее описывает условия: помещение (20–24 °C, влажность 40–60 %, шум, проходящие люди), не менее 15 минут покоя в положении записи, вещества, которых следует избегать накануне и в день записи (алкоголь, кофеин, никотин, шоколад, газированные напитки, энергетики), отсутствие тяжёлой нагрузки накануне, время суток, контроль частоты дыхания (при частоте ниже 8 вдохов в минуту диапазон LF Мощность LF Дисперсия интервалов между ударами в диапазоне 0,04–0,15 Гц, в мс²: ритмы, повторяющиеся каждые 7–25 секунд. Раньше её считали показателем симпатической активности, но теперь это оспаривают. Подробнее: Мощность LF интерпретировать нельзя), указанный критерий качества сигнала, число пропущенных и интерполированных ударов, а также то, задавали ли окна по времени или по числу ударов.

Рекомендации журнала Psychophysiology (Quigley и коллеги, 2024) добавляют минимальные требования к оборудованию: по умолчанию электрокардиограмма ЭКГ Запись электрической активности сердца с помощью электродов на коже. RR-интервалы, полученные по ЭКГ, — эталонный метод измерения ВСР. Подробнее: ЭКГ (ЭКГ) с частотой дискретизации 1000 Гц, фотоплетизмография ФПГ Оптический метод: свет направляют в кожу и улавливают пульсацию крови при каждом ударе сердца. Большинство часов, колец и приложений с камерой телефона оценивают по ней ВСР. Подробнее: ФПГ (ФПГ) и частота вплоть до 100 Гц — только там, где ЭКГ невозможна, число отредактированных ударов, любые пропуски дольше 5 секунд, проверка влияния редактирования, когда оно превышает 1 %, а для потребительских устройств — валидация относительно эталона.

Читайте статью с этими тремя списками перед глазами. Heathers (2014) изучил 97 исследований краткосрочной ВСР, опубликованных в 2012 году, и обнаружил, что перед записью только 45 % контролировали приём пищи, 23 % — воду и 4 % — полный мочевой пузырь. Исследование, которое сообщает о групповом различии HF без частоты дыхания, сообщает о различии, которое само же не может интерпретировать (урок 18).

Подробнее Вопросы о дизайне, которых нет в контрольных списках

Laborde и коллеги (2017) строят исследование вокруг трёх R: resting (уровень в покое), reactivity (изменение от исходного уровня к заданию) и recovery (изменение после него). Статья, измеряющая только одно из трёх, ответила лишь на треть вопроса, а эти три не взаимозаменяемы: более высокий вагусный тонус в покое обычно хорош, тогда как ослабление вагусного влияния под нагрузкой может быть адаптивным, а может и нет, в зависимости от задания. Спросите, какое R измеряли, совпадало ли исходное положение тела с положением во время задания и знал ли человек, очищавший удары, какое условие он очищает. Laborde рекомендует, чтобы он этого не знал.

Часть 2: четыре показателя согласия

Проверка устройства сравнивает новый прибор с эталоном по одним и тем же ударам. Появляются четыре числа, и они отвечают на разные вопросы.

Корреляция не измеряет согласие. r Пирсона говорит, растут ли и падают ли два показателя вместе. Часы, которые показывают ровно вдвое больше, чем ЭКГ, получают r = 1. Bland и Altman (1986) начали с этого замечания — «использование корреляции вводит в заблуждение» — и предложили метод, который её заменил.

Bland-Altman: смещение и границы согласия. Для каждой пары измерений возьмите разность dd и отложите её на графике против среднего значения этой пары. Средняя разность — это смещение, а границы согласия равны

dˉ±1.96 sd\bar{d} \pm 1.96\,s_d

Читайте так: средняя разность плюс-минус примерно два стандартных отклонения разностей. В работе 1986 года 1,96 округлено до 2. Если разности распределены примерно нормально, 95 % отдельных сравнений попадают внутрь этих границ. Смещение говорит, занижает или завышает устройство в среднем; границы говорят, насколько может ошибиться одно измерение. Валидация, в которой есть смещение без границ, сообщила среднее и скрыла разброс.

ICC: внутриклассовый коэффициент корреляции. Он рассматривает два прибора как двух оценщиков и спрашивает, какая доля общей дисперсии связана с настоящими различиями между людьми, а не с расхождением между оценщиками:

ICC=σbetween people2σbetween people2+σerror2\mathrm{ICC} = \frac{\sigma^2_{\text{between people}}}{\sigma^2_{\text{between people}} + \sigma^2_{\text{error}}}

Читайте так: дисперсия сигнала, делённая на сумму сигнала и шума. Koo и Li (2016) описывают десять форм ICC с разными допущениями и просят авторов указывать модель, тип и определение, которые они использовали, и приводить 95 % доверительный интервал. Их диапазоны: ниже 0,5 — плохо, от 0,5 до 0,75 — умеренно, от 0,75 до 0,9 — хорошо, выше 0,9 — отлично. ICC зависит от того, насколько различаются люди в выборке: широко разбросанная выборка делает устройство более надёжным на вид, чем то же устройство на похожих людях.

CCC: коэффициент корреляции конкордантности. Lin (1989) определил его как корреляцию между двумя измерениями вдоль линии в 45 градусов, проходящей через начало координат. Он умножает r Пирсона (прецизионность — насколько плотно облако) на поправочный множитель за смещение (точность — насколько облако удалено от линии полного согласия):

ρc=2ρ σxσyσx2+σy2+(μx−μy)2\rho_c = \frac{2\rho\,\sigma_x\sigma_y}{\sigma_x^2 + \sigma_y^2 + (\mu_x - \mu_y)^2}

Читайте так: удвоенная ковариация, делённая на две дисперсии плюс квадрат разности двух средних. Постоянный сдвиг или ошибка масштаба его снижают; плотное облако на линии даёт 1. Этот коэффициент использует исследование носимых устройств ниже.

MAPE: средняя абсолютная процентная ошибка. Среднее абсолютных разностей, каждая — в процентах от эталона:

MAPE=100n∑i=1n∣xi−riri∣\mathrm{MAPE} = \frac{100}{n}\sum_{i=1}^{n}\left|\frac{x_i - r_i}{r_i}\right|

Читайте так: на сколько процентов в среднем ошибается прибор. Её легко читать, но есть ловушка: для такого показателя, как RMSSD RMSSD Среднеквадратичное значение разностей соседних NN-интервалов: показатель ВСР в миллисекундах, построенный на том, насколько каждый интервал между ударами отличается от следующего. В основном отражает активность блуждающего нерва. Подробнее: RMSSD , который у одного человека может быть 15 мс, а у другого 90 мс, одна и та же ошибка в 5 мс составляет 33 % для первого и 6 % для второго.

Часть 3: разбор примера

Dial и коллеги (2025) попросили 13 здоровых взрослых (6 женщин, средний возраст 33 года) носить вместе с пятью носимыми устройствами нагрудный датчик Polar H10 (одноканальная ЭКГ с частотой дискретизации 1000 Гц), всего 536 ночей, и сравнили ночной RMSSD. Окна различались: кольца и Garmin усредняли 5-минутные фрагменты за ночь, WHOOP — среднее, в котором больший вес получает последний период медленного сна, а часы Polar — четыре часа после засыпания. Эталоном было значение датчика за всю ночь из Kubios, которое сохраняли только тогда, когда исправления требовали менее 5 % ударов.

УстройствоГраницы 95 % (мс)CCC (от 0 до 1)MAPE (%)
Oura Gen 4от −11,8 до +9,90,996,0
Oura Gen 3от −11,4 до +6,40,977,2
WHOOP 4.0от −12,5 до +10,90,948,2
Garmin Fenix 6от −15,2 до +11,60,8710,5
Polar Grit X Proот −14,3 до +23,6, как напечатано (см. ниже)0,8216,3

В статье границы приведены с двумя десятичными знаками; здесь они округлены до одного. Число ночей на устройство колебалось от 139 (Oura Gen 4) до 470 (Oura Gen 3), и все смещения в таблице статьи отрицательны: −1,0 мс для кольца Gen 4, −2,5 для Gen 3, −0,8 для WHOOP, −1,8 для Garmin и −4,7 для часов Polar. У четырёх устройств границы равны смещению плюс-минус 1,96 стандартного отклонения, как и должно быть. Строка Polar не такая: смещение −4,65 мс при стандартном отклонении 9,67 мс даёт границы от −23,6 до +14,3 мс, то есть напечатанные границы с переставленными знаками. Один из двух знаков в этой строке — опечатка, и в статье не сказано, какой именно. В любом случае ширина диапазона у часов Polar — около 38 мс, а у колец и WHOOP — примерно 18–23 мс.

Авторы оценивали CCC так: ниже 0,80 — неприемлемое согласие, 0,80–0,89 — низкое, 0,90–0,95 — умеренное, 0,95–0,99 — существенное, выше 0,99 — почти идеальное, и считали MAPE выше 10 % неприемлемой. Разберём таблицу по четырём показателям. У каждого устройства небольшое смещение: в среднем все они в пределах нескольких миллисекунд от датчика. Различают их границы: кольца и WHOOP удерживают 95 % ночей в пределах примерно 12 мс в обе стороны, тогда как ночи часов Polar разбросаны в диапазоне почти вдвое шире. CCC ранжирует их в этом порядке, потому что наказывает за разброс, а по MAPE два устройства выходят за порог, принятый авторами. Оговорки самого исследования: 13 человек, по всей видимости, здоровых, и фирменные алгоритмы, которые могут обновляться, поэтому авторы просят проводить повторную валидацию.

Часть 4: групповой результат и число одного человека

Таблица описывает, как устройства работают на многих людях. Ваш вопрос другой: достаточно ли точно это устройство, чтобы показать в ваших данных изменения, которые вам важны? Это зависит от числа, которое не даст ни одна валидация: размера вашего собственного разброса от дня ко дню.

Даже когда ничего не происходит, отдельные измерения колеблются от 37 до 57 мс. Среднее за 7 дней держится в пределах от 45 до 49 мс.

Добавьте событие
Смоделированный утренний RMSSD за последние 30 дн. Обычный диапазон — от 43 до 51 мс.
ДеньИзмерение (мс)Среднее за 7 дней (мс)Событие
29 дней назад4148Нет
28 дней назад4347Нет
27 дней назад4746Нет
26 дней назад4545Нет
25 дней назад4645Нет
24 дня назад5246Нет
23 дня назад5146Нет
22 дня назад4547Нет
21 день назад5248Нет
20 дней назад3747Нет
19 дней назад4246Нет
18 дней назад4045Нет
17 дней назад5245Нет
16 дней назад5746Нет
15 дней назад5347Нет
14 дней назад4246Нет
13 дней назад4146Нет
12 дней назад4647Нет
11 дней назад5349Нет
10 дней назад4648Нет
9 дней назад4346Нет
8 дней назад4545Нет
7 дней назад5146Нет
6 дней назад4647Нет
5 дней назад4747Нет
4 дня назад5247Нет
3 дня назад4146Нет
2 дня назад4046Нет
Вчера4546Нет
Сегодня5246Нет

Смоделированные измерения, а не реальные данные. Добавьте событие и сравните, насколько сдвигается одно утро и насколько — среднее за 7 дней.

График 30 смоделированных утренних измерений RMSSD, которые день ото дня скачут вверх и вниз, с более ровной линией среднего за 7 дней и затенённой полосой обычного диапазона. После графика идёт таблица со всеми измерениями.

Допустим, ваши утра разбросаны на ±20 мс вокруг среднего за 7 дней 60 мс. Устройство, чьи отдельные ночи лежат в пределах примерно ±12 мс от истины, добавляет меньше шума, чем у вас уже есть, и его тенденция пригодна. Допустим, ваш разброс ±5 мс вокруг 25 мс. Ошибка того же устройства теперь больше вашей собственной изменчивости, и изменение за одну ночь на 8 мс с равной вероятностью вызвано датчиком или вами. Средство то же, что в уроке 10: ориентируйтесь на среднее за 7 дней, а не на одну ночь. Смещение значит меньше, чем кажется, потому что устройство, всегда занижающее на 2,5 мс, всё равно верно показывает ваши изменения; оно важно только тогда, когда вы сравниваете своё значение с чужим или с таблицей норм, построенной на другом устройстве (урок 7).

Та же логика применима ко всем групповым выводам этого курса. Метаанализ, связывающий низкую ВСР с сердечно-сосудистыми событиями (урок 14), описывает популяцию. ICC 0,9 описывает, насколько хорошо устройство отличает людей друг от друга. Ни то ни другое ничего не говорит об одном измерении в одно утро.

Часть 5: как читать заголовок о риске

Исследования ВСР и здоровья заканчиваются числом другого рода: относительным риском. «У людей в группе с самой низкой ВСР риск был на 35 % выше» означает, что событие в этой группе было в 1,35 раза чаще, чем в группе сравнения. Но эта фраза не говорит, насколько часто событие встречается само по себе, а именно этот, абсолютный, риск определяет, о скольких людях идёт речь в заголовке.

Виджет ниже переводит относительный риск в число людей. Задайте исходный уровень: у скольких из 100 человек событие случилось бы и так; исходный уровень здесь придуман для примера, а не взят из чьего-то реального риска. Затем загрузите относительный риск из одного из трёх исследований, приведённых на странице ВСР и продолжительность жизни: метаанализ первых сердечно-сосудистых событий у людей без известной болезни сердца (Hillebrand 2013), метаанализ смерти от любой причины (Jarczok 2022) и исследование 1987 года о людях, переживших инфаркт (Kleiger). Посчитайте, сколько лишних случаев означает каждый из них.

Относительный риск в людях

Придуманный исходный риск, настоящие относительные риски

При исходном риске
4 из 100
Дополнительные случаи
1 из 100
  • Событие при исходном риске
  • Дополнительные при более высоком риске
  • Событие не происходит ни при каком риске

При исходном риске 4 из 100 относительный риск 1,35 означает около 5 из 100; дополнительных случаев — 1, а у 95 из 100 события нет ни при каком риске.

Ползунки работают при включённом JavaScript.

Сто человек. Задайте придуманный исходный риск и относительный риск из исследований выше и посмотрите, сколько дополнительных случаев стоит за громким числом. Относительные риски сравнивают группы; исходный риск для конкретного человека определяют инструменты, которыми пользуются врачи, а не ВСР.

Сетка 10 × 10 из 100 кружков. Закрашенные янтарные кружки — люди, у которых событие произошло бы при исходном риске, кружки с янтарным контуром — дополнительные случаи при более высоком относительном риске, серые кружки — люди, у которых события нет ни при каком риске. Один ползунок задаёт исходный риск от 1 до 30 из 100, другой — относительный риск от 1 до 6. Три кнопки загружают относительные риски из работ Hillebrand 2013, Jarczok 2022 и Kleiger 1987. Предложение подсчитывает людей в каждой группе.

Задавайте любому заголовку о риске три вопроса. Относительно кого: группы с самой высокой ВСР, остальной выборки или фиксированного порога? Из скольких: каков абсолютный риск в группе сравнения? И у кого: у здоровых взрослых или у пациентов после инфаркта, у которых ВСР измеряли медицинской ЭКГ за 24 часа, а не носимым устройством за ночь? Относительный риск 5 после инфаркта и 1,35 у здоровых взрослых — разные выводы о разных людях, и ни один из них не прогноз для одного человека.

Часть 6: что значит «валидировано» на коробке

«Клинически валидировано» — не регулируемая формулировка для потребительской ВСР. Встретив её, задайте шесть вопросов и ожидайте, что ответ на каждый находится в опубликованной статье, а не на странице продукта.

  1. Относительно чего? ЭКГ 1000 Гц или нагрудный датчик — эталон; другое носимое устройство — нет.
  2. На ком? Сколько людей, какого возраста, насколько здоровых и были ли среди них люди с нерегулярным ритмом. Исследование 2025 года — 13 взрослых со средним возрастом 33 года.
  3. Какой показатель и окно? Ночной RMSSD по 5-минутным фрагментам — не та же величина, что 5-минутное утреннее измерение или 24-часовой SDNN SDNN Стандартное отклонение всех NN-интервалов записи, в миллисекундах. Показывает общую вариабельность; его значение сильно зависит от длительности записи. Подробнее: SDNN (урок 17).
  4. Какой статистический показатель? Одна корреляция — не согласие. Ищите смещение с границами и ICC или CCC с интервалом.
  5. Какая прошивка? Алгоритмы меняются. Валидация прошивки прошлого года — исторический документ.
  6. Кто платил? Исследования на деньги производителя могут быть добротными, но о таком финансировании в них должно быть сказано.

Устройство, прошедшее все шесть проверок, валидировано для этого применения. Это всё равно не диагностический прибор, и тревожные признаки из урока 14 остаются в силе при любом его CCC.

Вы дошли до конца уровня 4 и всего курса. Страница сертификата проверяет все тесты, все двадцать, и выдаёт ваш сертификат, когда они пройдены.

Практическое задание

Найдите исследование валидации, на которое ссылается производитель носимого устройства, или самое новое для его модели, и ответьте на шесть вопросов: эталонное устройство, участники, показатель и окно, статистический показатель, дата прошивки и финансировал ли исследование производитель. Затем сравните границы согласия с величиной изменения, которое ваше приложение считает значимым.

Проверьте себя

Ответьте верно на 3 из 5 вопросов (60 %), чтобы пройти тест и отметить урок выполненным. Попыток может быть сколько угодно: варианты ответов каждый раз перемешиваются, а засчитывается лучший результат.

Вопросы ниже можно прочитать и без JavaScript, но для проверки ответов он нужен.

Вопрос 1 В проверке устройства для ночного RMSSD приведено только r = 0,96 между часами и ЭКГ. Что вы знаете о согласии?
Вопрос 2 В исследовании носимых устройств 2025 года у ночного RMSSD кольца Oura Gen 3 смещение составило около −2,5 мс, а границы согласия — примерно от −11 до +6 мс. Что это значит для одной ночи?
Вопрос 3 Какой диапазон ICC руководство Koo и Li называет хорошим?
Вопрос 4 На графике тенденции значения отдельных утр скачут вокруг более устойчивого среднего за 7 дней. Почему это важно при чтении проверки устройства?
Вопрос 5 На странице продукта сказано, что ВСР кольца «клинически валидирована». Какой вопрос наименее полезен?

Карта курса
  1. Что на самом деле делает ваше сердце10 мин
  2. Что такое ВСР — на одной картинке9 мин
  3. Нервная система, которая стоит за ВСР13 мин
  4. Почему более высокая ВСР не всегда лучше11 мин
  5. Как измеряют ВСР: от ударов к числу15 мин
  6. RMSSD, SDNN и lnRMSSD: простое объяснение13 мин
  7. ВСР и нормативные значения: возраст, пол, форма12 мин
  8. Почему устройства расходятся в оценке ВСР9 мин
  9. Как получить надёжное измерение ВСР11 мин
  10. Определение базового уровня ВСР и оценка тенденций12 мин
  11. Что меняет ВСР: самые важные факторы11 мин
  12. ВСР для тренировок и восстановления11 мин
  13. ВСР, стресс и дыхание15 мин
  14. Ограничения ВСР, мифы и тревожные признаки9 мин
  15. Итоговый тест по уровням 1–39 мин
  16. От электрода до числа: цепочка обработки сигнала18 мин
  17. Всё семейство показателей ВСР17 мин
  18. Частотный анализ без тумана17 мин
  19. Нелинейные показатели: DFA, энтропия и что они дают16 мин
  20. Как читать исследование ВСР или проверку устройства (вы здесь)19 мин

Источники

  1. Guidelines for Reporting Articles on Psychiatry and Heart rate variability (GRAPH): recommendations to advance research communication Translational Psychiatry, 2016
  2. Heart rate variability: are you using it properly? Standardisation checklist of procedures Brazilian Journal of Physical Therapy, 2020
  3. Publication guidelines for human heart rate and heart rate variability studies in psychophysiology, Part 1: physiological underpinnings and foundations of measurement Psychophysiology, 2024
  4. Statistical methods for assessing agreement between two methods of clinical measurement The Lancet, 1986
  5. A guideline of selecting and reporting intraclass correlation coefficients for reliability research Journal of Chiropractic Medicine, 2016
  6. A concordance correlation coefficient to evaluate reproducibility Biometrics, 1989
  7. Validation of nocturnal resting heart rate and heart rate variability in consumer wearables Physiological Reports, 2025
  8. Heart rate variability and cardiac vagal tone in psychophysiological research: recommendations for experiment planning, data analysis, and data reporting Frontiers in Psychology, 2017
  9. Everything Hertz: methodological issues in short-term frequency-domain HRV Frontiers in Physiology, 2014
  10. Heart rate variability and first cardiovascular event in populations without known cardiovascular disease: meta-analysis and dose-response meta-regression Europace, 2013
  11. Heart rate variability in the prediction of mortality: A systematic review and meta-analysis of healthy and patient populations Neuroscience & Biobehavioral Reviews, 2022
  12. Decreased heart rate variability and its association with increased mortality after acute myocardial infarction American Journal of Cardiology, 1987