diff options
Diffstat (limited to 'modules/statistics/macros/nanmax.sci')
| -rwxr-xr-x | modules/statistics/macros/nanmax.sci | 52 |
1 files changed, 52 insertions, 0 deletions
diff --git a/modules/statistics/macros/nanmax.sci b/modules/statistics/macros/nanmax.sci new file mode 100755 index 000000000..f357f92aa --- /dev/null +++ b/modules/statistics/macros/nanmax.sci @@ -0,0 +1,52 @@ + +// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab +// Copyright (C) 200 - 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 [s,index] = nanmax(x,orient) + // + //This function gives for a real or a numerical matrix x a his largest + //element m (but ignoring the NANs). + // + //For x, a numerical vector or matrix, m=nanmax(x) returns in scalar m + //the largest element of x (ignoring the NANs). The form [m,index] = + //nanmax(x,orient) gives in addition of the value of the largest element + //of x (ignoring the NANs) in scalar m, and the index of this element + //in x, as a 2-vector. + // + //m=nanmax(x,'r') gives in the 1xsize(x,2) matrix m the largest elements + //(ignoring the NANs) of each column of x. If the form + //[m,index]=nanmax(x,'r') is used, the elements of the 1xsize(x,2) + //matrix index are the indexes of the largest elements (ignoring the + //NANs) of each column of x in the corresponding column. + // + //m=nanmax(x,'c') gives in the size(x,2)x1 matrix m the largest elements + //(ignoring the NANs) of each row of x. If the form + //[m,index]=nanmax(x,'c') is used, the elements of the size(x,2)x1 + //matrix index are the indexes of the largest elements (ignoring the + //NANs) of each row of x in the corresponding row. + // + // + // correction of bug 1239 (Bruno) + + [lhs, rhs] = argn(0) + if rhs == 0 then + error(msprintf(gettext("%s: Wrong number of input arguments: %d to %d expected.\n"),"nanmax",1,2)) + end + if rhs==1 then orient="*",end + if orient==1 then orient="r",end + if orient==2 then orient="c",end + if x==[]|(size(x,"*")==1&isnan(x)) then s=[],index=[],return,end + isn=isnan(x) + x(isn)=-%inf + [s,index]=max(x,orient) + s(find(and(isn,orient)))=%nan + +endfunction |