		load lowfinput   
		X=lowf;
		%
		% Options for regem
		%
		OPTIONS.regress = 'ttls';
		OPTIONS.regpar = 3
		OPTIONS.stagtol = 1e-4;
		OPTIONS.maxit = 400;
		OPTIONS.inflation = 1; % this works pretty well all around, but it may vary based on your application if
									the calibration period is short and the network sparse. The only 
		                            way to set this is by cross validation. One could do a short calibration (1900-1980 in our case)
		                            to set this. For Experiment "b" 1.09, for "j" 1.5 and for "k" 1.7
		
		
		OPTIONS.disp = 1;
		%
		[X, M, C, Xerr] = regemlow(X, OPTIONS);
		save lowf.filled
		
