[R] spatstat output
rolf.turner at xtra.co.nz
Sun Jul 7 13:05:28 CEST 2013
On 07/07/13 22:12, catalin roibu wrote:
> Dear R users,
> Is there a possibility to extract only the r, CI's envelope and L function
> from the output of spatstat?
> I use this code
> E <- alltypes(df1, Kest, nsim = 100, envelope =
> And second question, is there a possibility to modify the margin of plot in
> plot(E, sqrt(./pi) - r ~ r, ylab = "L(r)-r",main=NULL,sub=NULL,las=1)
> Thank you very much for your help!
(1) This is a question about the spatstat package and as such should be
in the first instance directed to the authors of this package.
(2) Your example is *not reproducible* since we have no access to "df1".
(3) Your use of the expression "CI's envelope" indicates that you are
The abbreviation "CI" would appear to denote "confidence interval". The
do ***NOT*** consist of confidence intervals. They are ***critical
(Repeat after me, 50 times: "Critical envelope, critical envelope,
....."). See the help
(4) The use of nsim=100, while not incorrect, is bizarre. The
p-values equal to fractions whose denominator is (nsim+1). Hence nsim=99
(the default) is usually a much better choice than nsim=100.
(5) The object E returned by your code is of class "fasp" (function
spatial processes). It is thus a list whose first entry is names
"fns". In turn
"fns" is a list whose entries correspond to the marks of the pattern to
you applied alltypes(). Each entry is an object of class "fv". See the
fv.object. If you examine the names of these objects --- which are data
with extra attributes --- you should easily be able to see how to
items that you desire.
(6) In respect of changing the margins, see the help for plot.fasp. Does
the argument mar.panel do what you want?
More information about the R-help