13.07.2015 Views

The gstat Package - NexTag Supports Open Source Initiatives

The gstat Package - NexTag Supports Open Source Initiatives

The gstat Package - NexTag Supports Open Source Initiatives

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

16 juraNoteRock see book and belowCd see bookCo see bookCr see bookCu see bookNi see bookPb see bookZn see book<strong>The</strong> points data sets were obtained from http://home.comcast.net/~goovaerts/book.html, the grid data were kindly provided by Pierre Goovaerts.Rock Types: 1: Argovian, 2: Kimmeridgian, 3: Sequanian, 4: Portlandian, 5: Quaternary.Land uses: 1: Forest, 2: Pasture (Weide(land), Wiese, Grasland), 3: Meadow (Wiese, Flur, Matte,Anger), 4: Tillage (Ackerland, bestelltes Land)Points 22 and 100 in the validation set (validation.dat[c(22,100),]) seem not to lieexactly on the grid origininally intended, but are kept as such to be consistent with the book.Author(s)Data preparation by David Rossiter (rossiter@itc.nl) and Edzer PebesmaReferencesGoovaerts, P. 1997. Geostatistics for Natural Resources Evaluation. Oxford Univ. Press, New-York,483 p. Appendix C describes (and gives) the Jura data set.Atteia, O., Dubois, J.-P., Webster, R., 1994, Geostatistical analysis of soil contamination in theSwiss Jura: Environmental Pollution 86, 315-327Webster, R., Atteia, O., Dubois, J.-P., 1994, Coregionalization of trace metals in the soil in the SwissJura: European Journal of Soil Science 45, 205-218Examplesdata(jura)summary(prediction.dat)summary(validation.dat)summary(transect.dat)summary(juragrid.dat)# the commands to create the spatial objects:require(sp)jura.pred = prediction.datjura.val = validation.datjura.grid = juragrid.datjura.pred$Landuse = factor(prediction.dat$Landuse, labels=levels(juragrid.dat$Landuse))

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!