Измерение направленного информационного потока в данных fNIRS-гиперсканирования с использованием метода когерентности с частичным вейвлет-преобразованием

Этот протокол описывает когерентность частичного вейвлет-преобразования (pWTC) для вычисления запаздывающего во времени паттерна межличностной нейронной синхронизации (INS) для вывода направления и временного паттерна информационного потока во время социального взаимодействия. Эффективность pWTC в устранении путаницы автокорреляции сигнала на INS была доказана двумя экспериментами.

Социальное взаимодействие имеет жизненно важное значение для людей. В то время как подход гиперсканирования широко используется для изучения межличностной нейронной синхронизации (INS) во время социальных взаимодействий, функциональная ближняя инфракрасная спектроскопия (fNIRS) является одним из самых популярных методов гиперсканирования натуралистических социальных взаимодействий из-за его относительно высокого пространственного разрешения, звуковой анатомической локализации и исключительно высокой толерантности к артефактам движения. Предыдущие исследования гиперсканирования на основе fNIRS обычно вычисляют запаздывающую во времени INS с использованием когерентности вейвлет-преобразования (WTC) для описания направления и временной структуры информационного потока между людьми. Однако результаты этого метода могут быть сбиты с толку автокорреляционным эффектом сигнала fNIRS каждого человека. Для решения этой проблемы был введен метод, называемый когерентностью частичного вейвлет-преобразования (pWTC), который был направлен на устранение эффекта автокорреляции и поддержание высокого разрешения временного спектра сигнала fNIRS. В этом исследовании сначала был проведен симуляционный эксперимент, чтобы показать эффективность pWTC в устранении влияния автокорреляции на INS. Затем было предложено пошаговое руководство по работе pWTC на основе набора данных fNIRS из эксперимента по социальному взаимодействию. Кроме того, было проведено сравнение между методом pWTC и традиционным методом WTC, а также сравнение между методом pWTC и методом причинности Грейнджера (GC). Результаты показали, что pWTC может быть использован для определения разницы INS между различными экспериментальными условиями и направленной и временной картины INS между людьми во время натуралистических социальных взаимодействий. Кроме того, он обеспечивает лучшее временное и частотное разрешение, чем традиционный ЦМТ, и лучшую гибкость, чем метод GC. Таким образом, pWTC является сильным кандидатом для вывода направления и временной модели информационного потока между людьми во время натуралистических социальных взаимодействий.

Социальное взаимодействие имеет жизненно важное значение для людей ^1,2. Для понимания нейрокогнитивного механизма социального взаимодействия двойного мозга в последнее время широко используется подход гиперсканирования, показывающий, что паттерны межличностной нейронной синхронизации (INS) могут хорошо характеризовать процесс социального взаимодействия ^{3,4,5,6,7,8,9,10,11,12 ,13,14}. Среди недавних исследований интересным открытием является то, что разница ролей людей в диаде может привести к замедленной во времени картине INS, то есть INS возникает, когда активность мозга одного человека отстает от активности мозга другого человека на секунды, например, от слушателей к динамикам ^5,9, от лидеров к последователям⁴, от учителей до учеников⁸, от матерей до детей^13,15, а от женщин до мужчин в романтической паре⁶. Самое главное, существует хорошее соответствие между интервалом запаздывающих во времени INS и поведения социального взаимодействия, например, между учителями, задающими вопросы, и учащимися, отвечающими^{на 8}, или между родительским поведением матерей и поведением соответствия детей¹⁵. Таким образом, запаздывающая во времени ИНС может отражать направленный информационный поток от одного человека к другому, как это предлагается в недавней иерархической модели межличностной вербальной^{коммуникации 16}.

Ранее замедленная по времени ИНС в основном рассчитывалась по функциональному сигналу ближней инфракрасной спектроскопии (fNIRS) из-за его относительно высокого пространственного разрешения, анатомической локализации звука и исключительно высокой толерантности артефактов движения¹⁷ при изучении натуралистических социальных взаимодействий. Более того, чтобы точно охарактеризовать соответствие между нейронным временным лагом и поведенческим временным лагом во время социального взаимодействия, важно получить силу INS для каждого временного лага (например, от отсутствия временного лага до временного лага 10 с). С этой целью ранее широко применялась процедура когерентности вейвлет-преобразования (ВТЦ) после смещения сигнала мозга одного индивидуума вперед или назад относительно сигнала другого индивидуума ^5,6,18. При использовании этой традиционной процедуры ЦМТ для сигналов fNIRS существует потенциальная проблема, поскольку наблюдаемая задержка во времени INS может быть спутана эффектом автокорреляции сигнала fNIRS для отдельного 19,20,21. Например, во время диадического процесса социального взаимодействия сигнал участника A в точке времени t может быть синхронизирован с сигналом участника B в той же точке времени. Между тем, сигнал участника A в точке времени t может быть синхронизирован с сигналом участника A в более поздней точке времени t+1 из-за эффекта автокорреляции. Таким образом, между сигналом участника A в точке времени t и сигналом участника B в точке времени t и сигналом участника B в точке времени t+1 может возникнуть ложная задержка во времени.

