comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Circular statistics in IDL
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: Circular statistics in IDL [message #94035 is a reply to message #94019] Thu, 29 December 2016 17:16 Go to previous message
Dick Jackson is currently offline  Dick Jackson
Messages: 347
Registered: August 1998
Senior Member
On Wednesday, 21 December 2016 09:25:04 UTC-8, Craig Markwardt wrote:
> On Tuesday, December 20, 2016 at 11:05:03 AM UTC-5, Gadi Goelman wrote:

>> What I need is to calculate significance of complex numbers or phases. Do you know if correct calculations of mean, variance etc is sufficient? I guess that you cant assume normal distributions and different distributions are needed such circular Gaussian.
>
> I'm a big fan of the Rayleigh (Z) statistic. It's pretty easy to calculate. For an array of unit complex number U, it's
>
> Z = total(abs(U)^2)
>
> In the case of the null hypothesis (only noise), the Z statistic is distributed as a chi-square with 2 degrees of freedom. You you can use IDL's chi-square statistic to test statistical significance, or my equivalent MPCHITEST() function.
>
> Craig

In general, I trust anything Craig says :-), and other sources point in the same direction:
https://en.wikipedia.org/wiki/Directional_statistics#Goodnes s_of_fit_and_significance_testing

Some MATLAB functions written for this seem to follow this:
https://github.com/circstat/circstat-matlab/blob/master/circ _r.m
https://github.com/circstat/circstat-matlab/blob/master/circ _rtest.m

I hope this is helpful!

Cheers,
-Dick

Dick Jackson Software Consulting Inc.
Victoria, BC, Canada --- http://www.d-jackson.com
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Attempt to call undefined procedure/function -- error
Next Topic: help needed in timegen

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Wed Oct 08 11:00:19 PDT 2025

Total time taken to generate the page: 0.00442 seconds