[R-sig-Geo] simulation envelopes with Kmark(), for an inhomogeneous point pattern

2024-06-21 Thread Christopher W. Ryan via R-sig-Geo
It's been a while since I've worked with point patterns. I have a marked point pattern X consisting of household locations, with the marks being waiting time to get a lab test. I'm interested in the spatial distribution of those waiting times. Among other exploratory methods, I'm looking at

Re: [R-sig-Geo] Learning Resources Spatial Regression Models from the ground up

2024-04-24 Thread Christopher W. Ryan via R-sig-Geo
Josiah-- I've found the following very helpful over the years: Geographic Information Analysis, by David O'Sullivan and David Unwin Spatial Point Patterns, by Adrian Baddeley, Ege Rubak, and Rolf Turner Applied Spatial Data Analysis with R, by Roger Bivand, Edzer Pebesma, and Virgilio

Re: [R-sig-Geo] making sequential subsets of a ppp object, based on value of a mark

2021-01-06 Thread Christopher W. Ryan
Ah, it's split(), of course. That was silly of me. --Chris Ryan Christopher W. Ryan wrote: > I have a marked ppp object, about 5300 points. One of the marks is > called period and has values 1:13. What is an efficient way to make 13 > period-specific subset point patterns? All I've

[R-sig-Geo] making sequential subsets of a ppp object, based on value of a mark

2021-01-06 Thread Christopher W. Ryan
I have a marked ppp object, about 5300 points. One of the marks is called period and has values 1:13. What is an efficient way to make 13 period-specific subset point patterns? All I've been able to come up with, embarrasingly, is 13 lines like the following: per1.ppp <- subset(analytical.ppp,

Re: [R-sig-Geo] How do i plot a graph with x axis vertically labelled?

2018-08-13 Thread Christopher W Ryan
Vasana-- May I ask why you want to produce such a graph? What relationship(s) in your data are you trying to show? I can't be sure from your original post, but it seems to me that you are trying to do this the hard way. R can help you. What happens when you run the code you posted? What error

Re: [R-sig-Geo] how to limit an image or a funxy to an irregular polygonal window, instead of the whole enclosing rectangle?

2017-12-15 Thread Christopher W. Ryan
Thanks Rolf. I'm going to have to reflect more on my code and my data, to understand better what is going on. Obviously this won't help you much, without having access to all my data and preceeding code, but the error message that is tripping me up is: > rhohat(m12, pov.f) Error: the fitted

[R-sig-Geo] how to limit an image or a funxy to an irregular polygonal window, instead of the whole enclosing rectangle?

2017-12-15 Thread Christopher W. Ryan
Using R 3.3.3 and spatstat 1.50-0 on Windows 7. MWE below. I have a SpatialPolygonsDataFrame of census tract poverty levels in 3 contiguous counties in the US, called sremsPoverty. I want to use this as a predictor in a ppm model. The window for the point pattern is the three counties--so an

Re: [R-sig-Geo] simple online CSR pattern generator?

2017-09-05 Thread Christopher W Ryan
Unfortunately R-studio is not an option, but r-fiddle is perfect! Thanks. --Chris On Tue, Sep 5, 2017 at 1:58 PM, VIRGILIO GOMEZ RUBIO wrote: > Hi, > > You can use r-fiddle: > > http://www.r-fiddle.org/#/fiddle?id=kuUG5uW6 > > Click on run code and you will see a very

[R-sig-Geo] simple online CSR pattern generator?

2017-09-05 Thread Christopher W Ryan
I may have to do a very simple demo about random variation in spatial epidemiology, without access to R. I would have internet access. Is anyone aware of a website where one can generate repeated examples of a simple CSR point pattern, akin to rpoispp() in spatstat? Thanks. --Chris Ryan Broome

[R-sig-Geo] using spatialpolygonsdataframe in ppm (or, converting spatialpolygonsdataframe to pixel image or other object useful in ppm)

2017-09-01 Thread Christopher W. Ryan
Hello. What is the best way to use a spatialpolygonsdataframe, with a numerical variable of interest for each polygon (proportion of households in poverty for US census tracts in the region of interest) as a predictor in ppm() in spatstat? I don't think I can use it directly on the RHS of ppm(),

Re: [R-sig-Geo] For each point, distance to nearest point in second dataset

2016-05-20 Thread Christopher W Ryan
See also the spatstat package, and its nncross function. --Chris Ryan On Thu, May 19, 2016 at 6:48 PM, Nick Eubank wrote: > Last note for future searchers: also suggested was the SearchTree library: >

Re: [R-sig-Geo] removing internal boundaries from a spatstat window that was created from three adjacent counties

2016-03-01 Thread Christopher W. Ryan
Adrian, Rolf, et al-- Thanks so much; worked like a charm. Any value of r below 1 (meter, I assume?) seemed to serve well. Missing FAQ aside, your book is phenomenal. I'm working through it in concert with my current project and learning a ton. --Chris Christopher W. Ryan, MD, MS

[R-sig-Geo] removing internal boundaries from a spatstat window that was created from three adjacent counties

2016-02-15 Thread Christopher W. Ryan
never I plot it or plot subsequent point patterns that use the window. In essence, in plots it looks like 3 polygons instead of one big one. I'd prefer not to have the inter-county boundaries be visible--I'd rather have just one big polygon for the whole area. How can I remove them? Or should I crea

[R-sig-Geo] geocoding street addresses within R, not using any web-based services

2016-01-22 Thread Christopher W. Ryan
Hello. I'm a longtime R and R-help user, and recent listener on R-sig-geo. First post here. I use R on both Win 7 and Linux Mint. Is there a way to geocode street addresses in the US within R on my local machine, that is, without transmitting the addresses to any web-based service? I have about

[R-sig-Geo] density from spatstat less than zero

2015-10-15 Thread Christopher W. Ryan
. Whittling it down to 100 or so points still yields the same problem, but I think that changes the game and would not be a true reflection of my situation. Thanks. --Chris -- Christopher W. Ryan, MD, MS cryanatbinghamtondotedu https://www.linkedin.com/in/ryancw Early success is a terrible teacher. You’re