[R-sig-Geo] units for spatial range in stvariogram gstat

sara osama sarosama at feps.edu.eg
Wed May 17 20:56:02 CEST 2017


Dear All,
I have question concerning the unit of measurement for the range parameter
in the fit.stvariogram in gstat, i have fitted the following spatio
temporal variogram for data on monthly averages pf pm10 and the longitude
and latitude are in decimel format

> pm10_sp
SpatialPoints:
     longitude latitude
1     30.92500 29.95300
61    31.36000 30.26000
121   31.33300 29.84720
 and this is the commands for fitting the variogram
> varpm10 <-
variogramST(lnpm10only~1,data=lnpm10_stfdf,cutoff=0.8,assumeRegular=T,na.omit=T)
> plot(varpm10,map=F)
> sumMetric <- vgmST("sumMetric", space = vgm(psill=5,"Mat", range=50,
nugget=0),time = vgm(psill=500,"Mat", range=50, nugget=0,kappa=0.3), joint
= vgm(psill=10,"Mat", range=100, nugget=1), stAni=200)
> sumMetric_Vgm <- fit.StVariogram(varpm10, sumMetric,
method="L-BFGS-B",lower=pars.l,upper=pars.u)
> plot(varpm10,sumMetric_Vgm,map=F)

> extractPar(sumMetric_Vgm)
      sill.s      range.s     nugget.s       sill.t      range.t
nugget.t
  1.73598444  50.63845035   0.03212857   0.01570379  49.99811448
0.00000000
     sill.st     range.st    nugget.st         anis
  0.00000000 100.00004052   0.09151787 200.00000000

-- 
now i want to know the units for these parameters in order to be able to
interpret it for example is it in km or m.

I really appreciate any help that I deeply need

Best Regards
Sara

Faculty of Economics & Political Science
Cairo University

Tel:(202)35728055-(202)35728116-(202)35736608-(202)35736605
Fax:(202)35711020
Follow us on twitter:https://twitter.com/fepsnews
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20170517/35f109d4/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: variogram.png
Type: image/png
Size: 13032 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20170517/35f109d4/attachment.png>


More information about the R-sig-Geo mailing list