library(maptools)
library(gstat)
points5 <- as.data.frame">as.data.frame(matrix(c(2,3,9,6,5,2,7,9,5,3,3,4,2,4,6),nrow=5,ncol=3,dimnames=list(c("p1","p2","p3","p4","p5"),c("X","Y","Z"))))
coordinates(points5) <- ~X+Y
points5.plt <- bubble(points5, "Z", fill=F, scales=list(draw=T))
print(points5.plt)
grid10 = expand.grid(x=seq(0,10,10), y=seq(0,10,10))
gridded(grid10) <- ~x+y
newpoint <- as.data.frame">as.data.frame(matrix(c(5,5),nrow=1,ncol=2,dimnames=list(c("x1"),c("X","Y"))))
coordinates(newpoint) <- ~X+Y
krige(Z~1, points5, newpoint, vgm(nugget=2.5,"Sph",psill=7.5,range=10))