Calculate four standard summary functions of a point pattern.
Calculates the $F$, $G$, $J$, and $K$
summary functions for an unmarked point pattern.
Returns them as a function array (of class
allstats(pp, dataname=NULL, verb=FALSE)
- The observed point pattern, for which summary function
estimates are required. An object of class
"ppp". It must not be marked.
- A character string giving an optional (alternative) name for the point pattern.
- A logical value meaning ``verbose''. If
TRUE, progress reports are printed during calculation.
This computes four standard summary statistics for a
point pattern: the empty space function $F(r)$,
nearest neighbour distance distribution function $G(r)$,
van Lieshout-Baddeley function $J(r)$
and Ripley's function $K(r)$.
The real work is done by
Kest respectively. Consult the help files for these functions
for further information about the statistical interpretation
of $F$, $G$, $J$ and $K$.
TRUE, then ``progress reports''
(just indications of completion) are printed out when the
calculations are finished for each of the four function types.
The overall title of the array of four functions
(for plotting by
will be formed from the argument
If this is not given, it defaults to the expression
pp given in the call to
- A function array (an object of class
fasp.object). This can be plotted using
The function array has dimensions $2 \times 2$ with the entries at positions
[2,2]representing $F(r)$, $G(r)$, $J(r)$ and $K(r)$ respectively. Each function entry
fns[[i]]retains the format of the output of the relevant estimating routine
The default formulae for plotting these functions are
cbind(km,theo) ~ rfor F, G, and J, and
cbind(trans,theo) ~ rfor K.
data(swedishpines) a <- allstats(swedishpines,dataname="Swedish Pines") plot(a) plot(a, subset=list("r<=15","r<=15","r<=15","r<=50"))