The attached figure is a mark correlation function produced by the spatstat package in R. The summary returned from the object used to produce this figure is as follows:
r theo trans iso
Min. : 0.00 Min. :1 Min. :0.1663 Min. :0.1678
1st Qu.:15.27 1st Qu.:1 1st Qu.:0.1759 1st Qu.:0.1765
Median :30.53 Median :1 Median :0.5172 Median :0.4858
Mean :30.53 Mean :1 Mean :0.5695 Mean :0.5652
3rd Qu.:45.80 3rd Qu.:1 3rd Qu.:0.9435 3rd Qu.:0.9286
Max. :61.06 Max. :1 Max. :1.1325 Max. :1.2666
When I enter the name of the object on its own, I get the following:
lo k[mm][lo](r) lower pointwise envelope of k[mm](r) from simulations
hi k[mm][hi](r) upper pointwise envelope of k[mm](r) from simulations
I would like to return the actual values, as in the former output, for the simulations labelled "lo" and "hi" in the latter output.
In other words, I would like to be able to obtain a range of values at which the observed (black) line falls outside the grey envelope.
Is this possible?
I get the impression it may have something to do with the argument "savefuns", but I am at a loss to know how to implement it and return the values I need.
It would be easier to answer if you supplied the R code you're using, and sample data. In any case, it may be helpful to point out that the str() function will tell you about the contents of an object. When you name the object, you are calling the print() function on it.
In any case, you should probably consult the documentation on whatever function you're using to do the calculations (look for the subsection called "Value"). In most cases that will tell you what you want. R documentation tends to be of high quality, and you'll find it helpful, I think (or hope!).
hi Tom,
I hope the following code might be of some help to you.
###R-code###
library(spatstat)
X
Thanks to both! Using your answers, I have hit upon a solution: I used
str(Xenv$lo)
which returned the type, number of data points and the first ten points, then simply
Xenv$lo
which printed all 513 data points that form the lower envelope. From this, I can plot a subset of the data with:
plot(Xenv$lo[402:421]), to check where the observed data intersects with the lower envelope.
Which was my (initially unstated) intention. Thanks again.
- Badges
- Science topic
- Similar topics
- Mathematics
- Statistics
- Statistical Software
- R