comp.lang.idl-pvwave archive: archive » Re: constraining parameters in multi-Gaussian 1D fitting

Home » Public Forums » archive » Re: constraining parameters in multi-Gaussian 1D fitting

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Re: constraining parameters in multi-Gaussian 1D fitting [message #45360]

Mon, 05 September 2005 19:35

Craig Markwardt
Messages: 1869
Registered: November 1996

Senior Member

JD Smith <jdsmith@as.arizona.edu> writes:
> On Mon, 05 Sep 2005 12:07:42 -0500, Craig Markwardt wrote:
>
>>
>> "Jess" <jobrien@mso.anu.edu.au> writes:
>>> One constraint I am unable yet to do is: I = would like to be able to
>>> tie the peak flux of the Gaussians such that the peak flux of last
>>> Gaussian is always greater than that of the first Gaussian.
>>> I tried using
>>> parinfo((n_gauss-1)*3).tied = 'GT P[0]'
>>>
>>> However the tied structure of parinfo doesn't seem to be meant to
>>> accept operators like GT,LT, etc. ...
>>
>> True. MPFIT's TIED fields are limited to equality constraints only.
>
> What if you availed yourself of the ITERPROC procedure to enforce the
> constraint, dragging the fitter (kicking and screaming if necessary)
> back into line if it attempts to step out? Any reason this wouldn't
> work?

It might work, it might not. I suspect that in general the fitter
might get stuck. For example, if we *started* out by thinking the
tallest peak was on the left part of the curve - and made a model
function to match that - but the truth is that a peak on the right is
truly he tallest, then we can enforce all the constraints in the world
and won't come out with the best fit.

Now, with the additional info that the original poster provided, this
will probably not be the case, so s/he might be alright.

Craig

--
------------------------------------------------------------ --------------
Craig B. Markwardt, Ph.D. EMAIL: craigmnet@REMOVEcow.physics.wisc.edu
Astrophysics, IDL, Finance, Derivatives | Remove "net" for better response
------------------------------------------------------------ --------------

Report message to a moderator

[Message index]

		Re: constraining parameters in multi-Gaussian 1D fitting By: Craig Markwardt on Mon, 05 September 2005 19:35
		Re: constraining parameters in multi-Gaussian 1D fitting By: Craig Markwardt on Mon, 05 September 2005 18:58
		Re: constraining parameters in multi-Gaussian 1D fitting By: JD Smith on Mon, 05 September 2005 17:33
		Re: constraining parameters in multi-Gaussian 1D fitting By: Jess on Mon, 05 September 2005 17:02
		Re: constraining parameters in multi-Gaussian 1D fitting By: Craig Markwardt on Mon, 05 September 2005 10:07

Previous Topic:	Concatenating arrays across chosen dimension
Next Topic:	Another IDL 6.2 Project Problem

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Fri Oct 10 18:08:12 PDT 2025

Total time taken to generate the page: 1.52029 seconds