Миханович и его коллеги²² сначала представили метод, называемый когерентностью частичного вейвлет-преобразования (pWTC), а затем применили его в морской науке^23,24. Первоначальной целью этого метода было управление экзогенным смешанным шумом при оценке когерентности двух сигналов. Здесь, чтобы решить проблему автокорреляции в данных гиперсканирования fNIRS, метод pWTC был расширен для расчета ins с задержкой во времени на сигнале fNIRS. Точнее, запаздывающий во времени INS (и направленный поток информации) от участника A к участнику B может быть рассчитан с использованием уравнения ниже (уравнение 1)²³.

Здесь предполагается, что есть два сигнала, A и B, от участников A и B соответственно. Появление сигнала B всегда предшествует сигналу A с временным лагом n, где WTC (A_t, B_{t +n}) является традиционным WTC с задержкой во времени. ЦМТ (A_t, A_t+n) — автокоррелированный ЦМТ в участнике A. ЦМТ (A_t, B_t) — выровненный по времени ЦМТ в точке времени t между участником A и B. * — комплексный сопряженный оператор (рисунок 1A).

Рисунок 1: Обзор pWTC. (A) Логика pWTC. В диаде есть два сигнала А и В. Встречаемость A всегда следует за В с лагом n. Серый прямоугольник — это вейвлет-окно в определенной точке времени t или t+n. На основе уравнения pWTC (представленного на рисунке) необходимо рассчитать три WTC: ВТЦ с задержкой во времени A_t+n и B_t; автокоррелированный ЦМТ у участника A из A_t и A_t+n; и выровненный по времени ВТЦ в точке времени t, A_t и B_t. (Б) Компоновка комплектов оптодных зондов. CH11 был размещен в T3, а CH25 был размещен в T4 после международной системы 10-20^27,28. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Этот протокол впервые представил эксперимент по моделированию, чтобы продемонстрировать, насколько хорошо pWTC решает проблему автокорреляции. Затем он объяснил, как проводить pWTC шаг за шагом на основе эмпирического эксперимента натуралистических социальных взаимодействий. Здесь для введения метода использовался контекст коммуникации. Это связано с тем, что ранее задержка во времени INS обычно рассчитывалась в натуралистическом коммуникационном контексте ^{3,4,6,8,13,15,18}. Кроме того, было также проведено сравнение между pWTC и традиционным WTC и валидация с тестом причинности Грейнджера (GC).

Протокол эксперимента на людях был одобрен Советом по институциональному обзору и Комитетом по этике Государственной ключевой лаборатории когнитивной неврологии и обучения в Пекинском педагогическом университете. Все участники дали письменное информированное согласие до начала эксперимента.

1. Симуляционный эксперимент

1. Сгенерируйте два временных ряда сигналов, которые коррелируют друг с другом, причем один сигнал имеет автокорреляцию с временным лагом 4 с. Установите коэффициент корреляции r между двумя сигналами равным 0,4.
2. Кроме того, генерируйте два временных ряда сигналов без какой-либо корреляции, но с автокорреляцией в одном сигнале.
3. Рассчитайте значения традиционных 4-секундных INS с помощью уравнения 2 на основе сгенерированных сигналов с корреляцией или без нее, которые можно назвать INS_WTC с задержкой во времени с автокорреляцией и запаздывающей по времени базовой INS_WTC с автокорреляцией.
   ПРИМЕЧАНИЕ: Здесь традиционный ВТЦ с задержкой во времени выражается следующим уравнением (уравнение 2)²⁵
   
   где , C обозначает оператор непрерывного вейвлет-преобразования в разных масштабах i и временных точках t. S обозначает оператор сглаживания. * обозначает комплексный сопряженный оператор. W и M обозначают два отдельных временных ряда сигналов.
