From 322b3a4f40233b26425e2e42ac1b9923d062528d Mon Sep 17 00:00:00 2001
From: Norbert Marwan <norbert.marwan@gmail.com>
Date: Mon, 23 Sep 2024 18:23:10 +0200
Subject: [PATCH] remove dependency from statistics toolbox

---
 rqa.m | 11 ++++++-----
 1 file changed, 6 insertions(+), 5 deletions(-)

diff --git a/rqa.m b/rqa.m
index 12e26ae..bdffa71 100644
--- a/rqa.m
+++ b/rqa.m
@@ -165,9 +165,9 @@
 l_prob = l_hist(l_hist~=0)/sum(l_hist(l_hist~=0)); % get probability distribution from histogram
 ent_Sum = (l_prob .* log(l_prob));
 if l_classes > 1
-    y(5) = -nansum(ent_Sum)/log(N(1));
+    y(5) = -sum(ent_Sum(~isnan(ent_Sum)))/log(N(1));
 else
-    y(5) = -nansum(ent_Sum);
+    y(5) = -sum(ent_Sum(~isnan(ent_Sum)));
 end
 
 % histogram of vertical lines
@@ -244,15 +244,16 @@
 rt_prob = rt_hist(rt_hist~=0)/sum(rt_hist(rt_hist~=0)); % get probability distribution from histogram
 ent_Sum = (rt_prob .* log(rt_prob));
 if rt_classes > 1
-    y(11) = -nansum(ent_Sum)/log(N(1));
+    y(11) = -sum(ent_Sum(~isnan(ent_Sum)))/log(N(1));
 else
-    y(11) = -nansum(ent_Sum);
+    y(11) = -sum(ent_Sum(~isnan(ent_Sum)));
 end
 
 if netw
     % clustering
     kv = sum(x_theiler,1); % degree of nodes
-    y(12) = nanmean(diag(x_theiler*x_theiler*x_theiler)' ./ (kv .* (kv-1)));
+    clustering_coeff = diag(x_theiler*x_theiler*x_theiler)' ./ (kv .* (kv-1));
+    y(12) = mean(clustering_coeff(~isnan(clustering_coeff)));
 
     % transitivity
     denom = sum(sum(x_theiler * x_theiler));