changeset 27984:b09432b20a84

maint: Remove special cases of old version control keywords in code base. * common_size.m, cumtrapz.m, interp1.m, polyarea.m, trapz.m, inpolygon.m, summer.m, javamem.m, duplication_matrix.m, subspace.m, starting_stepsize.m, text.m, errorbar.m, loglogerr.m, peaks.m, __errplot__.m, semilogxerr.m, semilogyerr.m, stemleaf.m, streamtube.m, __pltopt__.m, printd.m, pchip.m, arch_fit.m, arch_rnd.m, arch_test.m, arma_rnd.m, autoreg_matrix.m, bartlett.m, blackman.m, diffpara.m, durbinlevinson.m, fractdiff.m, hamming.m, hanning.m, hurst.m, ifftshift.m, periodogram.m, rectangle_lw.m, rectangle_sw.m, triangle_lw.m, triangle_sw.m, sinetone.m, sinewave.m, spectral_adf.m, spectral_xdf.m, spencer.m, stft.m, synthesis.m, yulewalker.m, bicg.m, bicgstab.m, cgs.m, gmres.m, pcr.m, sprand.m, tfqmr.m, betaincinv.m, betaln.m, gammaincinv.m, hadamard.m, center.m, cov.m, discrete_inv.m, discrete_pdf.m, discrete_rnd.m, empirical_cdf.m, empirical_inv.m, empirical_pdf.m, empirical_rnd.m, iqr.m, kendall.m, mean.m, meansq.m, moment.m, prctile.m, quantile.m, range.m, run_count.m, spearman.m, statistics.m, var.m, zscore.m, datenum.m, datevec.m: Remove special cases of old version control keywords in code base.
author Rik <rik@octave.org>
date Tue, 21 Jan 2020 14:04:31 -0800
parents d3dc3f7f7999
children 9f9ac219896d
files scripts/general/common_size.m scripts/general/cumtrapz.m scripts/general/interp1.m scripts/general/polyarea.m scripts/general/trapz.m scripts/geometry/inpolygon.m scripts/image/summer.m scripts/java/javamem.m scripts/linear-algebra/duplication_matrix.m scripts/linear-algebra/subspace.m scripts/ode/private/starting_stepsize.m scripts/plot/appearance/text.m scripts/plot/draw/errorbar.m scripts/plot/draw/loglogerr.m scripts/plot/draw/peaks.m scripts/plot/draw/private/__errplot__.m scripts/plot/draw/semilogxerr.m scripts/plot/draw/semilogyerr.m scripts/plot/draw/stemleaf.m scripts/plot/draw/streamtube.m scripts/plot/util/__pltopt__.m scripts/plot/util/printd.m scripts/polynomial/pchip.m scripts/signal/arch_fit.m scripts/signal/arch_rnd.m scripts/signal/arch_test.m scripts/signal/arma_rnd.m scripts/signal/autoreg_matrix.m scripts/signal/bartlett.m scripts/signal/blackman.m scripts/signal/diffpara.m scripts/signal/durbinlevinson.m scripts/signal/fractdiff.m scripts/signal/hamming.m scripts/signal/hanning.m scripts/signal/hurst.m scripts/signal/ifftshift.m scripts/signal/periodogram.m scripts/signal/private/rectangle_lw.m scripts/signal/private/rectangle_sw.m scripts/signal/private/triangle_lw.m scripts/signal/private/triangle_sw.m scripts/signal/sinetone.m scripts/signal/sinewave.m scripts/signal/spectral_adf.m scripts/signal/spectral_xdf.m scripts/signal/spencer.m scripts/signal/stft.m scripts/signal/synthesis.m scripts/signal/yulewalker.m scripts/sparse/bicg.m scripts/sparse/bicgstab.m scripts/sparse/cgs.m scripts/sparse/gmres.m scripts/sparse/pcr.m scripts/sparse/sprand.m scripts/sparse/tfqmr.m scripts/specfun/betaincinv.m scripts/specfun/betaln.m scripts/specfun/gammaincinv.m scripts/special-matrix/hadamard.m scripts/statistics/center.m scripts/statistics/cov.m scripts/statistics/discrete_inv.m scripts/statistics/discrete_pdf.m scripts/statistics/discrete_rnd.m scripts/statistics/empirical_cdf.m scripts/statistics/empirical_inv.m scripts/statistics/empirical_pdf.m scripts/statistics/empirical_rnd.m scripts/statistics/iqr.m scripts/statistics/kendall.m scripts/statistics/mean.m scripts/statistics/meansq.m scripts/statistics/moment.m scripts/statistics/prctile.m scripts/statistics/quantile.m scripts/statistics/range.m scripts/statistics/run_count.m scripts/statistics/spearman.m scripts/statistics/statistics.m scripts/statistics/var.m scripts/statistics/zscore.m scripts/time/datenum.m scripts/time/datevec.m
diffstat 85 files changed, 30 insertions(+), 287 deletions(-) [+]
line wrap: on
line diff
--- a/scripts/general/common_size.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/general/common_size.m	Tue Jan 21 14:04:31 2020 -0800
@@ -47,9 +47,6 @@
 ## @seealso{size, size_equal, numel, ndims}
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Optimized-By: Jaroslav Hajek
-
 function [err, varargout] = common_size (varargin)
 
   if (nargin < 2)
--- a/scripts/general/cumtrapz.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/general/cumtrapz.m	Tue Jan 21 14:04:31 2020 -0800
@@ -49,11 +49,6 @@
 ## @seealso{trapz, cumsum}
 ## @end deftypefn
 
-## Author:      Kai Habel <kai.habel@gmx.de>
-##
-## also: June 2000 Paul Kienzle (fixes,suggestions)
-## 2006-05-12 David Bateman - Modified for NDArrays
-
 function z = cumtrapz (x, y, dim)
 
   if (nargin < 1) || (nargin > 3)
