|WikiProject Statistics||(Rated C-class, High-importance)|
Circular standard deviation | Factor of 2
I noticed an oddity in the definition of the circular variance, which probably has historical reasons: While it is not generally true that the estimator of the circular standard deviation squared yields the circular standard deviation, I would have expected this equality to hold for very small values of the circular standard deviation, which recovers the linear case. However, there is a factor of 2 missing! You can test this numerically yourself with the following line of code
I therefore added the following sentence to the article:
- Note that for small , we have .
Does anybody know a reference for this?
The final step of the formula given for the example is unclear; apparently it is written for result in radians. Wouldn't it be better written this way:
Mean angle = arctangent (mean sine / mean cosine). If mean cosine < 0, add 180 degrees to the result. If mean sine <0 and mean cosine >= 0, add 360 degrees to the result.
The modulus method only works in few specific cases (such as the example given).
Consider the same example but rotated further to the left so that the three angles are 330, 340 and 350 degrees. Taking the modulus 360 of the sum (1020) results in 2 remainder 280, which divided by 3 is clearly not 340.
As the remainder (in a modulus 360 operation) can only range from 0 to 360, the 'average' can only range from 0 to 360/n.
Unless I'm misunderstanding the procedure, I would suggest taking it out.
Fink3412 09:28, 20 August 2007 (UTC)
Agreed - I've removed it. Also, this article is about Directional statistics itself, not about various methods to calculate a circular average. Tomixdf 09:42, 20 August 2007 (UTC)
I've added the wrapped normal distribution as an example of the way a wrapped distribution can be made from the pdf of another distribution. I'm also redirecting the nonexistent wrapped normal distribution to this page. I hope this is considered useful! digfarenough (talk) 22:18, 11 November 2008 (UTC)