changelog shortlog tags changeset files revisions annotate raw

scripts/statistics/base/meansq.m

changeset 10289: 4b124317dc38
parent:22a7e4690742
author: John W. Eaton <jwe@octave.org>
date: Tue Feb 09 20:58:55 2010 -0500 (65 minutes ago)
permissions: -rw-r--r--
description: base_properties::set_children: account for hidden children 
     1## Copyright (C) 1995, 1996, 1997, 1998, 2000, 2002, 2004, 2005, 2006,

     2##               2007 Kurt Hornik

     3## Copyright (C) 2009 Jaroslav Hajek

     4##

     5## This file is part of Octave.

     6##

     7## Octave is free software; you can redistribute it and/or modify it

     8## under the terms of the GNU General Public License as published by

     9## the Free Software Foundation; either version 3 of the License, or (at

    10## your option) any later version.

    11##

    12## Octave is distributed in the hope that it will be useful, but

    13## WITHOUT ANY WARRANTY; without even the implied warranty of

    14## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU

    15## General Public License for more details.

    16##

    17## You should have received a copy of the GNU General Public License

    18## along with Octave; see the file COPYING.  If not, see

    19## <http://www.gnu.org/licenses/>.

    20

    21## -*- texinfo -*-

    22## @deftypefn {Function File} {} meansq (@var{x})

    23## @deftypefnx {Function File} {} meansq (@var{x}, @var{dim})

    24## For vector arguments, return the mean square of the values.

    25## For matrix arguments, return a row vector containing the mean square

    26## of each column.  With the optional @var{dim} argument, returns the

    27## mean squared of the values along this dimension.

    28## @end deftypefn

    29

    30## Author: KH <Kurt.Hornik@wu-wien.ac.at>

    31## Description: Compute mean square

    32

    33function y = meansq (x, dim)

    34

    35  if (nargin != 1 && nargin != 2)

    36    print_usage ();

    37  endif

    38

    39  if (nargin < 2)

    40    t = find (size (x) != 1);

    41    if (isempty (t))

    42      dim = 1;

    43    else

    44      dim = t(1);

    45    endif

    46  endif

    47

    48  y = sumsq (x, dim) / size (x, dim);

    49

    50endfunction