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
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
