Estimation
Every slope in the companion comes from one equation. This page writes it down, defines each term, and shows how the regime comparison, the trade-elasticity curve, and the Castilian timing extend it.
The equation
- \(y_{j\ell d}\) is the diversity or the distinctiveness of the titles printed in city \(j\) in language \(\ell\) in period \(d\). The measures are defined in the measures explainer.
- \(\log \mathrm{MA}_{j\ell d}\) is the log of the city's market access in that language and period, defined in the market access explainer. The city's own population is excluded and the sum runs over the cities that printed in the language, or over the university cities for Latin.
- \(\alpha_{j\ell}\) is a city-and-language fixed effect. A city that printed in three languages contributes three series.
- \(\delta_d\) is a period fixed effect. Where a pool includes Latin it becomes \(\delta_{\ell d}\), one for each language and period, because Latin market access sums over a different set of cities from any vernacular's.
The slope \(\beta\) measures how a city's print in a language changed when its market access changed relative to that city's own average in that language and to the change common to all cities in the period.
The unit of observation
An observation is a city, a language, and a 50-year period, 1450–1499, 1500–1549, 1550–1599, or 1600–1649. Each measure is first computed for every city, language, and decade with at least five titles. The period value is the average of the decade values weighted by their number of titles.
A city-and-language series observed in only one period has no variation left after \(\alpha_{j\ell}\) is removed and does not enter the slope. The doctrine diversity sample keeps 1,209 of 1,578 observations, in 246 cities. Diversity over all titles keeps 1,458 of 1,876, in 291 cities.
Weights
Each observation is weighted by \(\sqrt{n}\), where \(n\) is its number of titles. A cell average built from \(n\) titles has sampling noise that scales as \(1/\sqrt{n}\), so the square root is the inverse-standard-deviation weight (Wooldridge 2010, ch. 19). Solon, Haider, and Wooldridge (2015) discuss weighting in panel regressions. The weighted six-language doctrine diversity slope is 0.250 (SE 0.050). Unweighted it is 0.219 (SE 0.066).
Standard errors
Standard errors are clustered by city, so the errors of one city are allowed to be correlated across its languages and periods. A second standard error, following Conley (1999), allows the errors of cities within 250 km of one another to be correlated. Figures 6, 9, and 12 of the companion draw intervals with the larger of the two. Figures 8 and 10 use the city-clustered standard error. For German doctrine diversity the slope is 0.746, with a city-clustered standard error of 0.194 and a Conley standard error of 0.255.
Magnitudes in standard deviations
The companion reports most slopes as the change in the outcome, in its standard deviations, when log market access rises by one within-city standard deviation.
Here \(\sigma_w(\log \mathrm{MA})\) is the standard deviation of log market access around each city-and-language mean, among series observed in more than one period, and \(\sigma(y)\) is the standard deviation of the outcome in the sample. For six-language doctrine diversity, \(0.250 \times 0.161 / 0.104 = 0.39\). For German, \(0.746 \times 0.140 / 0.081 = 1.28\). A within-city standard deviation of 0.16 log units is about the rise in market access a typical city saw over half a century.
Three extensions
The trade elasticity
Market access is rebuilt at each value of the trade elasticity \(\theta\) and the equation is estimated again. Each slope is scaled by the within-city standard deviation of log market access at that \(\theta\), so slopes at different elasticities are comparable. These are the curves in Figure 8.
Free entry and licensing
German and Italian are pooled as the free-entry regime and French and Spanish as the licensing regime. Each pool is estimated with the equation above. The difference between the two slopes comes from one regression on both pools,
where \(L_\ell\) is one for French and Spanish. The period effects \(\delta_{dL}\) are specific to each regime, so \(\gamma\) equals the licensing slope minus the free-entry slope. For doctrine diversity the two slopes are 0.683 and 0.022 and \(\gamma\) is −0.66 (SE 0.14). For doctrine distinctiveness \(\gamma\) is 1.24 (SE 0.48). The paper reports the two slopes. The test of their difference is computed for the companion and is not in the paper.
The Castilian timing
For Spanish the equation is estimated on decade cells, with city and decade fixed effects, separately before and from a cut year. The doctrine diversity slope is 1.375 (SE 0.872) before 1560 and −0.221 (SE 0.082) from 1560. These are the estimates in Figure 10.
Novelty and bridging
The variety results use the same equation on the same doctrine cells, with \(y_{j\ell d}\) the mean novelty or the mean bridging of the cell's doctrine titles. The novelty and bridging explainers define the two measures.
What the specification does not absorb
The variation in a city's market access comes from population change in other cities, which its printers did not control. The period effects absorb shocks common to all cities in a period. Shocks specific to a region and a period remain. A war, a confessional change, or a regional boom can change the populations around a city and what its presses print at the same time. The paper reads \(\beta\) as an association for that reason.
References
- Timothy G. Conley (1999). GMM Estimation with Cross Sectional Dependence. Journal of Econometrics 92(1), 1–45.
- Gary Solon, Steven J. Haider, and Jeffrey M. Wooldridge (2015). What Are We Weighting For? Journal of Human Resources 50(2), 301–316.
- Jeffrey M. Wooldridge (2010). Econometric Analysis of Cross Section and Panel Data. Cambridge, MA, MIT Press.