[R-sig-Geo] how to fit 3d variogram in gstat

[Edzer Pebesma]

On 28/11/16 10:18, maogui.hu wrote:
> Dear list,
>     gstat has a powerful function in 2d and spatial-temporal variogram fitting. I'm wondering how can i fit a 3d automatically? Here is a simple snippet to create a 3d dataset:
> 

It doesn't; either the problem is not simple enough or we don't
understand it well enough in order to let you do that.

> library(sp) 
> library(gstat) 
> 
> data(DE_RB_2005, package = "gstat") 
> st <- as.data.frame(DE_RB_2005[,1:100]) 
> st <- st[,c(1,2,6,16)] 
> names(st) <- c("x", "y", "z", "v") 
> 
> st$x <- (st$x - mean(st$x)) / sd(st$x) * 100 
> st$y <- (st$y - mean(st$y)) / sd(st$y) * 100 
> st$z <- (st$z - mean(st$z)) / sd(st$z) * 100 
> 
> coordinates(st) <- ~x+y+z 
> v <- variogram(v~1, st)  
> print(v)
> 
> 
> Thanks very much!
>  
> 
> 
> Maogui Hu, Ph.D.
> State Key Laboratory of Resources and Environmental Information System
> Institute of Geographic Sciences and Natural Resources Research, CAS
> 
> 	[[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
Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83 33081
Journal of Statistical Software:   http://www.jstatsoft.org/
Computers & Geosciences:   http://elsevier.com/locate/cageo/

-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 473 bytes
Desc: OpenPGP digital signature
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20161128/e4481cbc/attachment.bin>