lognormal fitting
| pjfd | August 23, 2010 - 03:48 | ||
|---|---|---|---|
|
I have tried the email list, but i do not know what happens and it is not working for me. i come back to the forum. the problem: i have an xy data set that looks like a nice lognormal distribution. however, when i fit a get this ridiculously high chi-square value, which does not mean anything because i do not have a weighting wave -nor will i have one. the question: which is a general question for fitting - how can i evaluate goodness of fitting without a weighting wave? is there something like a kolmogorov test that i can do (or anything else)? how? i mean, the xy data set comes from two different waves, and cannot be compared to the fitted function, which is just a single wave. If i try to interpolate the xy set to a single wave, i get meaningless waves. what am i doing wrong? cheers P |
|||

Joined: 2007-06-29
Location: United States
John Weeks
WaveMetrics, Inc.
support@wavemetrics.com
Joined: 2007-06-21
Location: United States
But if the X values are not sorted, interpolate will give pretty weird results. Also, if you have NaN's (Not a Number, or a blank cell) in your data, you will get local patches within the interpolated data that are NaN's.
The Interpolate2 operation (Analysis->Interpolate) automatically sorts input data and removes NaNs before doing the interpolation.
I believe the cubic spline interpolator (Interpolate... at the bottom of the Analysis menu) will fail completely with NaN's under some cirumstances.
Only in the "X Coords From Dest Wave" mode when the destination is an XY pair and the X destination wave contains NaNs. This is very rare. The next release of Interpolate2 will tolerate the NaNs in this situation.
Joined: 2008-02-22
Location: United States
Then plot the Y wave of model values against the Y wave of data values, while is called a Q-Q plot (for Quantile-Quantile). The closer this line is to a straight line with slope=1, the better the fit. That gives you a visual assessment of fit. A Kolmorgov-Smirnov test (StatsKSTest) of the two waves (model values and data values) will assess this statistically.
Joined: 2010-08-20
Location: Spain
Hi guys,
thanks for the help. the newwave and kolmogorov-smirnov solved the problem. this forum is very helpful.
i am still troubled by the weighting function. how can i get one? actually, i do not even know what the residuals are.
the reduced chisqr looks ok, but how do a get the degrees of freedom? it is not going to be N-2 if my data set is over 1000 points.
P
Joined: 2007-06-29
Location: United States
the reduced chisqr looks ok, but how do a get the degrees of freedom? it is not going to be N-2 if my data set is over 1000 points.
John Weeks
WaveMetrics, Inc.
support@wavemetrics.com