Issue with shapefile and Stata

[Dario Maimone Ansaldo Patti]

Dear All,

I have a shapefile which contains information about some European regions (a picture is attached). Some areas have been merged together, such that some "macro-regions are created". I use:

spshape2dta "\NUTS_RG_01M_2013

to export information in Stata.

I have a dataset which contains information at grid level (lat and lon are provided), covering a larger number of regions/state than those I am working on. So I try to map those lat and lon into my regions using:

geoinpoly _Y _X using NUTS_RG_01M_2013_shp

According to my understanding, geoinpoly should capture only those coordinates that correspond to my regions, excluding the others. Morever, it should use the same level of aggregation I have in my dataset. But if so, why when I use geoinpoly it maps also regions/countries outside my shapefile (for instance Austria) and with a different level of aggregation (I have region BE1 in my shapefile, while after geopoly I find BE11, BE12 and so on)?

Am I missing something?

Thanks in advance for any help you may want to provide.

Best,

Dario

[Manh Hoang Ba]

Maybe you haven't saved the new shapefile with the merged regions.
I use QGIS to split/merge the features (regions) and save it as a new shapefile. It works. You can try it.

[Dario Maimone Ansaldo Patti]

Manh Hoang Ba The shapefile is saved correctly, but I am not sure about auxiliary files. I will double check this

[Dario Maimone Ansaldo Patti]

Solved apparently I did not change the auxiliary files. I have done this through QGIS and now it works