spatial clusters

country<-wrld_simpl[wrld_simpl$ISO3%in%c("AZE", "ALB", "ARM", "BIH",

ske<-spCbind(country, rclasif)
ske at data
de <- data.frame(scale(ske at data[13]))   #col 13 contains the health

nbe <- poly2nb(country)
ce <- nbcosts(nbe, de)
le <- nb2listw(nbe, ce, style="B", zero.policy=TRUE)

On Sat, 18 Dec 2010, Roman Lu?trik wrote:

Perhaps you could follow the instructions from this

see if you can extract the desired country (assuming the resolution
suits you).

I don't think that we know which countries Dorina needs, or at which 
administrative level. If there are only 30 country observations, and the

countries are not contiguous, then the whole exercise seems unnecessary.

Can Dorina use the simplified world map in the maptools package:

library(maptools)
?wrld_simpl

adding the data using the ISO country names? If these are countries at 
some administrative division level, then maybe www.gadm.org will have

is needed, but the countries will have to be joined together, or the

map subsetted. If these are say EU NUTS2 regions, then use can be made

GISCO shapefiles from Eurostat. But we don't know what boundaries are 
needed, and of course to use Skater for clustering, you need the 
boundaries.

Roger

Cheers,
Roman



On Sat, Dec 18, 2010 at 6:20 PM, Dorina Lazar <drnlazar at yahoo.com>

Georg,

Thanks for being very patient with me.

I have not a shapefile for countries.

I have tried to obtain it. I have loaded a shapefile for Finland, GIS
data,
from

and I tried to read the file FI_1K.shp
Fin<- readShapePoly(system.file("etc/shapes/ FI_1K.shp",
package="spdep"))
but...
Error in getinfo.shape(filen) : Error opening SHP file

How can I obtain, for each country,  the shapefile spatial polygon in

(SpatialPolygons-DataFrame)?

Thanks again,
Dorina


--- On Thu, 12/16/10, Georg Ru? <research at georgruss.de> wrote:

From: Georg Ru? <research at georgruss.de>
Subject: Re: [R-sig-Geo] spatial clusters
To: "Dorina Lazar" <dorina.lazar at econ.ubbcluj.ro>
Cc: "R-sig-geo mailing list" <r-sig-geo at r-project.org>
Date: Thursday, December 16, 2010, 10:41 AM
On 16/12/10 09:06:33, dorina.lazar
wrote:

How to create a shapefile including my data about some

socio-demographic

indicators (7),  for about 30 countries around

the world (to be read with

readShapePoly)?

Hi Dorina,

I'm not quite sure what you have and what you want.

I assume that you'd like to create a shapefile such that
for every country
(which is represented as a polygon) there's a vector of
seven attributes
attached to it:

country1:  attr1 attr2 attr3 attr4 attr5 attr6 attr7
country2:  attr1 attr2 attr3 attr4 attr5 attr6 attr7
...

Do you have your countries as spatial polygons in R
(SpatialPolygons-
DataFrame or similar)? If yes, you should be able to add
your socio-demo-
graphic indicators to this data frame (just as in normal
data frames). In
this way you'd join the data you have inside R.

If you have the shapefile for the countries available and
would like to
edit this directly to add the indicators, I'd suggest using
a GIS like
GRASS http://grass.fbk.eu/ Afterwards, this can be easily
read into R.

Would be useful to suggest some good introductory

lectures in spatial

statistics (links or books).

I think the following is what Roger would suggest:

 author = {Bivand, Roger S. and Pebesma, Edzer J. and
G?mez-Rubio, Virgilio},
 title = {Applied Spatial Data Analysis with R},
 series = {Use R},
 publisher = {Springer},

http://www.asdar-book.org/
http://www.springerlink.com/content/uw07v1/
(Depending on your login/institute you may have fulltext
access, but the
book is definitely worth buying anyway.)

Regards,
Georg.
--
Research Assistant
Otto-von-Guericke-Universit?t Magdeburg
research at georgruss.de
http://research.georgruss.de