Dear Patrick,
(a) Your attachment did not come through; only attachments with
extensions belonging to a limited class are passed through by the list
software.
Perhaps you could email your data set to me privately,
(b) I have only had a quick glance at your code. A few comments:
* In your first call to ppp(), "range(ColN), range(Row)" will be
ignored. They would be used to construct a (rectangular) window
except for the fact that you have specified "window" explicitly.
* the argument "maxsize" gets used only if the marks involved are
*numeric*; it'll be ignored when the marks are a factor.
* the argument "chars" should take values which would be acceptable
as values of the "pch" argument to plot(). The levels of
"Concept" do not meet this requirement.
* if I were me I wouldn't clutter things up by plotting the density
estimate initially. Just plot the pattern, until you get the
sort of result you want. *Then* overlay it on top of a density
plot.
(c) If you send me the data off-list and explain in words what it is
that you want to achieve ("Tell me what you want to do not how you want
to do it.") then I'm pretty sure that I can show you how to achieve it.
cheers,
Rolf
On 15/05/14 14:25, p_connolly wrote:
The spatstat package has hundreds of useful functions but I'm having
trouble understanding the intricate ways it does things. I've read lots of
Rspatialcourse_CMIS_PDF Standard.pdf from here:
http://www.csiro.au/resources/pf16h but can't find what I need to know. I'm
particularly interested in how to set the plotting character and its
size in density plots.
require(spatstat)
load("CharacterSizeRData") # (attached to this message)
ls()
[1] "all.win" "bb9"
head(bb9)
Col Row ColN Concept Pch VinesJune VinesAug New Event
1 A 10 3.333333 Gold 4x 4 42 41 1 TRUE
8 A 16 3.333333 Gold 4x 4 42 41 1 TRUE
10 A 18 3.333333 Gold 4x 4 42 41 1 TRUE
11 A 19 3.333333 Gold 4x 4 42 41 1 TRUE
14 A 21 3.333333 Gold 4x 4 42 41 1 TRUE
18 A 24 3.333333 Gold 4x 4 42 34 8 TRUE
## If I make a ppp class object from bb9, like so:
ppp9A <- with(bb9, ppp(ColN, Row, range(ColN), range(Row), window =
all.win,
marks = Concept))
## I can do a density plot:
plot(density(ppp9A, 5, edge = FALSE, weights = bb9$New), main = "ppp9A")
## Adding to show the points:
plot(ppp9A, add = TRUE, cols = "red", maxsize = max(bb9$New)/20)
## works except that maxsize is ignored.
##
## If I make the ppp object like this:
ppp9B <- with(bb9, ppp(ColN, Row, range(ColN), range(Row), window =
all.win,
marks = New))
plot(density(ppp9B, 5, edge = FALSE, weights = bb9$New), main = "ppp9B")
plot(ppp9B, add = TRUE, cols = "red", maxsize = max(bb9$New)/20)
## maxsize works how I intended.
##
## I tried to combine plotting character and size with a fancier use of
marks
ppp9C <- with(bb9, ppp(ColN, Row, range(ColN), range(Row), window =
all.win,
marks = data.frame(New, Concept)))
## Then try
plot(density(ppp9C, 5, edge = FALSE, weights = marks(ppp9C)$New),
main = "")
plot(ppp9C, add = TRUE, cols = "red", chars =
levels(marks(ppp9C)$Concept),
maxsize = max(marks(ppp9C)$New)/20)
## chars is ignored and no error message to complain or explain
What I wish to know is how to get plotting character and size
information displayed simultaneously. I tried ways of specifying it
from bb9 dataframe instead of the ppp object, or specifying a fixed
value but it's always ignored. Ideally, I'd like to customize the
plotting characters used, but getting the defaults to work would be a
start.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
