Является ли стандартное отклонение в круговой упаковке в R правильным? - PullRequest
0 голосов
/ 26 апреля 2019

Я использую циклический пакет в R, и некоторые функции довольно просты для моих целей. mean(vector1,na.rm = TRUE) бум, готово. sd(vector1, na.rm = TRUE) бум, готово.

У меня есть 42 подшипника, которые я анализирую, и они могут быть найдены практически в любом направлении, но я надеюсь, что они имеют тенденцию особым образом. Я думал, что доверяю среднему значению, но затем, когда я запустил SD, это дало мне результат в 2,2 градуса. Мое понимание SD - то, что 2.2 слишком мал, чтобы покрыть широкий диапазон в подшипниках, которые я имею.

Есть ли что-то другое в угловых SD? Или у меня неправильно установлены параметры? Точные ли результаты я получаю?

anglecir.nest
Circular Data: 
Type = angles 
Units = degrees 
Template = geographics 
Modulo = 2pi 
Zero = 1.570796 
Rotation = clock 
 [1] 318 306 334 120  38 323 219 140 129 148  41 176 301 261 230 118 211 357 268 203 180 218 308 236 176  65  33 187  48 127  50 28 348 191 208 190 344 278  24  14 298 348

Среднее 272,937784889062 SD 2.23099490661604

Итак, если вы хотите попробовать это сами в циклическом пакете, вот код, который вы можете использовать:

anglecir.nest <- c(318, 306, 334, 120,  38, 323, 219, 140, 129, 148,  41, 
176, 301, 261, 230, 118, 211, 357, 268, 203, 180, 218, 308, 236, 176,  65,  
33, 187,  48, 127,  50, 28, 348, 191, 208, 190, 344, 278,  24,  14, 298, 348)

mean(anglecir.nest, na.rm=TRUE)
sd(anglecir.nest, na.rm=TRUE)

1 Ответ

1 голос
/ 26 апреля 2019

Вывод sd.circular вероятен в радианах, даже если данные вводятся в градусах. На что обратить внимание, если вы используете этот пакет.

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...