[R-sig-Geo] readOGR problem in Windows

2014-12-23 Thread James Rooney
Dear List, I have some code I have developed on a Mac that I'm trying to transfer onto a higher spec windows 8.1 machine and I am running into a problem with the readOGR function in rgdal. My code is below. This works fine under Mac, but under windows readOGR gives me the error: Error in

Re: [R-sig-Geo] readOGR problem in Windows

2014-12-23 Thread James Rooney
/SA.basic.analysis.data.csv, header=TRUE) SA@data = data.frame(SA@data, df[match(SA@data$SA_CODE, df$SA_CODE),]) From: Alex Mandel [tech_...@wildintellect.com] Sent: 23 December 2014 16:43 To: James Rooney; r-sig-geo@r-project.org Subject: Re: [R-sig-Geo] readOGR

Re: [R-sig-Geo] readOGR problem in Windows

2014-12-23 Thread James Rooney
Fantastic! That fixed it. Thanks much appreciated! James From: Michael Sumner [mdsum...@gmail.com] Sent: 23 December 2014 21:00 To: James Rooney; r-sig-geo@r-project.org Subject: Re: [R-sig-Geo] readOGR problem in Windows Try removing the trailing slash

Re: [R-sig-Geo] rgdal writeOGR - clipping variable names

2014-07-07 Thread James Rooney
:12 To: James Rooney Cc: R-sig-Geo@r-project.org Subject: Re: [R-sig-Geo] rgdal writeOGR - clipping variable names On Sun, 6 Jul 2014, James Rooney wrote: Hi all, I am using writeOGR to save ESRI shapefiles after a ton of processing. Typical warning message I get is as follows: writeOGR(ED

Re: [R-sig-Geo] rgdal writeOGR - clipping variable names

2014-07-07 Thread James Rooney
...@lancaster.ac.uk] Sent: 07 July 2014 09:44 To: James Rooney Cc: R-sig-Geo@r-project.org Subject: Re: [R-sig-Geo] rgdal writeOGR - clipping variable names Ten characters in field names should be enough for everyone - Jack Dangermond (ESRI founder and billionaire) Okay, he didn't really say that (just as Bill

[R-sig-Geo] rgdal writeOGR - clipping variable names

2014-07-06 Thread James Rooney
Hi all, I am using writeOGR to save ESRI shapefiles after a ton of processing. Typical warning message I get is as follows: writeOGR(ED, /Users/jrooney/Documents/My Projects/ALS Spatial - Advanced/Irish Data/Processed Data/, layer=ED.cluster,driver=ESRI Shapefile,overwrite=T) Warning

Re: [R-sig-Geo] Choropleth map

2014-07-03 Thread James Rooney
From: Janka VANSCHOENWINKEL [janka.vanschoenwin...@uhasselt.be] Sent: 03 July 2014 09:35 To: James Rooney; roger.biv...@nhh.no Cc: r-sig-geo@r-project.org Subject: Re: [R-sig-Geo] Choropleth map Hi James and Roger, You both pointed out that maybe I should not use

Re: [R-sig-Geo] Choropleth map

2014-06-25 Thread James Rooney
Dear Janka, I can't answer your ggplot questions as I don't really use it much however I did spot something in your code that will cause problems. You used a merge to combine data between a dataframe and an SPDF: ValueMapDf - merge(nuts, Value, by.x=NUTS_ID, by.y=NUTS_ID) I've forgotten the

Re: [R-sig-Geo] BYM model in R-INLA - how to specify priors individually

2014-05-06 Thread James Rooney
of variation. What is if that controls the degree of smoothing for want of a better term! Many thanks, James From: VIRGILIO GOMEZ RUBIO [virgilio.go...@uclm.es] Sent: 05 May 2014 23:40 To: James Rooney Cc: r-sig-geo@r-project.org Subject: Re: [R-sig-Geo] BYM

[R-sig-Geo] BYM model in R-INLA - how to specify priors individually

2014-05-05 Thread James Rooney
Hi all, I am attempting to implement the Besag-York-Mollie model in R-INLA. Specifically to recreate the WinBUGS model specified in ASDAR fig 11.10. However I cannot figure out how to specify the spatial and random priors to be different as they are in the WinBUGS model fig11.10. i.e. from

