Re: [R-sig-Geo] Spatial Weights

2008-10-20 Thread Roger Bivand
On Mon, 20 Oct 2008, Raphael Saldanha wrote: Hi! Can I create queen or rook spatial weights with R? Yes, see ?poly2nb in spdep. I'm migrating from Geoda to R, is there any good introductory text for basic functions of spatial econometrics with R ? See Luc Anselin's tutorials: http://www

[R-sig-Geo] Spatial Weights

2008-10-20 Thread Raphael Saldanha
Hi! Can I create queen or rook spatial weights with R? I'm migrating from Geoda to R, is there any good introductory text for basic functions of spatial econometrics with R ? Best regards, Raphael Saldanha [[alternative HTML version deleted]] ___

Re: [R-sig-Geo] help to Improve "elegant" code

2008-10-20 Thread James Nylen
Follow the technique that your code already uses (a for() loop to iterate over a series of values, and a paste() function to generate the filename). To be more specific: for(method in 1:3) { ### here, insert the code you posted ## but change: # paste("DCM_1_power", i+14, ".sgrd", sep="") ## to #

[R-sig-Geo] help to Improve "elegant" code

2008-10-20 Thread Alessandro
Hi all, I wish improve this code. I have 3 raster DCM_1, DCM_2, and DCM_3 (where: DCM= Digital Crown Model) obtained from three methodology. To apply my analysis I wish to elevate every raster to power of 1.5 to 4.0 by 0.1. Respect to write three times the below formula I wish to understatin

Re: [R-sig-Geo] Res: Problem with points.geodata and plot.geodata

2008-10-20 Thread James Nylen
In that case, you could always try stacking the panels vertically instead of horizontally? If this is even possible with the function you are using. Or just plot each panel to an image and superimpose them using another method? I've done things like that where I had R create PNG images and then c

Re: [R-sig-Geo] Res: Problem with points.geodata and plot.geodata

2008-10-20 Thread Rubén Roa-Ureta
Gledson Picharski wrote: Hi Rubén, Try this, my.func <- set.coords.lims #if you want this geoR function latter. set.coords.lims <-function(coords)apply(coords,2,range) plot(z1) points(z1,pt.divide=c("data.proportional"),lambda=1,col="gray",yl="",xl="") ## to get set.coords.lims back. set.co

[R-sig-Geo] spsample and small sample size

2008-10-20 Thread Patrick Giraudoux
Dear listers, I met a trouble with spsample which has been solved by Roger Bivand off list. I though it was general enough to be of interest for other users. In brief I needed to get 50 random points on a SpatialLinesDataFrame called "hydntline" and got: spsample(hydntline,50,type="random")