4. Удалите автокорреляцию из генерируемых сигналов. Затем рассчитайте значения традиционных 4-секундных INS_WTC с помощью уравнения 2 на основе сгенерированных сигналов с корреляцией или без нее, которые можно назвать временно запаздывающими INS_WTC без автокорреляции и отстающими по времени базовыми INSWTC без автокорреляции.
5. Рассчитайте значения pWTC с задержкой во времени 4 с помощью уравнения 3 на основе сгенерированных сигналов с корреляцией или без нее, названных ins_pWTC с задержкой во времени и базовых значений INS_pWTC с задержкой во времени.
   ПРИМЕЧАНИЕ: pWTC может быть рассчитан на основе следующего уравнения (уравнение 3)²³
   
   где ЦМТ (W_t, M_t+n) — традиционный ЦМТ с задержкой во времени. ЦМТ (M_t, M_t+n) — автокоррелированный ЦМТ одного человека. WTC (W_t, M_t) - это выровненный по времени ВТЦ. * — комплексный сопряженный оператор.
6. Повторите вышеуказанные процедуры 1000 раз.
7. После вычитания базового INS сравните результаты ins_wwtc с задержкой во времени с автокорреляцией, INSWTC без автокорреляции и ins_pWTC с задержкой во времени с использованием метода дисперсионного анализа (ANOVA).
   ПРИМЕЧАНИЕ: Здесь ожидается, что запаздывающий по времени INS_WTC с автокорреляцией будет значительно выше, чем INS_WTC без автокорреляции и INS_pWTC с задержкой во времени, и не ожидается существенной разницы между ins_wtc без автокорреляции и ins_pWTC с задержкой во времени.

2. Эмпирический эксперимент

1. Участники и процедура
   1. Наберите соответствующих участников.
      ПРИМЕЧАНИЕ: В этом исследовании двадцать две пары близких друзей противоположного пола (средний возраст женщин = 20,95, стандартное отклонение (SD) = 1,86; средний возраст мужчин = 20,50, SD = 1,74) были набраны через рекламу от студентов университетов в Пекине. Все участники были правшами и имели нормальное или скорректированное до нормального зрение. Кроме того, ни у одного участника не было никаких языковых, неврологических или психических расстройств.
   2. Попросите каждую пару участников сесть лицом к лицу во время эксперимента. Попросите их свободно общаться на поддерживающую тему в одной сессии и на тему конфликта в другой сессии.
      ПРИМЕЧАНИЕ: Темы использовались для индуцирования предполагаемой положительной или отрицательной эмоциональной валентности. Каждый сеанс общения длился 10 минут, а порядок тем был уравновешен.
   3. Попросите участников сообщить о поддерживающих и конфликтующих темах в качестве стандартного правила настройки. Попросите каждого партнера оценить положительный или отрицательный уровень валентности, который мог быть вызван по определенной балльной шкале. Затем ранжируйте темы, о которых сообщается, в соответствии с рейтингом.
      ПРИМЕЧАНИЕ: В этой работе темы были выбраны следующими тремя шагами. Во-первых, для поддерживающих тем каждый участник должен был сообщить о 1-3 личных проблемах, связанных с тем, что он / она хочет улучшить в своей жизни. Каждый участник должен был сообщить о 1-3 случаях, которые вызвали или могли бы вызвать конфликт между ними или которые могли бы поставить под угрозу их отношения по темам конфликта. Во-вторых, каждый партнер должен был оценить уровень положительной или отрицательной валентности, которую каждая тема может вызвать по 7-балльной шкале (1 = совсем нет, а 7 = очень много). В-третьих, заявленные темы были ранжированы в соответствии с рейтингом. Были выбраны первые две темы в списке вспомогательных тем и тем конфликтов.
2. Сбор данных fNIRS
   1. Используйте 26-канальную топографическую систему fNIRS (см. Таблицу материалов) для сбора данных fNIRS.
      ПРИМЕЧАНИЕ: Два индивидуальных оптодных зонда покрывали двустороннюю лобную, височную и теменную кору (рисунок 1B).
   2. Точнее, попросите каждого участника надеть кепку с двумя индивидуальными наборами зондов (см. Таблицу материалов).
   3. Выровняйте мастоиды носа, иона и уха с Fpz, Opz, T7 и T8, которые являются типичными ориентирами 10-20 международной системы²⁶.
   4. Выровнять канал (CH) 11 к T3 и CH25 к T4 в соответствии с международной системой 10-20 для двух наборов зондов^27,28.
   5. Проверка анатомического местоположения наборов зондов путем сканирования данных магнитно-резонансной томографии (МРТ) от типичного участника с быстрой последовательностью градиент-эхо, подготовленной с высоким разрешением T1-взвешенной намагниченностью (TR = 2530 мс; TE = 3,39 мс; угол сальто = 7°; толщина среза = 1,3 мм; размер вокселя = 1,3 x 1 x 1,3 мм).
   6. Используйте статистическое параметрическое отображение 12 (SPM12) для нормализации изображения до стандартного пространства координат Монреальского института визуализации (координата MNI)²⁹. Затем используйте набор инструментов NIRS_SPM (см. Таблицу материалов) для проецирования координат MNI зондов в шаблон автоматизированной анатомической маркировки (AAL).
   7. Сбор данных о оптической плотности ближнего инфракрасного света на трех длинах волн (780, 805 и 830 нм) с частотой дискретизации 55,6 Гц (параметры оборудования по умолчанию).
   8. Проверьте качество сигнала с помощью встроенного программного обеспечения топографической системы fNIRS (см. Таблицу материалов).
   9. Начните запись сигнала.
      ПРИМЕЧАНИЕ: Некоторые опубликованные протоколы продемонстрировали, как собирать сигналы fNIRS с помощью различного оборудования и систем 30,31,32.
