https://github.com/cran/spatstat
Tip revision: 9f84c79450f33061219ca8cc5518816c31a9be64 authored by Adrian Baddeley on 21 November 2017, 07:39:44 UTC
version 1.54-0
version 1.54-0
Tip revision: 9f84c79
plot.rppm.Rd
\name{plot.rppm}
\alias{plot.rppm}
\title{
Plot a Recursively Partitioned Point Process Model
}
\description{
Given a model which has been fitted to point pattern data
by recursive partitioning, plot the partition tree
or the fitted intensity.
}
\usage{
\method{plot}{rppm}(x, \dots, what = c("tree", "spatial"), treeplot=NULL)
}
\arguments{
\item{x}{
Fitted point process model of class \code{"rppm"}
produced by the function \code{\link{rppm}}.
}
\item{what}{
Character string (partially matched) specifying whether to plot the
partition tree or the fitted intensity.
}
\item{\dots}{
Arguments passed to \code{\link[rpart]{plot.rpart}}
and \code{\link[rpart]{text.rpart}} (if \code{what="tree"})
or passed to \code{\link{plot.im}} (if \code{what="spatial"})
controlling the appearance of the plot.
}
\item{treeplot}{
Optional. A function to be used to plot and label the partition tree,
replacing the two functions \code{\link[rpart]{plot.rpart}}
and \code{\link[rpart]{text.rpart}}.
}
}
\details{
If \code{what="tree"} (the default), the partition tree will be plotted
using \code{\link[rpart]{plot.rpart}}, and labelled using
\code{\link[rpart]{text.rpart}}.
If the argument \code{treeplot} is
given, then plotting and labelling will be performed by
\code{treeplot} instead. A good choice is the function
\code{prp} in package \pkg{rpart.plot}.
If \code{what="spatial"}, the predicted intensity
will be computed using \code{\link{predict.rppm}}, and
this intensity will be plotted as an image using \code{\link{plot.im}}.
}
\value{
If \code{what="tree"}, a list containing \code{x} and \code{y}
coordinates of the plotted nodes of the tree.
If \code{what="spatial"}, the return value of \code{\link{plot.im}}.
}
\author{
\spatstatAuthors
}
\seealso{
\code{\link{rppm}}
}
\examples{
# Murchison gold data
mur <- solapply(murchison, rescale, s=1000, unitname="km")
mur$dfault <- distfun(mur$faults)
#
fit <- rppm(gold ~ dfault + greenstone, data=mur)
#
opa <- par(mfrow=c(1,2))
plot(fit)
plot(fit, what="spatial")
par(opa)
}
\keyword{spatial}
\keyword{hplot}
\keyword{models}