comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Re: How to determine the WEIGHTS in MPFITFUN for distribution function fit?
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: How to determine the WEIGHTS in MPFITFUN for distribution function fit? [message #55764] Thu, 06 September 2007 22:26 Go to previous message
duxiyu@gmail.com is currently offline  duxiyu@gmail.com
Messages: 88
Registered: March 2007
Member
On 9 7 , 2 26 , Craig Markwardt
<craigm...@REMOVEcow.physics.wisc.edu> wrote:
>
> If Y is measured counts, then the YERR should be ~ SQRT(Y), right?
>
> But overall I agree, the original poster needs to specify how the
> numbers were measured in order to estimate the measurement error.
>
> Craig
>

The effect of Yerr=SQRT(Y) is the same as the effect of WEIGHTS=1/y
for the procedure MPFITFUN, isn't it?

The array named D in IDL means the original mesured data.
Y = histogram ( D, locations = X )
I want use a specified function to fit this curve in the x-y plane.
The computed function is the distribution function for the measured
data D, isn't it?
If there are not the measurement errors for the data D, I should set
the Yerr = SQRT(Y) when I use the MPFITFUN to fit the curve.
If there ara the measurement errors Derr which have the same
dimensions as the data D, how do I caculate the Yerr?

The another question is how to determine the START_PARAMS in the
procedure MPFITFUN?
When I use the matlab to fit the curve, the procedure can select the
start_value automatically.
The START_PARAMS in <PFITFUN must be given by user.
The START_PARAMS is set to some values, and the procedure can return
the results without any error.
But when the START_PARAMS change a little, the procedurel gives the
warning "Program caused arithmetic error: Floating illegal operand"
and return the different results.
I do not know the criterion of its selection.

I am looking forward to your reply.

Best regards,
jdu
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Call IDL (Virtual Machine) Save file from Java
Next Topic: Random problem with Delaunay triangulation

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Fri Oct 10 23:51:25 PDT 2025

Total time taken to generate the page: 0.00278 seconds