3. Предварительная обработка данных fNIRS
   1. Экспортируйте файлы данных с оборудования.
      ПРИМЕЧАНИЕ: В текущем эксперименте встроенное программное обеспечение автоматически преобразовывало полностью оптические данные плотности в изменения концентрации оксигемоглобина (HbO) на основе модифицированного закона Бира-Ламберта.
   2. Удалите первые и последние 15 секунд данных для каждого сеанса, чтобы избежать временных ответов.
   3. Используйте встроенную функцию MATLAB decimate для понижения уровня данных с 55,6 Гц до 11,1 Гц.
      ПРИМЕЧАНИЕ: Схемы спектра мощности между 55,6 Гц и 11,1 Гц довольно похожи (дополнительный рисунок 1).
   4. Используйте встроенную функцию приложения MATLAB (Homer3, см. Таблицу материалов) с соответствующей функцией фильтрации, чтобы применить метод фильтра дискретного вейвлет-преобразования для коррекции артефактов движения.
   5. Используйте встроенную функцию MATLAB pca для удаления глобального физиологического шума. Удалите верхние 80% дисперсии из сигналов.
   6. Удалить физиологический шум на основе предыдущих исследований³³. Точно удалите полосы частот каждого сигнала выше 0,7 Гц, чтобы избежать сглаживания высокочастотного физиологического шума (например, сердечной деятельности).
   7. Затем удалите полосы частот каждого сигнала ниже 0,01 Гц, чтобы отфильтровать очень низкочастотные колебания.
   8. Наконец, удалите полосы частот каждого сигнала в пределах 0,15-0,3 Гц, чтобы исключить потенциальное влияние дыхательной активности.
