James Sethna > Sloppy Models > Fitting Polynomials: Where is sloppiness from?

Cube, fourth, and fifth power of x plotted
Sloppiness comes from parameters which affect the predictions in similar ways. Here the monomials x3, x4, and x5 all are flat near zero, and increase smoothly near one. We can trade the coefficient a4 for a suitable combination of a3 and a5, and get almost exactly the same final function.

Fitting Polynomials: Where is sloppiness from?

Multiparameter models are usually sloppy - their parameters can vary over huge ranges without changing the fits to data very much. Where does this sloppiness come from? We provide three answers to this question, all really the same answer in disguise.

Consider the standard problem of fitting polynomials to data. Our model is a sum of monomials xn, with coefficients an:

y(a0, a1, a2, ... a5) = a0 + a1 x + a2 x2 + a3 x3 + a4 x4 + a5 x5
If we view a0 ... a5 as the six parameters of our model, it is extremely sloppy (even for polynomials with a rather small number of terms). At right, you can see that the monomials xn all have roughly the same shape, and can be traded for one another without changing the fits.

1. Sloppiness comes from having groups of parameters whose effects can be traded for one another.

Indeed, if we are fitting a continuous function in the range between zero and one, we can explicitly calculate the Hessian matrix whose range of eigenvalues we use to decide whether the system is sloppy. This Hessian is the Hilbert matrix:

Hij = 1/(i+j+1) =
11/21/31/41/51/6
1/21/31/41/51/61/7
1/31/41/51/61/71/8
1/41/51/61/71/81/9
1/51/61/71/81/91/10
1/61/71/81/91/10 1/11
a matrix famous for being ill-conditioned (having a huge ratio between its largest and smallest eigenvalue).
Wiggly Legendre polynomials plotted
Changing the parameters can remove the sloppiness. When fitting polynomials to data, one is told to vary not the monomial coefficients an, but the coefficients bn of a suitable class of orthogonal polynomials: here the shifted Legendre polynomials Ln. Here are plotted L3, L4, and L5. By design, these polynomials can't be traded for one another. The same polynomial fit to the same data, expressed in terms of new variables bn, is not sloppy.

Those readers who actually have fit polynomials to data will remember, though, that one doesn't usually calculate the monomial coefficients an. Instead, one fits to a sum of orthogonal polynomials

Y(..., b3, b4, b5) = ... + b3 L3(x) + b4 L4(x) + b5 L5(x)
(figure at left), designed so that every polynomial is a different "shape". This means that if we use the coefficients b0 ... b5 as the parameters in our model, it should not be sloppy. Indeed, the Hessian for the bn is the identity matrix, and all the eigenvalues are equal to one. Notice that the space of functions is the same (polynomials up to fifth order), the best fit is the same, the ensemble of acceptable fits is the same - all we have changed is the way we've labeled the functions in our model.

2. Sloppiness comes from the way we choose the variables in our theory.

Notice that in the case of fitting polynomials, we often don't think that the monomial coefficients are very important: the b variables are almost as natural as the a variables. But in a biological model or an engineering system, the variables we use are usually natural and sensible, and using strange combinations of variables just to avoid sloppiness will usually not make scientific sense.
Three red arrows, pointing almost in the same direction
Skewness. The change of variables from the bare parameters (monomial coefficients an) to parameters natural to the model (coefficients of the orthogonal polynomials bn is not a rotation in parameter space. It is more of a tortured, skewed stretching. The red arrows represent the change in natural coordinates when different bare coordinates are varied. Notice that one can get to nearly the same final behavior (tip of the skewed cube) by going along any of the three bare coordinates (they all affect the predictions in similar ways). Notice also that the volume of the skewed cube is small; it is proportional to the determinant of the mapping from bare to skewed coordinates.

Why did changing from the monomial coefficients an to the orthogonal polynomial coefficients bn make such a big difference? Unlike many familiar transformations, this change is not just a rotation (in six-dimensional function space); it is instead a very skewed, shearing kind of transformation (figure at left). The meaning of "perpendicular" in parameter space is subtle: most scientifically natural choices for parameters are sloppy for this reason.
3. Sloppiness comes a severe skewness between the scientifically natural "bare" parameters and the parameters naturally governing the model behavior.

Can we understand the key properties of sloppy models mathematically? See Why sloppiness? The Sloppy Universality Class

References

More on sloppiness:



Last Modified: June 11, 2008

James P. Sethna, sethna@lassp.cornell.edu; This work supported by the Division of Materials Research of the U.S. National Science Foundation, through grant DMR-070167.

Statistical Mechanics: Entropy, Order Parameters, and Complexity, now available at Oxford University Press (USA, Europe).