Re: [R-sig-Geo] geoR Beta Estimates Not Equal to lm Coefficients: Why?

2007-09-02 Thread Edzer J. Pebesma
Joshua, My guess is that to estimate beta geoR uses generalized least squares and lm and SAS ordinary least squares. If you use a pure nugget model, or some other model with a range parameter sufficiently close to zero, i.e. model the observations as independent, the estimates should be the s

[R-sig-Geo] geoR Beta Estimates Not Equal to lm Coefficients: Why?

2007-09-02 Thread Joshua Palmer
Hello everyone, I am using the excellent geoR package to perform Kriging with External Drift. As to be expected, for any given set of covariates, trend coefficients I obtain performing least-squares regression in SAS equal those coefficients obtained by using R's lm function. However, those tren

Re: [R-sig-Geo] polygon buffers?

2007-09-02 Thread Agustin Lobo
Very interesting work, Paulo. Can you display from R on the TerraView window? I dream with a GIS tool (A GIS display tool for R? An R plugging for Qgis? etc) that could display spatial R objects. This is beyond current tools, based on exporting (hence duplicating space and troubles) R objects to

Re: [R-sig-Geo] polygon buffers?

2007-09-02 Thread Paulo Justiniano Ribeiro Jr
Peter There may be other alternatives, but this is one of the functionalities in an external software (terralib) that the aRT package tries to make usage. The aRT web page has one example: www.leg.ufpr.br/aRT Paulo Justiniano Ribeiro Jr LEG (Laboratório de Estatística e Geoinformação) Universida

Re: [R-sig-Geo] polygon buffers?

2007-09-02 Thread Tim Keitt
Hi Peter, I've been doing this in postgis. Once you get used to it, its not too difficult. Load data with shp2pgsql (or use rgdal + RODBC or similar if you want to stay in R) and then create a new table to hold the results. After that you can do "insert into restab select buffer(the_geom, dist) fr

[R-sig-Geo] polygon buffers?

2007-09-02 Thread Peter Adler
Dear R gurus, I am very excited about all the spatials tools in R because they have almost liberated me from ESRI products. One of the last things I have yet to figure out how to do in R is to create buffers around polygons (or just expand polygons). I searched the R-sig-geo archives (using Go

[R-sig-Geo] Calculate the angles for a spatial point patern?

2007-09-02 Thread zhijie zhang
Dear Rusers, I'd like to take the cases of cancer of the larynx in chorley(spatstat) to explain my question. I want to join the points of cancer of the larynx with the disused industrial incinerator to generate lines, and then calculate the angles of these line comparing the horizontal line?