An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20130707/814ac3eb/attachment.pl>
spatstat output
3 messages · catalin roibu, Rolf Turner, Hodgess, Erin
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 = TRUE,savepatterns=TRUE,correction="isotropic") And second question, is there a possibility to modify the margin of plot in spatstat? 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
confused.
The abbreviation "CI" would appear to denote "confidence interval". The
envelopes
do ***NOT*** consist of confidence intervals. They are ***critical
envelopes***.
(Repeat after me, 50 times: "Critical envelope, critical envelope,
....."). See the help
for envelope().
(4) The use of nsim=100, while not incorrect, is bizarre. The
envelopes yield
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
array for
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
which
you applied alltypes(). Each entry is an object of class "fv". See the
help for
fv.object. If you examine the names of these objects --- which are data
frames
with extra attributes --- you should easily be able to see how to
extract the
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?
cheers,
Rolf Turner
Hello! For the first question, here is an example:
E <- alltypes(bramblecanes, Kest, nsim = 100, envelope =
+ TRUE,savepatterns=TRUE,correction="isotropic")
E$fns[[1]]$r E$fns[[1]]$obs
From: r-sig-geo-bounces at r-project.org [r-sig-geo-bounces at r-project.org] on behalf of catalin roibu [catalinroibu at gmail.com]
Sent: Sunday, July 07, 2013 5:12 AM
To: r-sig-geo at stat.math.ethz.ch
Cc: r-help at r-project.org
Subject: [R-sig-Geo] spatstat output
Sent: Sunday, July 07, 2013 5:12 AM
To: r-sig-geo at stat.math.ethz.ch
Cc: r-help at r-project.org
Subject: [R-sig-Geo] spatstat output
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 =
TRUE,savepatterns=TRUE,correction="isotropic")
And second question, is there a possibility to modify the margin of plot in
spatstat?
plot(E, sqrt(./pi) - r ~ r, ylab = "L(r)-r",main=NULL,sub=NULL,las=1)
Thank you very much for your help!
CR
--
---
Catalin-Constantin ROIBU
Lecturer PhD, Forestry engineer
Forestry Faculty of Suceava
Str. Universitatii no. 13, Suceava, 720229, Romania
office phone +4 0230 52 29 78, ext. 531
mobile phone +4 0745 53 18 01
+4 0766 71 76 58
FAX: +4 0230 52 16 64
silvic.usv.ro
[[alternative HTML version deleted]]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo at r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo