Hi All,

I want to aggregate my NUTS-3 regions in bigger ones, for which, of course, I have corresponding codes.
I have two main problems:

1) In terms of map visualization, I used unionSpatialPolygons(), starting from my original shape file. At this point I would like to write the resulting object out to a shape file by writePolyShape(), but R gives me the following answer:
class(x) == "SpatialPolygonsDataFrame" is not TRUE

2) In order to use spatial econometric techniques, I need to obtain a new neighbours list. Because of islands problems that have already been solved, rather than creating a new spatial weights matrix, I would like to obtain this by aggregating the one based on the original regions, or by aggregating directly the neighbour's set (a listw object for example), if possible.

Can anyone help on these two problems?

Thanks and best regards,
Roberto Patuelli



********************
Roberto Patuelli, Ph.D.
Post-doc researcher
Institute for Economic Research (IRE)
University of Lugano (USI)
via Maderno 24, CP 4361
CH-6904 Lugano
Switzerland
Phone: +41-(0)58-666-4166
Email: roberto.patuelli@lu.unisi.ch
********************
	[[alternative HTML version deleted]]