--- a/scripts/general/interp1.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/general/interp1.m	Tue Jan 21 14:04:31 2020 -0800
@@ -116,15 +116,6 @@
 ## @seealso{pchip, spline, interpft, interp2, interp3, interpn}
 ## @end deftypefn
 
-## Author: Paul Kienzle
-## Date: 2000-03-25
-##    added 'nearest' as suggested by Kai Habel
-## 2000-07-17 Paul Kienzle
-##    added '*' methods and matrix y
-##    check for proper table lengths
-## 2002-01-23 Paul Kienzle
-##    fixed extrapolation
-
 function yi = interp1 (x, y, varargin)
 
   if (nargin < 2 || nargin > 6)
--- a/scripts/general/polyarea.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/general/polyarea.m	Tue Jan 21 14:04:31 2020 -0800
@@ -39,22 +39,13 @@
 ##
 ## @end deftypefn
 
-## todo:  Add moments for centroid, etc.
+## FIXME: Add moments for centroid, etc.
 ##
-## bugs and limitations:
+## Bugs and limitations:
 ##        Probably ought to be an optional check to make sure that
 ##        traversing the vertices doesn't make any sides cross
 ##        (Is simple closed curve the technical definition of this?).
 
-## Author: David M. Doolin <doolin@ce.berkeley.edu>
-## Date: 1999-04-14
-## Modified-by:
-##    2000-01-15 Paul Kienzle <pkienzle@kienzle.powernet.co.uk>
-##    * use matlab compatible interface
-##    * return absolute value of area so traversal order doesn't matter
-##    2005-10-13 Torsten Finke
-##    * optimization saving half the sums and multiplies
-
 function a = polyarea (x, y, dim)
 
   if (nargin != 2 && nargin != 3)
--- a/scripts/general/trapz.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/general/trapz.m	Tue Jan 21 14:04:31 2020 -0800
@@ -65,11 +65,6 @@
 ## @seealso{cumtrapz}
 ## @end deftypefn
 
-## Author:      Kai Habel <kai.habel@gmx.de>
-##
-## also: June 2000 - Paul Kienzle (fixes,suggestions)
-## 2006-05-12 David Bateman - Modified for NDArrays
-
 function z = trapz (x, y, dim)
 
   if (nargin < 1) || (nargin > 3)
--- a/scripts/geometry/inpolygon.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/geometry/inpolygon.m	Tue Jan 21 14:04:31 2020 -0800
@@ -38,11 +38,7 @@
 ## @seealso{delaunay}
 ## @end deftypefn
 
-## Author: Frederick (Rick) A Niles <niles@rickniles.com>
-
-## Vectorized by Søren Hauberg <soren@hauberg.org>
-
-## The method for determining if a point is in a polygon is based on
+## Algorithm: The method for determining if a point is in a polygon is based on
 ## the algorithm shown on
 ## http://local.wasp.uwa.edu.au/~pbourke/geometry/insidepoly/
 ## and is credited to Randolph Franklin.
--- a/scripts/image/summer.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/image/summer.m	Tue Jan 21 14:04:31 2020 -0800
@@ -33,9 +33,6 @@
 ## @seealso{colormap}
 ## @end deftypefn
 
-## Author:  Kai Habel <kai.habel@gmx.de>
-## Date:  06/03/2000
-
 function map = summer (n)
 
   if (nargin > 1)
--- a/scripts/java/javamem.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/java/javamem.m	Tue Jan 21 14:04:31 2020 -0800
@@ -63,12 +63,6 @@
 ## Java's own cleaning up and your operating system's memory management.
 ## @end deftypefn
 
-## Author: Philip Nienhuis
-## Updates:
-## 2010-03-26 Changed name to javamem & indentation to double spaces
-## 2010-08-25 Corrected text on java memory assignments
-## 2010-09-05 Further overhauled help text
-
 function jmem = javamem ()
 
   rt = javaMethod ("getRuntime", "java.lang.Runtime");
--- a/scripts/linear-algebra/duplication_matrix.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/linear-algebra/duplication_matrix.m	Tue Jan 21 14:04:31 2020 -0800
@@ -65,9 +65,6 @@
 ## Calculus with Applications in Statistics and Econometrics}.
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Adapged-By: jwe
-
 function d = duplication_matrix (n)
 
   if (nargin != 1)
--- a/scripts/linear-algebra/subspace.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/linear-algebra/subspace.m	Tue Jan 21 14:04:31 2020 -0800
@@ -29,8 +29,8 @@
 ## spanned by the columns of matrices @var{A} and @var{B}.
 ## @end deftypefn
 
-## reference:
-## [1]  Andrew V. Knyazev, Merico E. Argentati:
+## Reference:
+## Andrew V. Knyazev, Merico E. Argentati:
 ## Principal Angles between Subspaces in an A-Based Scalar Product:
 ## Algorithms and Perturbation Estimates.
 ## SIAM Journal on Scientific Computing, Vol. 23 no. 6, pp. 2008-2040
--- a/scripts/ode/private/starting_stepsize.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/ode/private/starting_stepsize.m	Tue Jan 21 14:04:31 2020 -0800
@@ -34,8 +34,8 @@
 ## condition.  @var{AbsTol} and @var{RelTol} are the absolute and relative
 ## tolerance on the ODE integration taken from an ode options structure.
 ##
-## References:
-## [1] E. Hairer, S.P. Norsett and G. Wanner,
+## Reference:
+## E. Hairer, S.P. Norsett and G. Wanner,
 ## @cite{Solving Ordinary Differential Equations I: Nonstiff Problems},
 ## Springer.
 ## @end deftypefn
--- a/scripts/plot/appearance/text.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/plot/appearance/text.m	Tue Jan 21 14:04:31 2020 -0800
@@ -85,8 +85,8 @@
 ## @seealso{gtext, title, xlabel, ylabel, zlabel}
 ## @end deftypefn
 
