JuliaDynamics / RecurrenceAnalysis.jl

Recurrence Quantification Analysis in Julia
https://juliadynamics.github.io/RecurrenceAnalysis.jl/stable/
MIT License
44 stars 12 forks source link

Hypothesis testing with Recurrence Matrices #154

Closed fargolo closed 1 year ago

fargolo commented 1 year ago

Hello,

I intend to calculate the probabilities associated with structures in Rec. Matrices (trapping time, vertical, diagonal).

I sketched a way of calculating the probabilities associated with the sequences as binary words. This allows one to calculate, for instance, the probability of a sequence (e.g. vertical run, 'trapping time') considering its section (column size with all entries a_i,j).

How would you include the information of the lagged sections (other columns)? Could you point me to any references?

Datseris commented 1 year ago

Hi, the question falls outside my expertise but maybe @pucicu has some idea.

pucicu commented 1 year ago

Yes, I was already thinking about it, but to be honest, I don’t completely understand the question. I also do not get the point how the provided link is related to the question.

Am 13.09.2023 um 16:20 schrieb George Datseris @.***>:

Hi, the question falls outside my expertise but maybe @pucicu https://github.com/pucicu has some idea.

— Reply to this email directly, view it on GitHub https://github.com/JuliaDynamics/RecurrenceAnalysis.jl/issues/154#issuecomment-1717730862, or unsubscribe https://github.com/notifications/unsubscribe-auth/AAHNNEG4HCYD5U5GTCRB2YTX2G6LVANCNFSM6AAAAAA4RPJTL4. You are receiving this because you were mentioned.

fargolo commented 1 year ago

Hello, @pucicu .

I was just reading your paper and the problem is related to section 2.3

"A fundamental finding was elaborated by [67], who mathematically developed the connection between correlation sum C(m) and the RQA measures recurrence rate RR(m), determinism DET(m), and average diagonal line length L(m). (...) These considerations have been further elaborated by [68, 69], which have further derived the analytical expressions for several RQA measures for certain stochastic processes, fractional Gaussian noise, and correlated noise (first analytical solutions were already given in [70]). Analytical and asymptotical expressions for RQA measures of specific random processes are important for defining baselines for benchmarking and testing."

I did not have the chance to check these references yet, The method that I linked calculates probabilities for consecutive runs in random binary sequences.

One could calculate the probability of observing Lmax or Lavg if the column was random.

Datseris commented 1 year ago

(I am closing this issue because it is not related with the code base RecurrenceAnalysis.jl, but rather it is a scientific question; please do feel free to continue the discussion as much as you would like here!)