|
Анализ данных хронобиологических
исследований Ю.В. Корягина, С.В. Нопин Сибирский государственный университет физической культуры
и спорта Омск, Россия Презентацию
“Анализ данных хронобиологических исследований” в формате MS Power Point можно скачать здесь. Для
анализа волновых процессов и обработки хронобиологических данных принято
использовать Косинор-анализ, предложенный в Так
как любой периодический процесс можно представить в виде бесконечной суммы
периодических слагаемых (ряд Фурье), Ф. Халберг предложил аппроксимировать
методом наименьших квадратов экспериментальные данные первым членом этого
ряда – синусоидой (гармоникой). При
помощи гармонического анализа можно однозначно представить при достаточно
общих предположениях любой процесс, то есть практически любую форму волны
исследуемого биоритма. Гармоники в данном случае есть синусоиды с
целочисленными периодами по отношению к периоду объекта исследования (волны).
Число гармоник подсчитывается тривиально на основании входной информации. Аппроксимация
индивидуальных суточных кривых гармониками с заданным периодом является
первым этапом построения модели биоритмов методом Косинор-анализа. Каждый
массив наблюдений за один день называется хронограммой, расчетные кривые – синусоидами.
Входной информацией для анализа ритмов
является массив хронограмм. Каждая хронограмма должна быть построена на
основе не менее трех измерений, и чем больше число измерений, тем точнее
окончательный результат. Измерения могут быть равноотстоящие или
разноотстоящие. Требования к получению хронограмм по измерениям в одни и те
же часы суток необязательны. Выходной
информацией косинор-анализа являются основные параметры ритмов: мезор, т.е.
величина среднего уровня синусоиды (h), амплитуда синусоиды (А) и акрофаза
(Phi), то есть время наступления максимума функции. Эти данные рассчитываются
как для отдельных синусоид, так и для средней синусоиды. Кроме того,
рассчитывается и ряд дополнительных показателей: x и y (центр эллипса), r
(коэффициент корреляции между координатами точки в прямоугольной системе
координат), a (большая полуось эллипса), b (малая полуось эллипса), Sx и Sy
(ошибки x и y) , Teta (угол наклона большой оси эллипса к оси абсцисс), Delta
(среднеквадратическое отклонение синусоиды от соответствующей хронограммы). Вторым
этапом является усреднение индивидуальных данных, определение математического
ожидаемого и доверительных интервалов. Построение эллипса ошибок необходимо
для определения достоверности существования ритмов на принятом доверительном
уровне (например на уровне 0,95). Синусоида изображается па плоскости точкой,
полярные координаты которой — амплитуда и акрофаза. Все полученные таким
образом точки в декартовых координатах рассматриваются как реализации
двумерной случайной величины с гипотетически нормальным законом
распределения, и строится эллипс рассеивания ошибок генерального среднего. Реализацию
Косинор-анализа на персональном компьютере можно осуществить с помощью
специальных программ. Расчетные данные Косинор-анализа получаем с помощью
компьютерной программы “Cosinor-Analisis 2.4 for Excel 2000/XP”,
разработанной С.Н. Шереметьевым (2005). Работа с программой Cosinor-Analisis 2.4
for Excel 2000/XP По
умолчанию программа инсталлируется в директорию “Documents and Settings\User Application
Data\Microsoft\AddIns”. После инсталляции
программы Вам необходимо открыть Microsoft Excel 2000/XP/2003. В меню “Сервис
> Надстройки” пометьте флажком пункт “Косинор-анализ”. Программа
автоматически создаст панель инструментов “ExcelStat” и кнопку “Cosinor”
(знак песочных часов).
Рис 1. Окно программы MS Excel c меню надстройки Нажатие кнопки “Cosinor” запускает программу Cosinor-Analisis 2.4 for Excel 2000/XP. В появившемся окне программы необходимо задать анализируемое время, при анализе данных суточных ритмов обычно берется 15 измерений (3 дня по 5 измерений). Величина пробелов между измерениями задается в фактических часах измерения. Например, измеряем 7, 11, 15, 19 и 23 ч. В следующие сутки измерения начинаются опять с 7 ч, но в столбце мы указываем 31 (7+24), чтобы не нарушать временного ряда и продолжаем 35, 39, 43, 47 и т.д... Полученный временной ряд и есть анализируемое время, т.е. входной диапазон времени. В настройках программы также необходимо указать массив хронограмм и диапазон для вывода результатов, периоды ритмов для анализа. Функция расчет хронограмм позволяет рассчитать синусоиды как для каждого испытуемого, так и средние по группе. Нажатие кнопки “Считать” запускает косинор-анализ.
Рис 2. Окно программы Cosinor-Analisis 2.4 for Excel 2000/XP с установкой параметров косинор-анализа В результате указанных действий получаются данные косинор-анализа. Например, расчет параметров 24 часовых ритмов индивидуальной минуты (табл.) и времени простой сенсомоторной реакции на звук (табл.). Таблица 1 Фрагмент данных косинор-анализа 24 часовой гармоники индивидуальной минуты
Примечание: x и y - центр эллипса, h - мезор, А – амплитуда, Phi – акрофаза, Delta - среднеквадратическое отклонение
синусоиды от соответствующей хронограммы, Sx и Sy - ошибки x и y,
r - коэффициент
корреляции между координатами точки в прямоугольной системе координат, a и b - большая и малая полуоси эллипса, Teta - угол наклона большой оси эллипса к оси
абсцисс. Таблица 2 Фрагмент данных косинор-анализа 24 часовой гармоники времени простой сенсомоторной реакции на звук
По данным расчета индивидуальных и групповых хронограмм с помощью средств MS Excel можно построить синусоиду или аппроксимированный ритм времени реакции.
х, ч Рис 3. Аппроксимированный 24 ч суточный ритм (синусоида) времени простой сенсомоторной реакции на звук Дальнейшие действия косинор-анализа заключаются в выяснении достоверности существования полученных ритмов (гармоник). Для этого строится доверительный интервал или эллипс ошибок (эллипс рассеивания). Графическое представление данных Косинор-анализа
с построением доверительных интервалов (эллипсов рассеивания) можно
осуществить с помощью компьютерной программы “Cosinor Ellipse Косинор
анализ с использованием программы Cosinor Ellipse 2006 Программа предназначена для операционных систем Windows 95/98/2000/XP и NT. Свидетельство об официальной регистрации программы представлено на рисунке. Для установки программы достаточно запустить файл установки “setup” и зарегистрироваться. Ярлык программы появится в меню ПУСК – Программы – Cosinor Ellipse 2006. Файлы с результатами автоматически сохраняются в директорию программы в папке Program Files.
Рис 4. Свидетельство об официальной регистрации программы “Cosinor
Ellipse Внешний вид программы “Cosinor Ellipse
Рис 5. Внешний вид программы Cosinor Ellipse 2006 с меню “Файл” Меню “Правка” позволяет выполнить необходимые настройки: задать диапазон поиска (в зависимости от количества анализируемых таблиц, масштаб (в зависимости от единиц измерения, чем больше единицы, тем меньше масштаб), установить необходимость анализировать все таблицы одновременно или по одной, средние по группе или данные каждого испытуемого и т.д.
Рис 6. Внешний вид программы Cosinor Ellipse 2006 с меню “Правка” Следующим этапом выбирается лист в открытом файле MS Excel на котором находятся необходимые данные, например листы с названиями “ИМ” и “реакция звук”. Программа получает из открытого файла данные на основании которых будет проведено графического представление результатов косинор-анализа хронограмм. Функция “Номер группы (таблицы)” позволяет пользователю выбирать номер таблицы результатов косинор-анализа хронограмм из таблиц сохраненных в оперативной памяти ПЭВМ.
Рис 7. Окно программы Cosinor Ellipse 2006 с меню выбора листа в файле MS Excel Указанные действия приводят к построению синусоид аппроксимирующих хронограммы, каждая из которых изображается на плоскости точкой, полярные координаты которой - амплитуда и акрофаза. Аппроксимирующая хронограммы усредненная синусоида изображается крестиком. Строится эллипс рассеивания ошибок генерального среднего на основании следующих величин: средние x и y (центр эллипса), полуоси эллипса a и b, угол наклона эллипса к основной системе координат. При анализе достоверности ритмов учитываем, что усредненная аппроксимирующая хронограммы синусоида (изображенная крестиком) должна входить в эллипс, а сам эллипс не должен проходить через центр координат (т.к. в данном случае акрофаза будет приходиться на весь 24 ч период). При соблюдении данных условий рассчитанный нами ритм является достоверным, в противном случае данный ритм не существует. Рассмотрим анализируемые в наших примерах данные. В первом случае рассчитан и изображен на рисунке эллипс ошибок 24 часовой гармоники индивидуальной минуты: средняя синусоида расположена в центре эллипса, но эллипс проходит через центр координат (разброс акрофаз от 0 до 24 ч), такая гармоника не является достоверной, т.е. 24 часового ритма индивидуальной минуты у исследуемой группы лиц не наблюдается.
Рис 8. Окно программы Cosinor Ellipse 2006 с результатами анализа 24 часовой гармоники индивидуальной минуты Во втором случае – строим эллипс рассеивания для 24 часовой гармоники времени простой сенсомоторной реакции на звук. Эллипс 24 часовой гармоники соблюдает оба условия, доверительный интервал акрофазы с 12,96 часов до 19,10 часов, таким образом данный суточный ритм является достоверным.
Рис 9. Окно программы Cosinor Ellipse 2006 с результатами анализа 24 часовой гармоники времени простой сенсомоторной реакции на звук (акрофазы приходятся на период с 12,96 часов до 19,10 часов) Помимо
графического изображения эллипса программа выводит следующие расчетные данные,
необходимые для анализа биоритмов: Amin и Amax (минимальная и
максимальная амплитуда), PhiMin и PhiMax (минимальная и
максимальная акрофаза – интервал акрофазы), MinH и MaxH (минимальный и максимальный мезор – разброс среднего значения), CkoH/n (ошибка среднего) (рис.). Таким образом, представленные программы позволяют проанализировать хронобиологические данные, выяснить наличие ритмичности различных функциональных показателей, а также рассчитать параметры установленных ритмов. Презентацию
“Анализ данных хронобиологических исследований” в формате MS Power Point можно скачать здесь. |