-## Note: The following code is rigged for Matlab compatibility and is
-##       full of hidden assumptions.  Be very wary when modifying.
+## Caution: The following code is rigged for Matlab compatibility and is
+##          full of hidden assumptions.  Be very wary when modifying.
 
 function h = text (varargin)
 
--- a/scripts/plot/draw/errorbar.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/plot/draw/errorbar.m	Tue Jan 21 14:04:31 2020 -0800
@@ -170,9 +170,6 @@
 ## @seealso{semilogxerr, semilogyerr, loglogerr, plot}
 ## @end deftypefn
 
-## Author: Teemu Ikonen <tpikonen@pcu.helsinki.fi>
-## Keywords: errorbar, plotting
-
 function h = errorbar (varargin)
 
   [hax, varargin] = __plt_get_axis_arg__ ("errorbar", varargin{:});
--- a/scripts/plot/draw/loglogerr.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/plot/draw/loglogerr.m	Tue Jan 21 14:04:31 2020 -0800
@@ -54,9 +54,6 @@
 ## @seealso{errorbar, semilogxerr, semilogyerr}
 ## @end deftypefn
 
-## Author: Teemu Ikonen <tpikonen@pcu.helsinki.fi>
-## Keywords: errorbar, plotting
-
 function h = loglogerr (varargin)
 
   [hax, varargin] = __plt_get_axis_arg__ ("loglogerr", varargin{:});
--- a/scripts/plot/draw/peaks.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/plot/draw/peaks.m	Tue Jan 21 14:04:31 2020 -0800
@@ -62,7 +62,7 @@
 ## @seealso{sombrero, meshgrid, mesh, surf}
 ## @end deftypefn
 
-## Expression for the peaks function was taken from the following paper:
+## Reference: Expression for peaks function was taken from the following paper:
 ## http://www.control.hut.fi/Kurssit/AS-74.115/Material/GENALGgoga.pdf
 
 function [X_out, Y_out, Z_out] = peaks (x, y)
--- a/scripts/plot/draw/private/__errplot__.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/plot/draw/private/__errplot__.m	Tue Jan 21 14:04:31 2020 -0800
@@ -28,8 +28,6 @@
 ## Undocumented internal function.
 ## @end deftypefn
 
-## Author: Teemu Ikonen <tpikonen@pcu.helsinki.fi>
-## Keywords: errorbar, plotting
 function retval = __errplot__ (caller, hax, varargin)
 
   if (nargin < 4)
--- a/scripts/plot/draw/semilogxerr.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/plot/draw/semilogxerr.m	Tue Jan 21 14:04:31 2020 -0800
@@ -56,9 +56,6 @@
 ## @seealso{errorbar, semilogyerr, loglogerr}
 ## @end deftypefn
 
-## Author: Teemu Ikonen <tpikonen@pcu.helsinki.fi>
-## Keywords: errorbar, plotting
-
 function h = semilogxerr (varargin)
 
   [hax, varargin] = __plt_get_axis_arg__ ("semilogxerr", varargin{:});
--- a/scripts/plot/draw/semilogyerr.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/plot/draw/semilogyerr.m	Tue Jan 21 14:04:31 2020 -0800
@@ -56,9 +56,6 @@
 ## @seealso{errorbar, semilogxerr, loglogerr}
 ## @end deftypefn
 
-## Author: Teemu Ikonen <tpikonen@pcu.helsinki.fi>
-## Keywords: errorbar, plotting
-
 function h = semilogyerr (varargin)
 
   [hax, varargin] = __plt_get_axis_arg__ ("semilogyerr", varargin{:});
--- a/scripts/plot/draw/stemleaf.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/plot/draw/stemleaf.m	Tue Jan 21 14:04:31 2020 -0800
@@ -80,9 +80,6 @@
 ## @seealso{hist, printd}
 ## @end deftypefn
 
-## Author: Michael D. Godfrey <michaeldgodfrey@gmail.com>
-## Description: Compute stem and leaf plot
-
 function plotstr = stemleaf (x, caption, stem_sz)
   ## Compute and display a stem and leaf plot of the vector x.  The x
   ## vector is converted to integer by x = fix(x).  If an output argument
--- a/scripts/plot/draw/streamtube.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/plot/draw/streamtube.m	Tue Jan 21 14:04:31 2020 -0800
@@ -74,7 +74,7 @@
 ## @seealso{stream3, streamline}
 ## @end deftypefn
 
-## References:
+## Reference:
 ##
 ## @inproceedings{
 ##    title = {Visualization of 3-D vector fields - Variations on a stream},
--- a/scripts/plot/util/__pltopt__.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/plot/util/__pltopt__.m	Tue Jan 21 14:04:31 2020 -0800
@@ -114,10 +114,6 @@
 endfunction
 
 ## Really decode plot option strings.
-
-## Author: Rick Niles <niles@axp745.gsfc.nasa.gov>
-## Maintainer: jwe
-
 function [options, valid] = __pltopt1__ (caller, opt, err_on_invalid)
 
   options = __default_plot_options__ ();
--- a/scripts/plot/util/printd.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/plot/util/printd.m	Tue Jan 21 14:04:31 2020 -0800
@@ -38,9 +38,6 @@
 ## @seealso{stemleaf}
 ## @end deftypefn
 
-## Author: Michael D. Godfrey <michaeldgodfrey@gmail.com>
-## Description: Convert objects into other file formats.
-
 function pr_out = printd (obj, filename)
   ## Convert any object acceptable to disp() into various display formats.
   ## obj is the input object.
--- a/scripts/polynomial/pchip.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/polynomial/pchip.m	Tue Jan 21 14:04:31 2020 -0800
@@ -64,10 +64,8 @@
 ## @seealso{spline, ppval, mkpp, unmkpp}
 ## @end deftypefn
 
