Re: [R-sig-Geo] Cluster identification using space-time permutations

2010-08-30 Thread Virgilio Gómez-Rubio
Dear Jesse, You can look at the functions in the DCluster package for the (purely) spatial scan statistic. You could start from there. I have also some other code for spatio-temporal clustering that I have not yet released (and I believe that the permutation algorithm may not be tuned) but I will

Re: [R-sig-Geo] Spatial data tower of babel

2010-08-30 Thread Virgilio Gómez-Rubio
Hi, > geoR has SpatialPointsDataFrame to geodata - but does it have the > other way round too? Or is that in sp? It doesn't matter too much, > since students will find them either way, but does > as.sp(as.geodata('meuse","zinc")) get you back where you started? But the problem is that this is no

Re: [R-sig-Geo] Adding miniature plots at specific locations on a map

2010-04-05 Thread Virgilio Gómez-Rubio
Hi, You could also check the book "R Graphics" by Paul Murrell. You can see an example of a map (of Australia) with some other graphs on it here: http://www.stat.auckland.ac.nz/~paul/RGraphics/chapter7.html He also provides the source code. Hope this helps, Virgilio El dom, 04-04-2010 a las 1

[R-sig-Geo] Call for papers: useR! 2010 conference

2010-01-20 Thread Virgilio Gómez-Rubio
Dear all, I believe that this conference will be of interest to some of the people on this list. We are hoping to organise a specific session(s) on spatial and spatio-temporal data analysis, depending on the number of submissions. There will also be a 3-hour tutorial on spatial data analysis the

Re: [R-sig-Geo] Converting nb list containing areas with no neighbours with nb2WB()

2009-11-25 Thread Virgilio Gómez-Rubio
Hi, > I checked the adjacency matrix of the spatial weights list output object, > and a zero integer was placed at the position corresponding to the area with > no neighbours. Yes. Probably we need to fix this so that nothing is added when an area has no neighbours. > When I re-ran the analysis

Re: [R-sig-Geo] Space-time cluster analysis

2009-10-26 Thread Virgilio Gómez-Rubio
Hi Joe, > I am using space-time cluster analysis to search for areas/times of high > disease incidence. My data is aggregated by geographic areas with known > background population numbers. So under the null hypothesis, the number of > cases in a given location/time period is Poisson distribute

Re: [R-sig-Geo] variogram estimation with adaptive sampling weights

2009-10-19 Thread Virgilio Gómez-Rubio
Hi John, Not sure how to address your initial question but it looks like what you are doing is similar to preferential sampling. If so, you may want to look at the paper Geostatistical Inference under Preferential Sampling. The authors provide provide code and data at http://www.lancs.ac.uk/~diggl

Re: [R-sig-Geo] Problem with Sweave and spplot

2009-09-01 Thread Virgilio Gómez-Rubio
Hi, > I'm facing problems using Sweave with spplot... there is no error > message but the pdf produced has errors and then latex compilation > fails... I really don't know where that come from, did anyone > experience something similar? Have you embedded the spplot() function inside a print() fun

Re: [R-sig-Geo] Question concerning the Kulldorff's method

2009-08-20 Thread Virgilio Gómez-Rubio
Dear Everton, > Regarding the rotine that identifies the possible cluster locations, I do > not know which value should I consider for the fractpop. The spatial > clustering is very associated to this value and, depending of the considered > value for the fractpop, I get a different spatial confi

Re: [R-sig-Geo] Bi- or Tri-variate maps for R

2009-08-11 Thread Virgilio Gómez-Rubio
Hi, > Does anyone know any examples of Bi- or Tri-variate maps > for R like this ArcVie example? I do not know of any code to do this. However, I am not sure if this is the best way of displaying bivariate data in one map. The map in Section IV is difficult to interpret (to me at least). Wouldn'

Re: [R-sig-Geo] PBSmapping to SpatialPolygonsDF

2007-03-26 Thread Virgilio Gómez-Rubio
Hi, I am also interested in this convertion function. In fact, yesterday night I wrote some code to convert from PolySet to SP, but Roger's code is much better than mine... If nobody is already doing it, can I take Roger's code, add a few checks and make a function to be included in, say, maptoo

Re: [R-sig-Geo] help on Polygon(sp)?

2007-03-21 Thread Virgilio Gómez-Rubio
Hi, > #errors occurs > env.lab <- Kenv.label(case_pts,control_pts,splancs_poly, nsim=999, seq > (0,23000,500)) > error:v error in erifyclass(X, "ppp") : argument 'X' is not of class > 'ppp' > > Q1: Now that khat(splancs) goes fine, why not Kenv.label(splancs)? > From the help, it didn't

Re: [R-sig-Geo] help on Polygon(sp)?

2007-03-21 Thread Virgilio Gómez-Rubio
Hi, El mié, 21-03-2007 a las 22:45 +0800, zhijie zhang escribió: > Dear friends, > I have a .shp file of region boundary, and i want to use it as the polygon > in the khat(splancs) function, but it goes wrong. > # khat(pts,poly,s,newstyle=FALSE) > > guichi<-readShapePoly("d:/d/guichi2.shp") #

[R-sig-Geo] Small Area Estimation with R (Request for Comments)

2007-02-02 Thread Virgilio Gómez-Rubio
Dear all, First of all, apologies if I am making a wrong use of this list, but we have been working on a package to implement EBLUP and Spatial EBLUP estimators (code provided by Nicola Salvati) for Small Area Estimation and we would like to announce it in this list. The package itself contains a

Re: [R-sig-Geo] CVS checkout instructions out of date at http://r-spatial.sourceforge.net/

2007-02-02 Thread Virgilio Gómez-Rubio
Hi David, My setting to connect to the CVS are: export CVSROOT=":ext:[EMAIL PROTECTED]:/cvsroot/r-spatial" export CVS_RSH=ssh You will need to change the server to r-spatial.cvs.sf.net. I do not know why, but it works. Best, Virgilio El vie, 02-02-2007 a las 14:45 -0600, David Forrest escrib

Re: [R-sig-Geo] Help with the Kmeasure function from the Spatstat package

2006-12-08 Thread Virgilio Gómez-Rubio
Hi, El vie, 08-12-2006 a las 09:59 +0800, Adrian Baddeley escribió: > Virgilio Gomez-Rubio writes: > > > Diggle (2003), for example, describes what you ask and most of the > > examples in the book can be reproduced using package Splancs. Check > > manual page of function mse2d and references (in

Re: [R-sig-Geo] Help with the Kmeasure function from the Spatstat package

2006-12-06 Thread Virgilio Gómez-Rubio
Dear Debarchana, I think that you could check any of the good books about point patterns that there are currently out there. Diggle (2003), for example, describes what you ask and most of the examples in the book can be reproduced using package Splancs. Check manual page of function mse2d and ref

Re: [R-sig-Geo] Error in 'Jest' and 'Fest' funcstions of Spatstat package

2006-12-04 Thread Virgilio Gómez-Rubio
Hi, I think that your problem is related to the scale of your data. Can you give us the bounding box of the region that contains your data. Perhaps you need to change that 0.01 to a more reasonable value. Best, Virgilio El lun, 04-12-2006 a las 13:35 -0600, Debarchana Ghosh escribió: > ppp_wn