help with "Tps" function to create a surface
Sorry this went to the wrong thread. I don't know. I can't reproduce your example, and a simple example on the meuse data set: library(sp) library(fields) data(meuse) Tps(data.frame(meuse$x,meuse$y), meuse$zinc) seemed to work. You could provide a reproducable example, or you could contact the author(s) of the fields package, I believe they're not listening to this list. -- Edzer
On 10/03/2010 05:51 PM, gianni lavaredo wrote:
Dear Researchers, I have a problem with Tps function. I tried several days and I am worked saturday night and I am working sunday, but I have always the same Error. I didn't find a solution in the web sorry. I have a "SpatialPointsDataFrame" with X, Y and H value
class( input.train)
[1] "SpatialPointsDataFrame"
attr(,"package")
[1] "sp"
# extract X and Y
x <- input.train at coords[,1]
y <- input.train at coords[,2]
xy <- data.frame(x,y)
# extract the value
v <- input.train at data[,c("H")]
Tps(xy,v)
Errore in rep(0, n1 * n2) : argomento 'times' non valido
Errore in out$W2 %d*% do.call(out$cov.function.name, c(out$args, list(x1 =
out$knots, :
error in evaluating the argument 'y' in selecting a method for function
'%d*%'
Errore in t(out$W2 %d*% do.call(out$cov.function.name, c(out$args, list(x1 =
out$knots, :
error in evaluating the argument 'x' in selecting a method for function
't'
Sorry and thanks for Help
Gianni
[[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Edzer Pebesma Institute for Geoinformatics (ifgi), University of M?nster Weseler Stra?e 253, 48151 M?nster, Germany. Phone: +49 251 8333081, Fax: +49 251 8339763 http://ifgi.uni-muenster.de http://www.52north.org/geostatistics e.pebesma at wwu.de