[R-sig-Geo] help with "Tps" function to create a surface

[Edzer Pebesma]

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