On 13/02/16 11:59, Christopher W. Ryan wrote:
In spatstat, I want to create a single window from three adjacent
counties in New York State, USA. My original data is a shapefile,
"cty036" showing all the New York State county boundaries. Here's what
I've done so far:
== begin code
Hi Dominik,
If you use the gdalUtils package there is no significant difference in
the results using CLI or R:
library(gdalUtils)
gdaltransform(s_srs="+proj=longlat +datum=WGS84", t_srs="+proj=lcc
+lat_1=28 +lat_2=50 +lat_0=39.70001220694445 +lon_0=-98 +x_0=0 +y_0=0
+ellps=sphere
Sorry, forgot to add the details. rgdal and gdal 1.11.3 were installed from
kyngchaos.com
> sessionInfo()
R version 3.2.3 (2015-12-10)
Platform: x86_64-apple-darwin14.5.0 (64-bit)
Running under: OS X 10.11.3 (El Capitan)
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
Hi,
I'm struggling to use a custom projection. I am seeing differences with
someone using python proj4 bindings and when I compared my R results with
my commandline results I got even more confused. the coordinate
transformation is different for the two different methods.
could someone explain to
Hello,
You could have a look at the ClustGeo package :
>ClustGeo: Clustering of Observations with Geographical Constraints
>Functions which allow to integrate geographical constraints in Ward
hierarchical clustering. Geographical maps of typologies obtained can be
displayed with the use of
Hi Maryia,
What function are you using to fit the models? If you're using the fitting
routines in spdep, (lagsarlm, errorsarlm, sacsarlm, spautolm) then the AIC
prints for each model using the summary() function. You can also get a pseudo
R^2 using the summary(fit, Nagelkerke=T) option
for
Dear list members,
In moran scatter plot, I checked
[1] Slope of lm((spatial lagged variable) ~ (variable)): moran coefficient
(spdep:::moran.plot),
[2] Slope of lm((spatial lagged standardized variable) ~ (standardized
variable)): moran coefficient, and
[3] Slope of lm((standardized spatial
In spatstat, I want to create a single window from three adjacent
counties in New York State, USA. My original data is a shapefile,
"cty036" showing all the New York State county boundaries. Here's what
I've done so far:
== begin code =
library(shapefiles)
library(spatstat)
Dear Roger,
thanks for your answer!
> > Dear list members,
> >
> > I'm currently estimating a fixed effects panel model and I want to
> > control for spatial dependence. Thus, I also estimated two spatial
> > fe-models, one with a spatial error term and one with a spatial error
> > term and
