# HG changeset patch # User jwe # Date 1069180206 0 # Node ID cae0b7c46842dc48e150f68743b50d2a21281b52 # Parent ca70857bdbd1fe2e712f1e11c888c9cca291c42d [project @ 2003-11-18 18:30:06 by jwe] diff -r ca70857bdbd1 -r cae0b7c46842 scripts/ChangeLog --- a/scripts/ChangeLog Tue Nov 18 18:20:11 2003 +0000 +++ b/scripts/ChangeLog Tue Nov 18 18:30:06 2003 +0000 @@ -1,3 +1,7 @@ +2003-11-18 Danilo Piazzalunga + + * statistics/base/iqr.m: Handle matrices. + 2003-11-18 John W. Eaton * general/issymmetric.m: Don't fail if norm (x) == 0. diff -r ca70857bdbd1 -r cae0b7c46842 scripts/statistics/base/iqr.m --- a/scripts/statistics/base/iqr.m Tue Nov 18 18:20:11 2003 +0000 +++ b/scripts/statistics/base/iqr.m Tue Nov 18 18:30:06 2003 +0000 @@ -34,6 +34,15 @@ usage ("iqr (x)"); endif - y = empirical_inv (3/4, x) - empirical_inv (1/4, x); +  if (rows (x) == 1) +    x = x.'; +  endif + +  [r, c] = size (x); +  y = zeros (1, c); + +  for i = 1:c; +    y(i) = empirical_inv (3/4, x(:,i)) - empirical_inv (1/4, x(:,i)); +  endfor endfunction