Corrected dof for demeaning

glimGls.cc

@@ -22,7 +22,7 @@ using namespace fslsurface_name;
 
 namespace FILM {
 
-  GlimGls::GlimGls(const int pnumTS, const int psizeTS, const int pnumParams, const int pnContrasts, const int pnFContrasts) : 
+  GlimGls::GlimGls(const int pnumTS, const int psizeTS, const int pnumParams, const int pnContrasts, const int pnFContrasts, const int dofAdjustment) : 
     numTS(pnumTS),
     sizeTS(psizeTS),
     nContrasts(pnContrasts),
@@ -33,7 +33,7 @@ namespace FILM {
     varcopes(nContrasts, numTS),  
     fstats(nFContrasts,numTS),
     sigmaSquareds(numTS),
-    dof(sizeTS - numParams)
+    dof(sizeTS - numParams + dofAdjustment)
     {
       I=IdentityMatrix(sizeTS);
       fDof.resize(nFContrasts,0);

glimGls.h

@@ -26,7 +26,7 @@ namespace FILM {
   class GlimGls
     {
     public:
-      GlimGls(const int pnumTS, const int psizeTS, const int pnumParams, const int pnContrasts, const int pnFContrasts);
+      GlimGls(const int pnumTS, const int psizeTS, const int pnumParams, const int pnContrasts, const int pnFContrasts, const int dofAdjustment);
      
       void setData(const ColumnVector& p_y, const Matrix& p_x, const int ind, const Matrix& tContrasts, const Matrix& fContrasts);       
       void Save(const NEWIMAGE::volume<float>& mask, NEWIMAGE::volume4D<float>& saveVolume,fslsurface_name::fslSurface<float, unsigned int>& saveSurface,const string& saveMode,const float reftdim=1.0);