https://github.com/cran/RandomFields
Raw File
Tip revision: f082dc8b0950aff830aab568d89a74af74f10e14 authored by Martin Schlather on 12 August 2014, 00:00:00 UTC
version 3.0.35
Tip revision: f082dc8
RMangle.Rd
\name{RMangle}
\alias{RMangle}
\title{Anisotropy matrix given by angle}
\description{
  \command{RMangle} delivers an anisotropy matrix for the
  argument \code{Aniso} in \command{\link{RMmodel}} in two dimensions.
  \command{RMangle} requires one two strechting values, passed by
  \code{ratio} or \code{diag}, and an \command{angle}.
  For \code{angle} \eqn{a} and code{diag} equal to \eqn{(d_1, d_2)} the
  anisotropy matrix \eqn{A} is
  
  \code{A = diag(d_1, d_2) \%*\% matrix(ncol=2, c(cos(a), sin(a),
    -sin(a), cos(a)))}
}
 

\usage{
RMangle(angle, ratio, diag) 
}

\arguments{
 \item{angle}{angle \code{a}}
 \item{ratio}{equivalent to \code{diag=c(1, 1/ratio)}}
 \item{diag}{the diagonal components of the matrix}
}

\value{
 \command{\link{RMangle}} returns an object of class \code{\link[=RMmodel-class]{RMmodel}}.
}

\author{Martin Schlather, \email{schlather@math.uni-mannheim.de}
}
\seealso{
 \command{\link{RMtrafo}}, \command{\link{RMmodel}}
  }

\examples{

RFoptions(seed=0) ## *ANY* simulation will have the random seed 0; set
##                   RFoptions(seed=NA) to make them all random again

model <- RMexp(Aniso=RMangle(angle=pi/4, ratio=3))
plot(model, dim=2)

x <- seq(0, 2, if (interactive()) 0.05 else 1)
z <- RFsimulate(x, x, model=model)
plot(z)


\dontshow{FinalizeExample()}

}


\keyword{spatial}
\keyword{models}




back to top