[R-sig-Geo] Trouble with universal kriging

[Edzer Pebesma]

Dominik, if you want to receive help, please explain why you are
frustrated, and instead of claims of "won't work", explain in detail
what you did and what happened, preferably with a well-prepared example
that people on the list can reproduce without downloading data first.


On 06/20/2013 09:57 PM, Dominik Schneider wrote:
> 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]]
> 
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> 

-- 
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster
Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251
83 33081, Fax: +49 251 8339763 http://ifgi.uni-muenster.de