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]]

