|
7220
|
1 ## Copyright (C) 2007 David Bateman |
|
|
2 ## Copyright (C) 2003 Alberto Terruzzi |
|
|
3 ## |
|
|
4 ## This file is part of Octave. |
|
|
5 ## |
|
|
6 ## Octave is free software; you can redistribute it and/or modify it |
|
|
7 ## under the terms of the GNU General Public License as published by |
|
|
8 ## the Free Software Foundation; either version 3 of the License, or (at |
|
|
9 ## your option) any later version. |
|
|
10 ## |
|
|
11 ## Octave is distributed in the hope that it will be useful, but |
|
|
12 ## WITHOUT ANY WARRANTY; without even the implied warranty of |
|
|
13 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
|
|
14 ## General Public License for more details. |
|
|
15 ## |
|
|
16 ## You should have received a copy of the GNU General Public License |
|
|
17 ## along with Octave; see the file COPYING. If not, see |
|
|
18 ## <http://www.gnu.org/licenses/>. |
|
|
19 |
|
|
20 ## -*- texinfo -*- |
|
|
21 ## @deftypefn {Function File} {} pareto (@var{x}) |
|
|
22 ## @deftypefnx {Function File} {} pareto (@var{x}, @var{y}) |
|
|
23 ## @deftypefnx {Function File} {} pareto (@var{h}, @dots{}) |
|
|
24 ## @deftypefnx {Function File} {@var{h} =} pareto (@dots{}) |
|
|
25 ## Draw a Pareto chart, also called ABC chart. A Pareto chart is a bar graph |
|
|
26 ## used to arrange information in such a way that priorities for process |
|
|
27 ## improvement can be established. It organizes and displays information |
|
|
28 ## to show the relative importance of data. The chart is similar to the |
|
|
29 ## histogram or bar chart, except that the bars are arranged in decreasing |
|
|
30 ## order from left to right along the abscissa. |
|
|
31 ## |
|
|
32 ## The fundamental idea (Pareto principle) behind the use of Pareto |
|
|
33 ## diagrams is that the majority of an effect is due to a small subset of the |
|
|
34 ## causes, so for quality improvement the first few (as presented on the |
|
|
35 ## diagram) contributing causes to a problem usually account for the majority |
|
|
36 ## of the result. Thus, targeting these "major causes" for elimination |
|
|
37 ## results in the most cost-effective improvement scheme. |
|
|
38 ## |
|
|
39 ## The data are passed as @var{x} and the abscissa as @var{y}. If @var{y} is |
|
|
40 ## absent, then the abscissa are assumed to be @code{1 : length (@var{x})}. |
|
|
41 ## @var{y} can be a string array, a cell array of strings or a numerical |
|
|
42 ## vector. |
|
|
43 ## |
|
|
44 ## An example of the use of @code{pareto} is |
|
|
45 ## |
|
|
46 ## @example |
|
|
47 ## @group |
|
|
48 ## Cheese = @{"Cheddar", "Swiss", "Camembert", "Munster", "Stilton", ... |
|
|
49 ## "Blue"@}; |
|
|
50 ## Sold = [105, 30, 70, 10, 15, 20]; |
|
|
51 ## pareto(Sold, Cheese); |
|
|
52 ## @end group |
|
|
53 ## @end example |
|
|
54 ## @end deftypefn |
|
|
55 |
|
|
56 function h = pareto (varargin) |
|
|
57 |
|
|
58 [ax, varargin, nargin] = __plt_get_axis_arg__ ("pareto", varargin{:}); |
|
|
59 |
|
|
60 if (nargin != 1 && nargin != 2) |
|
|
61 print_usage (); |
|
|
62 endif |
|
|
63 |
|
|
64 x = varargin {1}(:).'; |
|
|
65 if (nargin == 2) |
|
|
66 y = varargin {2}(:).'; |
|
|
67 if (! iscell (y)) |
|
|
68 if (ischar (y)) |
|
|
69 y = cellstr (y); |
|
|
70 else |
|
|
71 y = num2cell (y); |
|
|
72 endif |
|
|
73 endif |
|
|
74 else |
|
|
75 y = cellfun (@(x) int2str (x), num2cell (1 : numel(x)), |
|
|
76 "UniformOutput", false); |
|
|
77 endif |
|
|
78 |
|
|
79 [x, idx] = sort (x, "descend"); |
|
|
80 y = y (idx); |
|
|
81 cdf = cumsum (x); |
|
|
82 maxcdf = cdf(end); |
|
|
83 cdf = cdf ./ cdf (end); |
|
|
84 [dummy, idx95] = min (abs (cdf - .95)); |
|
|
85 idx95 - idx95(1); |
|
|
86 |
|
|
87 [ax, hbar, hline] = plotyy (ax, 1 : idx95, x (1 : idx95), |
|
|
88 1 : length(cdf), 100 .* cdf, |
|
|
89 @bar, @plot); |
|
|
90 |
|
|
91 axis (ax(1), [1 - 0.6, idx95 + 0.6, 0, maxcdf]); |
|
|
92 axis (ax(2), [1 - 0.6, idx95 + 0.6, 0, 100]); |
|
|
93 set (ax(2), "ytick", [0, 20, 40, 60, 80, 100], |
|
|
94 "yticklabel", {"0%", "20%", "40%", "60%", "80%", "100%"}); |
|
|
95 set (ax(1), "xtick", 1 : idx95, "xticklabel", y (1: idx95)); |
|
|
96 set (ax(2), "xtick", 1 : idx95, "xticklabel", y (1: idx95)); |
|
|
97 |
|
|
98 if (nargout > 0) |
|
|
99 h = [hbar; hline]; |
|
|
100 endif |
|
|
101 |
|
|
102 endfunction |
|
|
103 |
|
|
104 %!demo |
|
|
105 %! close |
|
|
106 %! Cheese = {"Cheddar", "Swiss", "Camembert", "Munster", "Stilton", "Blue"}; |
|
|
107 %! Sold = [105, 30, 70, 10, 15, 20]; |
|
|
108 %! pareto(Sold, Cheese); |
|
|
109 |
|
|
110 %!demo |
|
|
111 %! close |
|
|
112 %! % Suppose that we want establish which products makes 80 % of turnover. |
|
|
113 %! Codes = {"AB4","BD7","CF8","CC5","AD11","BB5","BB3","AD8","DF3","DE7"}; |
|
|
114 %! Value = [2.35 7.9 2.45 1.1 0.15 13.45 5.4 2.05 0.85 1.65]'; |
|
|
115 %! SoldUnits = [54723 41114 16939 1576091 168000 687197 120222 168195, ... |
|
|
116 %! 1084118 55576]'; |
|
|
117 %! pareto (Value.*SoldUnits, Codes); |
