Hi all,
I have a circular data set, in degrees, with a known standard deviation (SD) of between 60-80°— the range in SD is because there are different formulas.
Anyways, when I use the sd.circular()
function, the answer is clearly wrong. Doesn't even pass the eye test.
What am I doing wrong here?
library(tidyverse)
library(circular)
#>
#> Attaching package: 'circular'
#> The following objects are masked from 'package:stats':
#>
#> sd, var
SD_CircTest <- tibble(c(228,64,259,124,150,152,309,335,256,242,233,334,34,
198,229,225,105,203,260,134,231,143,165,166,156,227,
317,154,102,208,191,15,144,329,184,321,230,276,155,
73,289,204,201,5,246,179,301,208,202,55,223,159,112,
276,85,137,226,174,158,155))
SD__CircTest <- circular(SD_CircTest, units = "degrees", template = "geographics")
sd.circular(SD__CircTest)
#> c(...)
#> 2.257284