[R] interp function in package akima

[David Winsemius]

On Feb 8, 2010, at 5:06 PM, Sacha Viquerat wrote:

> hello! does anyone know how to get values out of an interp object  
> (from akima package)? ive constructed an elevation map and would  
> like to numerically see the calculated z values of data points,  
> which were not empirically sampled (basically: id like to see the  
> altitude at locations x,y which have been computed by the interp  
> function. is there any way? thanks in advance!

Have you read the help page for interp???  The interp function returns  
a list, the third element of which is a matrix named "z":

interpFn <- interp( < arguments > )

interpFn$z  # should return that matrix

If your pick out a particular element associated with the 10th x-value  
and the 20th y-value, you should be able get the x,y,z triplet by

c(interpFn$x[10] , interpFn$y[20] ,interpFn$z[10,20] )

If you have particular x,y values you want to access then you need to  
include them in the xo[] and yo[] vectors supplied to interp().  
interp() does not provide interpolation between points on the output  
grid, at least as as I read the help page.
There are strategies that could do so. Whenever I estimate a crossed- 
spline model in package Design, I could extract a 2d-spline function  
if I chose to.

-- 

David Winsemius, MD
Heritage Laboratories
West Hartford, CT