Computing Laminarity

Everything about quantification of recurrence plots and recurrence networks.
Post Reply
holistic
Senior
Posts: 16
Joined: Wed Feb 19, 2014 09:53
Affiliation (Univ., Inst., Dept.): Univ. Bochum
Location: Bochum, Germany
Research field: Neuroscience/Psychology

Computing Laminarity

Post by holistic »

Hi,

I have 10 timeseries and I want to compute the Lamniarity for all of those without using the gui-interface. Unfortunately I don't know how the commands are to save to compute the crp and the LAM and save it into an array. Can anybody help me out with which matlab code to use?
User avatar
Norbert
Expert
Posts: 183
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: Computing Laminarity

Post by Norbert »

There are different ways, depending whether the time series are in files or are already in your Matlab workspace. Let's start simple and assume your time series are the columns of an array:

x = [x1, x2, x3, … x10];
size(x) gives something like N, 10, where N is the length of your time series.

Now you can use

Code: Select all

for i = 1:10
   Rtemp = crqa(x(:,i),m,t,e,'silent');
   LAM(i) = Rtemp(6);
end
(where m and t are the embedding parameters and e the recurrence threshold which you have to chose)

In case you need the time resolved values, i.e., applying moving windows, simply use

Code: Select all

for i = 1:10
   Rtemp = crqa(x(:,i),m,t,e,w,ws,'silent');
   LAM(:,i) = Rtemp(:,6);
end
(where w and ws are window size and window step)
Post Reply