4. Обработка данных первого уровня fNIRS
   1. Во-первых, рассчитайте INS, используя традиционный WTC (INS_WTC).
      ПРИМЕЧАНИЕ: Здесь было предсказано, что запаздывающая во времени картина INS_WTC под руководством женщин возникает между мозговой активностью женщин и мужчин, потому что предыдущие исследования предполагали разные роли женщин и мужчин во время разговора^34,35. Традиционный ЦМТ рассчитал эту модель INS_WTC, сместив активность мозга мужчин назад по сравнению с активностью женщин (см. Уравнение 2).
   2. Рассчитайте женское значение 2 s-laged INS_WTC после удаления начальных 2 s данных от женщин и последних 2 s данных от мужчин с помощью уравнения 2. Аналогичным образом, после удаления начальных 2 с данных от мужчин и последних 2 с данных от женщин, рассчитайте мужское значение INS_WTC 2 s с уравнением 4.
      ПРИМЕЧАНИЕ: Здесь использовалась функция wcoherence, которая является встроенной функцией вейвлет-инструментария MATLAB (см. Таблицу материалов).
   3. Повторите эту процедуру с различными временными лагами n, т.е. n = 2 с, 4 с, 6 с, 8 с во всех потенциальных парах CH (например, CH2 у женщин и CH10 у мужчин, всего 676 пар). Кроме того, вычислите силу запаздывающего во времени INS_WTC под руководством мужчин таким же образом (уравнение 4).
      
   4. Во-вторых, рассчитайте INS с помощью pWTC (INS_pWTC).
      ПРИМЕЧАНИЕ: pWTC был рассчитан на основе уравнения 3. Расчет INS_pWTC повторялся с различными временными лагами n, т.е. n = 2 с, 4 с, 6 с, 8 с по всем потенциальным парам каналов (например, CH2 у женщин и CH10 у мужчин, всего 676 пар). Кроме того, сила управляемого мужчинами INS_pWTC была рассчитана таким же образом (уравнение 5).
      
   5. Генерация временных рядов сигналов fNIRS с задержкой во времени при различных временных лагах.
   6. Рассчитайте значения ЦМТ с задержкой во времени при разных временных лагах.
   7. Генерация автокоррелированных временных рядов сигналов fNIRS с различными временными лагами. Чтобы рассчитать значение 2 s-автокорреляции для мужчин, удалите первые 2 s данных у мужчин и последние 2 s данных у мужчин.
   8. Вычисление автокоррелированных значений ЦМТ при разных временных лагах.
   9. Генерация выровненных по времени временных рядов сигналов fNIRS с различными временными лагами. Чтобы рассчитать 2 с выровненного по времени ЦМТ, удалите первые 2 с данных из данных мужчин и первые 2 с данных женщин.
   10. Рассчитайте выровненные по времени значения ЦМТ.
   11. Введите выровненные по времени значения WTC, WTC с задержкой по времени и автокоррелированные значения WTC при разном временном лаге в уравнение 3 и уравнение 5 - уравнение pWTC, генерируя INS_pWTC.
   12. Наконец, рассчитайте INS с использованием метода GC (INS_GC).
       ПРИМЕЧАНИЕ: Для дальнейшей проверки метода pWTC и оценки его преимуществ и недостатков INS на основе GC был рассчитан с использованием метода GC (INS_GC).
   13. Основываясь на результате pWTC, полоса пропускания фильтрует сигнал HbO каждого человека на SMC (т.е. 0,4-0,6 Гц, см. Репрезентативные результаты).
   14. Проведите тест GC (Econometric toolbox, MATLAB) внутри каждой диады в поддерживающей и конфликтной темах отдельно.
       ПРИМЕЧАНИЕ: Для INS_GC получены четыре группы F-значений: (1) от женщин к мужчинам по поддерживающей теме (W2M_supp); 2) от мужчин к женщинам по вспомогательной теме (M2W_supp); 3) от женщин к мужчинам по теме конфликта (W2M_conf); и 4) от мужчин к женщинам по теме конфликта (М2 W_conf). F-значения используются для индексации_{INS GC}.
5. Обработка данных fNIRS второго уровня
   1. Преобразуйте INS с преобразованием Фишера-z, а затем усредните INS во временной размерности.
      ПРИМЕЧАНИЕ: Здесь преобразование Fisher-z проводилось с помощью пользовательского скрипта MATLAB с equation 6³⁶:
      
      где r — значение WTC или pWTC, а z — преобразованное значение Fisher-z WTC или pWTC.
   2. Для усредненных INS при каждом временном лаге проведите парный t-тест с двумя выборками (поддерживающий против конфликта) на каждой паре CH в диапазоне частот. Затем определите все значимые частотные кластеры (P < 0,05).
   3. Проведите кластерный тест перестановки, чтобы установить пороговое значение для результатов.
      1. Перераспределить диадические отношения, случайным образом распределив участников по новым парам из двух членов, то есть участникам диады, которые никогда не общались друг с другом. Пересчитывайте INS при каждом временном лаге, снова выполняйте парные t-тесты в новой выборке и снова идентифицируйте значительные частотные кластеры.
      2. Выберите кластер с наибольшим суммированным t-значением. Повторите вышеуказанные процедуры 1000 раз, чтобы сгенерировать нулевое распределение максимальных ложноположительных t-значений.
         ПРИМЕЧАНИЕ: Распределение служит в качестве уровня шанса. Частота ошибок семейства (FWER) контролируется при q = 0,05, что означает, что только верхние 5% нулевого распределения ложноположительных t-значений превышают порог (R*).
      3. Сравните суммированное t-значение каждого идентифицированного частотного кластера в исходной выборке с нулевым распределением для получения значимых статистических результатов.
   4. Проведите контекстный (поддерживающий, конфликтный) x направление (женщины к мужчинам, мужчины к женщинам) анализ дисперсии (ANOVA) для проверки разницы в направлении INS между различными условиями (т.е. темами) (p < 0,05).
   5. Провести парный двухвыборный двуххвостый t-тест между результатами WTC (Wt, Mt + n) и WTC (Mt, Mt + n) для проверки потенциального влияния автокорреляции на INS.
      ПРИМЕЧАНИЕ: INS ЦМТ (Mt, Mt + n) отражает автокорреляцию.

