4.102 sddssplinefit

• description: sddssplinefit fits splines to column data using the gsl library. The fits are of the form y = Σ_i{A[i] * B(x - xOffset,i)}, where B(x,i) is the ith basis spline function evaluated at x. sddssplinefit internally computes the A[i], writes the y values in the output file and estimates of the errors in the values. One can specify the order of the spline and the number of breakpoints (or alternatively the number of coefficients). The options are very similar to those for sddsmpfit.
• examples:
  sddssplinefit data.sdds fit.sdds -independent=x -dependent=y -order=3 -breakpoints=10  
  sddssplinefit input.sdds output.sdds -independent=t -dependent=z -coefficients=15 -normalize=1  
      

• synopsis:
  sddssplinefit [-pipe=[input][,output]] [inputfile] [outputfile]  
    -independent=xName -dependent=yname1-wildcard[,yname2-wildcard...]  
    [-sigmaIndependent=xSigma] [-sigmaDependent=ySigmaFormatString]  
    [-order=number] [-coefficients=number] [-breakpoints=number]  
    [-xOffset=value] [-xFactor=value]  
    [-sigmas=value,{absolute | fractional}]  
    [-modifySigmas] [-generateSigmas[={keepLargest, keepSmallest}]]  
    [-sparse=interval] [-range=lower,upper[,fitOnly]]  
    [-normalize[=termNumber]] [-verbose]  
    [-evaluate=filename[,begin=value][,end=value][,number=integer]]  
    [-fitLabelFormat=sprintf-string] [-infoFile=filename]  
    [-copyParameters]  
      

• files: inputFile is an SDDS file containing columns of data to be fit. If it contains multiple pages, they are processed separately. outputFile is an SDDS file containing one page for each page of inputFile. It contains columns of the independent and dependent variable data, plus columns for error bars (“sigmas”) as appropriate. The values of the fit and of the residuals are in columns named yNameFit and yNameResidual. In addition, various parameters having names beginning with yName are created that give reduced chi-squared, slope, intercept, and so on.
• switches:
  • -pipe[=input][,output] — The standard SDDS Toolkit pipe option.
  • -evaluate=filename[,begin=value][,end=value][,number=integer] — Specifies creation of an SDDS file called filename containing points from evaluation of the fit. The fit is normally evaluated over the range of the input data; this may be changed using the begin and end qualifiers, though the spline routines do not allow exceeding the range of the fit. Normally, the number of points at which the fit is evaluated is the number of points in the input data; this may be changed using the number qualifier.
  • -infoFile=filename — Specifies creation of an SDDS file containing results of the fits in columns. Under construction.
  • -order=number — Specifies the order of the spline, which is one more than the order of the local polynomials.
  • -breakpoints=number — The number of spline pieces the data will be split into for fitting.
  • -coefficients=number — The number of coefficients used in the spline fitting, which essentially controls the number of splines since the number of coefficients equals the number of splines minus two plus the order.
  • -xOffset=value, -xFactor=value — Specify offsetting and scaling of the independent data prior to fitting. The transformation is x → (x - xOffset)∕Factor. This feature can be used to make a fit about a point other than x = 0, or to scale the data to make high-order fits more accurate. The benefits for spline fitting are unknown, as it is left over from polynomial fitting.
  • -sparse=interval — Specifies sparsing of the input data prior to fitting. This can greatly speed computations when the number of data points is large.
  • -range=lower,upper — Specifies the range of independent variable over which to do fitting.
  • -verbose — Specifies that the results of the fit be printed to the standard error output.
• see also:
  • Data for Examples
  • sddspfit()
  • sddsmpfit()
  • sddsoutlier()
• author: M. Borland, ANL/APS.