[R] possible overflow/underflow error causing NaN in .trval when invoking trmat with degree 4

2013-12-15 Thread Nerd of Darkness
Hello all, I have two copies of the same data; I want to compare trend matrices for degrees 3 and 4, so I name the copies "triplo" and "quatro". > triploSurface<-surf.ls(3,triplo) > triploMat<-trmat(triploSurface,-150.,150.,-150.,150., 300) [This produces the expected triploSurface

Re: [R] What is the easiest way to interpolate vertical values on a square section of a nearly-planar 3D surface

2013-12-03 Thread Nerd of Darkness
On 2013年12月04日 08:21, David Winsemius wrote: > library(car) > # will also need rgl > scatter3d(dat$X, dat$Y, dat$Value) > > library(akima) > akima.li <- interp(dat$X, dat$Y, dat$Value, >xo=seq(min(dat$X), max(dat$X), length = 100), >yo=seq(min(dat$Y), max