[R-sig-Geo] areal data with point covariates

2014-04-17 Thread James Rooney
Dear List, This is not strictly an R question so I hope it is allowed. I am a little unsure of what type of analysis I need to use next so hoping someone can point me in the right direction. I have point data representing cases of disease. I also have population data attached to polygons. I

Re: [R-sig-Geo] algorthirm to join polygons based on population properties

2014-01-09 Thread James Rooney
islands and peninsulas etc. Now that I know what the problem is I should be able to fix it. Many thanks! James From: Roger Bivand [roger.biv...@nhh.no] Sent: 08 January 2014 19:58 To: James Rooney Cc: r-sig-geo@r-project.org Subject: RE: [R-sig-Geo] algorthirm

Re: [R-sig-Geo] algorthirm to join polygons based on population properties

2014-01-07 Thread James Rooney
that will create me a join list based on such rules ? I can of course write some code to do it but I was hoping not to have to spend the time on it! James From: Roger Bivand [roger.biv...@nhh.no] Sent: 07 January 2014 08:28 To: James Rooney Cc: r-sig-geo@r

Re: [R-sig-Geo] algorthirm to join polygons based on population properties

2014-01-07 Thread James Rooney
Ok thanks Roger I'll read up on that! Many thanks! James From: Roger Bivand [roger.biv...@nhh.no] Sent: 07 January 2014 10:12 To: James Rooney Cc: r-sig-geo@r-project.org Subject: Re: [R-sig-Geo] algorthirm to join polygons based on population properties

Re: [R-sig-Geo] algorthirm to join polygons based on population properties

2014-01-07 Thread James Rooney
to be called James! James From: Kuroda [kuroda@gmail.com] Sent: 07 January 2014 11:19 To: James Rooney Cc: r-sig-geo@r-project.org Subject: Re: [R-sig-Geo] algorthirm to join polygons based on population properties Professor Rooney: As professor

[R-sig-Geo] algorthirm to join polygons based on population properties

2014-01-06 Thread James Rooney
Dear all, I have dataset with very many more polygons than cases. I wish to apply Bayesian smoothing to areal disease rates, however I have too many polygons and need a smart way to combine them so that there are less overall polygons. Bascially I need to only combine polygons of similar

[R-sig-Geo] spatial polygon - max id length

2013-12-31 Thread James Rooney
Hi, I'm wondering does anyone know the maximum allowed length for the ID slot of a spatial polygons ? I'm having some troubles I think may be caused by trying to assign slot names that are too long. Can change it of course but it would be useful to know the max length to determine how much I

[R-sig-Geo] Resolved - RE: spatial polygon - max id length

2013-12-31 Thread James Rooney
...@r-project.org] On Behalf Of James Rooney [roone...@tcd.ie] Sent: 31 December 2013 14:55 To: r-sig-geo@r-project.org Subject: [R-sig-Geo] spatial polygon - max id length Hi, I'm wondering does anyone know the maximum allowed length for the ID slot of a spatial polygons ? I'm having some

Re: [R-sig-Geo] add image to plot

2013-12-11 Thread James Rooney
Anyone got experience with grImport even ? I'm at my wits end with this. Think I'm going to have to just do it in photoshop or something after the fact. James From: r-sig-geo-boun...@r-project.org [r-sig-geo-boun...@r-project.org] On Behalf Of James

Re: [R-sig-Geo] add image to plot

2013-12-11 Thread James Rooney
have it. James From: Agustin Lobo [alobolis...@gmail.com] Sent: 11 December 2013 12:42 To: James Rooney Cc: r-sig-geo@r-project.org Subject: Re: [R-sig-Geo] add image to plot James, Are you sure this is the right list for this question? I would rather use

[R-sig-Geo] add image to plot

2013-12-09 Thread James Rooney
Hi all, I need to add a logo to a plot. I have the logo in eps and jpg formats though can convert it to something else if needed. My preference is to use the eps version since I am plotting my graphs to an eps file. However I'm having a problem - with grImport - it is not detecting all the

