normalisation of design

fsl_glm.cc

@@ -51,6 +51,9 @@ using namespace std;
 	Option<int> dofset(string("--dof"),0,
 		string("        set degrees-of-freedom explicitly"),
 		false, requires_argument);
+	Option<bool> normdes(string("--des_norm"),FALSE,
+		string("        switch on normalisation of the design matrix columns to unit std. deviation"),
+		false, no_argument);
 	Option<bool> perfvn(string("--vn"),FALSE,
 		string("        perfrom MELODIC variance-normalisation on data"),
 		false, no_argument);
@@ -179,6 +182,9 @@ int setup(){
 	meanR=mean(data,1);
 	data = remmean(data,1);
 	design = remmean(design,1);
+	
+	if(normdes.value())
+		design =  SP(design,ones(design.Nrows(),1)*pow(stdev(design,1),-1));
 
 	if(fncontrasts.value()>""){//read contrast		
 		contrasts = read_ascii_matrix(fncontrasts.value());
@@ -244,6 +250,7 @@ int main(int argc,char *argv[]){
 			options.add(fnmask);
 			options.add(fnftest);
 			options.add(dofset);
+			options.add(normdes);
 			options.add(perfvn);
 			options.add(help);
 			options.add(outcope);