
Last chance! 50% off unlimited learning
Sale ends in
"fasp"
.
A method for plot
.plot(x)
plot(x, formula=formula)
plot(x, formula=formula, subset)
plot(x, formula=formula, subset, lty, col, title, ...)
"fasp"
representing a
function array.formula
is a list, its $k^{th}$ component
should be applicable to the $(i,j)^{th}$
plot where x$which[i,j]=k
. Isubset
is a list, its $k^{th}$ component
should beplot.default
to control other
features of the individual plot panels.subset
argument may be a
logical vector (of the same length as the vectors of data which
are extracted from x
), or a vector of indices, or an
expression such as expression(r<=0.2)< code="">, or a text string,
such as "r<=0.2"< code="">. Attempting a syntax such as subset = r<=0.2< code=""> (without
wrapping r<=0.2< code=""> either in quote marks or in expression()
)
will cause this function to fall over.
Variables referred to in any formula must exist in the data frames
stored in x
. What the names of these variables are will
of course depend upon the nature of x
.
"fasp"
represents
an array of summary functions, usually associated with a point
pattern. See fasp.object
for details.
Such an object might be created, for example, by alltypes
or allstats
. The function plot.fasp
is
a method for plot
. It calls plot.fv
to plot the
individual panels.
For information about the interpretation of the
arguments formula
, subset
, lty
and col
,
see plot.fv
.
The argument title
, if present, will determine the
overall title of the plot. If it is absent, it defaults to x$title
.
Titles for the individual plot panels will be taken from
x$titles
.
alltypes
,
allstats
,
plot.fv
,
fasp.object
library(spatstat)
# Bramble Canes data.
data(bramblecanes)
X.G <- alltypes(bramblecanes,type="G",dataname="Bramblecanes",verb=TRUE)
plot(X.G)
plot(X.G,subset="r<=0.2")
plot(X.G,formula=cbind(asin(sqrt(km)),
asin(sqrt(theo))) ~ asin(sqrt(theo)))
plot(X.G,fo=cbind(km-theo,0)~r,"r<=0.2")
# Swedish pines.
data(swedishpines)
X <- allstats(swedishpines,dataname="Swedish Pines")
plot(X,subset=list("r<=20","r<=20","r<=20","r<=50"))
# Simulated data.
pp <- runifpoint(350, owin(c(0,1),c(0,1)))
pp$marks <- factor(c(rep(1,50),rep(2,100),rep(3,200)))
X.K <- alltypes(pp,type="K",verb=TRUE,dataname="Fake Data")
plot(X.K,fo=cbind(border,theo)~theo,"theo<=0.75")
Run the code above in your browser using DataLab