Commandline Recurrence Plots

Available software for computation of recurrence plots and their quantification. Announcements of new releases.
Post Reply
User avatar
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

Commandline Recurrence Plots

Post by Norbert »

The bunch of free available software for computing RPs and their quantifications was extended by a commandline version, which is fast and can handle long data series. It includes the most recent developed RQA measures and can be used in batch processes.

It can be found at http://www.agnld.uni-potsdam.de/~marwan ... oad/rp.php.
User avatar
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

Order pattern recurrence plots

Post by Norbert »

Order patterns recurrence plots are now included in the commandline version and in the plugin of the Matlab toolbox (rather fast).
User avatar
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: Commandline Recurrence Plots

Post by Norbert »

New location of the commandline RQA:
http://tocsy.pik-potsdam.de/commandline-rp.php
viniciusjdv
Junior
Posts: 4
Joined: Sun Feb 10, 2013 21:29
Affiliation (Univ., Inst., Dept.): Federal University of Campina Grande
Location: Campina Grande, Brazil
Research field: Engeneering

Re: Commandline Recurrence Plots

Post by viniciusjdv »

Professor Marwan,

I have some doubts about the use of command line reccurrence plots:

Let's see the use of rp like the official page:

dos('rp -i temp_data -s -m 3 -t 5 -e 0.7 -o rqa.dat');

Doubts:

1) Is this threshold parameter an absolute value of the radius or a percentage value of the maximum diameter of the phase space?

2) Is there a way to vary the parameter m and t? I mean, instead of using "-m 3 -t 7", Could I to use something like that (adapted from the official page):

Code: Select all

% import some data
x=load('soi3.dat'); 

% data length
N = length(x); 
% window length and window step size
w = 100; s = 10; 

for i = 1:s:N-w
    m=i; tau=round(N/i);
    % make sub-string of data
    sub_x = x(i:i+w-1); 
    % save sub-string of data
    save('temp_data','-ascii','sub_x') 
    % call extern rp programme
    dos('./rp -i temp_data -s -m m(i) -t tau(i) -e 0.7 -o rqa.dat'); 
end
delete('temp_data')
Instead of use numbers, could I to use m(i) and tau(i)? I have tested like this and it did not work.

3) The similar use of command line rp in CRP toolbox is like that: (?)

Code: Select all

X=temp_data;
l_min=2; v_min=2;
crqa(X,3,5, 0.7, [], [], l_min, v_min, 1, ...
              'euclidean', 'nonormalize', 'silent');
Thanks in advance!

Best regards,

Vinicius Veira.
Vinícius Jefferson Dias Vieira, M.Sc.
Electrical Engineering Department,
Federal University of Campina Grande, Brazil
User avatar
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: Commandline Recurrence Plots

Post by Norbert »

viniciusjdv wrote:1) Is this threshold parameter an absolute value of the radius or a percentage value of the maximum diameter of the phase space?
Yes, it is an absolute value.
viniciusjdv wrote:2) Is there a way to vary the parameter m and t? I mean, instead of using "-m 3 -t 7", …
Yes, you can do it. You have to use either num2str or sprintf:

Code: Select all

m = 1:5; tau = 1:10;
for i = 1:length(m)
   for j = 1:length(tau)
       % call extern rp programme
       dos(['./rp -i temp_data -s -m ',num2str(m(i)),' -t ',num2str(tau(j)),' -e 0.7 -o rqa.dat']); 
   end
end
viniciusjdv wrote:3) The similar use of command line rp in CRP toolbox is like that: (?)
I don't understand your question. Sorry.
Post Reply