bug fix

gnu-octave · May 27, 2022 · 810758d · 810758d
1 parent 3cdeebb
commit 810758d
Show file tree

Hide file tree

Showing 3 changed files with 5 additions and 7 deletions.
diff --git a/DESCRIPTION b/DESCRIPTION
@@ -1,5 +1,5 @@
 name: statistics-bootstrap
-version: 3.7.9
+version: 3.7.10
 date: 2022-05-27
 author: Andrew Penn <andy.c.penn@gmail.com>
 maintainer: Andrew Penn <andy.c.penn@gmail.com>

diff --git a/inst/helper/maxstat.m b/inst/helper/maxstat.m
@@ -44,9 +44,8 @@
         % Quick calculation for the standard error of the mean
         SE(j) = std(Y(g==gk(j),:),0) / sqrt(nk(j));
       elseif strcmp (func2str (bootfun), 'smoothmedian')
-        theta(j) = smoothmedian(Y(g==gk(j),:));
-        % Quick calculation for the smoothmedian and it's standard error
-        [theta(j), SE(j)] = feval(bootfunY(g==gk(j),:));
+        % Quick calculation for the smoothed median and it's standard error
+        [theta(j), SE(j)] = feval(bootfun, Y(g==gk(j),:));
       else
         theta(j) = feval(bootfun,Y(g==gk(j),:));
         % If requested, compute unbiased estimates of the standard error using jackknife resampling

diff --git a/inst/ibootnhst.m b/inst/ibootnhst.m
@@ -970,9 +970,8 @@
         % Quick calculation for the standard error of the mean
         SE(j) = std(data(g==gk(j),:),0) / sqrt(nk(j));
       elseif strcmp (func2str (bootfun), 'smoothmedian')
-        theta(j) = smoothmedian(data(g==gk(j),:));
-        % Quick calculation for the smoothmedian and it's standard error
-        [theta(j), SE(j)] = feval(bootfun,data(g==gk(j),:));
+        % Quick calculation for the smoothed median and it's standard error
+        [theta(j), SE(j)] = feval(bootfun, data(g==gk(j),:));
       else
         theta(j) = feval(bootfun,data(g==gk(j),:));
         % If requested, compute unbiased estimates of the standard error using jackknife resampling