pss function and other methods to find proper radius for CRQA
Posted: Fri Oct 1, 2021 14:10
Dear all community members,
I'm a clinical psychologist and a beginner in the usage of these measures so I kindly ask for your support. Going through the literature I found different methods to find the proper radius for computing cross recurrence measures. So I tried to compute some of them.
- Firstly, after having found the appropriate time (mi function) and dimension (fnn function), I tried to compute cross recurrence
measures on the whole series, fixing the recurrence rate at 0.02 and 0.05, in the following way:
crqa(timeseries1, timeseries2,m,t,0.02,[],'rr','normalize')
crqa(timeseries1, timeseries2,m,t,0.05,[],'rr','normalize')
- Secondly, I was trying to input as radius the 5% of the maximal phase space diameter computed with pss function. But i'm not sure to have
properly used the function. Following the line of code I used:
mpsd = 0.05*pss(timeseries1,timeseries2,m,t,'euclidean').
- Lastly I tried to find the radius computing iteratively the recurrence rate percentages for different radii searching for a linear scaling
region in the loglog plot (attached an example) of radii (on x-axis) and recurrence rates % (y-axis), but i'm not able to recognize such
linear scaling region in the plot.
Any hint or help would be of great support for me.
Thank you very much
Best regards
I'm a clinical psychologist and a beginner in the usage of these measures so I kindly ask for your support. Going through the literature I found different methods to find the proper radius for computing cross recurrence measures. So I tried to compute some of them.
- Firstly, after having found the appropriate time (mi function) and dimension (fnn function), I tried to compute cross recurrence
measures on the whole series, fixing the recurrence rate at 0.02 and 0.05, in the following way:
crqa(timeseries1, timeseries2,m,t,0.02,[],'rr','normalize')
crqa(timeseries1, timeseries2,m,t,0.05,[],'rr','normalize')
- Secondly, I was trying to input as radius the 5% of the maximal phase space diameter computed with pss function. But i'm not sure to have
properly used the function. Following the line of code I used:
mpsd = 0.05*pss(timeseries1,timeseries2,m,t,'euclidean').
- Lastly I tried to find the radius computing iteratively the recurrence rate percentages for different radii searching for a linear scaling
region in the loglog plot (attached an example) of radii (on x-axis) and recurrence rates % (y-axis), but i'm not able to recognize such
linear scaling region in the plot.
Any hint or help would be of great support for me.
Thank you very much
Best regards