-## Author:  Kai Habel <kai.habel@gmx.de>
-## Date: 9. mar 2001
-##
-## S_k = a_k + b_k*x + c_k*x^2 + d_k*x^3; (spline polynom)
+## Algorithm:
+## S_k = a_k + b_k*x + c_k*x^2 + d_k*x^3; (spline polynomial)
 ##
 ## 4 conditions:
 ## S_k(x_k) = y_k;
--- a/scripts/signal/arch_fit.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/arch_fit.m	Tue Jan 21 14:04:31 2020 -0800
@@ -56,9 +56,6 @@
 ## for the scoring algorithm.
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: Fit an ARCH regression model
-
 function [a, b] = arch_fit (y, x, p, iter, gamma, a0, b0)
 
   if (nargin < 3 || nargin == 6 || nargin > 7)
--- a/scripts/signal/arch_rnd.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/arch_rnd.m	Tue Jan 21 14:04:31 2020 -0800
@@ -45,9 +45,6 @@
 ## @end smallexample
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: Simulate an ARCH process
-
 function y = arch_rnd (a, b, t)
 
   if (nargin != 3)
--- a/scripts/signal/arch_test.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/arch_test.m	Tue Jan 21 14:04:31 2020 -0800
@@ -69,8 +69,7 @@
 ## If no output argument is given, the @math{p}-value is displayed.
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: Test for conditional heteroscedascity
+## Description: Test for conditional heteroscedascity.
 
 function [pval, lm] = arch_test (y, x, p)
 
--- a/scripts/signal/arma_rnd.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/arma_rnd.m	Tue Jan 21 14:04:31 2020 -0800
@@ -47,9 +47,6 @@
 ## is omitted, @var{n} = 100 is used.
 ## @end deftypefn
 
-## Author: FL <Friedrich.Leisch@ci.tuwien.ac.at>
-## Description: Simulate an ARMA process
-
 function x = arma_rnd (a, b, v, t, n)
 
   if (nargin == 4)
--- a/scripts/signal/autoreg_matrix.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/autoreg_matrix.m	Tue Jan 21 14:04:31 2020 -0800
@@ -35,9 +35,6 @@
 ## The resulting matrix may be used as a regressor matrix in autoregressions.
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: Design matrix for autoregressions
-
 function X = autoreg_matrix (y, k)
 
   if (nargin != 2)
--- a/scripts/signal/bartlett.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/bartlett.m	Tue Jan 21 14:04:31 2020 -0800
@@ -33,9 +33,6 @@
 ## @cite{Discrete-Time Signal Processing}.
 ## @end deftypefn
 
-## Author: AW <Andreas.Weingessel@ci.tuwien.ac.at>
-## Description: Coefficients of the Bartlett (triangular) window
-
 function c = bartlett (m)
 
   if (nargin != 1)
--- a/scripts/signal/blackman.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/blackman.m	Tue Jan 21 14:04:31 2020 -0800
@@ -40,9 +40,6 @@
 ## @cite{Discrete-Time Signal Processing}.
 ## @end deftypefn
 
-## Author: AW <Andreas.Weingessel@ci.tuwien.ac.at>
-## Description: Coefficients of the Blackman window
-
 function c = blackman (m, opt)
 
   if (nargin < 1 || nargin > 2)
--- a/scripts/signal/diffpara.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/diffpara.m	Tue Jan 21 14:04:31 2020 -0800
@@ -44,9 +44,6 @@
 ## Theory and Methods}. @nospell{Springer} 1987.
 ## @end deftypefn
 
-## Author: FL <Friedrich.Leisch@ci.tuwien.ac.at>
-## Description: Estimate the fractional differencing parameter
-
 function [d, dd] = diffpara (x, a, b)
 
   if (nargin < 1 || nargin > 3)
--- a/scripts/signal/durbinlevinson.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/durbinlevinson.m	Tue Jan 21 14:04:31 2020 -0800
@@ -36,9 +36,6 @@
 ## the algorithm are performed.
 ## @end deftypefn
 
-## Author: FL <Friedrich.Leisch@ci.tuwien.ac.at>
-## Description: Perform one step of the Durbin-Levinson algorithm
-
 function [newphi, newv] = durbinlevinson (c, oldphi, oldv)
 
   if (nargin != 1 && nargin != 3)
--- a/scripts/signal/fractdiff.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/fractdiff.m	Tue Jan 21 14:04:31 2020 -0800
@@ -29,9 +29,6 @@
 ## denotes the lag-operator and @math{d} is greater than -1.
 ## @end deftypefn
 
-## Author: FL <Friedrich.Leisch@ci.tuwien.ac.at>
-## Description: Compute fractional differences
-
 function retval = fractdiff (x, d)
 
   if (nargin != 2)
--- a/scripts/signal/hamming.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/hamming.m	Tue Jan 21 14:04:31 2020 -0800
@@ -39,9 +39,6 @@
 ## @cite{Discrete-Time Signal Processing}.
 ## @end deftypefn
 
-## Author: AW <Andreas.Weingessel@ci.tuwien.ac.at>
-## Description: Coefficients of the Hamming window
-
 function c = hamming (m, opt)
 
   if (nargin < 1 || nargin > 2)
--- a/scripts/signal/hanning.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/hanning.m	Tue Jan 21 14:04:31 2020 -0800
@@ -39,9 +39,6 @@
 ## @cite{Discrete-Time Signal Processing}.
 ## @end deftypefn
 
-## Author: AW <Andreas.Weingessel@ci.tuwien.ac.at>
-## Description: Coefficients of the Hanning window
-
 function c = hanning (m, opt)
 
   if (nargin < 1 || nargin > 2)