Результаты моделирования
Результаты показали, что запаздывающий по времени INS_WTC с автокорреляцией был значительно выше, чем INS_WTC без автокорреляции (t(1998) = 4,696, p < 0,001) и ins_pWTC (t(1998) = 5,098, p < 0,001). Кроме того, не было существенной разницы между запаздываемым во времени INS_WTC без автокорреляции и INS_pWTC (t(1998) = 1,573, p = 0,114, рисунок 2A). Эти результаты показывают, что pWTC может эффективно устранить влияние эффекта автокорреляции на INS. Кроме того, когда значение WTC было установлено близким к 0 или 1, замедленный по времени INS_pWTC по-прежнему показывал надежные результаты, когда значение WTC было далеко от 0 или 1 (дополнительный рисунок 2).

Результаты эмпирического эксперимента
Шаблон INS с использованием традиционного метода ЦМТ
Результаты показали, что при 0,04-0,09 Гц INS_WTCв сенсомоторной коре (SMC, CH20) как женщин, так и мужчин был значительно выше в поддерживающей теме, чем в конфликтной теме, когда активность мозга мужчин отставала от активности мозга женщин на 2 с, 4 с и 6 с (2 с: t(21) = 3,551, p = 0,0019; лаг 4 с: t(21) = 3,837, p = 0,0009; лаг 6 с: t(21) = 3,725, p = 0,0013). Кроме того, при 0,4-0,6 Гц INS_WTC в SMC был значительно выше в теме конфликта, чем в поддерживающей теме, когда мозговая активность мужчин отставала от женской на 4 с (t(21) = 2,828, p = 0,01, рисунок 2B).

Кроме того, для сравнения направления INS_WTC в разных темах, тема (поддерживающая, конфликтная) x направление (женщины к мужчинам, мужчины к женщинам) ANOVA была впервые проведена на INS_WTC SMC с временным лагом 2-6 с. Результаты 0,04-0,09 Гц не показали каких-либо значительных эффектов взаимодействия при любом временном лаге (ps > 0,05). Для диапазона частот 0,4-0,6 Гц результаты показали, что эффект взаимодействия был незначительно значимым (F(1, 21) = 3,23, p = 0,086). Парные сравнения показали, что INS_WTC от женщин к мужчинам был значительно выше в теме конфликта, чем в поддерживающей теме (M.D. = 0,014, S.E. = 0,005, p = 0,015), тогда как INS_WTC от мужчин к женщинам существенно не отличался между темами (M.D. = 0,002, S.E. = 0,006, p = 0,695).

Наконец, чтобы проверить влияние автокорреляции на результаты традиционных запаздывающих по времени INS_WTC, INS_WTC сравнивали между WTC(W_t, M_t+4) и WTC(M_t, M_t+4) при 0,04-0,09 Гц и 0,4-0,6 Гц соответственно. Обратите внимание, что INS_WTC WTC (M_t, M_t+4) отражает автокорреляцию. Результаты показали, что при 0,4-0,6 Гц не было существенной разницы между INS_WTC (W_t, M_t+4) и WTC(M_t, M_t+4) (t(21) = 0,336, p = 0,740). При 0,04-0,09 Гц INS_{WTC WTC} (Mt, Mt+4) был значительно выше, чем у WTC (W_t, M_t+4) (t(21) = 4,064, p < 0,001). Также было проведено сравнение диапазонов частот 0,04-0,09 Гц и 0,4-0,6 Гц относительно INS_WTC WTC (M_t, M_t+4). Результаты показали, что INS_WTC WTC (M_t, M_t+4) был значительно выше при 0,04-0,09 Гц, чем при 0,4-0,6 Гц (t(21) = 5,421, p < 0,001). Эти результаты показывают, что задержка во времени INS_WTC была затронута автокорреляцией как в низкочастотном, так и в высокочастотном диапазонах, но воздействие было больше для диапазона более низких частот, чем для диапазона более высоких частот.

Шаблон INS с использованием метода pWTC
Результаты показали, что разница в INS_pWTC между конфликтными и поддерживающими темами достигла значимости в SMC как женщин, так и мужчин при 0,4-0,6 Гц, когда мужская мозговая активность отставала от активности мозга женщин на 4 с (t(21) = 4,224, p = 0,0003). при 0,04-0,09 Гц; однако не было найдено никаких существенных результатов, равно как и их эффективных результатов в других частотных диапазонах (Ps > 0,05, рисунок 2C).

Дополнительный тест ANOVA был проведен на INS_pWTC SMC при 0,4-0,6 Гц. Результаты показали, что взаимодействие между темой и направлением было незначительно значимым (F(1,21) = 3,48, p = 0,076). Дальнейшие попарные сравнения показали, что INS_pWTC от женщин к мужчинам был значительно выше в теме конфликта, чем в поддерживающей теме (M.D. = 0,016, S.E. = 0,004, p = 0,002), тогда как INS_pWTC от мужчин к женщинам существенно не различался между темами (M.D. = 0,0007, S.E. = 0,006, p = 0,907, рисунок 2D).

