Dear all,
For a noisy time series, it is mentioned in the literature that the threshold should be around 5*sigma. I wanted to ask if this sigma is the standard deviation of the time series?
recurrence threshold
-
- Junior
- Posts: 5
- Joined: Sat Feb 27, 2021 10:30
- Affiliation (Univ., Inst., Dept.): IIT Ropar
- Location: Ropar, Punjab, India
- Research field: Engine combustion
- Norbert
- Expert
- Posts: 194
- Joined: Wed Jan 4, 2006 11:03
- Affiliation (Univ., Inst., Dept.): Potsdam Institute for Climate Impact Research, Germany
- Location: Potsdam, Germany
- Location: Potsdam Institute for Climate Impact Research, Germany
Re: recurrence threshold
No, it is the standard deviation of the observation noise. This threshold value (5σ) is actually much too high. The study which you refer to is more a theoretical consideration. This choice for the threshold means that the structures in the RP will be preserved. But it will also include additional structures that are related to noise. These new noise-related structures will strongly bias the final results.ajay singh wrote: ↑Tue Mar 2, 2021 14:07 Dear all,
For a noisy time series, it is mentioned in the literature that the threshold should be around 5*sigma. I wanted to ask if this sigma is the standard deviation of the time series?
Optimal theshold selection depends on the research question. This was summarized in several papers, e.g., in http://dx.doi.org/10.1016/j.physrep.2006.11.001 or http://dx.doi.org/10.1142/S0218127411029008.
-
- Junior
- Posts: 3
- Joined: Wed Mar 10, 2021 04:21
- Affiliation (Univ., Inst., Dept.): zhengzhou university Henan China
- Location: china Henan zhengzhou
- Research field: Bearing fault diagnosis
Re: recurrence threshold
Dir sir:Norbert wrote: ↑Wed Mar 3, 2021 15:24No, it is the standard deviation of the observation noise. This threshold value (5σ) is actually much too high. The study which you refer to is more a theoretical consideration. This choice for the threshold means that the structures in the RP will be preserved. But it will also include additional structures that are related to noise. These new noise-related structures will strongly bias the final results.ajay singh wrote: ↑Tue Mar 2, 2021 14:07 Dear all,
For a noisy time series, it is mentioned in the literature that the threshold should be around 5*sigma. I wanted to ask if this sigma is the standard deviation of the time series?
Optimal theshold selection depends on the research question. This was summarized in several papers, e.g., in http://dx.doi.org/10.1016/j.physrep.2006.11.001 or http://dx.doi.org/10.1142/S0218127411029008.
I use CRP toolbox ,i have some puzzle about the function ‘crp(x,m,t,e)’ and ‘pss’; e represent the threshold,but when i choose the norm
‘euclidean’,and set the e =2,like 'crp(x,3,1,2,'euclidean')',I noticed that there is a title on the picture ‘Dimention:3 Delay:1 threshold: 2σ’ ,i dont konw what ‘2σ’ mans,what is the σ?? Is σ equal to 1??(x is normalized when use crp??) Or σ is the origin sequence's Standard deviation before normalize??‘2’ means.
i know the 'pss' function can calculate the max phase space radius , and set it to 0.1*pss(x,m,t,'euclidean'),but the result can be imported directly into crp (Do I need use (0.1*maxdist)/std(x) )????
tks sir , i am really need some help
And if the time sequemce have pulse,threshold calculate by 10% phase radiu is big ,it make the rp almost black what should i do??
Last edited by Yuzhi233 on Fri Mar 12, 2021 14:01, edited 1 time in total.
-
- Junior
- Posts: 5
- Joined: Sat Feb 27, 2021 10:30
- Affiliation (Univ., Inst., Dept.): IIT Ropar
- Location: Ropar, Punjab, India
- Research field: Engine combustion
Re: recurrence threshold
Thanks sir for replying
- Norbert
- Expert
- Posts: 194
- Joined: Wed Jan 4, 2006 11:03
- Affiliation (Univ., Inst., Dept.): Potsdam Institute for Climate Impact Research, Germany
- Location: Potsdam, Germany
- Location: Potsdam Institute for Climate Impact Research, Germany
Re: recurrence threshold
So many questions …Yuzhi233 wrote: ↑Thu Mar 11, 2021 23:28
Dir sir:
I use CRP toolbox ,i have some puzzle about the function ‘crp(x,m,t,e)’ and ‘pss’; e represent the threshold,but when i choose the norm
‘euclidean’,and set the e =2,like 'crp(x,3,1,2,'euclidean')',I noticed that there is a title on the picture ‘Dimention:3 Delay:1 threshold: 2σ’ ,i dont konw what ‘2σ’ mans,what is the σ?? Is σ equal to 1??(x is normalized when use crp??) Or σ is the origin sequence's Standard deviation before normalize??‘2’ means.
i know the 'pss' function can calculate the max phase space radius , and set it to 0.1*pss(x,m,t,'euclidean'),but the result can be imported directly into crp (Do I need use (0.1*maxdist)/std(x) )????
tks sir , i am really need some help
And if the time sequemce have pulse,threshold calculate by 10% phase radiu is big ,it make the rp almost black what should i do??
- The default is that the time series will be automatically normalized to have standard deviation σ=1. Therefore, the threshold is in units of σ.
- If you disable the normalisation by the argument 'nonorm', you will see that this σ will vanish.
- I do not understand clearly your question on 'pss'. This function also normalises the time series to have σ=1. If you prefer to work with non-normalised time series, simply use the 'nonorm' flag in these functions.
-
- Junior
- Posts: 3
- Joined: Wed Mar 10, 2021 04:21
- Affiliation (Univ., Inst., Dept.): zhengzhou university Henan China
- Location: china Henan zhengzhou
- Research field: Bearing fault diagnosis
Re: recurrence threshold
thank U sir .I finally figured it out, thank you very muchNorbert wrote: ↑Wed May 5, 2021 16:15So many questions …Yuzhi233 wrote: ↑Thu Mar 11, 2021 23:28
Dir sir:
I use CRP toolbox ,i have some puzzle about the function ‘crp(x,m,t,e)’ and ‘pss’; e represent the threshold,but when i choose the norm
‘euclidean’,and set the e =2,like 'crp(x,3,1,2,'euclidean')',I noticed that there is a title on the picture ‘Dimention:3 Delay:1 threshold: 2σ’ ,i dont konw what ‘2σ’ mans,what is the σ?? Is σ equal to 1??(x is normalized when use crp??) Or σ is the origin sequence's Standard deviation before normalize??‘2’ means.
i know the 'pss' function can calculate the max phase space radius , and set it to 0.1*pss(x,m,t,'euclidean'),but the result can be imported directly into crp (Do I need use (0.1*maxdist)/std(x) )????
tks sir , i am really need some help
And if the time sequemce have pulse,threshold calculate by 10% phase radiu is big ,it make the rp almost black what should i do??
- The default is that the time series will be automatically normalized to have standard deviation σ=1. Therefore, the threshold is in units of σ.
- If you disable the normalisation by the argument 'nonorm', you will see that this σ will vanish.
- I do not understand clearly your question on 'pss'. This function also normalises the time series to have σ=1. If you prefer to work with non-normalised time series, simply use the 'nonorm' flag in these functions.