utils/matlab/stats.m

function [varargout] = stats(A)
% stats(A) writes the basic statistics of the A to the terminal which
% are i) the minimum finite value
%    ii) the maximum finite value
%   iii) the mean of the finite values
%    iv) the S.D. of the finite values ( RMS of [A-mean] )
%     v) the fraction of non-finite elements excluded from calculations
%
% e.g.
%  >> stats(topo)
% Min -4555     Max 0         Mean -2331.07  SD 1207.3441   N-Z 1024/1024
%
% [Min Max Mean SD]=stats(topo); returns the statistics in Min, Max,
% Mean and SD and does not write to the terminal.

% $Header: /u/gcmpack/MITgcm/utils/matlab/stats.m,v 1.4 2004/06/04 17:03:50 adcroft Exp $
% $Name:  $

A=A(:);

%ii=find(A~=0);
ii=find(isfinite(A));
if isempty(ii)
 ii=1;
end
sZ=prod(size(A));
minA=min(A(ii));
maxA=max(A(ii));
meanA=mean(A(ii));
sdA=sqrt( mean( (A(ii)-meanA).^2 ) );
rmsA=sqrt( mean( A(ii).^2 ) );
nZ=sum(~isfinite(A));
switch max(nargout)
 case {0}
% disp( ...
%  sprintf('Min %0.5g  Max %0.5g  Mean %0.5g  SD %0.5g   NaN %i/%i',...
%        minA,maxA,meanA,sdA,nZ,sZ) );
  disp( ...
   sprintf('Min %0.5g  Max %0.5g  Mean %0.5g  RMS %0.5g   NaN %i/%i',...
         minA,maxA,meanA,rmsA,nZ,sZ) );
 case {1}
  varargout(1)={minA};
 case {2}
  varargout(1)={minA};
  varargout(2)={maxA};
 case {3}
  varargout(1)={minA};
  varargout(2)={maxA};
  varargout(3)={meanA};
 case {4}
  varargout(1)={minA};
  varargout(2)={maxA};
  varargout(3)={meanA};
  varargout(4)={sdA};
 otherwise
  error('Too many return arguments requested')
end
1	adcroft	1.2	function [varargout] = stats(A)
2			% stats(A) writes the basic statistics of the A to the terminal which
3			% are i) the minimum finite value
4			% ii) the maximum finite value
5			% iii) the mean of the finite values
6			% iv) the S.D. of the finite values ( RMS of [A-mean] )
7			% v) the fraction of non-finite elements excluded from calculations
8			%
9			% e.g.
10			% >> stats(topo)
11			% Min -4555 Max 0 Mean -2331.07 SD 1207.3441 N-Z 1024/1024
12			%
13			% [Min Max Mean SD]=stats(topo); returns the statistics in Min, Max,
14			% Mean and SD and does not write to the terminal.
15	jmc	1.5
16			% $Header: /u/gcmpack/MITgcm/utils/matlab/stats.m,v 1.4 2004/06/04 17:03:50 adcroft Exp $
17			% $Name: $
18	adcroft	1.3
19			A=A(:);
20	adcroft	1.2
21			%ii=find(A~=0);
22			ii=find(isfinite(A));
23			if isempty(ii)
24			ii=1;
25			end
26			sZ=prod(size(A));
27			minA=min(A(ii));
28			maxA=max(A(ii));
29			meanA=mean(A(ii));
30			sdA=sqrt( mean( (A(ii)-meanA).^2 ) );
31	adcroft	1.4	rmsA=sqrt( mean( A(ii).^2 ) );
32	adcroft	1.3	nZ=sum(~isfinite(A));
33	adcroft	1.2	switch max(nargout)
34			case {0}
35	adcroft	1.4	% disp( ...
36			% sprintf('Min %0.5g Max %0.5g Mean %0.5g SD %0.5g NaN %i/%i',...
37			% minA,maxA,meanA,sdA,nZ,sZ) );
38	adcroft	1.2	disp( ...
39	adcroft	1.4	sprintf('Min %0.5g Max %0.5g Mean %0.5g RMS %0.5g NaN %i/%i',...
40			minA,maxA,meanA,rmsA,nZ,sZ) );
41	adcroft	1.2	case {1}
42			varargout(1)={minA};
43			case {2}
44			varargout(1)={minA};
45			varargout(2)={maxA};
46			case {3}
47			varargout(1)={minA};
48			varargout(2)={maxA};
49			varargout(3)={meanA};
50			case {4}
51			varargout(1)={minA};
52			varargout(2)={maxA};
53			varargout(3)={meanA};
54			varargout(4)={sdA};
55			otherwise
56			error('Too many return arguments requested')
57			end