[R-sig-Geo] About GeoBUGS and R

2007-11-19 Thread            
Hi, all I use WinBUGS by way of R (R2WinBUGS). So,I want to use GeoBUGS like this or another method (batch / script etc.. not GUI. I would like to use R,if possible). Do you know what to do? Thanks. --- Minag ([EMAIL PROTECTED]) -- New Design Yahoo! JAPAN

Re: [R-sig-Geo] Transforming from cartographic to arbitrary local coordinate system

2007-11-19 Thread Roger Bivand
On Mon, 19 Nov 2007, Don MacQueen wrote: > I have virtually no experience with cartographic coordinate systems, > and transforming them, so I would like to ask for some help. > > Someone has supplied me with a shapefile. When I open it in Qgis, > Qgis tells me its "spatial reference system" is: >

[R-sig-Geo] finding points outside but very near a polygon

2007-11-19 Thread Medina, Hernan
I used the code below to find points in a North American Regional Reanalysis data set that were within the 48 contiguous states. # The 48 contiguous states m <- map("state",fill=TRUE,plot=FALSE) indx1 <- in.polygon(m,list(x=my_lon,y=my_lat)) st <- map.where("state",my_lon,my_lat) Unfortunately,

Re: [R-sig-Geo] Transforming from cartographic to arbitrary local coordinate system

2007-11-19 Thread Dylan Beaudette
On Monday 19 November 2007, Don MacQueen wrote: > I have virtually no experience with cartographic coordinate systems, > and transforming them, so I would like to ask for some help. > > Someone has supplied me with a shapefile. When I open it in Qgis, > Qgis tells me its "spatial reference system"

[R-sig-Geo] Transforming from cartographic to arbitrary local coordinate system

2007-11-19 Thread Don MacQueen
I have virtually no experience with cartographic coordinate systems, and transforming them, so I would like to ask for some help. Someone has supplied me with a shapefile. When I open it in Qgis, Qgis tells me its "spatial reference system" is: +proj=lcc +lat_1=35.46667 +lat_2=34.033

Re: [R-sig-Geo] spplot(), SpatialGridDataFrame

2007-11-19 Thread Edzer J. Pebesma
Ingo Holz wrote: > Hi, > > I have a SpatialGridDataFrame and want to plot it with spplot(). > > Is it possible to plot only the grids that have a special value (eg. 312)? > fullgrid(x) = FALSE spplot(x[x$value == 312,"value"] Note that if fullgrid(x) is TRUE, x[rows,cols] would select on row

[R-sig-Geo] spplot(), SpatialGridDataFrame

2007-11-19 Thread Ingo Holz
Hi, I have a SpatialGridDataFrame and want to plot it with spplot(). Is it possible to plot only the grids that have a special value (eg. 312)? If it is not possible to do this with spplot() how would it be done with plot()? Thanks, Ingo ___ R-si