MathGroup Archive 2010

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: restricting interpolating functions to be positive

On Jan 17, 5:11 am, Ray Koopman <koop... at> wrote:
> On Jan 15, 12:18 am, schochet123 <schochet... at> wrote:
> > On Jan 14, 12:46 pm, DrMajorBob <btre... at> wrote:
> >> For some data, that works pretty well; for other samples
> >> it has HUGE  peaks, reaching far above any of the data:
> > This problem can be avoided by using a transformation that
> > is close to the identity for positive values:
> > trans[x_] = (Sqrt[1 + x^2] + x)/2
> > inv[y_] = (4 y^2 - 1)/(4 y)
> There's a whole family of transformations here:
> trans[c_,x_] = (Sqrt[4c + x^2) + x)/2
> inv[c_,y_] = y - c/y
> Extremely small values of c give results that approach those obtained
> by chopping the interpolation at 0.

These transformations have the property that the interpolation behaves
differently for small and large numbers, with the transition at a
scale given by Sqrt[c]. The Log/Exp pair I suggested is scaleless.
{Sqrt[x],y^2} is another scaleless transformation.

All such transformations will exhibit artifacts from some point of
view given some data set: there is no flawless universal interpolation

The  {Sqrt[x],y^2} pair is useful for data that's proportional to
counts of independent events (like electric current over a potential
barrier), as it approximately gives each point its proper "weight" in
the analysis, given the correlation of the shot noise variance with
the expectation. However, the interpolation may "bounce" off the x-

  • Prev by Date: Re: Re: AstronomicalData - solar ecliptic longitude
  • Next by Date: Re: Testing Mathematica Expressions?
  • Previous by thread: Re: restricting interpolating functions to be positive
  • Next by thread: Re: restricting interpolating functions to be positive