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]]
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@r-project.org [r-sig-geo-bounces@r-project.org] on behalf of catalin roibu [catalinroibu@gmail.com] Sent: Sunday, July 07, 2013 5:12 AM To: r-sig-geo@stat.math.ethz.ch Cc: r-help@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@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
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 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