## Computing Laminarity

Everything about quantification of recurrence plots and recurrence networks.
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

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?
Norbert
Expert
Posts: 195
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

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)