Hi, I am trying to follow the examples in  Bivand et al 'Applied Spatial
Data Analysis' and the examples at
http://casoilresource.lawr.ucdavis.edu/drupal/node/442 but I am frustrated.
I finally managed an ordinary kriging example to work but with little
understanding of why/when it should work.
I posted some data: https://www.dropbox.com/s/rtg7gh67ziai3ec/mar2011.RData
Another item I've noticed, if I  define the projection with
projection(object)='+proj=lonlat +datum=NAD83'  for any of my spatial
objects it won't work.  It only works with an defined projection. any help
would be appreciated. Thanks.

load('mar2011.RData')

vt=variogram(snotel~recon,data=mar2011)

vt.fit=fit.variogram(vt,vgm(psill=0.025,nugget=0.015,model='Sph',range=2))

#plot(vt,vt.fit)

ok<-krige(snotel~1,locations=mar2011,newdata=swe.grid,model=vt.fit)

uk<-krige(snotel~mar2011$recon,locations=mar2011,newdata=swe.grid,model=
vt.fit)

#spplot(ok)
Dominik

	[[alternative HTML version deleted]]

