Hi all. This question relates to the previous one I had about plotting
kriging variance and function of grid spacing. I found what I think is
an example in "gstat's User Manual". This example is given below...
data(meuse) #HERE I WOULD SUBSTITUTE MY DATA "BRI2000"
coordinates(meuse) = ~x+y
data(meuse.grid) #HERE I WOULD SUBSTITUTE MY GRID WITH RESOLUTION 100 X 143 "SD100"
gridded(meuse.grid) = ~x+y
m<- vgm(.59, "Sph", 874, .04)
# ordinary kriging:
x<- krige(log(zinc)~1, meuse, meuse.grid, model = m)
spplot(x["var1.pred"], main = "ordinary kriging predictions")
spplot(x["var1.var"], main = "ordinary kriging variance")
My questions: 1)How do I convert a spatial points data frames, BRI2000
and SD100 to an S4 data frame that can be used with the above code?
2) This seems where I get my kriging variance as a function of grid
spacing ( ie. SUBSTITUTE BRI3000, BRI2000, BRI1000,BRI500, BRI250), correct?
Thanks,
Kerry
--
**********************
Kerry Ritter, Ph.D.
statistician
Southern California Coastal Water Research Project
3535 Harbor Blvd., Suite 110
work: 714-755-3210
cell: 714-420-3346
fax: 714-755-3299
email: kerryr@sccwrp.org
[[alternative HTML version deleted]]