Commandline Recurrence Plots

[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.

[Norbert]

Order patterns recurrence plots are now included in the commandline version and in the plugin of the Matlab toolbox (rather fast).

[Norbert]

New location of the commandline RQA:
http://tocsy.pik-potsdam.de/commandline-rp.php

[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.

[Norbert]

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.

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

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

I don't understand your question. Sorry.