[R-sig-Geo] Resolved - RE: Help - calling OpenBugs through R2OpenBugs problem

2013-12-01 Thread James Rooney
From: r-sig-geo-boun...@r-project.org [r-sig-geo-boun...@r-project.org] On Behalf Of James Rooney [roone...@tcd.ie] Sent: 27 November 2013 21:13 Cc: r-sig-geo@r-project.org Subject: Re: [R-sig-Geo] Help - calling OpenBugs through R2OpenBugs problem Dear Virgilio, Many thanks

Re: [R-sig-Geo] Help - calling OpenBugs through R2OpenBugs problem

2013-11-28 Thread James Rooney
...@gmail.com] Sent: 28 November 2013 09:11 To: James Rooney Cc: r-sig-geo@r-project.org Subject: Re: [R-sig-Geo] Help - calling OpenBugs through R2OpenBugs problem Hi All I try to follow the above example but i got stuck in reading Spatial Data Frame any Idea or Help? Here are my R codes with error. I am

Re: [R-sig-Geo] Help - calling OpenBugs through R2OpenBugs problem

2013-11-28 Thread James Rooney
I don't know then Yuster. Is your rgdal package up to date ? From: Yuster Ronoh [yust...@gmail.com] Sent: 28 November 2013 11:08 To: James Rooney Cc: r-sig-geo@r-project.org Subject: Re: [R-sig-Geo] Help - calling OpenBugs through R2OpenBugs problem Yes

Re: [R-sig-Geo] Help - calling OpenBugs through R2OpenBugs problem

2013-11-27 Thread James Rooney
Hi everyone, Has anyone looked at this problem below ? Even if nobody has an answer I'd be keen to know if the problem replicated for others. Many thanks, James From: r-sig-geo-boun...@r-project.org [r-sig-geo-boun...@r-project.org] On Behalf Of James

Re: [R-sig-Geo] Help - calling OpenBugs through R2OpenBugs problem

2013-11-27 Thread James Rooney
' figure out how or why. I wonder is there a file size limit for passing things into OpenBugs ? Appreciate any help! James From: Virgilio Gómez-Rubio [virgilio.go...@uclm.es] Sent: 27 November 2013 20:29 To: James Rooney Cc: r-sig-geo@r-project.org Subject: Re

[R-sig-Geo] Help - caling OpenBugs thourgh R2OpenBugs problem

2013-11-23 Thread James Rooney
Hi all, I am working on a Bayesian spatial smoothing project - and have implemented the BYM model, similar to the description in Chapter 11 of the ASDAR book. However, I have been having a huge problem calling OpenBugs from within R and I cannot figure out why. I would veyr much liek to do

Re: [R-sig-Geo] DCluster questions

2013-09-08 Thread James Rooney
Thanks for the reply Virgilio, 2) Will opgam/kn.iscluster only detect hotspots or will they also detect coldspots i.e. areas of statistically unlikely lower incidence rates ? You get a p-value, so you could take the areas with very large p-values as cold-spots. But detection of cold-spots is

[R-sig-Geo] DCluster questions

2013-09-06 Thread James Rooney
Hi all, I'm wondering if I can get a little advice on using DCluster. I have produced a map of areal incidence rates and I'd like to try and detect clusters. I have also implemented Bayesian smoothing and have therefore pre and post smoothing maps. I have some code as follow for detecting

Re: [R-sig-Geo] Openbugs doesn't seen to want to thin.. ?

2013-08-26 Thread James Rooney
these solves the problem that is ok. Thanks! James From: Virgilio Gómez-Rubio [virgilio.go...@uclm.es] Sent: 26 August 2013 11:51 To: James Rooney Cc: r-sig-geo@r-project.org Subject: Re: [R-sig-Geo] Openbugs doesn't seen to want to thin.. ? Hi James, How many

[R-sig-Geo] Openbugs doesn't seen to want to thin.. ?

2013-08-26 Thread James Rooney
Hi all, Doing some spatial smoothing using the BYM model from the ASDAR book. I was using Winbugs but had lots of troubles so switched to Openbugs. Unfortunately having a new problem running out of memory. I have tried to fix it with thinning - however Openbugs does not seem to want to thin. My