Шаблон INS с использованием метода GC
Тест ANOVA был проведен на INS_GC в SMC только в пределах 0,4-0,6 Гц. Результаты показали значительное взаимодействие между темой и направлением (F(1,21) = 8,116, p = 0,010). Парный анализ показал, что INS_GC от женщин к мужчинам был значительно выше в конфликтной теме, чем в поддерживающей теме (MD = 5,50, SE = 2,61, p = 0,043). Напротив, INS_GC от мужчин к женщинам существенно не отличался между темами (MD = 1,42, SE = 2,61, p = 0,591, рисунок 2E).

Рисунок 2: Результаты моделирования и эмпирического эксперимента. (А) Результаты моделирования трех смоделированных образцов. Задержка во времени INS_WTC с автокорреляцией была значительно выше, чем у INS_WTC без автокорреляции и INS_pWTC. Не было существенной разницы между запаздывающей по времени INS_WTC без автокорреляции и pWTC. (B) T-карта INS_WTC в эмпирическом эксперименте, показывающая значительные контекстные эффекты в пределах 0,04-0,09 Гц, когда активность SMC мужчин отставала от активности женщин на 2-6 с. Наблюдался также незначительно значительный контекстный эффект в пределах 0,4-0,6 Гц, когда активность SMC мужчин отставала от активности женщин на 4 с. (C) T-карта INS_pWTC, показывающая значительный контекстный эффект в пределах 0,4-0,6 Гц, когда активность SMC мужчин отставала от активности женщин на 4 с. (D) Сравнение направленного INS_pWTC по различным темам pWTC. Направленная ИНС от женщин к мужчинам значительно выше в конфликтных контекстах, чем в поддерживающих. Е) Валидация направленных ИНС с помощью испытания ГК (ИНС_ГК). Результирующая картина INS_GC аналогична INS_pWTC. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Дополнительный рисунок 1: График спектра мощности для частоты дискретизации при 11,1 Гц (синяя линия) и 55,6 Гц (красная линия). Схема энергетического спектра для этих двух довольно похожа. Пожалуйста, нажмите здесь, чтобы загрузить этот файл.

Дополнительный рисунок 2: Карты pWTC пола и потолка ВТЦ. (A) Левая панель: карта ВТЦ с задержкой во времени, генерируемая двумя одинаковыми сигналами, ось X - это точка времени, а ось Y - полоса частот. Среднее значение ЦМТ во всех точках составляет ~1. Правая панель: pWTC карта двух похожих сигналов. Карта pWTC очень похожа на карту WTC. (B) Левая панель: карта ВТЦ с задержкой во времени, генерируемая двумя случайными сигналами, ось X - это точка времени, а ось Y - полоса частот. Среднее значение ЦМТ во всех точках равно ~0. Правая панель: pWTC карта двух похожих сигналов. Карта pWTC очень похожа на карту WTC. Пожалуйста, нажмите здесь, чтобы загрузить этот файл.

В исследованиях гиперсканирования обычно важно описать направленные и временные закономерности информационного потока между людьми. Большинство предыдущих исследований гиперсканирования fNIRS использовали традиционный WTC²⁵ для вывода этих характеристик путем расчета INS с задержкой во времени. Однако, как одна из неотъемлемых особенностей сигнала fNIRS^20,21, эффект автокорреляции может сбить с толку задержку во времени INS. Для решения этой проблемы в настоящем протоколе был введен способ, называемый pWTC²². Этот метод оценивает запаздывающую по времени INS после частичной автокорреляции и сохраняет преимущества метода WTC. Этот протокол предлагает пошаговое руководство о том, как проводить pWTC и проверяет результаты pWTC, сравнивая его результаты с результатами традиционных тестов WTC и GC.

В этом протоколе демонстрируются критические этапы применения pWTC в данных гиперсканирования на основе fNIRS. В частности, во-первых, чтобы рассчитать запаздывающий по времени ЦМТ, автокоррелированный ЦМТ и выровненный по времени ЦМТ должны быть рассчитаны на основе временных рядов fNIRS с задержкой во времени. Далее pWTC вычисляются на разных временных лагах в соответствии с уравнением 1. Результаты pWTC возвращают матрицу частот времени x, а значения в матрице находятся в диапазоне от 0 до 1. Таким образом, на этих значениях могут быть проведены дальнейшие статистические тесты.

