added checks for datatype

Author: acp29

inst/bootbayes.m

@@ -197,13 +197,21 @@
     error ('bootbayes: DATA must be provided')
   end
   [n, q] = size (Y);
+  % Check that Y contains floating point numbers
+  if (~ any (strcmpi (class (Y), {'single', 'double'})))
+    error ('bootwild: Y must contain single or double precision numbers.');
+  end
 
   % Evaluate the design matrix
   if ( (nargin < 2) || (isempty (X)) )
     X = ones (n, 1);
   elseif (size (X, 1) ~= n)
     error ('bootbayes: X must have the same number of rows as y')
   end
+  % Check that Y contains floating point numbers
+  if (~ any (strcmpi (class (X), {'single', 'double'})))
+    error ('bootwild: Y must contain single or double precision numbers.');
+  end
 
   % Remove rows of the data whose outcome or value of any predictor is NaN or Inf
   excl = any ([isnan([Y, X]), isinf([Y, X])], 2);

inst/bootridge.m

@@ -551,12 +551,20 @@
   % Get dimensions of the data
   [m, n] = size (X);
   q = size (Y, 2);
+  % Check that Y contains floating point numbers
+  if (~ any (strcmpi (class (Y), {'single', 'double'})))
+    error ('bootwild: Y must contain single or double precision numbers.');
+  end
 
   % Check that the first column is X are all equal to 1, if not create one
   if ( ~all (X(:, 1) == 1) )
     X = cat (2, ones (m, 1), X);
     n = n + 1;
   end
+  % Check that X contains floating point numbers
+  if (~ any (strcmpi (class (X), {'single', 'double'})))
+    error ('bootwild: X must contain single or double precision numbers.');
+  end
 
   % If categor is not provided, set it to empty
   if ( (nargin < 3) || isempty (categor) )
@@ -629,7 +637,7 @@
     seed = 1;
   else
     if ( isinf (seed) || isnan (seed) || (numel (seed) > 1) || ...
-         seed ~= fix(seed))
+         seed ~= fix (seed))
       error ('bootridge: The seed must be a finite integer');
     end
   end
@@ -1408,6 +1416,19 @@
 
 %-------------------------------------------------------------------------------
 
+%!demo
+%!
+%! % Analysis of Rohwer's dataset
+%! [group, SES, SAT, PPVT, Raven, n, s, ns, na, ss] = ...
+%!      textread ('./rohwer_data.csv', '%f %s %f %f %f %f %f %f %f %f', ...
+%!                'Delimiter', ',', 'HeaderLines', 1);
+%! MAT = bootlm (SAT, {SES, n, s, ns, na, ss}, 'model', 'linear', ...
+%!                     'nboot', 0, 'display', 'off', 'continuous', [2:6], ...
+%!                     'contrasts', 'helmert');
+%! 
+%! % Analysis using bootridge
+%! bootridge([SAT, PPVT, Raven], MAT.X, 2, 1999, .05);
+
 %!demo
 %!
 %! % Simple linear regression. The data represents salaries of employees and

inst/bootwild.m

@@ -164,7 +164,7 @@
     end
   end
 
-  % Calculate the length of y
+  % Calculate the length and datatype of y
   if (nargin < 1)
     error ('bootwild: DATA must be provided')
   end
@@ -173,13 +173,21 @@
     error ('bootwild: y must be a column vector')
   end
   n = numel (y);
+  % Check that y contains floating point numbers
+  if (~ any (strcmpi (class (y), {'single', 'double'})))
+    error ('bootwild: y must contain single or double precision numbers.');
+  end
 
   % Evaluate the design matrix
   if ( (nargin < 2) || (isempty (X)) )
     X = ones (n, 1);
   elseif (size (X, 1) ~= n)
     error ('bootwild: X must have the same number of rows as y')
   end
+  % Check that X contains floating point numbers
+  if (~ any (strcmpi (class (X), {'single', 'double'})))
+    error ('bootwild: X must contain single or double precision numbers.');
+  end
 
   % Remove rows of the data whose outcome or value of any predictor is NaN or Inf
   excl = any ([isnan([y, X]), isinf([y, X])], 2);