Mercurial > forge

changeset 3301:fb8523a8c81c octave-forge

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix incorrect interpretation of prctile([0,1],[25,75])
author pkienzle
date Tue, 03 Apr 2007 02:02:19 +0000
parents cac692d87f7e
children 262165a1bc20
files main/statistics/inst/prctile.m
diffstat 1 files changed, 4 insertions(+), 1 deletions(-) [+]
line wrap: on
line diff
--- a/main/statistics/inst/prctile.m	Sun Apr 01 07:30:30 2007 +0000
+++ b/main/statistics/inst/prctile.m	Tue Apr 03 02:02:19 2007 +0000
@@ -35,5 +35,8 @@
   if size (y,1) == 1, y = y(:); endif
   trim = 1 + (size(y,1)-1)*p(:)*0.01;
   delta = (trim - floor(trim))*ones(1,size(y,2));
-  a = y(floor(trim), :) .* delta + y(ceil(trim), :) .* (1-delta);
+  a = y(floor(trim), :) .* (1-delta) + y(ceil(trim), :) .* delta;
 endfunction
+
+%!assert(prctile([0,1],[25,75]),[.25;.75],eps)
+%!assert(prctile([0,0;1,2],[25,75]),[.25,.5;.75,1.5],eps)