A function to compute the second derivative of a function (of several real variables) using a quadratic approximation on a grid of points defined by the list argRanges. Also returns the local maximum.
QuadApprox(fun, npts, argRanges, plot = FALSE, ...)a 2 by 2 matrix containing the curvature at the maximum and the (x,y) value at which the maximum occurs
a function
integer number of points in each direction
a list of ranges on which to construct the grid for each parameter
whether to plot the quadratic approximation of the posterior (for two-dimensional parameters only)
other arguments to be passed to fun