Skip to content

Commit 322b3a4

Browse files
committed
remove dependency from statistics toolbox
1 parent cb96f0a commit 322b3a4

File tree

1 file changed

+6
-5
lines changed

1 file changed

+6
-5
lines changed

rqa.m

+6-5
Original file line numberDiff line numberDiff line change
@@ -165,9 +165,9 @@
165165
l_prob = l_hist(l_hist~=0)/sum(l_hist(l_hist~=0)); % get probability distribution from histogram
166166
ent_Sum = (l_prob .* log(l_prob));
167167
if l_classes > 1
168-
y(5) = -nansum(ent_Sum)/log(N(1));
168+
y(5) = -sum(ent_Sum(~isnan(ent_Sum)))/log(N(1));
169169
else
170-
y(5) = -nansum(ent_Sum);
170+
y(5) = -sum(ent_Sum(~isnan(ent_Sum)));
171171
end
172172

173173
% histogram of vertical lines
@@ -244,15 +244,16 @@
244244
rt_prob = rt_hist(rt_hist~=0)/sum(rt_hist(rt_hist~=0)); % get probability distribution from histogram
245245
ent_Sum = (rt_prob .* log(rt_prob));
246246
if rt_classes > 1
247-
y(11) = -nansum(ent_Sum)/log(N(1));
247+
y(11) = -sum(ent_Sum(~isnan(ent_Sum)))/log(N(1));
248248
else
249-
y(11) = -nansum(ent_Sum);
249+
y(11) = -sum(ent_Sum(~isnan(ent_Sum)));
250250
end
251251

252252
if netw
253253
% clustering
254254
kv = sum(x_theiler,1); % degree of nodes
255-
y(12) = nanmean(diag(x_theiler*x_theiler*x_theiler)' ./ (kv .* (kv-1)));
255+
clustering_coeff = diag(x_theiler*x_theiler*x_theiler)' ./ (kv .* (kv-1));
256+
y(12) = mean(clustering_coeff(~isnan(clustering_coeff)));
256257

257258
% transitivity
258259
denom = sum(sum(x_theiler * x_theiler));

0 commit comments

Comments
 (0)