Refined data analysis for local marriages in Spain.

Panel a: log-log scale plot of LMαi,y versus Γαi,y, where different regions are denoted in different colours and symbols, as explained in the legend. These data are properly binned (green squares) and best-fitted (solid line) by , in agreement with the theoretical result (2). The best-fit coefficients are a = 1.852 and b = −5.16 ⋅ 10−4. Notice that the parameter b accounts for the standard deviation associated to binned data. It is clear, even by visual inspection, that a few regions slightly deviate from the best-fit curve. In the insets we show, as examples, the data pertaining to two selected regions, namely Castilla La Mancha (upper inset) and Galicia (lower inset). Notice that for both the best-fit previously obtained for the whole data set (solid line) still provides a proper fit. Panel b: For each region we calculate ρα, as defined in the text and deepened in the Theoretical Protocol Section. The horizontal line is drawn as a reference for the unitary value. Notice that the largest deviation from the unitary value is for Melila. From this plot we can distinguish regions displaying a relatively large number of marriages (i.e., ρα > 1) and regions displaying a relatively small number of marriages (i.e., ρα c. Interestingly, regions exhibiting analogous deviations share a certain degree of geographical proximity: regions with ρα > 1 (dark shading) correspond to the Southern part of Spain and to Basque region (with some adjacent regions), while regions with ρα d: the two clusters of regions highlighted are analyzed separately. For each we bin the related raw data and get a best fit, still according to the function , obtaining aup = 1.95 and bup = −1.79 ⋅ 10−4 (R2 = 0.97) for the set of regions with ρα > 1, and adown = 1.84 and bdown = −8.37 ⋅ 10−4 (R2 = 0.99) for the set of regions with ρα aup/adown ≈ 1.1. Binned data for the former set (triangles) and for the latter set (square) are shown in the main panel, together with the related best fits, in a log-log scale plot. These fits are not significantly better that the one obtained at the country level, suggesting that possible internal heterogeneities may be rather limited. In the insets we compare these best fits with raw data for two regions (Andalucia and Comunitat Valenciana) with ρα > 1 (upper inset) and two regions (Castilla y Leon and Cataluna) with ρα