[Homer-users] (no subject)

[Huppert, Ted]

In matlab, you can always set a debug break point and stop the program on the line to investigate the issue.    Its really hard to diagnose problems without your exact dataset and the exact processing pipeline.

In matlab type:
"dbstop in hmrDeconvTB_SS3rd at 71”

This will set a debug break right at the line you are having troubles with.  You can then try to look at the variables and diagnose the problem.  Alternatively, you can use the command “dbstop if error”  which will stop anywhere in the program that an error occurred.

My guess from looking at that function is that your measurement list variable in the NIRS file is not sorted by wavelength (e.g. The two wavelength indexes in the 4th column are mixed).  E.g. The line "lst = find(ml(:,4)==1);” finds entries that are in the second 1/2 of the ml variable.

David-
The quick solution is to change lines 58,59 and 73 to read (e.g.) “rhoSD( lst(iML)) =  “ instead of "rhoSD(iML) = "



----------------
Theodore Huppert, PhD
Associate Professor
Departments of Radiology and Bioengineering
University of Pittsburgh
Email: huppertt at upmc.edu
Phone: (412) 647-8459



"Insanity: doing the same thing over and over and expecting different results"-  Einstein

From: <Chiu>, "Bryan (PHTH)" <bryan.chiu at ubc.ca<mailto:bryan.chiu at ubc.ca>>
Reply-To: "homer-users at nmr.mgh.harvard.edu<mailto:homer-users at nmr.mgh.harvard.edu>" <homer-users at nmr.mgh.harvard.edu<mailto:homer-users at nmr.mgh.harvard.edu>>
Date: Tuesday, August 19, 2014 at 2:04 PM
To: "homer-users at nmr.mgh.harvard.edu<mailto:homer-users at nmr.mgh.harvard.edu>" <homer-users at nmr.mgh.harvard.edu<mailto:homer-users at nmr.mgh.harvard.edu>>
Subject: Re: [Homer-users] (no subject)

hmrDeconvTB_SS3rd