Refined data analysis for local marriages in Germany.

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.87 and b = 8.62 ⋅ 10−5. Notice that the parameter b accounts for the standard deviation associated to binned data. In general, the various regions seem to be homogeneously scattered around the best-fit curve. In the insets we show, as examples, the data pertaining to two selected regions, namely Bayern (upper inset) and Sachsen (lower inset). Notice that for both the insets the best-fit previously obtained for the whole data set (solid line) still provides a very good 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 for most regions the deviation with respect the unitary value is within the error bar. The largest deviation from the unitary value is for two city-states, namely Hamburg and Berlin. Given that all the regions (a few cases apart) are compatible with the overall best fit no further analysis on territorial homogeneity is performed. The colormap presented in panel c highlights the deviation of ρα with respect to the unitary value.