Общую форму белка можно параметризовать как последовательность точек на единичной сфере . Показаны два вида сферической гистограммы таких точек для большой коллекции белковых структур. Статистическая обработка таких данных относится к области направленной статистики. [1]
Тот факт, что 0 градусов и 360 градусов являются идентичными углами , так что, например, 180 градусов не являются разумным средним значением 2 градусов и 358 градусов, служит иллюстрацией того, что для анализа некоторых типов данных требуются специальные статистические методы (в данном случае случай, угловые данные). Другие примеры данных, которые можно рассматривать как направленные, включают статистику, включающую временные периоды (например, время суток, неделю, месяц, год и т. д.), направления по компасу, двугранные углы в молекулах, ориентации, вращения и т. д.
Эту концепцию можно распространить на многомерный контекст путем расширения простой суммы до ряда сумм, которые охватывают все измерения в пространстве признаков:
В следующих разделах показаны некоторые соответствующие циклические распределения.
Круговое распределение фон Мизеса
Распределение фон Мизеса — это круговое распределение, которое, как и любое другое круговое распределение, можно рассматривать как обертку определенного линейного распределения вероятностей вокруг круга. Основное линейное распределение вероятностей для распределения фон Мизеса математически неразрешимо; однако для статистических целей нет необходимости иметь дело с лежащим в основе линейным распределением. Полезность распределения фон Мизеса двояка: оно является наиболее математически понятным из всех круговых распределений, что позволяет упростить статистический анализ, и оно является близким приближением к завернутому нормальному распределению, которое, аналогично линейному нормальному распределению, важно, потому что это предельный случай суммы большого числа малых угловых отклонений. Фактически, распределение фон Мизеса часто называют «круговым нормальным» распределением из-за его простоты использования и его тесной связи с завернутым нормальным распределением (Fisher, 1993).
Проецируемое нормальное распределение представляет собой круговое распределение, представляющее направление случайной величины с многомерным нормальным распределением, полученное путем радиальной проекции переменной на единичную (n-1)-сферу. Из-за этого, в отличие от других часто используемых круговых распределений, оно не является ни симметричным, ни унимодальным.
Распределения на многомерных многообразиях
Три набора точек, выбранные из разных распределений Кента на сфере.
Распределение Бингама — это распределение по осям в N измерениях или, что то же самое, по точкам на ( N — 1)-мерной сфере с идентифицированными антиподами. [8] Например, если N = 2, оси представляют собой ненаправленные линии, проходящие через начало координат на плоскости. В этом случае каждая ось разрезает единичный круг на плоскости (которая является одномерной сферой) в двух точках, которые являются антиподами друг друга. Для N = 4 распределение Бингама представляет собой распределение по пространству единичных кватернионов ( версоров ). Поскольку версор соответствует матрице вращения, распределение Бингема для N = 4 можно использовать для построения распределений вероятностей в пространстве вращений, как и распределение матрицы-фон Мизеса-Фишера.
Необработанные векторные (или тригонометрические) моменты кругового распределения определяются как
где – любой интервал длины , – PDF кругового распределения и . Поскольку интеграл равен единице, а интервал интегрирования конечен, отсюда следует, что моменты любого кругового распределения всегда конечны и четко определены.
Аналогично определяются выборочные моменты:
Результирующий вектор совокупности, длина и средний угол определяются аналогично соответствующим параметрам выборки.
Кроме того, длины высших моментов определяются как:
в то время как угловые части высших моментов равны всего . Длины всех моментов будут лежать между 0 и 1.
Наиболее распространенной мерой местоположения является среднее круговое. Круговое среднее населения — это просто первый момент распределения, а среднее выборочное — это первый момент выборки. Выборочное среднее будет служить несмещенной оценкой среднего значения генеральной совокупности.
Когда данные сконцентрированы, медиану и моду можно определить по аналогии с линейным случаем, но для более рассредоточенных или мультимодальных данных эти концепции бесполезны.
Дисперсия
Наиболее распространенными мерами кругового распространения являются:
The Круговая дисперсия . Для выборки круговая дисперсия определяется как:
и для населения
Оба будут иметь значения от 0 до 1.
The круговое стандартное отклонение
со значениями от 0 до бесконечности. Это определение стандартного отклонения (а не квадратного корня из дисперсии) полезно, поскольку для завернутого нормального распределения оно является оценкой стандартного отклонения основного нормального распределения. Таким образом, это позволит стандартизировать круговое распределение, как и в линейном случае, для небольших значений стандартного отклонения. Это также относится к распределению фон Мизеса, которое близко приближается к завернутому нормальному распределению. Обратите внимание, что для small мы имеем .
The круговая дисперсия
со значениями от 0 до бесконечности. Эта мера разброса оказывается полезной при статистическом анализе дисперсии.
Распределение среднего значения
Учитывая набор из N измерений , среднее значение z определяется как:
что может быть выражено как
где
или, альтернативно, как:
где
Распределение среднего угла ( ) для кругового PDF P ( θ ) будет определяться следующим образом:
где находится в любом интервале длины, а интеграл подлежит ограничению, что и постоянны или, альтернативно, и постоянны.
Расчет распределения среднего значения для большинства круговых распределений аналитически невозможен, и для проведения дисперсионного анализа необходимы численные или математические аппроксимации. [14]
^ Бингхэм, К. (1974). «Антиподально-симметричное распределение на сфере». Анна. Стат . 2 (6): 1201–1225. дои : 10.1214/aos/1176342874 .
^ Пил, Д.; Уайтен, штат Вашингтон; Маклахлан, Дж.Дж. (2001). «Подбор смесей распределений Кента для помощи в идентификации совместного набора» (PDF) . Варенье. Стат. доц . 96 (453): 56–63. дои : 10.1198/016214501750332974. S2CID 11667311.
^ Кригер Лассен, Северная Каролина; Юул Йенсен, Д.; Конрадсен, К. (1994). «О статистическом анализе ориентировочных данных». Акта Кристаллогр . А50 (6): 741–748. дои : 10.1107/S010876739400437X.
^ Кент, Дж. Т., Хамелрик, Т. (2005). Использование распределения Фишера-Бингама в стохастических моделях структуры белка. В С. Барбере, П.Д. Бакстере, К.В.Мардиа и Р.Э. Уоллсе (ред.), «Количественная биология, анализ формы и вейвлеты», стр. 57–60. Лидс, Издательство Университета Лидса
^ Бусмма, Воутер; Мардия, Канти В.; Тейлор, Чарльз С.; Феркингхофф-Борг, Йеспер; Крог, Андерс; Хамелрик, Томас (2008). «Генеративная вероятностная модель локальной структуры белка». Труды Национальной академии наук . 105 (26): 8932–8937. Бибкод : 2008PNAS..105.8932B. дои : 10.1073/pnas.0801715105 . ПМК 2440424 . ПМИД 18579771.
^ Фишер, штат Нью-Йорк, Статистический анализ круговых данных , Cambridge University Press, 1993. ISBN 0-521-35018-2.
^ аб Джаммаламадака, С. Рао; Сенгупта, А. (2001). Темы круговой статистики. Нью-Джерси: World Scientific. ISBN978-981-02-3778-3. Проверено 15 мая 2011 г.