/usr/share/octave/packages/3.2/statistics-1.0.10/histfit.m is in octave-statistics 1.0.10-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 | ## Copyright (C) 2003 Alberto Terruzzi
##
## This program is free software; you can redistribute it and/or modify
## it under the terms of the GNU General Public License as published by
## the Free Software Foundation; either version 2 of the License, or
## (at your option) any later version.
##
## This program is distributed in the hope that it will be useful,
## but WITHOUT ANY WARRANTY; without even the implied warranty of
## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
## GNU General Public License for more details.
##
## You should have received a copy of the GNU General Public License
## along with this program; If not, see <http://www.gnu.org/licenses/>.
## -*- texinfo -*-
## @deftypefn {Function File} histfit (@var{data}, @var{nbins})
##
## Plot histogram with superimposed fitted normal density.
##
## @code{histfit (@var{data}, @var{nbins})} plots a histogram of the values in
## the vector @var{data} using @var{nbins} bars in the histogram. With one input
## argument, @var{nbins} is set to the square root of the number of elements in
## data.
##
## Example
##
## @example
## histfit (randn (100, 1))
## @end example
##
## @seealso{bar,hist, pareto}
## @end deftypefn
## Author: Alberto Terruzzi <t-albert@libero.it>
## Version: 1.0
## Created: 3 March 2004
function histfit (data,nbins)
if nargin < 1 || nargin > 2
usage("histfit (data,nbins)")
endif
if isvector (data) != 1
error ("data must be a vector.");
endif
row = sum(~isnan(data));
if nargin < 2
nbins = ceil(sqrt(row));
endif
[n,xbin]=hist(data,nbins);
if any(abs(diff(xbin,2)) > 10*max(abs(xbin))*eps)
error("histfit bins must be uniform width");
endif
mr = nanmean(data); ## Estimates the parameter, MU, of the normal distribution.
sr = nanstd(data); ## Estimates the parameter, SIGMA, of the normal distribution.
x=(-3*sr+mr:0.1*sr:3*sr+mr)';## Evenly spaced samples of the expected data range.
[xb,yb] = bar(xbin,n);
y = normal_pdf(x,mr,sr.^2);
binwidth = xbin(2)-xbin(1);
y = row*y*binwidth; ## Normalization necessary to overplot the histogram.
plot(xb,yb,";;b",x,y,";;r-"); ## Plots density line over histogram.
endfunction
|