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.tr...@coords[,1]
> y <- input.tr...@coords[,2]
> xy <- data.frame(x,y)
>
> # extract the value
> v <- input.tr...@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@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.pebe...@wwu.de
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
