https://github.com/cran/RandomFields
Tip revision: bd298816a60ec4ca975f1289dc6ad3475d6247bf authored by Martin Schlather on 09 January 2016, 13:56:44 UTC
version 3.1.8
version 3.1.8
Tip revision: bd29881
RRdistr.Rd
\name{RRdistr}
\alias{RRdistr}
\title{RRdistr}
\description{
\command{RRdistr} defines of distribution family given by \code{fct}.
It is used to introduce \link[=RR]{Random parameters} based on
distributions defined on R.
}
\usage{
RRdistr(fct, nrow, ncol, envir)
}
\arguments{
\item{fct}{an arbitrary family of distribution. E.g.
\code{norm()} for the family \command{dnorm}, \command{pnorm},
\command{qnorm}, \command{rnorm}.
}
\item{nrow, ncol}{The matrix size (or vector if \code{ncol=1})
the family returns. Except for very advanced modelling we always have
\code{nrow=ncol=1}, which is the default.}
\item{envir}{an environment; defaults to \code{\link[base]{new.env}()}.
}
}
\note{
\command{\link{RRdistr}} is the generic model introduced
automatically when distribution families in R are used in the model
definition. See the examples below.
}
\details{
\command{\link{RRdistr}} returns an object of class \code{\link[=RMmodel-class]{RMmodel}}.
}
\note{
See \link{Bayesian Modelling} for some less technical introduction to
hierarchical modelling.
The use of \command{RRdistr} is completely on the risk of the user. There is no
way to check whether the expressions of the user are mathematically
correct.
Further, \command{\link{RRdistr}} may not be used in connection of obsolete
commands of RandomFields.
}
\author{Martin Schlather, \email{schlather@math.uni-mannheim.de} \url{http://ms.math.uni-mannheim.de/de/publications/software}
}
\seealso{
\command{\link{RMmodel}},
\link{RR},
\command{\link{RFsimulate}},
\command{\link{RFdistr}}
}
\examples{
RFoptions(seed=0) ## *ANY* simulation will have the random seed 0; set
## RFoptions(seed=NA) to make them all random again
\dontshow{StartExample()}
## here a model with random scale parameter
model <- RMgauss(scale=exp(rate=1))
x <- seq(0,10,0.02)
n <- 10
\dontshow{if(RFoptions()$internal$examples_reduced){warning("reduced 'n'"); n<-1}}
for (i in 1:n) {
readline(paste("Simulation no.", i, ": press return", sep=""))
plot(RFsimulate(model, x=x, seed=i))
}
## another possibility to define exactly the same model above is
## model <- RMgauss(scale=exp())
## note that however, the following two definitions lead
## to covariance models with fixed scale parameter:
## model <- RMgauss(scale=exp(1)) # fixed to 2.7181
## model <- RMgauss(scale=exp(x=1)) # fixed to 2.7181
## here, just two other examples:
## fst
model <- RMmatern(nu=unif(min=0.1, max=2)) # random
for (i in 1:n) {
readline(paste("Simulation no.", i, ": press return", sep=""))
plot(RFsimulate(model, x=x, seed=i))
}
## snd
## note that the fist 'exp' refers to the exponential function,
## the second to the exponential distribution.
(model <- RMgauss(var=exp(3), scale=exp(rate=1)))
plot(z <- RFsimulate(model=model, x=1:100/10))
\dontshow{FinalizeExample()}
}
\keyword{spatial}
\keyword{models}