--- a/scripts/signal/hurst.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/hurst.m	Tue Jan 21 14:04:31 2020 -0800
@@ -31,9 +31,6 @@
 ## If @var{x} is a matrix, the parameter is estimated for every column.
 ## @end deftypefn
 
-## Author: FL <Friedrich.Leisch@ci.tuwien.ac.at>
-## Description: Estimate the Hurst parameter
-
 function H = hurst (x)
 
   if (nargin != 1)
--- a/scripts/signal/ifftshift.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/ifftshift.m	Tue Jan 21 14:04:31 2020 -0800
@@ -33,10 +33,6 @@
 ## @seealso{fftshift}
 ## @end deftypefn
 
-## Author: Vincent Cautaerts <vincent@comf5.comm.eng.osaka-u.ac.jp>
-## Modified-By: Paul Kienzle, converted from fftshift
-## Modified-By: David Bateman, add NDArray capability and option dim arg
-
 function retval = ifftshift (x, dim)
 
   if (nargin != 1 && nargin != 2)
--- a/scripts/signal/periodogram.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/periodogram.m	Tue Jan 21 14:04:31 2020 -0800
@@ -79,9 +79,6 @@
 ## @seealso{fft}
 ## @end deftypefn
 
-## Author: FL <Friedrich.Leisch@ci.tuwien.ac.at>
-## Description: Compute the periodogram
-
 function [pxx, f] = periodogram (x, varargin)
 
   ## check input arguments
--- a/scripts/signal/private/rectangle_lw.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/private/rectangle_lw.m	Tue Jan 21 14:04:31 2020 -0800
@@ -29,9 +29,6 @@
 ## estimation.
 ## @end deftypefn
 
-## Author: FL <Friedrich.Leisch@ci.tuwien.ac.at>
-## Description: Rectangular lag window
-
 function retval = rectangle_lw (n, b)
 
   retval = zeros (n, 1);
--- a/scripts/signal/private/rectangle_sw.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/private/rectangle_sw.m	Tue Jan 21 14:04:31 2020 -0800
@@ -29,9 +29,6 @@
 ## estimation.
 ## @end deftypefn
 
-## Author: FL <Friedrich.Leisch@ci.tuwien.ac.at>
-## Description: Rectangular spectral window
-
 function retval = rectangle_sw (n, b)
 
   retval = zeros (n, 1);
--- a/scripts/signal/private/triangle_lw.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/private/triangle_lw.m	Tue Jan 21 14:04:31 2020 -0800
@@ -29,9 +29,6 @@
 ## estimation.
 ## @end deftypefn
 
-## Author: FL <Friedrich.Leisch@ci.tuwien.ac.at>
-## Description: Triangular lag window
-
 function retval = triangle_lw (n, b)
 
   retval = 1 - (0 : n-1)' * b;
--- a/scripts/signal/private/triangle_sw.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/private/triangle_sw.m	Tue Jan 21 14:04:31 2020 -0800
@@ -29,9 +29,6 @@
 ## estimation.
 ## @end deftypefn
 
-## Author: FL <Friedrich.Leisch@ci.tuwien.ac.at>
-## Description: Triangular spectral window
-
 function retval = triangle_sw (n, b)
 
   retval = zeros (n,1);
--- a/scripts/signal/sinetone.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/sinetone.m	Tue Jan 21 14:04:31 2020 -0800
@@ -34,9 +34,6 @@
 ## @seealso{sinewave}
 ## @end deftypefn
 
-## Author: FL <Friedrich.Leisch@ci.tuwien.ac.at>
-## Description: Compute a sine tone
-
 function retval = sinetone (freq, rate = 8000, sec = 1, ampl = 64)
 
   if (nargin < 1 || nargin > 4)
--- a/scripts/signal/sinewave.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/sinewave.m	Tue Jan 21 14:04:31 2020 -0800
@@ -33,9 +33,6 @@
 ## @seealso{sinetone}
 ## @end deftypefn
 
-## Author: AW <Andreas.Weingessel@ci.tuwien.ac.at>
-## Description: Compute a sine wave
-
 function x = sinewave (m, n, d)
 
   if (nargin > 0 && nargin < 4)
--- a/scripts/signal/spectral_adf.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/spectral_adf.m	Tue Jan 21 14:04:31 2020 -0800
@@ -39,9 +39,6 @@
 ## @seealso{spectral_xdf}
 ## @end deftypefn
 
-## Author: FL <Friedrich.Leisch@ci.tuwien.ac.at>
-## Description: Spectral density estimation
-
 function retval = spectral_adf (c, win, b)
 
   if (nargin < 1 || nargin > 3)
--- a/scripts/signal/spectral_xdf.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/spectral_xdf.m	Tue Jan 21 14:04:31 2020 -0800
@@ -39,9 +39,6 @@
 ## @seealso{spectral_adf}
 ## @end deftypefn
 
-## Author: FL <Friedrich.Leisch@ci.tuwien.ac.at>
-## Description: Spectral density estimation
-
 function retval = spectral_xdf (x, win, b)
 
   if (nargin < 1 || nargin > 3)
--- a/scripts/signal/spencer.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/spencer.m	Tue Jan 21 14:04:31 2020 -0800
@@ -25,13 +25,10 @@
 
 ## -*- texinfo -*-
 ## @deftypefn {} {} spencer (@var{x})
-## Return @nospell{Spencer's} 15 point moving average of each column of
+## Return @nospell{Spencer's} 15-point moving average of each column of
 ## @var{x}.
 ## @end deftypefn
 
-## Author: FL <Friedrich.Leisch@ci.tuwien.ac.at>
-## Description: Apply Spencer's 15-point MA filter
-
 function retval = spencer (x)
 
   if (nargin != 1)
--- a/scripts/signal/stft.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/stft.m	Tue Jan 21 14:04:31 2020 -0800
@@ -64,9 +64,6 @@
 ## @seealso{synthesis}
 ## @end deftypefn
 