[R-sig-Geo] output data from R to SaTScan

2013-08-26 Thread James Rooney
Hi all, So I have created some maps in R, smoothed then using the BYM model in Openbugs and opened the CODA files in R once again to analyze results etc. I would like to now test for clusters using SaTScan and the post smoothing data - but I can find very little info anywhere on how to output

Re: [R-sig-Geo] output data from R to SaTScan

2013-08-26 Thread James Rooney
, col.regions=c(gray(1), gray(.5), gray(.8 Thanks, James From: Hugh Sturrock [hughsturr...@hotmail.com] Sent: 26 August 2013 17:10 To: James Rooney; r-sig-geo@r-project.org Subject: RE: [R-sig-Geo] output data from R to SaTScan Hi James

Re: [R-sig-Geo] R2Winbugs - problem with workign directory

2013-08-06 Thread James Rooney
Hi Virgilio, Many thanks. As it happens today I switched to OpenBugs and don't have this particular problem anymore. But thanks again, James From: Virgilio Gomez Rubio [virgilio.go...@uclm.es] Sent: 06 August 2013 17:21 To: James Rooney Cc: r-sig-geo@r

[R-sig-Geo] R2Winbugs - problem with workign directory

2013-08-05 Thread James Rooney
Hi all, Having a problem with the R2Winbugs working.directory option. My call to winbugs looks like this: MCMCres-bugs(debug=T,data=d,inits=inits,parameters.to.save=c(theta,alpha,tau,precv), n.chains=2,n.iter=2000,n.burnin=500,n.thin=10,

[R-sig-Geo] how to code for neighbour list editing

2013-07-15 Thread James Rooney
Hi all, Ok so I have a spatial data set with over 3000 polygons. I used poly2nb() to identify neighbour relationships and make an nb object. Of the polygons, about 30 or 40 are islands and have no neighbours according to polynb. Prior to doing Bayesian smoothing, I wish to assign neighbours to

Re: [R-sig-Geo] how to code for neighbour list editing

2013-07-15 Thread James Rooney
Many thanks Roger, I'll study that and see how I get on! James From: Roger Bivand [roger.biv...@nhh.no] Sent: 15 July 2013 19:51 To: James Rooney Cc: r-sig-geo@r-project.org Subject: Re: [R-sig-Geo] how to code for neighbour list editing On Mon, 15 Jul

Re: [R-sig-Geo] R2WinBugs - getting TRAP errors in Winbugs

2013-07-10 Thread James Rooney
': Error in chol.default(Q.current) : the leading minor of order 8 is not positive definite Thanks for the other options I will look into them. James From: Roger Bivand [roger.biv...@nhh.no] Sent: 10 July 2013 10:25 To: James Rooney; r-sig-geo@r-project.org

[R-sig-Geo] R2WinBugs - getting TRAP errors in Winbugs

2013-07-09 Thread James Rooney
Hi all, I am trying to implement the WinBUGS spatial conditional autoregression model on page 11.11 of Applied Spatial Data Analysis in R. I have about 3300 polygons with about 17,000 adjacencies in the set. I have changed the model only by removing the beta term as, for now, I don't have any

Re: [R-sig-Geo] Problem with poly2nb after unionSpatialPoly

2013-07-08 Thread James Rooney
16:46 To: James Rooney Cc: r-sig-geo@r-project.org Subject: Re: [R-sig-Geo] Problem with poly2nb after unionSpatialPoly On Mon, 8 Jul 2013, Roger Bivand wrote: On Mon, 8 Jul 2013, James Rooney wrote: Hi all, I am trying to setup a geospatial CAR model to implement areal spatial smoothing

Re: [R-sig-Geo] Problem with poly2nb after unionSpatialPoly

2013-07-08 Thread James Rooney
I've plotted look perfect. Thank you very much Roger for your help! James From: r-sig-geo-boun...@r-project.org [r-sig-geo-boun...@r-project.org] On Behalf Of James Rooney [roone...@tcd.ie] Sent: 08 July 2013 17:11 To: r-sig-geo@r-project.org Subject