Программа для спектрального анализа максимальной мощности на ЭЭГ
Спектральный анализ Мощность ЭЭГ Обработка Визуализация
ЛИСТИНГ
library(tidyverse)
## Warning: пакет 'tidyverse' был собран под R версии 4.1.3
## -- Attaching packages --------------------------------------- tidyverse 1.3.1 --
## v ggplot2 3.3.5     v purrr   0.3.4
## v tibble  3.1.6     v dplyr   1.0.8
## v tidyr   1.2.0     v stringr 1.4.0
## v readr   2.1.2     v forcats 0.5.1
## Warning: пакет 'ggplot2' был собран под R версии 4.1.3
## Warning: пакет 'tibble' был собран под R версии 4.1.3
## Warning: пакет 'tidyr' был собран под R версии 4.1.3
## Warning: пакет 'readr' был собран под R версии 4.1.3
## Warning: пакет 'purrr' был собран под R версии 4.1.3
## Warning: пакет 'dplyr' был собран под R версии 4.1.3
## Warning: пакет 'stringr' был собран под R версии 4.1.3
## Warning: пакет 'forcats' был собран под R версии 4.1.3
## -- Conflicts ------------------------------------------ tidyverse_conflicts() --
## x dplyr::filter() masks stats::filter()
## x dplyr::lag()    masks stats::lag()
library(readxl)
## Warning: пакет 'readxl' был собран под R версии 4.1.3
data<- read_excel("Эксель_мощности_спектров.xlsx")
data <- select(data,-"full range")
data <- data[-c(19:24), ]
data$name <- factor(data$name)
data <- data %>% 
    mutate(name = fct_relevel(name,
                           "Шанин",
                           "Дорофеев",
                           "Горин",
                           "Ширяев",
                           "Романов",
                           "Муратов",
                           "Сысоев",
                           "Аракелян С"
         )) 
levels(data$name) <- c("1","2","3","4","5","6","7","8")

dataplot<-filter(data, state == "Покой" | state == "Доска")
names(dataplot)[1]<- "Испытуемый" 
names(dataplot)[4:8]<- c("deltaritm", "tetaritm", "alpharitm","bettalf","bettahf")
dataplot$new_var <- paste(dataplot$state, dataplot$names, sep="-")
dataplot$new_var <- factor(dataplot$new_var)
levels(dataplot$new_var)
## [1] "Доска-EEG FP1-A1 [FP1]" "Доска-EEG FP2-A2 [FP2]" "Покой-EEG FP1-A1 [FP1]"
## [4] "Покой-EEG FP2-A2 [FP2]"
dataplot <- dataplot %>%
    mutate(new_var = fct_relevel(new_var,
                           "Доска-EEG FP2-A2 [FP2]",
                           "Покой-EEG FP2-A2 [FP2]",
                           "Доска-EEG FP1-A1 [FP1]",
                           "Покой-EEG FP1-A1 [FP1]",
         ))
library(ggridges)
## Warning: пакет 'ggridges' был собран под R версии 4.1.3
ggplot(data=dataplot,aes(x=deltaritm, y = new_var,fill = stat(x))) +
  geom_density_ridges_gradient(from = 0, rel_min_height = 0.01, gradient_lwd = 1.) +
  scale_fill_viridis_c(name = "Значение", option = "C") +
   labs(title="Дельта Ритм",x="МкВт",y="Состояния") +
    scale_x_continuous(expand = c(0.2, 0)) +
  theme_bw()
## Picking joint bandwidth of 356
ggplot(data=dataplot,aes(x=tetaritm, y = new_var,fill = stat(x))) +
  geom_density_ridges_gradient(from = 0,rel_min_height = 0.01, gradient_lwd = 1.) +
  scale_fill_viridis_c(name = "Значение", option = "A") +
         scale_x_continuous(expand = c(0.2, 0)) +
   labs(title="Тетта Ритм",x="МкВт",y="Состояния") +
  theme_bw()
## Picking joint bandwidth of 6.72
ggplot(data=dataplot,aes(x=alpharitm, y = new_var,fill = stat(x))) +
  geom_density_ridges_gradient(from = 0,rel_min_height = 0.01, gradient_lwd = 1.) +
  scale_fill_viridis_c(name = "Значение", option = "В") +
      scale_x_continuous(expand = c(0.2, 0)) +
   labs(title="Альфа Ритм",x="МкВт",y="Состояния") +
  theme_bw()
## Warning in viridisLite::viridis(n, alpha, begin, end, direction, option): Option
## 'В' does not exist. Defaulting to 'viridis'.
## Picking joint bandwidth of 1.26
ggplot(data=dataplot,aes(x=bettalf, y = new_var,fill = stat(x))) +
  geom_density_ridges_gradient(from = 0,rel_min_height = 0.01, gradient_lwd = 1.) +
  scale_fill_viridis_c(name = "Значение", option = "E") +
       scale_x_continuous(expand = c(0.2, 0)) +
   labs(title="Бетта НЧ Ритм",x="МкВт",y="Состояния") +
  theme_bw()
## Picking joint bandwidth of 0.305

Описание программы:
Спектральный анализ является одним из стандартных методов, используемых для количественной оценки ЭЭС. Спектральная плотность мощности (спектр мощности) исследуемых ЭЭС в пределах отдельных ритмов отражает «частотный состав» сигнала или распределение мощности сигнала по частоте.
Основной частью алгоритма является обработка таблиц, содержащих максимальные мощности ЭЭС в различных частотных диапазонах. Программа позволяет вывести данные в графическом виде, что упрощает последующий анализ данных человеком.
Проводится спектральный анализ максимальной мощности ЭЭС в различных частотных диапазонах: альфа (8-12 Гц), бета (12-35 Гц), тета (4-8 Гц), дельта (0,5-4 Гц). Анализ произведится в среде статистического программирования R версии 4.1.2 с использованием свободно распространяемой интегрированной среды для разработки (IDE –integrated development environment) на языке R – RStudio 2021.09.0+351 . Для обработки и визуализации данных использованы пакеты для языка R: easyalluvial, GGally, ggplot2, parcats, readxl, scales, tidyverse
Основной частью алгоритма является обработка таблиц, содержащих максимальные мощности ЭЭС в различных частотных диапазонах. Программа позволяет вывести данные в графическом виде, что упрощает последующий анализ данных человеком.
Денис
Чернышов
Сергеевич
-
Отсутствует
Отсутствует
-
-
РОССИЯ
Денис
Чернышов
Сергеевич
-
РОССИЯ
28.21.15
Теория сигналов
03.03.PT
MEDICAL INFORMATICS
Денис
Чернышов
Сергеевич
РОССИЯ
Александр
Тычков
Юрьевич
РОССИЯ
Валерий
Горбунов
Николаевич
РОССИЯ
Наталия
Бофанова
Сергеевна
РОССИЯ
Дмитрий
Овчинников
Леонидович
РОССИЯ
ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ УЧРЕЖДЕНИЕ "ФОНД СОДЕЙСТВИЯ РАЗВИТИЮ МАЛЫХ ФОРМ ПРЕДПРИЯТИЙ В НАУЧНО-ТЕХНИЧЕСКОЙ СФЕРЕ"
1027739479734
СЕРГЕЙ
ПОЛЯКОВ
ГЕННАДЬЕВИЧ
Организация
Федеральное государственное бюджетное учреждение "Фонд содействия развитию малых форм предприятий в научно-технической сфере"
Федеральные государственные бюджетные учреждения
Денис
Чернышов
Сергеевич
Физическое лицо
РОССИЯ