Re: [R-sig-Geo] Problem Setting Projection of Rasters upon Import

2013-07-15 Thread Rainer M Krug
Robert J. Hijmans r.hijm...@gmail.com writes: Michael, I did not try it, but I think the below accomplishes the same. If at all possible, you should avoid using ascii format files (that leads to slow processing). library(raster) reference - raster(D:/GIS/California/Hab Suitability

Re: [R-sig-Geo] Spatial panel Durbin model

2013-07-15 Thread Juan Tomas Sayago Gomez
Hi, package spdep, when you use the command lagsarlm select the option type mixed. install.packages(spdep) library(spdep) ?lagsarlm Best, Juan On Sun, Jul 14, 2013 at 11:54 PM, Chen, Xi Prof. xi.c...@quinnipiac.eduwrote: Hi, Is there a package that can run spatial panel durbin model in r?

Re: [R-sig-Geo] Spatial panel Durbin model

2013-07-15 Thread Roger Bivand
On Mon, 15 Jul 2013, Juan Tomas Sayago Gomez wrote: Hi, package spdep, when you use the command lagsarlm select the option type mixed. install.packages(spdep) library(spdep) ?lagsarlm Right for cross-section models, but the question was about spatial panel models. In that case. it may be

Re: [R-sig-Geo] Spatial panel Durbin model

2013-07-15 Thread Juan Tomas Sayago Gomez
Hello, Yes, to add something to what prof Bivand wrote. to manually create the WX remember to use the Kronecker product of the time Identity and the Weight matrix. and it should be straight forward. Best, Juan. On Mon, Jul 15, 2013 at 6:45 AM, Roger Bivand roger.biv...@nhh.no wrote: On Mon, 15

[R-sig-Geo] Plotting together Lucas County observations and map

2013-07-15 Thread ROMAN MINGUEZ SALIDO
Hi everyone: I would like to plot the observations of housing prices included in the Lucas County data of spdep package. The problem is that I want to plot these observations in a map of Lucas County. I have found a map of Lucas County in an shp polygon file but the spatial coordinates are in

Re: [R-sig-Geo] Plotting together Lucas County observations and map

2013-07-15 Thread Hodgess, Erin
Hi! Use spTransform to get the spdep data into long/lat From: r-sig-geo-boun...@r-project.org [r-sig-geo-boun...@r-project.org] on behalf of ROMAN MINGUEZ SALIDO [roman.ming...@uclm.es] Sent: Monday, July 15, 2013 6:02 AM To: r-sig-geo@r-project.org

Re: [R-sig-Geo] assign a slot to a name (character) in loop

2013-07-15 Thread ipampa
after searching and searching I couldn´t find a way to that. instead, I simply came back in the loops and inserted my data.frame in the assignation of the spatialpolygonsdataframe. -- View this message in context:

[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

[R-sig-Geo] Citation for spgwr

2013-07-15 Thread Milan Sharma
Hi Roger,  Please let me know if the following way of citing spgwr is incorect: In references:  Bivand, R., Yu, D. (2013) spgwr:  Geographically weighted regression. R package version 3.0 and within the text:  We usedthe spgwr package (Bivand and Yu 2013) for GWR in the statistical software

[R-sig-Geo] Package for fitting spatial model on large dataset

2013-07-15 Thread Moshood Agba Bakare
Dear All, Please which of these packages - gstat and geoR is better for fitting spatial covariance structure model for spatial data with 20,000 observations. Kindly recommend any other package you think can fit such model. Thank you. Moshood [[alternative HTML version deleted]]

[R-sig-Geo] number of triples

2013-07-15 Thread Amen Alyaari
Dear All, I have three datasets and would like co know how much N was used in my calculations. I read the data into a multi-dimensional array with dimensions (nx, ny, ntsteps, ndatasets), e.g. with a smaller example dataset: # nx ny nsteps ndatasets dat =

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

2013-07-15 Thread Roger Bivand
On Mon, 15 Jul 2013, James Rooney wrote: 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

Re: [R-sig-Geo] Package for fitting spatial model on large dataset

2013-07-15 Thread Paulo Justiniano Ribeiro Jr
likelihood based methods and full kriging neighbourhood in some (or most) of geoR functions are not suitable for such data size Paulo Justiniano Ribeiro Jr LEG (Laboratorio de Estatistica e Geoinformacao) Universidade Federal do Parana Caixa Postal 19.081 CEP 81.531-990 Curitiba, PR - Brasil

Re: [R-sig-Geo] Citation for spgwr

2013-07-15 Thread Roger Bivand
On Mon, 15 Jul 2013, Milan Sharma wrote: Hi Roger,  Please let me know if the following way of citing spgwr is incorect: For citing, use the citation() function: citation() # for R itself citation(package) # for package Roger In references:  Bivand, R., Yu, D. (2013) spgwr: 

Re: [R-sig-Geo] Spatial panel Durbin model

2013-07-15 Thread Chen, Xi Prof.
Juan and Roger, thanks for the replies. If I construct the spatially lagged independent variables for each year period, add them to the panel dataset, and then use ML estimation from splm package, will it work? Xi From: Juan Tomas Sayago Gomez [mailto:juantomas.say...@gmail.com] Sent: Monday,

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] Package for fitting spatial model on large dataset

2013-07-15 Thread Moshood Agba Bakare
Thanks. What about gstat package? Moshood On Mon, Jul 15, 2013 at 12:51 PM, Paulo Justiniano Ribeiro Jr paulo...@c3sl.ufpr.br wrote: likelihood based methods and full kriging neighbourhood in some (or most) of geoR functions are not suitable for such data size Paulo Justiniano Ribeiro Jr

Re: [R-sig-Geo] Package for fitting spatial model on large dataset

2013-07-15 Thread Colin Rundel
You can look into using the spBayes package using gaussian predictive process models to handle the large number of observations. -Colin On Jul 15, 2013, at 1:08 PM, Moshood Agba Bakare bak...@ualberta.ca wrote: Dear All, Please which of these packages - gstat and geoR is better for fitting

[R-sig-Geo] Directions in northing-easting coordinates of spatial data

2013-07-15 Thread Moshood Agba Bakare
Dear All, I am always confused to identify north-west direction or south-west direction whenever I visualize my spatial data in 2-dimensional easting-northing coordinate space. Which point is actually the reference point? If the reference point is the intersection of easting and northing axes, am