-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab\r
-// Copyright (C) INRIA\r
-// Copyright (C) Bruno Pincon\r
-// Copyright (C) 2010 - Samuel Gougeon\r
-// Copyright (C) 2012 - Scilab Enterprises - Adeline CARNIS\r
-//\r
-// This file must be used under the terms of the CeCILL.\r
-// This source file is licensed as described in the file COPYING, which\r
-// you should have received as part of this distribution. The terms\r
-// are also available at\r
-// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt\r
-\r
-\r
-function histplot(n,data,style,strf,leg,rect,nax,logflag,frameflag,axesflag,normalization)\r
- // histplot(n,data,<opt_arg_seq>)\r
- // draws histogram of entries in data put into n classes\r
- //\r
- // histplot(xi,data,<opt_arg_seq>)\r
- // generates the histogram of entries in data put into classes\r
- // [xi(1),xi(2)], (xi(k) xi(k+1)], k=2,..,n.\r
- // xi's are assumed st. increasing (this point is verified now).\r
- //\r
- // optionnal args:\r
- // 1/ the same than for a plot2d:\r
- // style,strf,leg,rect,nax,logflag,frameflag,axesflag\r
- // 2/ normalization flag (default value %t). When true the\r
- // histogram is normalized so that to approach a density:\r
- // xmax\r
- // /\r
- // | h(x) dx = 1 (true if xmin <= min(data) and max(data) <= xmax)\r
- // /\r
- // xmin\r
- //\r
- // Example : enter histplot()\r
- //\r
- // modifs to use dsearch (Bruno Pincon 10/12/2001)\r
- // others modifs from Bruno (feb 2005):\r
- // - may be a cleaner and shorter way to deal with optionnal arg ?\r
- // - now the histgram is drawn using only one polyline\r
- // (so properties (color, thickness,...) are easier\r
- // to change with new graphics).\r
- // - removed computation of nax and rect if they are not\r
- // passed (let plot2d doing it)\r
- // - modify a little the demo\r
- // - add some checking on n|x and data\r
- //\r
- [lhs,rhs]=argn()\r
-\r
- if rhs == 0 then // demo\r
- histplot([-4.5:0.25:4.5],rand(1,20000,'n'),style=2,axesflag=1,..\r
- frameflag=1,rect=[-4.5 0 4.5 0.47]);\r
- deff("[y]=f(x)","y=exp(-x.*x/2)/sqrt(2*%pi);");\r
- x=-4.5:0.125:4.5;\r
- x=x';\r
- plot2d(x,f(x),26,"000");\r
- titre= gettext("histplot() : (normalized) histogram plot");\r
- xtitle(titre,"C (Classes)","N(C) / (Nmax length(C))"); // Not clear\r
- legend(gettext("Gaussian random sample histogram"), ..\r
- gettext("Exact gaussian density"));\r
- return\r
- end\r
-\r
- if rhs < 2\r
- error(msprintf(gettext("%s: Wrong number of input argument(s): At least %d expected.\n"),"histplot",2));\r
- end\r
-\r
- if type(n) ~= 1 | ~isreal(n)\r
- error(msprintf(gettext("%s: Wrong type for input argument #%d: Real expected.\n"),"histplot",1));\r
- elseif type(data) ~= 1 | ~isreal(data)\r
- error(msprintf(gettext("%s: Wrong type for input argument #%d: Real expected.\n"),"histplot",2));\r
- end\r
-\r
- // this is the only specific optionnal argument for histplot\r
- if ~exists('normalization','local') then, normalization=%t,end\r
-\r
- // now parse optionnal arguments to be sent to plot2d\r
- opt_arg_seq = []\r
- opt_arg_list = ["style","strf","leg","rect","nax","logflag","frameflag","axesflag"]\r
- for opt_arg = opt_arg_list\r
- if exists(opt_arg,"local") then\r
- opt_arg_seq = opt_arg_seq +","+ opt_arg + "=" + opt_arg\r
- end\r
- end\r
-\r
- p = length(data)\r
- if length(n) == 1 then // the number of classes is provided\r
- if n < 1\r
- error(msprintf(gettext("%s: Wrong value for input argument #%d: Must be in the interval %s.\n"),"histplot",1,"[1, oo)"));\r
- elseif n~=floor(n)\r
- error(msprintf(gettext("%s: Wrong type for input argument #%d: Integer expected.\n"),"histplot",1))\r
- end\r
- mind = min(data);\r
- maxd = max(data);\r
- if (mind == maxd) then\r
- mind = mind - floor(n/2); \r
- maxd = maxd + ceil(n/2);\r
- end\r
- x = linspace(mind, maxd, n+1);\r
- else // the classes are provided\r
- x = matrix(n,1,-1) // force row form\r
- if min(diff(x)) <= 0 then\r
- error(msprintf(gettext("%s: Wrong values for input argument #%d: Elements must be in increasing order.\n"),"histplot",1))\r
- end\r
- n = length(x)-1\r
- end\r
-\r
- [ind , y] = dsearch(data, x)\r
-\r
- if normalization then y=y ./ (p *(x(2:$)-x(1:$-1))),end\r
-\r
- // now form the polyline\r
- // X = [x1 x1 x2 x2 x2 x3 x3 x3 x4 ... xn xn+1 xn+1]'\r
- // Y = [0 y1 y1 0 y2 y2 0 y3 y3 ... 0 yn yn 0 ]'\r
- X = [x(1);x(1);matrix([1;1;1]*x(2:n),-1,1);x(n+1);x(n+1)]\r
- // BUG#1885\r
- // We start the histplot line to %eps rather than 0\r
- // So when switching to logarithmic mode we do not fall\r
- // in log(0) special behaviour.\r
- Y = [matrix([%eps;1;1]*y,-1,1);%eps]\r
-\r
- if opt_arg_seq == [] then\r
- plot2d(X,Y)\r
- else\r
- execstr('plot2d(X,Y'+opt_arg_seq+")")\r
- end\r
-\r
-endfunction\r
-\r
-\r
-\r
-\r
-\r\r
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA
+// Copyright (C) Bruno Pincon
+// Copyright (C) 2010 - Samuel Gougeon
+// Copyright (C) 2012 - Scilab Enterprises - Adeline CARNIS
+//
+// 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-en.txt
+
+
+function histplot(n,data,style,strf,leg,rect,nax,logflag,frameflag,axesflag,normalization)
+ // histplot(n,data,<opt_arg_seq>)
+ // draws histogram of entries in data put into n classes
+ //
+ // histplot(xi,data,<opt_arg_seq>)
+ // generates the histogram of entries in data put into classes
+ // [xi(1),xi(2)], (xi(k) xi(k+1)], k=2,..,n.
+ // xi's are assumed st. increasing (this point is verified now).
+ //
+ // optionnal args:
+ // 1/ the same than for a plot2d:
+ // style,strf,leg,rect,nax,logflag,frameflag,axesflag
+ // 2/ normalization flag (default value %t). When true the
+ // histogram is normalized so that to approach a density:
+ // xmax
+ // /
+ // | h(x) dx = 1 (true if xmin <= min(data) and max(data) <= xmax)
+ // /
+ // xmin
+ //
+ // Example : enter histplot()
+ //
+ // modifs to use dsearch (Bruno Pincon 10/12/2001)
+ // others modifs from Bruno (feb 2005):
+ // - may be a cleaner and shorter way to deal with optionnal arg ?
+ // - now the histgram is drawn using only one polyline
+ // (so properties (color, thickness,...) are easier
+ // to change with new graphics).
+ // - removed computation of nax and rect if they are not
+ // passed (let plot2d doing it)
+ // - modify a little the demo
+ // - add some checking on n|x and data
+ //
+ [lhs,rhs]=argn()
+
+ if rhs == 0 then // demo
+ histplot([-4.5:0.25:4.5],rand(1,20000,'n'),style=2,axesflag=1,..
+ frameflag=1,rect=[-4.5 0 4.5 0.47]);
+ deff("[y]=f(x)","y=exp(-x.*x/2)/sqrt(2*%pi);");
+ x=-4.5:0.125:4.5;
+ x=x';
+ plot2d(x,f(x),26,"000");
+ titre= gettext("histplot() : (normalized) histogram plot");
+ xtitle(titre,"C (Classes)","N(C) / (Nmax length(C))"); // Not clear
+ legend(gettext("Gaussian random sample histogram"), ..
+ gettext("Exact gaussian density"));
+ return
+ end
+
+ if rhs < 2
+ error(msprintf(gettext("%s: Wrong number of input argument(s): At least %d expected.\n"),"histplot",2));
+ end
+
+ if type(n) ~= 1 | ~isreal(n)
+ error(msprintf(gettext("%s: Wrong type for input argument #%d: Real expected.\n"),"histplot",1));
+ elseif type(data) ~= 1 | ~isreal(data)
+ error(msprintf(gettext("%s: Wrong type for input argument #%d: Real expected.\n"),"histplot",2));
+ end
+
+ // this is the only specific optionnal argument for histplot
+ if ~exists('normalization','local') then, normalization=%t,end
+
+ // now parse optionnal arguments to be sent to plot2d
+ opt_arg_seq = []
+ opt_arg_list = ["style","strf","leg","rect","nax","logflag","frameflag","axesflag"]
+ for opt_arg = opt_arg_list
+ if exists(opt_arg,"local") then
+ opt_arg_seq = opt_arg_seq +","+ opt_arg + "=" + opt_arg
+ end
+ end
+
+ p = length(data)
+ if length(n) == 1 then // the number of classes is provided
+ if n < 1
+ error(msprintf(gettext("%s: Wrong value for input argument #%d: Must be in the interval %s.\n"),"histplot",1,"[1, oo)"));
+ elseif n~=floor(n)
+ error(msprintf(gettext("%s: Wrong type for input argument #%d: Integer expected.\n"),"histplot",1))
+ end
+ mind = min(data);
+ maxd = max(data);
+ if (mind == maxd) then
+ mind = mind - floor(n/2);
+ maxd = maxd + ceil(n/2);
+ end
+ x = linspace(mind, maxd, n+1);
+ else // the classes are provided
+ x = matrix(n,1,-1) // force row form
+ if min(diff(x)) <= 0 then
+ error(msprintf(gettext("%s: Wrong values for input argument #%d: Elements must be in increasing order.\n"),"histplot",1))
+ end
+ n = length(x)-1
+ end
+
+ [ind , y] = dsearch(data, x)
+
+ if normalization then y=y ./ (p *(x(2:$)-x(1:$-1))),end
+
+ // now form the polyline
+ // X = [x1 x1 x2 x2 x2 x3 x3 x3 x4 ... xn xn+1 xn+1]'
+ // Y = [0 y1 y1 0 y2 y2 0 y3 y3 ... 0 yn yn 0 ]'
+ X = [x(1);x(1);matrix([1;1;1]*x(2:n),-1,1);x(n+1);x(n+1)]
+ // BUG#1885
+ // We start the histplot line to %eps rather than 0
+ // So when switching to logarithmic mode we do not fall
+ // in log(0) special behaviour.
+ Y = [matrix([%eps;1;1]*y,-1,1);%eps]
+
+ if opt_arg_seq == [] then
+ plot2d(X,Y)
+ else
+ execstr('plot2d(X,Y'+opt_arg_seq+")")
+ end
+
+endfunction
+
+
+
+
+
projects / scilab.git / commitdiff