diff --git a/code/paradigms/ParadigmSpecCSP.m b/code/paradigms/ParadigmSpecCSP.m
index 042df8ad..f607be9f 100644
--- a/code/paradigms/ParadigmSpecCSP.m
+++ b/code/paradigms/ParadigmSpecCSP.m
@@ -155,8 +155,8 @@
                             for t = 1:size(s{c},1)
                                 s{c}(t,k) = w'*F{c}(:,:,k,t)*w; end
                         end
-                        mu_s{c} = mean(s{c});
-                        var_s{c} = var(s{c});
+                        mu_s{c} = mean(s{c},1);
+                        var_s{c} = var(s{c},0,1);
                     end
                     % 10. update alpha{j} according to Eqs. (4) and (5)
                     for c=1:2