variog: Compute Empirical Variograms and Permutation Envelopes

In geostatsp: Geostatistical Modelling with Likelihood and Bayes

Compute Empirical Variograms and Permutation Envelopes

Description

These are wrappers for \Sexpr[results=rd]{c( '\\\code{variog} in the \\\code{geoR} package', '\\\command{\\\\link[geoR]{variog}}' )[1+requireNamespace('geoR', quietly=TRUE)]} and \Sexpr[results=rd]{c( '\\\code{variog.mc.env} in the \\\code{geoR} package', '\\\command{\\\\link[geoR]{variog.mc.env}}' )[1+requireNamespace('geoR', quietly=TRUE)]}.

Usage

variog(geodata, ...)
## S3 method for class 'SpatVector'
variog(geodata, formula, ...)
## Default S3 method:
variogMcEnv(geodata, ...)
## S3 method for class 'SpatVector'
variogMcEnv(geodata, formula, ...)

Arguments

geodata	An object of class SpatVector or of a class suitable for \Sexpr[results=rd]{c( '\\\code{variog} in the \\\code{geoR} package', '\\\command{\\\\link[geoR]{variog}}' )[1+requireNamespace('geoR', quietly=TRUE)]}
formula	A formula specifying the response variable and fixed effects portion of the model. The variogram is performed on the residuals.
...	additional arguments passed to \Sexpr[results=rd]{c( '\\\code{variog} in the \\\code{geoR} package', '\\\command{\\\\link[geoR]{variog}}' )[1+requireNamespace('geoR', quietly=TRUE)]}

Value

As \Sexpr[results=rd]{c( '\\\code{variog} in the \\\code{geoR} package', '\\\command{\\\\link[geoR]{variog}}' )[1+requireNamespace('geoR', quietly=TRUE)]} and \Sexpr[results=rd]{c( '\\\code{variog.mc.env} in the \\\code{geoR} package', '\\\command{\\\\link[geoR]{variog.mc.env}}' )[1+requireNamespace('geoR', quietly=TRUE)]}

See Also

\Sexpr[results=rd]{c( '\\\code{variog} in the \\\code{geoR} package', '\\\command{\\\\link[geoR]{variog}}' )[1+requireNamespace('geoR', quietly=TRUE)] }

and \Sexpr[results=rd]{c( '\\\code{variog.mc.env} in the \\\code{geoR} package', '\\\command{\\\\link[geoR]{variog.mc.env}}' )[1+requireNamespace('geoR', quietly=TRUE)]}

Examples

data("swissRain")
swissRain = unwrap(swissRain)
swissRain$lograin = log(swissRain$rain)
swissv= variog(swissRain, formula=lograin ~ 1,option="bin")
swissEnv = variogMcEnv(swissRain, lograin ~ 1, obj.var=swissv,nsim=9) 
if(!is.null(swissv)){
	plot(swissv, env=swissEnv, main = "Swiss variogram")
}

geostatsp documentation built on Dec. 24, 2024, 3 a.m.