Grokbase Groups R r-help April 2002
FAQ
Hi !

I'm performing kriging on a squared grid 12cm x 12cm. Each cell of the
grid is 0.5cm x 0.5cm.

I'm using sgeostat package "krige" funtion. When I plot the kriging
predictions the image occupies only the upper corner of my grid. I don't
understand my error.

The code of my script is

grid <- list(x=seq(-6,6,by=0.5),y=seq(-6,6,by=0.5))

grid$xr <- range(grid$x)
grid$xs <- grid$xr[2] - grid$xr[1]
grid$yr <- range(grid$y)
grid$ys <- grid$yr[2] - grid$yr[1]
grid$max <- max(grid$xs, grid$ys)
grid$xy <- data.frame(cbind(c(matrix(grid$x, length(grid$x),
length(grid$y))), c(matrix(grid$y, length(grid$x), length(grid$y),
byrow=T))))

colnames(grid$xy) <- c("x", "y")
grid$point <- point(grid$xy)

grid$krige <-krige(grid$point, dbpoint, variomd, at="Counts")

Thanks,
Paola.

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

Search Discussions

Related Discussions

Discussion Navigation
viewthread | post
Discussion Overview
groupr-help @
categoriesr
postedApr 5, '02 at 9:02a
activeApr 5, '02 at 9:02a
posts1
users1
websiter-project.org
irc#r

1 user in discussion

Lecca Paola: 1 post

People

Translate

site design / logo © 2022 Grokbase