bugfix in estimation black and white vertical lines

crqa.m

@@ -149,6 +149,10 @@ function xout=crqa(varargin)
 % $Revision$
 %
 % $Log$
+% Revision 5.46  2013/08/21 09:09:59  marwan
+% bugfix for mixed output RTE and T2
+% bugfix for recurrence time estimation
+%
 % Revision 5.45  2012/11/14 13:07:03  marwan
 % bugfix for RTE for very short time series
 %
@@ -1268,7 +1272,7 @@ for i=1:wstep:Nx-w;
          warning off 
          errcode=271;
          b(b<lmin)=[];
-         [c d dw]=tt(X);
+         [c d dw]=tt(X_theiler);
          hist_v = [hist_v; d];
          hist_w = [hist_w; dw];         
          N_hist_v(i) = length(d);

tt.m

@@ -27,6 +27,9 @@ function [a_out, b_out, c_out]=tt(x)
 % $Revision$
 %
 % $Log$
+% Revision 3.7  2010/06/30 12:02:52  marwan
+% Help text modified
+%
 % Revision 3.6  2009/03/24 08:33:47  marwan
 % copyright address changed
 %
@@ -85,6 +88,12 @@ try
 
         %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
         % for white vertical lines
+        
+        % if LOI is not present then put it in!
+        if ~all(diag(x))
+           x = x + eye(size(x));
+        end
+        
         for i = 1:size(x,2)
             i_ = find(x(:,i));
             i1 = min(i_); i2 = max(i_);