В демонстрационном протоколе репрезентативные результаты традиционного ЦМТ показали два значительных эффекта в двух полосах частот: 0,4-0,6 Гц. Однако воздействие в пределах 0,04-0,09 Гц не выдержало порога в результатах pWTC, предполагая, что этот эффект может быть сомкрапирован автокорреляционным эффектом сигнала fNIRS. С другой стороны, результаты в диапазоне 0,4-0,6 Гц были хорошо воспроизведены методом pWTC. Эти результаты показывают, что после устранения эффекта автокорреляции pWTC обеспечивает более чувствительные и специфические изменения в выводе направленных и временных паттернов INS между людьми. Другая возможность, однако, заключается в том, что pWTC не восприимчив к направленным и временным паттернам INS в более низких частотных диапазонах, чем в более высоких частотных диапазонах, что приводит к недооценке эффекта INS. Для дальнейшего уточнения этих возможностей необходимы будущие исследования.

Сравнение с тестом GC еще раз подтверждает этот вывод. Результаты теста GC были очень похожи на результаты pWTC, показывая важный поток информации от женщин к мужчинам, но не от мужчин к женщинам. Была небольшая разница между результатами теста GC и pWTC, т.е. эффект взаимодействия между темой и направлением был незначительно значительным в результатах pWTC, но достиг значения в тесте GC. Это различие может быть связано с тем, что pWTC рассчитывается в более тонком масштабе времени, чем тест GC. Таким образом, хотя тесты pWTC и GC могут обеспечить надежные результаты при контроле эффекта автокорреляции, pWTC является преимуществом, поскольку нет необходимости делать стационарные предположения и имеет структуру с высоким временным спектром.

Метод pWTC также имеет свои ограничения. Подобно тесту GC, причинность, выведенная из pWTC, не является реальной причинностью^37,38. Вместо этого он указывает только на временную связь между сигналами А и В. Эту проблему следует иметь в виду при применении метода pWTC. Во-вторых, pWTC лишь частично устраняет эффект автокорреляции. Таким образом, другие потенциальные параллельные переменные, такие как общие среды или аналогичные действия, все еще могут влиять на результаты. Следовательно, выводы о направлении и временной структуре информационного потока должны быть сделаны после контроля этих смешанных факторов.

Кроме того, существовали некоторые сложные проблемы, связанные с предварительной обработкой данных fNIRS. Хотя fNIRS имеет высокую устойчивость к движениям головы, артефакты движения по-прежнему являются наиболее значительным источником шума³⁹. Большие движения головы по-прежнему приводят к смещению положения оптодов, создавая артефакты движения, такие как резкие шипы и смещения базовой линии. Для решения этих проблем были разработаны многие подходы к коррекции артефактов, такие как сплайн-интерполяция⁴⁰, фильтрация на основе вейвлетов³⁹, анализ основных^{компонентов 41} и улучшение сигнала на основе корреляции⁴² и т. Д. Купер и его коллеги⁴³ сравнили эти подходы, основанные на реальных данных fNIRS состояния покоя, и обнаружили, что фильтрация на основе вейвлетов обеспечивает наибольшее увеличение отношения контрастности к шуму. Кроме того, Бригадой и ее коллеги⁴⁴ также сравнили эти подходы в реальных данных лингвистических задач, а также обнаружили, что фильтрация на основе вейвлетов была наиболее эффективным подходом в коррекции артефактов движения. Таким образом, в этом исследовании была применена вейвлет-фильтрация, а также рекомендована для будущих исследований гиперсканирования fNIRS.

В целом, pWTC является ценным подходом в оценке направленных и временных закономерностей информационного потока во время социального взаимодействия. Что еще более важно, считается, что метод pWTC также подходит для псевдо-гиперсканирующих исследований (т.е. сигналы двух или нескольких мозгов не собираются одновременно^45,46). В таких экспериментах, хотя направление информационного потока фиксировано, также представляет интерес изучение длительности временного лага между входом сигнала и процессом сигнала. Таким образом, автокорреляция также может сбить с толку результаты запаздывающих по времени INS. В будущем этот метод может ответить на многие вопросы в гиперсканировании и других исследованиях межмозгов. Например, определить доминирующую роль в различных социальных отношениях, таких как учителя и студенты, врачи и пациенты, а также исполнители и аудитории. Кроме того, поскольку pWTC поддерживает временные структуры INS, также можно проверить динамическую картину INS, такую как конвергенция групповых отношений.

Авторы заявляют об отсутствии конкурирующих финансовых интересов.

Эта работа была поддержана Национальным фондом естественных наук Китая (61977008) и программой Young Top Notch Talents of Ten Thousand Talent Program.