From 5b0718ea69a590cbe366ab7aa23a5c33c1dde81c Mon Sep 17 00:00:00 2001
From: Mark Woolrich <woolrich@fmrib.ox.ac.uk>
Date: Wed, 7 Jul 2004 16:22:20 +0000
Subject: [PATCH] *** empty log message ***
---
glimGls.cc | 6 ++++--
1 file changed, 4 insertions(+), 2 deletions(-)
diff --git a/glimGls.cc b/glimGls.cc
index 07b0c75..8971889 100644
--- a/glimGls.cc
+++ b/glimGls.cc
@@ -36,10 +36,12 @@ namespace FILM {
// compute r
Matrix I(sizeTS, sizeTS);
Identity(I);
- r = (I-x*inv_xx*x.t())*y;
+ Matrix R = I-x*inv_xx*x.t();
+ r = R*y;
// compute sigma squareds
- sigmaSquareds(ind) = (r.t()*r/sizeTS).AsScalar();
+ // sigmaSquareds(ind) = (r.t()*r/sizeTS).AsScalar();
+ sigmaSquareds(ind) = (r.t()*r).AsScalar()/R.Trace();
// set corrections
SetCorrection(inv_xx, ind);
--
GitLab