CMSCOPE changed

1 files changed, 53 insertions, 0 deletions

diff --git a/modules/statistics/macros/nanmedian.sci b/modules/statistics/macros/nanmedian.sci
new file mode 100755
index 000000000..95d108cb3
--- /dev/null
+++ b/modules/statistics/macros/nanmedian.sci
@@ -0,0 +1,53 @@
+
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) 2000 - INRIA - Carlos Klimann
+//
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at
+// http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt
+//
+
+
+function [m]=nanmedian(x,orient)
+    //
+    //For a vector or a matrix x, [m]=median(x) returns in  the vector p the
+    //median of the values (ignoring the  NANs) of vector x.
+    //
+    //[m]=median(x,'r')  (or, equivalently, [m]=pctl(x,1))  are the  rowwise
+    //medians. It returns in each  position of the row  vector m the medians
+    //of data (ignoring the NANs) in the corresponding column of x.
+    //
+    //[m]=median(x,'c')   (or,  equivalently,   [m]=median(x,2))    are  the
+    //columnwise   percentiles.  It returns in  each  position of the column
+    //vector m the medians of data (ignoring the  NANs) in the corresponding
+    //row of x.
+    //
+    //
+    [lhs,rhs]=argn(0)
+    if rhs<1|rhs>2 then  error(msprintf(gettext("%s: Wrong number of input argument: %d to %d expected.\n"),"nanmedian",1,2)), end
+    if type(x)<>1 then error(msprintf(gettext("%s: Wrong type for input arguments #%d: Numerical expected.\n"),"nanmedian",1)) , end
+    if x==[] then m=[], return,end
+    if rhs==1 then
+        p=perctl(x(~isnan(x)),50)
+        if p==[] then p=%nan,end
+        m=p(1)
+    elseif orient=="r"|orient==1 then
+        m=[]
+        for i=x
+            p=perctl(i(~isnan(i)),50)
+            if p==[] then p=%nan,end
+            m=[m p(1)]
+        end
+    elseif orient=="c"|orient==2 then
+        m=[]
+        for i=x' do
+            p=perctl(i(~isnan(i)),50)
+            if p==[] then p=%nan,end
+            m=[m;p(1)];
+        end
+    else
+        error(msprintf(gettext("%s: Wrong value for input argument #%d: ''%s'', ''%s'', %d or %d expected.\n"),"nanmedian",3,"r","c",1,2))
+    end
+endfunction