-## Author: AW <Andreas.Weingessel@ci.tuwien.ac.at>
-## Description: Short-Time Fourier Transform
-
 function [y, c] = stft (x, win_size = 80, inc = 24, num_coef = 64, win_type = 1)
 
   if (nargin < 1 || nargin > 5)
--- a/scripts/signal/synthesis.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/synthesis.m	Tue Jan 21 14:04:31 2020 -0800
@@ -36,9 +36,6 @@
 ## @seealso{stft}
 ## @end deftypefn
 
-## Author: AW <Andreas.Weingessel@ci.tuwien.ac.at>
-## Description: Recover a signal from its short-term Fourier transform
-
 function x = synthesis (y, c)
 
   if (nargin != 2)
--- a/scripts/signal/yulewalker.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/signal/yulewalker.m	Tue Jan 21 14:04:31 2020 -0800
@@ -32,9 +32,6 @@
 ## @var{v}.
 ## @end deftypefn
 
-## Author: FL <Friedrich.Leisch@ci.tuwien.ac.at>
-## Description: Fit AR model by Yule-Walker method
-
 function [a, v] = yulewalker (c)
 
   if (nargin != 1)
--- a/scripts/sparse/bicg.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/sparse/bicg.m	Tue Jan 21 14:04:31 2020 -0800
@@ -212,15 +212,11 @@
 ## @end group
 ## @end example
 ##
-## References:
-##
-## @enumerate
+## Reference:
 ##
-## @item @nospell{Y. Saad}, @cite{Iterative Methods for Sparse Linear Systems},
+## @nospell{Y. Saad}, @cite{Iterative Methods for Sparse Linear Systems},
 ## Second edition, 2003, SIAM.
 ##
-## @end enumerate
-##
 ## @seealso{bicgstab, cgs, gmres, pcg, qmr, tfqmr}
 ## @end deftypefn
 
--- a/scripts/sparse/bicgstab.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/sparse/bicgstab.m	Tue Jan 21 14:04:31 2020 -0800
@@ -197,15 +197,11 @@
 ## @end group
 ## @end example
 ##
-## References:
-##
-## @enumerate
+## Reference:
 ##
-## @item @nospell{Y. Saad}, @cite{Iterative Methods for Sparse Linear
+## @nospell{Y. Saad}, @cite{Iterative Methods for Sparse Linear
 ## Systems}, Second edition, 2003, SIAM
 ##
-## @end enumerate
-##
 ## @seealso{bicg, cgs, gmres, pcg, qmr, tfqmr}
 ##
 ## @end deftypefn
--- a/scripts/sparse/cgs.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/sparse/cgs.m	Tue Jan 21 14:04:31 2020 -0800
@@ -187,13 +187,9 @@
 ##
 ## References:
 ##
-## @enumerate
-##
-## @item @nospell{Y. Saad}, @cite{Iterative Methods for Sparse Linear Systems},
+## @nospell{Y. Saad}, @cite{Iterative Methods for Sparse Linear Systems},
 ## Second edition, 2003, SIAM
 ##
-## @end enumerate
-##
 ## @seealso{pcg, bicgstab, bicg, gmres, qmr, tfqmr}
 ## @end deftypefn
 
--- a/scripts/sparse/gmres.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/sparse/gmres.m	Tue Jan 21 14:04:31 2020 -0800
@@ -229,14 +229,11 @@
 ## @end group
 ## @end example
 ##
-## References:
+## Reference:
 ##
-## @enumerate
-##
-## @item @nospell{Y. Saad}, @cite{Iterative Methods for Sparse Linear
+## @nospell{Y. Saad}, @cite{Iterative Methods for Sparse Linear
 ## Systems}, Second edition, 2003, SIAM
 ##
-## @end enumerate
 ## @seealso{bicg, bicgstab, cgs, pcg, pcr, qmr, tfqmr}
 ## @end deftypefn
 
--- a/scripts/sparse/pcr.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/sparse/pcr.m	Tue Jan 21 14:04:31 2020 -0800
@@ -169,9 +169,9 @@
 ## @end group
 ## @end example
 ##
-## References:
+## Reference:
 ##
-## [1] @nospell{W. Hackbusch}, @cite{Iterative Solution of Large Sparse
+## @nospell{W. Hackbusch}, @cite{Iterative Solution of Large Sparse
 ## Systems of Equations}, section 9.5.4; @nospell{Springer}, 1994
 ##
 ## @seealso{sparse, pcg}
--- a/scripts/sparse/sprand.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/sparse/sprand.m	Tue Jan 21 14:04:31 2020 -0800
@@ -47,16 +47,6 @@
 ## @seealso{sprandn, sprandsym, rand}
 ## @end deftypefn
 
-## Author: Paul Kienzle <pkienzle@users.sf.net>
-##
-## Changelog:
-##
-## Piotr Krzyzanowski <przykry2004@users.sf.net>
-##      2004-09-27      use Paul's hint to allow larger random matrices
-##                      at the price of sometimes lower density than desired
-## David Bateman
-##      2004-10-20      Texinfo help and copyright message
-
 function s = sprand (m, n, d, rc)
 
   if (nargin == 1 )
--- a/scripts/sparse/tfqmr.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/sparse/tfqmr.m	Tue Jan 21 14:04:31 2020 -0800
@@ -193,15 +193,11 @@
 ## @end group
 ## @end example
 ##
-## References:
-##
-## @enumerate
+## Reference:
 ##
-## @item @nospell{Y. Saad}, @cite{Iterative Methods for Sparse Linear Systems},
+## @nospell{Y. Saad}, @cite{Iterative Methods for Sparse Linear Systems},
 ## Second edition, 2003, SIAM
 ##
