[R] Variogram fitted by Cauchy

damien landais damien.landais at tdf.fr
Thu Jun 11 09:49:38 CEST 2009


Hi
I'm using R(2.9.0)  and gstat package under Windows  to  plot sample variograms.
When I want to fit them, I don't have the possibility to choose Cauchy functions. Under Mathematica for example, we have it and it seems to be the best one to fit the variograms I have
I tested the others functions, they can fit but I would,if possible, use the Cauchy one. If anybody knows something about it, I would know, even with another package

To show what I want to do, here is my code

mydata=read.table("97MHz_15cm_V3.txt")
colnames(mydata)=c("x","y","z","E","taille")
g=gstat(id="zinc",formula=mydata[,4]~1,locations=~x+y+z,data=mydata)
g.var=variogram(g,cutoff=2500,width=100)

modele=vgm(model="Exp",range=2000)
vario=fit.variogram(g.var,modele)


Thanks!

Cordialement
Damien Landais




More information about the R-help mailing list