Mercurial > hg > octave-lyh
annotate scripts/plot/pareto.m @ 9040:dbd0c77e575e
Cleanup documentation file plot.texi
Spellcheck
Stylecheck (Mostly double spaces after periods)
| author | Rik <rdrider0-list@yahoo.com> |
|---|---|
| date | Sun, 22 Mar 2009 14:40:24 -0700 |
| parents | eb63fbe60fab |
| children | 5247e89688e1 |
| rev | line source |
|---|---|
| 8920 | 1 ## Copyright (C) 2007, 2008, 2009 David Bateman |
| 7220 | 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{}) | |
|
9040
dbd0c77e575e
Cleanup documentation file plot.texi
Rik <rdrider0-list@yahoo.com>
parents:
8920
diff
changeset
|
25 ## Draw a Pareto chart, also called ABC chart. A Pareto chart is a bar graph |
| 7220 | 26 ## used to arrange information in such a way that priorities for process |
|
9040
dbd0c77e575e
Cleanup documentation file plot.texi
Rik <rdrider0-list@yahoo.com>
parents:
8920
diff
changeset
|
27 ## improvement can be established. It organizes and displays information |
|
dbd0c77e575e
Cleanup documentation file plot.texi
Rik <rdrider0-list@yahoo.com>
parents:
8920
diff
changeset
|
28 ## to show the relative importance of data. The chart is similar to the |
| 7220 | 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 | |
|
9040
dbd0c77e575e
Cleanup documentation file plot.texi
Rik <rdrider0-list@yahoo.com>
parents:
8920
diff
changeset
|
36 ## of the result. Thus, targeting these "major causes" for elimination |
| 7220 | 37 ## results in the most cost-effective improvement scheme. |
| 38 ## | |
|
9040
dbd0c77e575e
Cleanup documentation file plot.texi
Rik <rdrider0-list@yahoo.com>
parents:
8920
diff
changeset
|
39 ## The data are passed as @var{x} and the abscissa as @var{y}. If @var{y} is |
| 7220 | 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 | |
|
9040
dbd0c77e575e
Cleanup documentation file plot.texi
Rik <rdrider0-list@yahoo.com>
parents:
8920
diff
changeset
|
48 ## Cheese = @{"Cheddar", "Swiss", "Camembert", "Munster", "Stilton", @dots{} |
| 7220 | 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 if (nargin != 1 && nargin != 2) | |
| 59 print_usage (); | |
| 60 endif | |
| 61 | |
| 62 x = varargin {1}(:).'; | |
| 63 if (nargin == 2) | |
| 64 y = varargin {2}(:).'; | |
| 65 if (! iscell (y)) | |
| 66 if (ischar (y)) | |
| 67 y = cellstr (y); | |
| 68 else | |
| 7228 | 69 y = cellfun (@(x) num2str (x), num2cell (y), "UniformOutput", false); |
| 7220 | 70 endif |
| 71 endif | |
| 72 else | |
| 73 y = cellfun (@(x) int2str (x), num2cell (1 : numel(x)), | |
| 74 "UniformOutput", false); | |
| 75 endif | |
| 76 | |
| 77 [x, idx] = sort (x, "descend"); | |
| 78 y = y (idx); | |
| 79 cdf = cumsum (x); | |
| 7221 | 80 maxcdf = max(cdf); |
| 81 cdf = cdf ./ maxcdf; | |
| 7228 | 82 cdf95 = cdf - 0.95; |
| 83 idx95 = find(sign(cdf95(1:end-1)) != sign(cdf95(2:end)))(1); | |
| 7220 | 84 |
| 8208 | 85 [ax, hbar, hline] = plotyy (1 : idx95, x (1 : idx95), |
| 7220 | 86 1 : length(cdf), 100 .* cdf, |
| 87 @bar, @plot); | |
| 88 | |
| 89 axis (ax(1), [1 - 0.6, idx95 + 0.6, 0, maxcdf]); | |
| 90 axis (ax(2), [1 - 0.6, idx95 + 0.6, 0, 100]); | |
| 91 set (ax(2), "ytick", [0, 20, 40, 60, 80, 100], | |
| 92 "yticklabel", {"0%", "20%", "40%", "60%", "80%", "100%"}); | |
| 93 set (ax(1), "xtick", 1 : idx95, "xticklabel", y (1: idx95)); | |
| 94 set (ax(2), "xtick", 1 : idx95, "xticklabel", y (1: idx95)); | |
| 95 | |
| 96 if (nargout > 0) | |
| 97 h = [hbar; hline]; | |
| 98 endif | |
| 99 | |
| 100 endfunction | |
| 101 | |
| 102 %!demo | |
|
8790
a013ff655ca4
Trivial changes to demos to produce a more pleasant output for octave+gnuplot+aquaterm.
Ben Abbott <bpabbott@mac.com>
parents:
8208
diff
changeset
|
103 %! clf |
|
a013ff655ca4
Trivial changes to demos to produce a more pleasant output for octave+gnuplot+aquaterm.
Ben Abbott <bpabbott@mac.com>
parents:
8208
diff
changeset
|
104 %! colormap (jet (64)) |
| 7220 | 105 %! Cheese = {"Cheddar", "Swiss", "Camembert", "Munster", "Stilton", "Blue"}; |
| 106 %! Sold = [105, 30, 70, 10, 15, 20]; | |
| 107 %! pareto(Sold, Cheese); | |
| 108 | |
| 109 %!demo | |
|
8790
a013ff655ca4
Trivial changes to demos to produce a more pleasant output for octave+gnuplot+aquaterm.
Ben Abbott <bpabbott@mac.com>
parents:
8208
diff
changeset
|
110 %! clf |
| 7220 | 111 %! % Suppose that we want establish which products makes 80 % of turnover. |
| 112 %! Codes = {"AB4","BD7","CF8","CC5","AD11","BB5","BB3","AD8","DF3","DE7"}; | |
| 113 %! Value = [2.35 7.9 2.45 1.1 0.15 13.45 5.4 2.05 0.85 1.65]'; | |
| 114 %! SoldUnits = [54723 41114 16939 1576091 168000 687197 120222 168195, ... | |
| 115 %! 1084118 55576]'; | |
| 116 %! pareto (Value.*SoldUnits, Codes); |