-## @end enumerate
-##
 ## @seealso{bicg, bicgstab, cgs, gmres, pcg, qmr, pcr}
 ##
 ## @end deftypefn
--- a/scripts/specfun/betaincinv.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/specfun/betaincinv.m	Tue Jan 21 14:04:31 2020 -0800
@@ -205,7 +205,7 @@
 endfunction
 
 
-## Subfunctions: Bisection and Newton Methods
+## subfunctions: Bisection and Newton Methods
 function xc = bisection_method (F, xl, xr, a, b, y, maxit)
   F_l = F (xl, a, b, y);
   F_r = F (xr, a, b, y);
--- a/scripts/specfun/betaln.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/specfun/betaln.m	Tue Jan 21 14:04:31 2020 -0800
@@ -48,9 +48,6 @@
 ## @seealso{beta, betainc, betaincinv, gammaln}
 ## @end deftypefn
 
-## Author:   Nicol N. Schraudolph <nic@idsia.ch>
-## Keywords: log beta special function
-
 function retval = betaln (a, b)
 
   if (nargin != 2)
--- a/scripts/specfun/gammaincinv.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/specfun/gammaincinv.m	Tue Jan 21 14:04:31 2020 -0800
@@ -254,7 +254,7 @@
 
 endfunction
 
-## Subfunction: Newton's Method
+## subfunction: Newton's Method
 function x = newton_method (F, JF, y, a, x0, tol, maxit);
   l = numel (y);
   res = -F (y, a, x0) ./ JF (a, x0);
--- a/scripts/special-matrix/hadamard.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/special-matrix/hadamard.m	Tue Jan 21 14:04:31 2020 -0800
@@ -57,12 +57,11 @@
 ## @seealso{compan, hankel, toeplitz}
 ## @end deftypefn
 
-
 ## Reference [1] contains a list of Hadamard matrices up to n=256.
 ## See code for h28 in hadamard.m for an example of how to extend
 ## this function for additional p.
 ##
-## References:
+## Reference:
 ## [1] A Library of Hadamard Matrices, N. J. A. Sloane
 ##     http://www.research.att.com/~njas/hadamard/
 
--- a/scripts/statistics/center.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/statistics/center.m	Tue Jan 21 14:04:31 2020 -0800
@@ -42,9 +42,6 @@
 ## @seealso{zscore}
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: Center by subtracting means
-
 function retval = center (x, dim)
 
   if (nargin != 1 && nargin != 2)
--- a/scripts/statistics/cov.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/statistics/cov.m	Tue Jan 21 14:04:31 2020 -0800
@@ -77,9 +77,6 @@
 ## @seealso{corr}
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: Compute covariances
-
 function c = cov (x, y = [], opt = 0)
 
   if (nargin < 1 || nargin > 3)
--- a/scripts/statistics/discrete_inv.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/statistics/discrete_inv.m	Tue Jan 21 14:04:31 2020 -0800
@@ -30,9 +30,6 @@
 ## @var{v} with probabilities @var{p}.
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: Quantile function of a discrete distribution
-
 function inv = discrete_inv (x, v, p)
 
   if (nargin != 3)
--- a/scripts/statistics/discrete_pdf.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/statistics/discrete_pdf.m	Tue Jan 21 14:04:31 2020 -0800
@@ -30,9 +30,6 @@
 ## in @var{v} with probabilities @var{p}.
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: PDF of a discrete distribution
-
 function pdf = discrete_pdf (x, v, p)
 
   if (nargin != 3)
--- a/scripts/statistics/discrete_rnd.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/statistics/discrete_rnd.m	Tue Jan 21 14:04:31 2020 -0800
@@ -41,9 +41,6 @@
 ## @var{v} and @var{p}.
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: Random deviates from a discrete distribution
-
 function rnd = discrete_rnd (v, p, varargin)
 
   if (nargin < 2)
--- a/scripts/statistics/empirical_cdf.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/statistics/empirical_cdf.m	Tue Jan 21 14:04:31 2020 -0800
@@ -30,9 +30,6 @@
 ## the univariate sample @var{data}.
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: CDF of the empirical distribution
-
 function cdf = empirical_cdf (x, data)
 
   if (nargin != 2)
--- a/scripts/statistics/empirical_inv.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/statistics/empirical_inv.m	Tue Jan 21 14:04:31 2020 -0800
@@ -30,9 +30,6 @@
 ## univariate sample @var{data}.
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: Quantile function of the empirical distribution
-
 function inv = empirical_inv (x, data)
 
   if (nargin != 2)
--- a/scripts/statistics/empirical_pdf.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/statistics/empirical_pdf.m	Tue Jan 21 14:04:31 2020 -0800
@@ -30,9 +30,6 @@
 ## univariate sample @var{data}.
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: PDF of the empirical distribution
-
 function pdf = empirical_pdf (x, data)
 
   if (nargin != 2)
--- a/scripts/statistics/empirical_rnd.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/statistics/empirical_rnd.m	Tue Jan 21 14:04:31 2020 -0800
@@ -41,9 +41,6 @@
 ## of the sample @var{data}.
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: Bootstrap samples from the empirical distribution
-
 function rnd = empirical_rnd (data, varargin)
 
   if (nargin < 1)
--- a/scripts/statistics/iqr.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/statistics/iqr.m	Tue Jan 21 14:04:31 2020 -0800
@@ -39,9 +39,6 @@
 ## @seealso{bounds, mad, range, std}
 ## @end deftypefn
 
-## Author KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: Interquartile range
-
 function y = iqr (x, dim)
 
   if (nargin != 1 && nargin != 2)
--- a/scripts/statistics/kendall.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/statistics/kendall.m	Tue Jan 21 14:04:31 2020 -0800
@@ -91,9 +91,6 @@
 ## @seealso{ranks, spearman}
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: Kendall's rank correlation tau
-
 function tau = kendall (x, y = [])
 
   if (nargin < 1 || nargin > 2)
