diff scripts/statistics/tests/f_test_regression.m @ 3457:e031284eea27

[project @ 2000-01-19 08:49:56 by jwe]
author jwe
date Wed, 19 Jan 2000 08:50:14 +0000
parents 434790acb067
children 3e3e14ad5149
line wrap: on
line diff
--- a/scripts/statistics/tests/f_test_regression.m
+++ b/scripts/statistics/tests/f_test_regression.m
@@ -40,22 +40,20 @@
   endif
 
   [T, k] = size (X);
-  if !( is_vector (y) && (length (y) == T) )
+  if (! (is_vector (y) && (length (y) == T)))
     error ("f_test_regression: y must be a vector of length rows (X)");
   endif
   y = reshape (y, T, 1);
 
   [q, c_R ] = size (R);
   if (c_R != k)
-    error (["f_test_regression: ", ...
-            "R must have as many columns as X."]);
+    error ("f_test_regression: R must have as many columns as X");
   endif
 
   if (nargin == 4)
     s_r = size (r);
     if ((min (s_r) != 1) || (max (s_r) != q))
-      error (["f_test_regression: ", ...
-              "r must be a vector of length rows (R)."]);
+      error ("f_test_regression: r must be a vector of length rows (R)");
     endif
     r = reshape (r, q, 1);
   else
@@ -67,7 +65,7 @@
 
   [b, v] = ols (y, X);
   diff   = R * b - r;
-  f      = diff' * inv (R * inv (X' * X) * R') * diff / ( q * v );
+  f      = diff' * inv (R * inv (X' * X) * R') * diff / (q * v);
   pval  = 1 - f_cdf (f, df_num, df_den);
 
   if (nargout == 0)