Re: fitting, step

From: Maurice Leutenegger <Maurice.A.Leutenegger_at_email.domain.hidden>
Date: Wed, 7 Jul 2010 20:36:27 -0400
  John,

I think I have tracked down the problem. I had set the fit parameters to 
be very narrowly constrained:

set_par ("egauss(1).center", 3.0325, 0, 3.03, 3.035);
set_par ("egauss(1).sigma", 0.002, 0, 0.001, 0.003);

because I knew that the answer would be close to those numbers and 
didn't want to allow the fitting algorithm to waste time on other parts 
of parameter space.

When I relaxed the constraints:

set_par ("egauss(1).center", 3.0325, 0, 2.7, 3.3);
set_par ("egauss(1).sigma", 0.002, 0, 0.0, 0.01);

and then fit with mpfit / chisqr:sigma=gehrels, everything works 
beautifully.

So now my question is: suppose I really want to do a fit with the 
original parameter constraints. Do I need to use subplex? Is there a way 
to make mpfit work?

Thanks,

Maurice
----
You received this message because you are
subscribed to the isis-users list.
To unsubscribe, send a message to
isis-users-request_at_email.domain.hiddenwith the first line of the message as:
unsubscribe
Received on Wed Jul 07 2010 - 20:36:42 EDT

This archive was generated by hypermail 2.3.0 : Fri May 02 2014 - 08:35:46 EDT