--- a/scripts/statistics/mean.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/statistics/mean.m	Tue Jan 21 14:04:31 2020 -0800
@@ -89,9 +89,6 @@
 ## @seealso{median, mode}
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: Compute arithmetic, geometric, and harmonic mean
-
 function y = mean (x, varargin)
 
   if (nargin < 1 || nargin > 4)
--- a/scripts/statistics/meansq.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/statistics/meansq.m	Tue Jan 21 14:04:31 2020 -0800
@@ -55,9 +55,6 @@
 ## @seealso{var, std, moment}
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: Compute mean square
-
 function y = meansq (x, dim)
 
   if (nargin != 1 && nargin != 2)
--- a/scripts/statistics/moment.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/statistics/moment.m	Tue Jan 21 14:04:31 2020 -0800
@@ -128,9 +128,6 @@
 ## Can easily be made to work for continuous distributions (using quad)
 ## as well, but how does the general case work?
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: Compute moments
-
 function m = moment (x, p, opt1, opt2)
 
   if (nargin < 2 || nargin > 4)
--- a/scripts/statistics/prctile.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/statistics/prctile.m	Tue Jan 21 14:04:31 2020 -0800
@@ -44,9 +44,6 @@
 ## @seealso{quantile}
 ## @end deftypefn
 
-## Author: Ben Abbott <bpabbott@mac.com>
-## Description: Matlab style prctile function.
-
 function q = prctile (x, p = [], dim)
 
   if (nargin < 1 || nargin > 3)
--- a/scripts/statistics/quantile.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/statistics/quantile.m	Tue Jan 21 14:04:31 2020 -0800
@@ -154,10 +154,6 @@
 ## @seealso{prctile}
 ## @end deftypefn
 
-## Author: Ben Abbott <bpabbott@mac.com>
-## Description: Matlab style quantile function of a discrete/continuous
-##              distribution.
-
 function q = quantile (x, p = [], dim, method = 5)
 
   if (nargin < 1 || nargin > 4)
@@ -406,10 +402,7 @@
 ## (https://www.r-project.org/).  The default is @var{method} = 7.
 ## @seealso{prctile, quantile, statistics}
 
-## Author: Ben Abbott <bpabbott@mac.com>
-## Vectorized version: Jaroslav Hajek <highegg@gmail.com>
 ## Description: Quantile function of empirical samples
-
 function inv = __quantile__ (x, p, method = 5)
 
   if (nargin < 2 || nargin > 3)
--- a/scripts/statistics/range.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/statistics/range.m	Tue Jan 21 14:04:31 2020 -0800
@@ -40,9 +40,6 @@
 ## @seealso{bounds, iqr, mad, std}
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: Compute range
-
 function y = range (x, dim)
 
   if (nargin < 1 || nargin > 2)
--- a/scripts/statistics/run_count.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/statistics/run_count.m	Tue Jan 21 14:04:31 2020 -0800
@@ -34,9 +34,6 @@
 ## @seealso{runlength}
 ## @end deftypefn
 
-## Author: FL <Friedrich.Leisch@ci.tuwien.ac.at>
-## Description: Count upward runs
-
 function retval = run_count (x, n, dim)
 
   if (nargin != 2 && nargin != 3)
--- a/scripts/statistics/spearman.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/statistics/spearman.m	Tue Jan 21 14:04:31 2020 -0800
@@ -66,9 +66,6 @@
 ## @seealso{ranks, kendall}
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: Spearman's rank correlation rho
-
 function rho = spearman (x, y = [])
 
   if (nargin < 1 || nargin > 2)
--- a/scripts/statistics/statistics.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/statistics/statistics.m	Tue Jan 21 14:04:31 2020 -0800
@@ -37,9 +37,6 @@
 ## @seealso{min, max, median, mean, std, skewness, kurtosis}
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: Compute basic statistics
-
 function stats = statistics (x, dim)
 
   if (nargin != 1 && nargin != 2)
--- a/scripts/statistics/var.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/statistics/var.m	Tue Jan 21 14:04:31 2020 -0800
@@ -72,9 +72,6 @@
 ## @seealso{cov, std, skewness, kurtosis, moment}
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: Compute variance
-
 function retval = var (x, opt = 0, dim)
 
   if (nargin < 1 || nargin > 3)
--- a/scripts/statistics/zscore.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/statistics/zscore.m	Tue Jan 21 14:04:31 2020 -0800
@@ -47,9 +47,6 @@
 ## @seealso{mean, std, center}
 ## @end deftypefn
 
-## Author: KH <Kurt.Hornik@wu-wien.ac.at>
-## Description: Subtract mean and divide by standard deviation
-
 function [z, mu, sigma] = zscore (x, opt = 0, dim)
 
   if (nargin < 1 || nargin > 3 )
--- a/scripts/time/datenum.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/time/datenum.m	Tue Jan 21 14:04:31 2020 -0800
@@ -102,6 +102,7 @@
 ## @end deftypefn
 
 ## Algorithm: Peter Baum (http://vsg.cape.com/~pbaum/date/date0.htm)
+
 function [days, secs] = datenum (year, month = [], day = [], hour = 0, minute = 0, second = 0)
 
   ## Days until start of month assuming year starts March 1.
--- a/scripts/time/datevec.m	Tue Jan 21 13:20:01 2020 -0800
+++ b/scripts/time/datevec.m	Tue Jan 21 14:04:31 2020 -0800
@@ -51,8 +51,6 @@
 
 ## Algorithm: Peter Baum (http://vsg.cape.com/~pbaum/date/date0.htm)
 
-## Author: pkienzle <pkienzle@users.sf.net>
-
 ## The function __date_str2vec__ is based on datesplit by Bill Denney.
 
 function [y, m, d, h, mi, s] = datevec (date, f = [], p = [])