maximizeInterpolant: Maximize a function given a table of values by spline interpolation.
Description
Maximize a function given a table of values by spline interpolation.
Usage
maximizeInterpolant(x, y)
Arguments
x
numeric vector of the inputs of the function.
y
numeric matrix of function values at the values of x.
Columns correspond to x values and each row corresponds to a different function to be maximized.
Value
numeric vector of input values at which the function maximums occur.
Details
Calculates the cubic spline interpolant for each row the method of Forsythe et al (1977) using the function fmm_spline from splines.c in the stats package).
Then calculates the derivatives of the spline segments adjacant to the input with the maximum function value.
This allows identification of the maximum of the interpolating spline.
References
Forsythe, G. E., Malcolm, M. A. and Moler, C. B. (1977).
Computer Methods for Mathematical Computations, Prentice-Hall.