{ "cells": [ { "cell_type": "markdown", "id": "8f1abf17", "metadata": {}, "source": [ "# Feature importance per signature type\n", "\n", "This notebooks analyses which characters are more important for each individual signature type. In other words, what makes each cluster unique compared to all the other." ] }, { "cell_type": "code", "execution_count": 110, "id": "d86c027b", "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import pandas as pd\n", "import geopandas as gpd\n", "import dask.dataframe\n", "import matplotlib.pyplot as plt\n", "import urbangrammar_graphics as ugg\n", "import seaborn as sns\n", "\n", "from matplotlib.lines import Line2D\n", "from sklearn.ensemble import RandomForestClassifier" ] }, { "cell_type": "code", "execution_count": 3, "id": "3ba4607d", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "CPU times: user 19.2 s, sys: 3.26 s, total: 22.5 s\n", "Wall time: 21.4 s\n", "CPU times: user 72.3 ms, sys: 4.09 ms, total: 76.4 ms\n", "Wall time: 116 ms\n", "CPU times: user 39.6 ms, sys: 7.9 ms, total: 47.5 ms\n", "Wall time: 41.4 ms\n", "CPU times: user 18.9 ms, sys: 0 ns, total: 18.9 ms\n", "Wall time: 18.7 ms\n", "CPU times: user 2min 41s, sys: 1min 27s, total: 4min 8s\n", "Wall time: 2min 45s\n" ] } ], "source": [ "%time standardized_form = dask.dataframe.read_parquet(\"../../urbangrammar_samba/spatial_signatures/clustering_data/form/standardized/\").set_index('hindex')\n", "%time stand_fn = dask.dataframe.read_parquet(\"../../urbangrammar_samba/spatial_signatures/clustering_data/function/standardized/\")\n", "%time data = dask.dataframe.multi.concat([standardized_form, stand_fn], axis=1).replace([np.inf, -np.inf], np.nan).fillna(0)\n", "%time data = data.drop(columns=[\"keep_q1\", \"keep_q2\", \"keep_q3\"])\n", "%time data = data.compute()" ] }, { "cell_type": "code", "execution_count": 4, "id": "8806704b", "metadata": {}, "outputs": [], "source": [ "labels_l1 = pd.read_parquet(\"../../urbangrammar_samba/spatial_signatures/clustering_data/KMeans10GB.pq\")\n", "labels_l2_9 = pd.read_parquet(\"../../urbangrammar_samba/spatial_signatures/clustering_data/clustergram_cl9_labels.pq\")\n", "labels_l2_2 = pd.read_parquet(\"../../urbangrammar_samba/spatial_signatures/clustering_data/subclustering_cluster2_k3.pq\")\n", "labels = labels_l1.copy()\n", "labels.loc[labels.kmeans10gb == 9, 'kmeans10gb'] = labels_l2_9['9'].values + 90\n", "labels.loc[labels.kmeans10gb == 2, 'kmeans10gb'] = labels_l2_2['subclustering_cluster2_k3'].values + 20\n", "\n", "outliers = [98, 93, 96, 97]\n", "mask = ~labels.kmeans10gb.isin(outliers)" ] }, { "cell_type": "markdown", "id": "9d549248", "metadata": {}, "source": [ "## Feature importance per cluster" ] }, { "cell_type": "code", "execution_count": 15, "id": "49e9d0f6", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "array([ 4, 0, 6, 1, 21, 7, 3, 5, 90, 20, 8, 22, 98, 92, 94, 91, 95,\n", " 96, 93, 97], dtype=int32)" ] }, "execution_count": 15, "metadata": {}, "output_type": "execute_result" } ], "source": [ "labels.kmeans10gb.unique()" ] }, { "cell_type": "code", "execution_count": 19, "id": "900223ad", "metadata": {}, "outputs": [], "source": [ "imps = pd.DataFrame()" ] }, { "cell_type": "code", "execution_count": 21, "id": "76b85d70", "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "[Parallel(n_jobs=-1)]: Using backend ThreadingBackend with 16 concurrent workers.\n", "[Parallel(n_jobs=-1)]: Done 10 out of 10 | elapsed: 13.6min finished\n", "[Parallel(n_jobs=-1)]: Using backend ThreadingBackend with 16 concurrent workers.\n", "[Parallel(n_jobs=-1)]: Done 10 out of 10 | elapsed: 10.4min finished\n", "[Parallel(n_jobs=-1)]: Using backend ThreadingBackend with 16 concurrent workers.\n", "[Parallel(n_jobs=-1)]: Done 10 out of 10 | elapsed: 5.8min finished\n", "[Parallel(n_jobs=-1)]: Using backend ThreadingBackend with 16 concurrent workers.\n", "[Parallel(n_jobs=-1)]: Done 10 out of 10 | elapsed: 6.6min finished\n", "[Parallel(n_jobs=-1)]: Using backend ThreadingBackend with 16 concurrent workers.\n", "[Parallel(n_jobs=-1)]: Done 10 out of 10 | elapsed: 6.0min finished\n", "[Parallel(n_jobs=-1)]: Using backend ThreadingBackend with 16 concurrent workers.\n", "[Parallel(n_jobs=-1)]: Done 10 out of 10 | elapsed: 7.9min finished\n", "[Parallel(n_jobs=-1)]: Using backend ThreadingBackend with 16 concurrent workers.\n", "[Parallel(n_jobs=-1)]: Done 10 out of 10 | elapsed: 6.5min finished\n", "[Parallel(n_jobs=-1)]: Using backend ThreadingBackend with 16 concurrent workers.\n", "[Parallel(n_jobs=-1)]: Done 10 out of 10 | elapsed: 6.6min finished\n", "[Parallel(n_jobs=-1)]: Using backend ThreadingBackend with 16 concurrent workers.\n", "[Parallel(n_jobs=-1)]: Done 10 out of 10 | elapsed: 6.4min finished\n", "[Parallel(n_jobs=-1)]: Using backend ThreadingBackend with 16 concurrent workers.\n", "[Parallel(n_jobs=-1)]: Done 10 out of 10 | elapsed: 7.3min finished\n", "[Parallel(n_jobs=-1)]: Using backend ThreadingBackend with 16 concurrent workers.\n", "[Parallel(n_jobs=-1)]: Done 10 out of 10 | elapsed: 5.9min finished\n", "[Parallel(n_jobs=-1)]: Using backend ThreadingBackend with 16 concurrent workers.\n", "[Parallel(n_jobs=-1)]: Done 10 out of 10 | elapsed: 7.5min finished\n", "[Parallel(n_jobs=-1)]: Using backend ThreadingBackend with 16 concurrent workers.\n", "[Parallel(n_jobs=-1)]: Done 10 out of 10 | elapsed: 6.3min finished\n", "[Parallel(n_jobs=-1)]: Using backend ThreadingBackend with 16 concurrent workers.\n", "[Parallel(n_jobs=-1)]: Done 10 out of 10 | elapsed: 3.4min finished\n", "[Parallel(n_jobs=-1)]: Using backend ThreadingBackend with 16 concurrent workers.\n", "[Parallel(n_jobs=-1)]: Done 10 out of 10 | elapsed: 1.9min finished\n", "[Parallel(n_jobs=-1)]: Using backend ThreadingBackend with 16 concurrent workers.\n", "[Parallel(n_jobs=-1)]: Done 10 out of 10 | elapsed: 2.0min finished\n" ] } ], "source": [ "for cluster in labels.kmeans10gb.unique():\n", " if cluster not in outliers:\n", " cluster_bool = labels.loc[mask]['kmeans10gb'].apply(lambda x: 1 if x == cluster else 0)\n", "\n", " clf = RandomForestClassifier(n_estimators=10, n_jobs=-1, random_state=42, verbose=1)\n", " clf = clf.fit(data.loc[mask].values, cluster_bool.values)\n", "\n", " importances = pd.Series(clf.feature_importances_, index=data.columns).sort_values(ascending=False)\n", "\n", " imps[f'cluster_{cluster}'] = importances.head(50).index.values\n", " imps[f'cluster_{cluster}_vals'] = importances.head(50).values" ] }, { "cell_type": "code", "execution_count": 26, "id": "2387324c", "metadata": { "tags": [] }, "outputs": [ { "data": { "text/html": [ "
\n", " | cluster_0 | \n", "cluster_1 | \n", "cluster_20 | \n", "cluster_21 | \n", "cluster_22 | \n", "cluster_3 | \n", "cluster_4 | \n", "cluster_5 | \n", "cluster_6 | \n", "cluster_7 | \n", "cluster_8 | \n", "cluster_90 | \n", "cluster_91 | \n", "cluster_92 | \n", "cluster_94 | \n", "cluster_95 | \n", "
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
0 | \n", "sicCAR_q1 | \n", "ltcWRE_q3 | \n", "ssbCCM_q2 | \n", "stbCeA_q1 | \n", "sdbPer_q2 | \n", "ldsAre_q1 | \n", "sdcLAL_q1 | \n", "ssbElo_q1 | \n", "lcnClo_q3 | \n", "mdcAre_q2 | \n", "linPDE_q2 | \n", "sdbPer_q2 | \n", "sdbAre_q1 | \n", "ssbCCD_q2 | \n", "ssbERI_q2 | \n", "sicCAR_q2 | \n", "
1 | \n", "sicCAR_q2 | \n", "ltcRea_q3 | \n", "ssbCCM_q3 | \n", "linP4W_q2 | \n", "ssbCCM_q2 | \n", "mdsAre_q1 | \n", "sicCAR_q2 | \n", "ssbCCM_q3 | \n", "lcnClo_q2 | \n", "sicCAR_q2 | \n", "lcdMes_q3 | \n", "ssbERI_q1 | \n", "G, I. Distribution, hotels and restaurants_q2 | \n", "ssbCCM_q2 | \n", "ssbCCM_q2 | \n", "C. Manufacturing_q2 | \n", "
2 | \n", "ltbIBD_q2 | \n", "ltcAre_q2 | \n", "sdbAre_q3 | \n", "stbCeA_q2 | \n", "sdbPer_q3 | \n", "sddAre_q2 | \n", "sicCAR_q1 | \n", "ssbElo_q2 | \n", "ldePer_q1 | \n", "mtdMDi_q2 | \n", "lcdMes_q2 | \n", "ssbCCM_q2 | \n", "K, L, M, N. Financial, real estate, profession... | \n", "ssbSqu_q3 | \n", "ssbCCD_q2 | \n", "R, S, T, U. Other_q2 | \n", "
3 | \n", "sdcAre_q2 | \n", "sdcAre_q2 | \n", "population_q3 | \n", "ldeAre_q2 | \n", "sdbAre_q2 | \n", "sicCAR_q2 | \n", "sdcAre_q2 | \n", "ssbCCo_q1 | \n", "ldeAre_q2 | \n", "sicCAR_q1 | \n", "ssbERI_q1 | \n", "ssbSqu_q3 | \n", "R, S, T, U. Other_q2 | \n", "K, L, M, N. Financial, real estate, profession... | \n", "ssbCor_q2 | \n", "G, I. Distribution, hotels and restaurants_q2 | \n", "
4 | \n", "sddAre_q2 | \n", "misCel_q1 | \n", "sdbPer_q2 | \n", "stcOri_q2 | \n", "population_q3 | \n", "lddNDe_q3 | \n", "mtcWNe_q3 | \n", "ssbCCD_q2 | \n", "lcnClo_q1 | \n", "mdsAre_q1 | \n", "ssbCCo_q1 | \n", "sdbAre_q2 | \n", "G, I. Distribution, hotels and restaurants_q3 | \n", "sdbPer_q2 | \n", "K, L, M, N. Financial, real estate, profession... | \n", "sicCAR_q1 | \n", "
5 | \n", "mtdMDi_q2 | \n", "misCel_q2 | \n", "sdbAre_q2 | \n", "ssbERI_q1 | \n", "ssbSqu_q3 | \n", "mdsAre_q2 | \n", "mdsAre_q1 | \n", "sdbPer_q3 | \n", "ltcWRE_q3 | \n", "sscCCo_q2 | \n", "population_q1 | \n", "ssbCCD_q3 | \n", "K, L, M, N. Financial, real estate, profession... | \n", "sdbPer_q3 | \n", "G, I. Distribution, hotels and restaurants_q2 | \n", "C. Manufacturing_q3 | \n", "
6 | \n", "mdsAre_q1 | \n", "lisCel_q2 | \n", "ldePer_q1 | \n", "linP4W_q1 | \n", "ssbCCD_q3 | \n", "sicCAR_q1 | \n", "ltcAre_q1 | \n", "sdsSPW_q2 | \n", "linP4W_q2 | \n", "mdcAre_q1 | \n", "ssbElo_q2 | \n", "K, L, M, N. Financial, real estate, profession... | \n", "C. Manufacturing_q2 | \n", "sdbAre_q2 | \n", "sdbPer_q2 | \n", "ssbCCM_q2 | \n", "
7 | \n", "Code_18_112_q2 | \n", "mtcWNe_q1 | \n", "lteOri_q2 | \n", "ldePer_q1 | \n", "K, L, M, N. Financial, real estate, profession... | \n", "ldeAre_q2 | \n", "sdcAre_q3 | \n", "ssbCCo_q2 | \n", "sddAre_q1 | \n", "sisBpM_q2 | \n", "mdsAre_q2 | \n", "G, I. Distribution, hotels and restaurants_q2 | \n", "sdbPer_q2 | \n", "G, I. Distribution, hotels and restaurants_q2 | \n", "ssbSqu_q3 | \n", "sdbPer_q2 | \n", "
8 | \n", "mtcWNe_q2 | \n", "ltcAre_q1 | \n", "sdbPer_q3 | \n", "linPDE_q1 | \n", "ssbERI_q1 | \n", "lseCWA_q3 | \n", "mtbNDi_q2 | \n", "ltcAre_q1 | \n", "sddAre_q2 | \n", "ltcAre_q1 | \n", "sdsAre_q3 | \n", "sdbPer_q3 | \n", "ssbCCM_q2 | \n", "ssbCor_q3 | \n", "K, L, M, N. Financial, real estate, profession... | \n", "sdsSPO_q2 | \n", "
9 | \n", "sdcLAL_q2 | \n", "lisCel_q1 | \n", "ldeAre_q1 | \n", "lteOri_q1 | \n", "R, S, T, U. Other_q2 | \n", "sdcAre_q2 | \n", "ltbIBD_q2 | \n", "sdbPer_q2 | \n", "ltcWRE_q2 | \n", "sddAre_q3 | \n", "ssbCCo_q2 | \n", "sdbAre_q1 | \n", "Code_18_211_q1 | \n", "ssbCCD_q3 | \n", "ssbCCM_q1 | \n", "mean_q3 | \n", "
10 | \n", "mtbNDi_q2 | \n", "sicCAR_q1 | \n", "stcOri_q3 | \n", "sdbPer_q3 | \n", "sdbPer_q1 | \n", "lisCel_q1 | \n", "ltcWRE_q1 | \n", "mtbNDi_q3 | \n", "lcdMes_q3 | \n", "mdsAre_q3 | \n", "linPDE_q1 | \n", "R, S, T, U. Other_q3 | \n", "ssbCCD_q1 | \n", "ssbERI_q1 | \n", "ssbCCD_q3 | \n", "ssbSqu_q1 | \n", "
11 | \n", "ltcAre_q1 | \n", "ssbERI_q1 | \n", "stbCeA_q2 | \n", "stbOri_q2 | \n", "sdbAre_q3 | \n", "ssbElo_q2 | \n", "mtcWNe_q2 | \n", "ltbIBD_q3 | \n", "ldeAre_q1 | \n", "sdsAre_q3 | \n", "stcOri_q2 | \n", "ssbCCD_q2 | \n", "G, I. Distribution, hotels and restaurants_q1 | \n", "ssbSqu_q2 | \n", "F. Construction_q3 | \n", "sdbAre_q1 | \n", "
12 | \n", "mdcAre_q2 | \n", "sdcLAL_q3 | \n", "ssbCCo_q1 | \n", "ssbCCo_q3 | \n", "O,P,Q. Public administration, education and he... | \n", "ltcWRE_q1 | \n", "sdcLAL_q2 | \n", "ltcWRE_q2 | \n", "linP4W_q1 | \n", "mdsAre_q2 | \n", "stbOri_q2 | \n", "H, J. Transport and communication_q1 | \n", "sscCCo_q3 | \n", "O,P,Q. Public administration, education and he... | \n", "F. Construction_q2 | \n", "H, J. Transport and communication_q3 | \n", "
13 | \n", "sdsAre_q2 | \n", "ltcRea_q1 | \n", "stcOri_q2 | \n", "ssbCCD_q2 | \n", "ssbCCM_q3 | \n", "stbSAl_q2 | \n", "mtbNDi_q3 | \n", "population_q3 | \n", "sdsAre_q2 | \n", "mtbNDi_q3 | \n", "ssbCCD_q2 | \n", "ssbERI_q2 | \n", "F. Construction_q3 | \n", "F. Construction_q3 | \n", "K, L, M, N. Financial, real estate, profession... | \n", "R, S, T, U. Other_q1 | \n", "
14 | \n", "mtcWNe_q3 | \n", "sisBpM_q3 | \n", "stbOri_q2 | \n", "sdbAre_q2 | \n", "mean_q1 | \n", "stbCeA_q1 | \n", "mdcAre_q2 | \n", "ltcWRE_q1 | \n", "stbCeA_q2 | \n", "sdcAre_q2 | \n", "ssbElo_q1 | \n", "ssbCor_q3 | \n", "H, J. Transport and communication_q2 | \n", "ssbCCM_q3 | \n", "sicCAR_q3 | \n", "ssbCCD_q1 | \n", "
15 | \n", "sisBpM_q3 | \n", "ldsMSL_q1 | \n", "stbOri_q1 | \n", "population_q2 | \n", "ssbCCM_q1 | \n", "sddAre_q3 | \n", "sisBpM_q3 | \n", "sdsSPW_q3 | \n", "ldsCDL_q3 | \n", "sscERI_q2 | \n", "mdsAre_q3 | \n", "ssbSqu_q2 | \n", "sicCAR_q3 | \n", "ssbERI_q2 | \n", "mtbAli_q2 | \n", "sisBpM_q3 | \n", "
16 | \n", "sicCAR_q3 | \n", "sdsSWD_q3 | \n", "lcdMes_q2 | \n", "ssbCCM_q3 | \n", "ssbERI_q2 | \n", "ltcWRE_q2 | \n", "sisBpM_q1 | \n", "ltbIBD_q1 | \n", "mtdMDi_q1 | \n", "sscERI_q1 | \n", "lteOri_q1 | \n", "G, I. Distribution, hotels and restaurants_q3 | \n", "ssbERI_q2 | \n", "ssbCor_q2 | \n", "R, S, T, U. Other_q3 | \n", "sdcAre_q1 | \n", "
17 | \n", "sdcLAL_q1 | \n", "ssbCCo_q1 | \n", "sicCAR_q3 | \n", "lteWNB_q2 | \n", "ssbCCD_q2 | \n", "lisCel_q2 | \n", "ltcWRE_q2 | \n", "sicCAR_q2 | \n", "mtcWNe_q3 | \n", "linWID_q1 | \n", "mtcWNe_q1 | \n", "misCel_q3 | \n", "F. Construction_q2 | \n", "R, S, T, U. Other_q3 | \n", "sdcAre_q1 | \n", "O,P,Q. Public administration, education and he... | \n", "
18 | \n", "mtdMDi_q1 | \n", "ssbCCM_q2 | \n", "stcOri_q1 | \n", "stbOri_q1 | \n", "K, L, M, N. Financial, real estate, profession... | \n", "lddNDe_q2 | \n", "ltcAre_q2 | \n", "night_lights_q1 | \n", "lcdMes_q2 | \n", "population_q1 | \n", "sscCCo_q2 | \n", "sdbAre_q3 | \n", "A, B, D, E. Agriculture, energy and water_q3 | \n", "sdbAre_q3 | \n", "mean_q1 | \n", "H, J. Transport and communication_q2 | \n", "
19 | \n", "ldsAre_q2 | \n", "sdcLAL_q1 | \n", "stbCeA_q1 | \n", "ssbElo_q3 | \n", "sdbAre_q1 | \n", "stbCeA_q2 | \n", "sscERI_q2 | \n", "mtcWNe_q2 | \n", "sdcLAL_q1 | \n", "sdcAre_q3 | \n", "stcOri_q1 | \n", "O,P,Q. Public administration, education and he... | \n", "lseCWA_q1 | \n", "mean_q1 | \n", "ssbSqu_q1 | \n", "A, B, D, E. Agriculture, energy and water_q2 | \n", "
20 | \n", "sdsSPO_q1 | \n", "linPDE_q3 | \n", "mean_q1 | \n", "ssbCCM_q2 | \n", "sicCAR_q3 | \n", "linWID_q1 | \n", "mtdMDi_q3 | \n", "linP3W_q3 | \n", "mtdMDi_q2 | \n", "population_q3 | \n", "stcOri_q3 | \n", "misCel_q2 | \n", "linP4W_q3 | \n", "K, L, M, N. Financial, real estate, profession... | \n", "H, J. Transport and communication_q2 | \n", "ltcAre_q3 | \n", "
21 | \n", "ltbIBD_q1 | \n", "misCel_q3 | \n", "stbOri_q3 | \n", "lcdMes_q3 | \n", "R, S, T, U. Other_q3 | \n", "sssLin_q2 | \n", "sdsSPO_q2 | \n", "sdsSWD_q3 | \n", "ltcWRE_q1 | \n", "mtdMDi_q3 | \n", "population_q2 | \n", "mean_q1 | \n", "R, S, T, U. Other_q1 | \n", "G, I. Distribution, hotels and restaurants_q3 | \n", "ssbSqu_q2 | \n", "ldsCDL_q2 | \n", "
22 | \n", "sdsLen_q2 | \n", "sdcAre_q3 | \n", "ldeAre_q2 | \n", "lcdMes_q2 | \n", "G, I. Distribution, hotels and restaurants_q2 | \n", "ssbCCo_q1 | \n", "sdcAre_q1 | \n", "ssbERI_q1 | \n", "ldsAre_q2 | \n", "sdsSPO_q2 | \n", "lcdMes_q1 | \n", "ldePer_q2 | \n", "K, L, M, N. Financial, real estate, profession... | \n", "sicCAR_q3 | \n", "R, S, T, U. Other_q1 | \n", "mdcAre_q2 | \n", "
23 | \n", "mtbNDi_q3 | \n", "lisCel_q3 | \n", "sddAre_q3 | \n", "ltcWRE_q3 | \n", "night_lights_q3 | \n", "mtdMDi_q3 | \n", "sisBpM_q2 | \n", "ssbCCo_q3 | \n", "sdcAre_q1 | \n", "stbCeA_q1 | \n", "stbOri_q3 | \n", "K, L, M, N. Financial, real estate, profession... | \n", "A, B, D, E. Agriculture, energy and water_q2 | \n", "ssbCCD_q1 | \n", "C. Manufacturing_q2 | \n", "H, J. Transport and communication_q1 | \n", "
24 | \n", "lddNDe_q2 | \n", "sdsSWD_q2 | \n", "lteOri_q1 | \n", "lteOri_q2 | \n", "ssbCor_q3 | \n", "ldePer_q1 | \n", "Code_18_231_q1 | \n", "linP3W_q2 | \n", "population_q2 | \n", "sdcLAL_q2 | \n", "ssbCCo_q3 | \n", "ssbCCM_q1 | \n", "linP4W_q1 | \n", "misCel_q3 | \n", "ssbCCM_q3 | \n", "G, I. Distribution, hotels and restaurants_q1 | \n", "
25 | \n", "sisBpM_q2 | \n", "stbSAl_q2 | \n", "K, L, M, N. Financial, real estate, profession... | \n", "G, I. Distribution, hotels and restaurants_q3 | \n", "G, I. Distribution, hotels and restaurants_q3 | \n", "stbSAl_q3 | \n", "ltcAre_q3 | \n", "sicCAR_q3 | \n", "lseCWA_q1 | \n", "linWID_q2 | \n", "ldsCDL_q1 | \n", "ssbCCM_q3 | \n", "H, J. Transport and communication_q3 | \n", "sscERI_q3 | \n", "sdbAre_q1 | \n", "ldsAre_q1 | \n", "
26 | \n", "ldsAre_q1 | \n", "ssbElo_q1 | \n", "ssbCCD_q3 | \n", "linPDE_q3 | \n", "population_q2 | \n", "misCel_q2 | \n", "ldsCDL_q1 | \n", "mtcWNe_q3 | \n", "linWID_q2 | \n", "ltcWRE_q2 | \n", "ssbCCM_q3 | \n", "H, J. Transport and communication_q3 | \n", "linP3W_q3 | \n", "R, S, T, U. Other_q2 | \n", "lseCWA_q2 | \n", "ldeAre_q1 | \n", "
27 | \n", "sdsAre_q1 | \n", "ssbCCo_q2 | \n", "ltcAre_q1 | \n", "population_q3 | \n", "night_lights_q2 | \n", "ssbERI_q1 | \n", "mdsAre_q3 | \n", "stbCeA_q1 | \n", "ssbCCM_q1 | \n", "sscCCo_q1 | \n", "linPDE_q3 | \n", "K, L, M, N. Financial, real estate, profession... | \n", "R, S, T, U. Other_q3 | \n", "H, J. Transport and communication_q3 | \n", "Code_18_111_q2 | \n", "listed_counts | \n", "
28 | \n", "sdcAre_q3 | \n", "stbCeA_q1 | \n", "linP4W_q2 | \n", "stcOri_q3 | \n", "O,P,Q. Public administration, education and he... | \n", "ldePer_q2 | \n", "mean_q3 | \n", "linP4W_q2 | \n", "lteWNB_q1 | \n", "sddAre_q2 | \n", "sscCCo_q1 | \n", "ssbElo_q1 | \n", "sicCAR_q2 | \n", "sdcAre_q3 | \n", "sdbAre_q2 | \n", "K, L, M, N. Financial, real estate, profession... | \n", "
29 | \n", "ldsMSL_q2 | \n", "mdcAre_q2 | \n", "linP4W_q3 | \n", "lcdMes_q1 | \n", "ltcWRE_q2 | \n", "ldsMSL_q2 | \n", "Code_18_231_q2 | \n", "ltcWRE_q3 | \n", "linP4W_q3 | \n", "ltcWRE_q1 | \n", "stbOri_q1 | \n", "Code_18_111_q3 | \n", "ssbCCM_q1 | \n", "Code_18_111_q2 | \n", "G, I. Distribution, hotels and restaurants_q3 | \n", "mtbNDi_q3 | \n", "
30 | \n", "mdsAre_q2 | \n", "ldeAre_q2 | \n", "ssbSqu_q1 | \n", "misCel_q2 | \n", "sicCAR_q2 | \n", "linPDE_q3 | \n", "ldsCDL_q2 | \n", "linPDE_q2 | \n", "linPDE_q2 | \n", "ssbCCM_q3 | \n", "linP4W_q2 | \n", "listed_counts | \n", "sdsSWD_q1 | \n", "ssbSqu_q1 | \n", "ssbCCD_q1 | \n", "mtcWNe_q3 | \n", "
31 | \n", "ltcAre_q3 | \n", "linPDE_q2 | \n", "lcnClo_q2 | \n", "stbOri_q3 | \n", "lcnClo_q3 | \n", "linWID_q2 | \n", "Code_18_211_q3 | \n", "ldePer_q1 | \n", "ldePer_q2 | \n", "sdsSPO_q1 | \n", "ssbElo_q3 | \n", "H, J. Transport and communication_q2 | \n", "ssbSqu_q1 | \n", "H, J. Transport and communication_q2 | \n", "R, S, T, U. Other_q2 | \n", "stcOri_q2 | \n", "
32 | \n", "sdsSPO_q2 | \n", "ssbCCo_q3 | \n", "R, S, T, U. Other_q3 | \n", "stcOri_q1 | \n", "ssbCCo_q3 | \n", "ssbCCo_q2 | \n", "sicCAR_q3 | \n", "stbCeA_q3 | \n", "sdsSPW_q3 | \n", "night_lights_q1 | \n", "lteOri_q2 | \n", "R, S, T, U. Other_q2 | \n", "mdcAre_q3 | \n", "Code_18_111_q3 | \n", "lisCel_q2 | \n", "sisBpM_q2 | \n", "
33 | \n", "mdcAre_q1 | \n", "sisBpM_q2 | \n", "ssbERI_q1 | \n", "linP4W_q3 | \n", "ssbCCo_q2 | \n", "lseERI_q2 | \n", "night_lights_q2 | \n", "ldsCDL_q1 | \n", "linPDE_q3 | \n", "ssbCCo_q1 | \n", "ssbSqu_q3 | \n", "O,P,Q. Public administration, education and he... | \n", "ltcAre_q1 | \n", "lseCWA_q1 | \n", "H, J. Transport and communication_q3 | \n", "ssbCCD_q3 | \n", "
34 | \n", "lddNDe_q3 | \n", "stbCeA_q3 | \n", "sdcAre_q2 | \n", "ssbSqu_q1 | \n", "linP4W_q2 | \n", "sdsSPO_q3 | \n", "Code_18_211_q2 | \n", "mdcAre_q2 | \n", "ssbERI_q1 | \n", "sdsLen_q1 | \n", "lddNDe_q2 | \n", "population_q3 | \n", "O,P,Q. Public administration, education and he... | \n", "ssbCCM_q1 | \n", "mtbNDi_q3 | \n", "sscCCo_q3 | \n", "
35 | \n", "ltcAre_q2 | \n", "lcdMes_q2 | \n", "linPDE_q3 | \n", "ldeAre_q1 | \n", "ldePer_q1 | \n", "stbCeA_q3 | \n", "mtbNDi_q1 | \n", "listed_counts | \n", "sdsSWD_q3 | \n", "ltbIBD_q2 | \n", "linWID_q1 | \n", "night_lights_q3 | \n", "sdcLAL_q3 | \n", "ltcAre_q1 | \n", "sicCAR_q1 | \n", "sssLin_q3 | \n", "
36 | \n", "mean_q1 | \n", "stcOri_q3 | \n", "O,P,Q. Public administration, education and he... | \n", "lcnClo_q3 | \n", "population_q1 | \n", "ssbCCD_q2 | \n", "ldsCDL_q3 | \n", "night_lights_q2 | \n", "ssbCCo_q2 | \n", "night_lights_q2 | \n", "sisBpM_q2 | \n", "population_q2 | \n", "ldsCDL_q1 | \n", "G, I. Distribution, hotels and restaurants_q1 | \n", "A, B, D, E. Agriculture, energy and water_q2 | \n", "sdsSWD_q1 | \n", "
37 | \n", "sscERI_q1 | \n", "ssbCCM_q3 | \n", "sicCAR_q2 | \n", "sdbAre_q3 | \n", "ldeAre_q1 | \n", "sdsAre_q2 | \n", "mean_q2 | \n", "ldeAre_q1 | \n", "sdbAre_q1 | \n", "ltbIBD_q3 | \n", "ssbSqu_q2 | \n", "R, S, T, U. Other_q1 | \n", "sdsSPW_q3 | \n", "ssbElo_q2 | \n", "O,P,Q. Public administration, education and he... | \n", "K, L, M, N. Financial, real estate, profession... | \n", "
38 | \n", "night_lights_q2 | \n", "sisBpM_q1 | \n", "sdsSPW_q2 | \n", "ssbCCo_q2 | \n", "night_lights_q1 | \n", "stbSAl_q1 | \n", "sdsAre_q3 | \n", "sdsSPW_q1 | \n", "mtdMDi_q3 | \n", "stbCeA_q2 | \n", "mtdMDi_q2 | \n", "Code_18_111_q2 | \n", "ssbERI_q1 | \n", "K, L, M, N. Financial, real estate, profession... | \n", "H, J. Transport and communication_q1 | \n", "lseERI_q3 | \n", "
39 | \n", "sscCCo_q3 | \n", "ssbElo_q3 | \n", "stbCeA_q3 | \n", "sdbPer_q2 | \n", "ssbSqu_q2 | \n", "ssbElo_q1 | \n", "mdcAre_q1 | \n", "sdsAre_q3 | \n", "sdcLAL_q2 | \n", "sdsAre_q2 | \n", "F. Construction_q1 | \n", "ssbSqu_q1 | \n", "lisCel_q2 | \n", "sdbAre_q1 | \n", "sdbAre_q3 | \n", "G, I. Distribution, hotels and restaurants_q3 | \n", "
40 | \n", "ltbIBD_q3 | \n", "ltcRea_q2 | \n", "linP4W_q1 | \n", "lisCel_q2 | \n", "G, I. Distribution, hotels and restaurants_q1 | \n", "mdcAre_q3 | \n", "sdsSPO_q1 | \n", "stbOri_q2 | \n", "linWID_q1 | \n", "mtbNDi_q2 | \n", "ssbERI_q2 | \n", "ssbCCo_q1 | \n", "ssbSqu_q3 | \n", "ssbERI_q3 | \n", "sdbPer_q1 | \n", "R, S, T, U. Other_q3 | \n", "
41 | \n", "stbCeA_q1 | \n", "ssbElo_q2 | \n", "night_lights_q1 | \n", "ssbCCo_q1 | \n", "lcnClo_q2 | \n", "lseCWA_q2 | \n", "ltbIBD_q1 | \n", "stcOri_q3 | \n", "misCel_q2 | \n", "mean_q2 | \n", "lddNDe_q3 | \n", "sscERI_q2 | \n", "mtcWNe_q2 | \n", "listed_counts | \n", "lcnClo_q3 | \n", "ssbElo_q2 | \n", "
42 | \n", "mtbNDi_q1 | \n", "ltcWRE_q1 | \n", "stbSAl_q2 | \n", "ssbElo_q2 | \n", "ssbCCo_q1 | \n", "ssbSqu_q2 | \n", "sdcLAL_q3 | \n", "sddAre_q3 | \n", "stcOri_q3 | \n", "sscCCo_q3 | \n", "lddNDe_q1 | \n", "mean_q3 | \n", "sdbAre_q3 | \n", "lcnClo_q1 | \n", "sicCAR_q2 | \n", "F. Construction_q2 | \n", "
43 | \n", "sisBpM_q1 | \n", "stbSAl_q1 | \n", "G, I. Distribution, hotels and restaurants_q2 | \n", "mdcAre_q2 | \n", "F. Construction_q2 | \n", "ssbElo_q3 | \n", "Code_18_231_q3 | \n", "ssbElo_q3 | \n", "sssLin_q2 | \n", "sdbPer_q3 | \n", "mdsAre_q1 | \n", "night_lights_q1 | \n", "lseCCo_q3 | \n", "sdbPer_q1 | \n", "mtcWNe_q3 | \n", "lieWCe_q1 | \n", "
44 | \n", "mean_q2 | \n", "mtdMDi_q1 | \n", "night_lights_q2 | \n", "night_lights_q1 | \n", "K, L, M, N. Financial, real estate, profession... | \n", "mdsAre_q3 | \n", "K, L, M, N. Financial, real estate, profession... | \n", "sdsSWD_q2 | \n", "misCel_q1 | \n", "lddNDe_q1 | \n", "sdbPer_q3 | \n", "ssbCCo_q3 | \n", "misCel_q1 | \n", "sisBpM_q3 | \n", "ssbERI_q1 | \n", "O,P,Q. Public administration, education and he... | \n", "
45 | \n", "linWID_q2 | \n", "stbOri_q3 | \n", "ssbCCo_q2 | \n", "ssbElo_q1 | \n", "R, S, T, U. Other_q1 | \n", "population_q3 | \n", "lcnClo_q1 | \n", "linWID_q1 | \n", "stcOri_q2 | \n", "sdbPer_q2 | \n", "sscCCo_q3 | \n", "culture_nearest | \n", "nearest_water | \n", "lcdMes_q1 | \n", "sssLin_q3 | \n", "ldsAre_q3 | \n", "
46 | \n", "Code_18_211_q2 | \n", "ltcWRE_q2 | \n", "ldePer_q2 | \n", "lteOri_q3 | \n", "H, J. Transport and communication_q2 | \n", "ssbCCo_q3 | \n", "mdcAre_q3 | \n", "linP4W_q3 | \n", "sdbPer_q1 | \n", "ldsAre_q3 | \n", "population_q3 | \n", "sdsSWD_q2 | \n", "ldsCDL_q3 | \n", "night_lights_q2 | \n", "ldsAre_q1 | \n", "F. Construction_q3 | \n", "
47 | \n", "ltcWRE_q2 | \n", "sssLin_q2 | \n", "lcnClo_q3 | \n", "lcnClo_q2 | \n", "C. Manufacturing_q2 | \n", "ldsMSL_q1 | \n", "ldeAre_q3 | \n", "stcOri_q2 | \n", "lisCel_q1 | \n", "ssbElo_q1 | \n", "ssbSqu_q1 | \n", "F. Construction_q2 | \n", "stcOri_q2 | \n", "sicCAR_q2 | \n", "mtbNDi_q2 | \n", "ltbIBD_q3 | \n", "
48 | \n", "mtcWNe_q1 | \n", "ldeAre_q1 | \n", "sdcLAL_q3 | \n", "sisBpM_q1 | \n", "lcnClo_q1 | \n", "lddNDe_q1 | \n", "ldePer_q3 | \n", "sdsSWD_q1 | \n", "sisBpM_q2 | \n", "mtcWNe_q1 | \n", "sscERI_q3 | \n", "G, I. Distribution, hotels and restaurants_q1 | \n", "ssbCCD_q2 | \n", "O,P,Q. Public administration, education and he... | \n", "misCel_q2 | \n", "sicCAR_q3 | \n", "
49 | \n", "ldsAre_q3 | \n", "lcdMes_q1 | \n", "ltcWRE_q2 | \n", "stbSAl_q2 | \n", "ldeAre_q2 | \n", "misCel_q1 | \n", "ldeAre_q2 | \n", "linPDE_q3 | \n", "stbOri_q2 | \n", "sdsAre_q1 | \n", "lteOri_q3 | \n", "mtbAli_q2 | \n", "C. Manufacturing_q3 | \n", "population_q2 | \n", "sdbCoA_q3 | \n", "sdbCoA_q3 | \n", "