World and National Wheat Market Integration in the 19th Century

Size: px
Start display at page:

Download "World and National Wheat Market Integration in the 19th Century"

Transcription

1 Chapter 5 World and National Wheat Market Integration in the 19th Century 5.1 Introduction Although market integration is one of the subjects in 19th century economic history that has drawn much attention, there is still room for improvement in terms of analytical tools. As more prices become available, demand for dynamic methods that can accommodate large cross-sections of price data increases. The origin of market integration local, national or international becomes a matter of interest as data increases in the cross-section. 1 This chapter draws from the literature on international business cycles and uses its tools for market integration research. Bayesian dynamic factor models are especially promising, since they measure comovement of many time series and go beyond bivariate comparisons. 2 Their complexity therefore grows only proportionally in the cross-section and along the time axis, but not exponentially as that of bivariate models including almost all cointegration frameworks. The 1 GAUSS code for the model used here is available from Chris Otrok s website at http: //people.virginia.edu/~cmo3h/research/wfac3b.prg. The code with the necessary modifications is available from the author upon request. 2 The dynamic factor is referred to as common component in this chapter. 99

2 5.2 Motivation dynamic setup allows for quantifying the share of price fluctuations due to world price movements, changes in national market conditions and local shocks. Applying comovement analysis to 19th century European and American wheat prices shows that the U.S. experience after 1870 was perhaps not particularly revolutionary to world wheat trade in contrast to what the established convergence literature à la O Rourke & Williamson (1999) suggests. It seems to be fair instead to speak of a major producer accessing the world s biggest market for wheat Western Europe, including the U.K. The results also call for reconsidering the relation of national and international market development as their respective timing turns out to be diverse. The largest boost toward global wheat market integration occurred before 1860, according to the results obtained here, before the railroad or the steamship could have had substantial effects. In the last quarter of the 19th century world wheat market integration accelerated further, but at a slower pace than before The next section motivates comovement analysis taking into account historical and technical considerations. Section 5.2 provides a discussion of the relevant methodological literature. An intuitive introduction into the method followed by a formal description can be found in Section 5.4. Section 5.5 describes the data, while Section 5.6 explains the findings. The last section concludes. The appendix contains details about the estimation and discusses some data issues. 5.2 Motivation The main focus of market integration research in the past three decades or so was on transatlantic market integration after 1870 (Harley, 1980; O Rourke, 1997). The main argument stated that markets on both sides of the Atlantic merged because of lower transport costs, which resulted in declining price differences. In the past few years scholars have begun to gather data from more markets spanning longer periods and applied improved econometric tools, mainly based on cointegration or price dispersion (Federico & Persson, 2006; Jacks, 2005; Persson, 1999; Sharp, 2008). The methodological arms race was aimed at accommodating the increasing amounts of data in a meaningful way and resolving the question 100

3 5.2 Motivation as to how market integration should actually be measured. 1 However, the curse of dimensionality has not yet been overcome, at least in the case when dynamic relations are to be accounted for. Most methods rely on bivariate comparisons cause the amount of parameters to increase exponentially in the number of markets. The complexity of the method proposed here grows linear in data size and allows for studying a large number of markets over a long period incorporating dynamic relationships between prices. Federico (2008) proposes to use the coefficient of variation as a measurement tool for market integration. While this method can handle many markets, it does not incorporate dynamic relationships. The obvious patterns of covariance in commodity prices are not exploited. Comovement analysis uses these patterns to estimate an unobserved common component that is used as a benchmark against which each single price is compared. This corresponds directly with Federico s (2008) critique of bivariate price comparisons. He calls for comparing prices against a hypothetical world price. Interestingly, this is exactly in the spirit of Veblen (1893) who discussed wheat prices after the American Civil War in the very first issue of the Journal of Political Economy. Thus, there appears to be a long tradition of analyzing local wheat prices on the basis of a latent world price. With this chapter I aim to abandon the purely technical need to use bivariate comparisons and contribute to reviving a more intuitive approach of analyzing market prices across time and space. As methods based on the coefficient of variation go to a certain extent in that direction, they cannot track each market s development over time, which is a natural feature of comovement analysis. Ejrnaes et al. (2007) use a multivariate error correction model that explicitly incorporates dynamic relationships. However, it is restricted in the cross-sectional dimension as its complexity grows exponentially with the the number of time series. Since it is based on cointegration, this and related classes of models rely on assumptions about the asymptotic properties of unit root tests. The model I propose uses Bayesian econometrics enabling my results to hold even if the sample is not representative for the whole population or the unit root properties of the data are not guaranteed (Uhlig, 1994). 1 Federico (2008) discusses the most popular methods so far. 101

4 5.2 Motivation There has been a recent revival of an older discussion about the reason for decreasing price gaps in the Atlantic trade between Knick Harley (1988) and Douglass North (1958, 1968): North repeatedly rejected the advantage of steam power and metal hulls in decreasing freight rates, and claimed that organizational improvements played a more important role in lowering transport costs and spurring change in international market integration in the first half the 1800s. According to recent studies important steps toward integration have indeed occurred before railroad, intercontinental steam ships and the telegraph could have had substantial effects. Wars and trade policy are instead suggested as important driving forces of market integration (Federico & Persson, 2006; Jacks, 2005; Persson, 1999). Brautaset & Grafe (2005) present another argument that is based on constant transport technology. They propose scale economies in market efficiency as an explanation for market integration. That is, costs per unit go down as volumes of trade go up, holding technology constant. Some of my results favor this explanation as will be shown in this chapter. However, this debate only centers on the supply side of the market for trade services. What about the demand side? In the absence of transport cost changes, trade may still increase and price gaps decrease if the relationship between supply and demand changes. Sharp (2008) claims that the main reason for declining price gaps between the U.K. and the U.S. was the increase of American wheat supply. This issue needs to be discussed at the national level as well. Kopsidis (1998) argues that industrialization creates urban demand for agricultural goods, and leads to regionally integrated markets holding agricultural and transport productivity constant. This theory may explain national differences in the timing of the relative development of national and international market integration. The conventional view using transport costs as the main argument is that national markets integrate first, since relatively short distances imply low transport costs, and then international trade links are created as technology reduces the cost of distance. However, this order could be reversed. In developing economies some cities may already be linked quite well to international wheat trade due to, for example, a strategic geographical location, while land-locked rural areas are separated from national and international wheat trade. As industrialization 102

5 5.3 Related Literature sets in, urban demand increases and nationwide specialization begins, fostering national markets on the basis of international integration. The feature distinguishing this study from others is that it truly differentiates between national and international integration. It can focus on single markets while keeping an eye on the aggregate development. Section 5.6 therefore delivers a number of new results about national integration in an international context. 5.3 Related Literature Various uses of the common component approach in market integration research can be found. Qin et al. (2006) used a dynamic factor in a vector error correction model (VEC) as an aggregate of all foreign prices and compare it to an observed home price. By doing so, they augmented the VEC model to the multivariate case. My model is simpler. It conceptualizes the common component or factor as a manifestation of the comovement of prices in different markets and therefore as a manifestation of the law of one price. Common factors can also be used in panel cointegration frameworks to increase the power of multivariate unit root tests, which goes back to Bai & Ng (2004) and Pesaran (2007). Applying this method in a recent study on German regional prices levels, Dreger & Kosfeld (2007) found a persistent lack of price convergence among German regions in the period from 1995 to Principal component analysis has been used by Sánchez-Albornoz (1974) in a truly pioneering study. He analyzed annual Spanish wheat and barley prices between 1856 and 1889 and focused on the causal effects of wheat trade between regions. He identified regions along geographical and agricultural borders, but could not model the dynamic relationships contained in the data. Thus, in later studies he began to use univariate time series analysis and abandoned the common component approach. Moreover, there is a strand in the international finance literature similar to what I am pursuing. Bekaert & Harvey (1995) proposed a time-varying measure for integration of national markets into the world market for capital assets. The set of country price returns in their model is explained jointly by a world benchmark portfolio and idiosyncratic country risk. The varying degree to which each 103

6 5.3 Related Literature of the two factors explained the returns was interpreted as a measure of world capital market integration, which is comparable to the latent factor of all prices in my chapter. Technically, this study is closest to Kose et al. (2003a), although their research interest is not directly related to mine. They estimated the common component of output, consumption and investment between 1960 and 2001 for the G7 countries, and identified a world component and national components. I work with practically the same model but my focus is on price data and market integration. The relation between international business cycle transmission and world market integration is obvious. Both describe different strata of globalization. While I measure market integration as it is manifested in the price comovement of a traded and crucial commodity, business cycles represent integrated markets subject to common output variations. 1 Thus, it is appropriate to cite a recent paper that applies dynamic factor models to identify long run business cycle comovement between four Latin American economies starting in the last quarter of the 19th century, i.e., Argentina, Brazil, Chile and Mexico (Aiolfi et al., 2005). It finds strong exogenous shocks on the cyclical activities of Argentina, Chile and Mexico. Brazil was obviously better insulated from external influences. It remains a field of future research if this result is reflected in the behavior of traded goods in these countries. A large body of literature utilises mixed sets of time series, i.e., prices and volumes, where either the common component is interpreted as the business cycle, or nominal and real common components are identified (Stock & Watson, 1998). An application of a single component model to historical time series is described in Sarferaz & Uebele (2007). 2 One further application to price series is Reis & Watson (2006). They interpret the common component of a consumer basket of prices as the part of price fluctuations whereas relative prices stay the same (Bryan & Cecchetti, 1993). They use it to analyze the degree of money neutrality and find that prices were not neutral in the U.S. after The theoretical literature is not entirely conclusive about the correlation of trade flows and output variations, and I do not attempt to contribute to this question. 2 See Chapter 3 of this thesis. 104

7 5.4 The Model 5.4 The Model Intuition Comovement measures synchronous price movements in large cross-sections. It has a certain similarity to correlation, the main difference being that correlation is defined over pairs. Another important difference is that comovement measures linear dependence not only in a given period but across time. It represents the whole spectral matrix of leading and lagging correlations (Kose et al., 2003a, p. 1218). While correlation can be understood as a simple bivariate counterpart of comovement, convergence captures a different aspect of relative prices. Take for example Harley s (1980) classic paper about convergence among U.K. and U.S. wheat prices in the second half of the 19th century, which begins with a graph showing a shrinking price gap between Chicago wheat and the British Gazette price (Figure 5.1). Harley, as well as other scholars succeeding him, refers to this closing gap when defining market integration (O Rourke & Williamson, 1999). However, Figure 5.1 s second striking element which Harley does not discuss is the degree of correlation between the two prices. This element is another important feature and thus used as the main argument in this chapter. This bivariate comparison can be carried out for more than two markets by comparing all single series with a benchmark series that represents maximum comovement of all observed series. Consider a simple steady state example. Imagine splitting up each single price p i into a part c that is the same in all markets and a part u i that is the deviation from c: p i = c + u i (5.1) The common part c could be any rational number, but optimally it should be the one that minimizes the sum of the deviations u i. For the sake of the argument, the mean over all prices p i could be chosen, and the absolute deviations u i could be expressed in percentage of p i to show how much the common price element explains in each market i. An example is given in Table

8 5.4 The Model Figure 5.1: Price convergence between Chicago and Britain (Harley, 1980). In the left Column of Table 5.1 the observed prices are given. An arbitrarily chosen common component is given in Column 2 and the respective deviations in Column 3. Column 4 contains a normalized measure of how well the common component explains single prices in Column 1. The lower the percentage number, the more the common component explains single prices. This is an example in the steady state, but what if we include dynamics? How would the common component change if the prices in Column 1 changed? These questions are answered in the following formal discussion. It starts with a dynamic extension of Equation Single Common Component Given the relation between N prices and their common component in time t, consider what happens if all prices changed in the same direction and to the same degree. In this case their variation should be only due to the common part c but not the remaining part u i, i = 1,... N. If some prices changed to a different degree or in different direction from the others, the common part will explain only 106

9 5.4 The Model Table 5.1: Numerical example of common component. Price Common Deviation Absolute Component Percentage Deviation p i c u i u i /p i % % % % % % % % a fraction of the price variation and the rest will be due to the specific component u i,t. Thus the dynamic formulation of Equation 5.1 for each p i,t is p i,t = a i + λ i c t + u i,t (5.2) Here c t represents the common component, which is the same for all markets and therefore not indexed by i. There is a constant a i and a weight λ i that links the common price component to the i-th variable. u i,t, the idiosyncratic or specific component, accounts for local, market specific influences, e.g. local crop failures or temporary demand fluctuations. However, the idiosyncratic parts may experience their individual dynamic processes, i.e., they may be serially correlated, which is expressed as an AR(p)- process: u i,t = θ i,1 u i,t θ i,p u i,t p + χ i,t (5.3) Equation 5.2 resembles a linear regression, only that we do not observe the regressor c t. We can instead describe c t s dynamics by an AR(q)-process and 107

10 5.4 The Model treat it (together with Equation 5.3) as the transition equation in a state space model: c t = ϕ 1 c t ϕ q c t q + ν t (5.4) These three equations describe the basic setup. However, since many, not only one, common components are strived for, Section extends the model to the case with K common components. In Section I will show how the model parameters and the common component can be estimated. The error term assumptions will be discussed next Error Term Assumptions The local market shocks u i,t are assumed to be normal and uncorrelated in the cross-section: E [u i,t u j,t s ] = σ 2 u i i = j s = 0, 0 otherwise. The error term χ i,t in the local market shock s process is likewise normal, and serially and cross-sectionally uncorrelated: E [χ i,t χ j,t s ] = σ 2 χ i i = j s = 0, 0 otherwise. The common component s error term ν t is normal with E [ν t ν t s ] = σ 2 ν t for s = 0, 0 otherwise. The error of the common component ν t is uncorrelated with the error of the local component χ i,t : E [χ i,t ν t s ] = 0 i, s. 108

11 5.4 The Model Multiple Common Components So far I have explained the estimation of only one common price component. If this represents comovement of all prices in the sample this can be referred to as the world component. Each local price series is thus explained by its comovement with the world price and local shocks. However, additional shocks may arise from the national level. Those shocks may be different from global shocks, for example if there are strong border effects that do not allow the transmission of national demand shocks to other countries. In the framework proposed here it is possible to estimate both global and national common components in one model and assess their relative explanatory power. Essentially, the world component explains the variance in all price series and therefore its corresponding weights are all nonzero, while the national components explain only the variance of some price series identified by nationality. For example, the national component of Spain is identified by setting all weights that belong to cities outside Spain to zero. Identifying national components ex ante is opposed to obtaining multiple orthogonal common components endogenously and identifying them ex post. It implies that the national common components do not need to be orthogonal to each other. What is orthogonal, however, is each national component relative to the world component. In this setup, Equation 5.2 can be formulated as: P t = ΛC t + U t, (5.5) where P t is an N 1 vector of N price series, C t is a (K + 1) 1 vector of common components (the number of a column of ones plus common components), Λ is a N (K + 1) matrix of weights and U t is a N 1 vector of idiosyncratic components. In the case of an international and several national common components, there is one international common component, and R < N national components, with R being the number of countries in the model. For example, in the case of N = 10 markets and R = 2 countries, each being represented by half of the sample, K = 1 + R = 3. The 10 (K + 1) matrix Λ t contains a column of constants, one column of N elements for the world component and then R 109

12 5.4 The Model columns of weights for each country, which are only nonzero for the observations for the respective country. Equation 5.5 for this case could then be formulated as: P t = A + λ w c w t + λ 1 c 1 t + λ 2 c 2 t + U t, (5.6) where P t consists of 10 price observations in period t, A = [a 1, a 2,..., a 10 ] is a vector of constants, λ w = [λ w 1, λ w 2,..., λ w 10] is a vector of weights that are nonzero for all i, c w t is the value of the world price component in t, λ 1 = [λ 1 1, λ 1 2,..., λ 1 10] is the national component for country 1, where only those λ 1 i are nonzero that correspond to cities of country 1. All other λ 1 i are set to zero. Accordingly, the elements contained in λ 2 are only nonzero if corresponding to cities in country 2. U t is a 10 1 vector and contains price elements not explained by either the world component or the respective national component. The transposed matrix of weights Λ then looks like: world component national component 1 national component 2 a 1 a 2 a 3 a 4 a 5 a 6 a 7 a 8 a 9 a 10 λ w 1 λ w 2 λ w 3 λ w 4 λ w 5 λ w 6 λ w 7 λ w 8 λ w 9 λ w 10 λ 1 1 λ 1 2 λ 1 3 λ 1 4 λ λ 2 6 λ 2 7 λ 2 8 λ 2 9 λ 2 10 = Λ Accordingly, a single price observation p i,t is composed of the following elements: p i,t = a i + λ w i c w t + λ 1 i c 1 t + λ 2 i c 2 t + u i,t, (5.7) where one of the national weights is zero. I follow Kose et al. (2003a) in estimating the multiple common components. They apply a sequence of single common component models. c 1 t and c 2 t are estimated for the variance unexplained by c w t. The rest of the model is the same as 110

13 5.4 The Model in the single component model above, only that there are K = 1 + R common component AR-processes now: where k {w, 1,..., R} and c k t = ϕ k 1c k t ϕ k qc k t q + ν k t, (5.8) E [ ] νt k ν j t s = σ 2 ν k = j s = 0; 0 k, s. k Estimation The classical estimation of a state space system is standard in multivariate time series econometrics (Geweke, 1977; Hamilton, 1994; Stock & Watson, 1989). However, I follow the Bayesian way of estimation, in part because the dimensionality is much less of a problem as I can use Gibbs sampling that reduces the curse of dimensionality. Another reason is that it is a convenient way to deal with uncertainty about the unit root properties of the variables (Uhlig, 1994). The single component model is described in detail in Otrok & Whiteman (1998). In the next subsection I relate mostly to Kose et al. (2003a), which features the estimation of K common orthogonal components. Kose et al. (2003a, p. 1220f) explain why Bayesian estimation is convenient for this class of models. The challenging aspect of the estimation is that both a linear regression with serially correlated errors and the AR-coefficients of c t have to be determined simultaneously. This is done in classical statistics by utilizing the linearity of the model in the observables. A Kalman filter-smoother procedure leads to the unobserved parameters likelihood function. It is Gaussian and can be estimated by maximum likelihood (Stock & Watson, 1998). However, if the the model becomes large in the cross-section it is difficult to estimate, since dimensionality increases exponentially (Kose et al., 2003a, p. 1220f). Bayesian methods allow for estimating the common component c t and the other parameters (φ, θ i, σ u, etc.) of the model separately. In Bayesian statistics the unknowns are treated as random variables, as opposed to in classical statistics where they are treated as constants. Treating the model parameters and the 111

14 5.4 The Model common component as random implies determining their probability distribution. Unfortunately, for the model above, the joint distribution of the parameters and the common price component is nonstandard. This problem can be solved by decomposing the joint distribution of the parameters and the common component into conditional distributions. One is the distribution of the parameters conditional on c t, and the other is the distribution of c t conditional on the model parameters. These conditional distributions have standard forms and are therefore computable. Moreover, the optimization problem for the variable specific parameters is done separately for each observable p i,t and does not increase exponentially with the number of variables N, since the covariance matrix of the u i,t is diagonal, i.e., all cross-sectional correlation is contained in the common components (Kose et al., 2003a, p. 1220) and not in u i,t. The sampling, i.e., making random draws from posterior distributions derived from the model, is done by a Markov chain Monte Carlo (MCMC) procedure. 1 Upon iterating on sampling, the Markov property of the marginal distributions, which is to converge to an asymptotic distribution, is utilized. To begin, a vector of arbitrary starting values is chosen for the common component. The distribution of the parameters conditional on that value is then determined and a vector of values for the parameters is sampled, which finishes the first iteration. In the second iteration, a new value for the common component is drawn conditional on the draw for the parameters from the previous iteration. Then, new values for the parameters are sampled conditional on the new common component draw. The procedure is repeated until convergence is achieved. It can be shown that the conditional posterior distributions converge to the true desired marginal posterior distributions as the number of iteration steps goes to infinity (Geman & Geman, 1984) Here the number of draws is 24,000 of which I use 20,000 for inference, and discard the first The latter is done in case the starting value was chosen inferiorly. As a convergence check, I repeat the procedure several times with different starting values and compare the respective results. The AR-order for the world common components is chosen as q = 8, which reflects business cycle frequency with annual data. For the variable 1 Monte Carlo, because artificial data are generated, and Markov chain, because the distribution conditions only on the last iteration and not on the whole history. 112

15 5.4 The Model specific processes, p = 3 is chosen following Kose et al. (2003a). I have estimated several variations of this setup and found that the results are robust to the choice of the AR-orders. Two approaches are proposed in the literature for sampling the common component conditional on the parameter draw. One is Otrok & Whiteman (1998) and the other Kim & Nelson (1999b). I have followed the former, having compared the outcomes of both methods. The results are basically the same. In the appendix, Section C.1, I formally describe the sequence of draws Identification The identification issues here have been discussed extensively in Otrok & Whiteman (1998). The principal problem is that the weights λ i and the latent variable c t are determined jointly in Equation 5.2 and the following two cases are observationally equivalent: λ i c t and ( λ i )( c t ). This problem can be solved by pinning down one (and only one, since this in turn pins down c t ) λ i to be positive. In the example chosen above identification is achieved through setting λ w 1 greater than zero, as well as λ 1 1 and λ 2 6, the first national factor loadings in rows 3 and 4, respectively (Kose et al., 2003a, p. 1219). Here I choose the weight corresponding to the price of wheat in London to be positive; i.e., to be positively correlated with the world price, which does not seem to be a very strong restriction. The cities whose prices are assumed to be positively correlated with their respective national common component are Paris for France, Berlin for Germany, Stockholm for Sweden, Vienna for Austria-Hungary, Brussels for Belgium, New York for the U.S., Oslo for Norway, and Santander for Spain. A similar problem arises if for example c t is measured in centimeters and λ i in inches or vice versa. The scale of the common component is undetermined, which is due to the fact that the variance of the common components error term ν t is not identified. Following, among others, Sargent & Sims (1977) it is set to one, but it could be set to any other constant likewise. 113

16 5.4 The Model Priors The priors I use are the same as those in Kose et al. (2003a, p. 1221). Five prior distributions must be chosen. The first two are the distributions of the AR-parameters for the common component (Equation 5.8) and the local shocks (Equation 5.3). Next is the prior distribution of the factor loadings λ followed by the prior distributions of the variances of the local shocks and the common components error terms, σ 2. The latter is the easiest, since for identification ν k purposes explained above it must be set to a constant and thus has no distribution (Section 5.4.6). The variance of the local shock s error term has an inverted gamma prior distribution: σ prior χ IG (6, 0.001) which implies a fairly loose prior. Thus I do not claim to have important prior knowledge about the idiosyncratic error variance and leave the setting of its value mostly to the data. The AR-parameters of both the common component and the local shocks have normally distributed prior distributions with zero mean, implying the assumption that they are not serially correlated. The more distant the lag is, the more certain this assumption becomes, and thus the variance around zero decreases exponentially: ϕ prior N(0 q 1, Σ) where with x = p or q. Σ = x, 114

17 5.4 The Model Presentation of Results This section explains how the results from the estimation are presented. It especially describes how the variance decomposition is carried out. In order to capture changes of market integration I choose subperiods. It means that the model is estimated separately for subsequent time periods. In each subperiod, a new world price component and new national price components are derived from the data. Section 5.6 starts with showing the importance of working in subperiods. I have experimented with different subperiods, and finally decided to choose periods of about 25 years starting from 1806 in order to divide 19th century into even quarters: This choice is supported by major historical events over these 102 years. The first quarter captures a period of disintegrated world markets, as a result of the Napoleonic Wars and British import tariffs of those years. The next period up to the mid-1850s possibly exhibits increasing market integration as steam ship technology, the railroad and organizational improvements started to proliferate, fewer wars occurred on the European continent, and liberal trade politics became more widespread. The following quarter should continue that development although it includes the American Civil War, which is likely to have had a severe impact on world wheat trade. In the same subperiod, tariffs were reduced due to the treaties induced by Cobden-Chevalier, which however seemed to have little effect on wheat trade (Lampe, 2007). The last subperiod starting in 1880 is likely to exhibit a strong drive toward Atlantic market integration according to O Rourke & Williamson (1999). On the other hand, some countries increased tariffs that had been lowered or abolished earlier in the century. 115

18 5.5 Data In order to assess the relative explanatory power of the components for each price series I follow Kose et al. (2003a) who employ variance decomposition for each price series i of the form var(p i ) = (λ w i ) 2 var(c w ) + (λ 1 i ) 2 var(c 1 ) + (λ 2 i ) 2 var(c 2 ) + var(u i ) (5.9) resulting in the fraction of volatility explained by the world component: (λ w i ) 2 var(c w ). (5.10) var(p i ) Since sampling from conditional distributions yields sampling error, the orthogonality of the common components is not automatically given, although they are uncorrelated. Thus, at each step of the Markov chain the national components are orthogonalized relative to the world component. Numerically, this does not change the results in any relevant way, but ensures that the volatility shares add up to 1 (Kose et al., 2003a, p. 1226). In order to give valuable insights into the relative explanatory power of each component, I present arithmetic averages of the volatility shares Data The data set is taken from Jacks (2005), Jörberg (1972, Sweden), and Jacobs & Richter (1935, Germany). 2 I do not use all series, because some start too late or end too early. For several reasons, I work with annual data here: first, it increases data coverage, second, the problem of seasonality does not arise, and third, it is interesting for economic historians if the proposed method is applicable even to low frequency data. The data set with which I work contains 60 annual wheat price series ranging from 1806 to The Norwegian and the Spanish series as 1 I carried out the same for the standard errors of the decomposed variances, which can be found in the appendix. This was not done at every step of the Markov chain, but still provides a good view of the average accuracy of results of the variance decomposition. 2 All series are annual or seasonal variations have been controlled for in another appropriate way (Jörberg, 1972). 116

19 5.5 Data well as the observations for Cincinnati only start in In the results section, I discuss how I included them in the estimation. Wheat prices are observed in the following markets: Austria-Hungary (4): Vienna, Lwow, Krakow, Ljubljana (Krakow did not belong to the Hapsburg monarchy for the whole period) Germany (4): Königsberg/Kaliningrad, Hamburg, Berlin, Munich Belgium (3): Ghent, Bruges, Brussels U.K. (12): London, Manchester, Liverpool, Exeter, Carmarthen, Dover, Gloucester, Worcester, Cambridge, Norwich, Leeds, Newcastle France (11): Bayeux, Saint-Brieuc, Toulouse, Bordeaux, Chateauroux, Mende, Barleduc, Arras, Pau, Lyon, Paris U.S. (4): New York, Alexandria, Philadelphia, Cincinnati Sweden (11): Stockholm, Uppsala, Södermanland, Östergötland, Kalmar, Halland, Skaraborg, Örebro, Västmanland, Gästrikland, Hälsingland Spain (9): Cordoba, Gerona, Granada, Lerida, Oviedo, Segovia, Zaragoza, Santander, Burgos Norway (2): Bergen, Christiania/Oslo Figure 5.2 shows the geographical extent of the European markets in the sample. Figure 5.3 presents the same for the U.S. markets. Although scattered data exists for Italy, Odessa and more German cities, the coverage is not sufficient. There are 19th century wheat prices from many German cities, but only up to the 1860s (Fremdling & Hohorst, 1979). Italian data is plentiful, too, but it either starts in the 1860s or ends in The reason why there is sometimes better data coverage for the first half of the 1900 may 1 Earlier prices for Cincinnati exist, and have been used to form wholesale price indexes for the Ohio delta (Berry, 1935). Future research will complete the data set using the sources given there. 117

20 5.6 Results be that administrations tried to control prices to preserve domestic peace, but progressing political and economic liberalization led the states to abandon those attempts after the middle of the 19th century. 1 For the empirical model employed here it is not necessary to convert the coins and weights and volumes as long as they remain constant over time. Unit differences only represent permanently different means that do not affect comovement. The means of all data series are therefore normalized to 0. Similarly, the variance of each series is normalized to 1. The price data provided by Jacks (2005) is converted to American dollars per 100kg. I ensure that all prices are expressed in gold dollars but not in greenbacks. There are large relative price variations during the 1860s, which make such an exercise advisable. I used gold denominated benchmark price series from independent sources to compare them with the original prices. I also directly look at the relevant exchange rates found in the Global Financial Database. The appendix documents this in detail. I find that Spanish, Austria-Hungarian, and English prices are converted to gold dollars, while the others are provided in greenbacks. As a consequence, instead of recalculating all currency conversions, I continue to work with Jacks data, deflating the series given in greenbacks such that all series are denominated in gold dollars per 100kg. Although overall inflation in the 19th century was small, I take out long run trends from the data by applying a Hodrick-Prescott filter with the Ravn-Uhlig lambda of Results Evidence for Structural Change In order to give an impression of the changing comovement in the 19th century Figure 5.4 plots a world price component spanning 102 year from 1806 to I thank Michael Kopsidis for this remark. 2 Alternative filters like Baxter-King and Christiano-Fitzgerald yield very similar results (Baxter & King, 1999; Christiano & Fitzgerald, 2003). 118

21 5.6 Results Figure 5.2: European markets. 119

22 5.6 Results Figure 5.3: American markets. against subperiod world price components, each spanning about 25 years. 1 The standard deviation of all series price components is normalized to 0.1, which is the average standard deviation of English wheat prices in the 19th century. The fluctuations are percentage deviations around a smooth trend as discussed in the data section. The vertical lines divide subperiods. We can clearly observe the price peak in 1847, which was induced by bad harvests throughout Europe in 1846 (Berger & Spoerer, 2001). The common component confirms also the notion that wheat did not strongly take part in a worldwide speculative inflation of Veblen (1893, p. 20). If structural change did not matter, estimating the model for the whole period or breaking it up would make no difference to the common price component. What Figure 5.4 demonstrates instead is that, in three out of four periods, the subsample world price (solid line) is different from the world price estimated from 1 Kose et al. (2003b) estimate a similar model for subperiods in order to capture structural breaks in world business cycle comovement. 120

23 5.6 Results Figure 5.4: World price component during 19th century. Whole sample estimate vs subsample estimates World Price Component 48 Markets Whole Sample Subsamples the full sample (dotted line). Before 1830 the subperiod world component does not exhibit the strong peak around 1812 as the full sample component does. Conversely, the M-shaped price deviations centered on 1895 are much stronger when taking only information after 1880 into account. In the 1850s and 1860s, the restricted sample world price fluctuates less than its full sample counterpart. Apart from misspecification of the common component, taking the period as a whole prohibits the relative weights λ i to change over time. They represent each market s sensitivity to world price fluctuations. Without breaking the sample up, it would be impossible to capture changes in the degree of how much single markets take part in world wheat trade. 121

24 5.6 Results Antebellum International Market Integration The trend obtained for worldwide comovement is shown in Figure 5.5 as a cross country average. It reveals the integrating forces of world trade in the 19th century. Globalization according to Figure 5.5 manifested itself in a strong increase between the first and the second quarter of the 19th century: the world component explains on average only 34% price variance between 1806 and 1829, but 58% between 1830 and This first shift to integrated markets is followed by a second, albeit smaller one, as the flatter shape of the upper curve in Figure 5.5 shows. At the same time national market integration, representing solely national comovement, steadily declined in the 1800s, first quickly, then at a slower pace. These two developments add up to decreasing market separation; i.e., on average 19% of price variance was subject to local shocks before 1830, compared to only 8% after 1880 (not shown). These results emphasize the role of market integration in the first half of the 19th century, as opposed to the First Wave of Globalization story famously put forward by O Rourke (1997), among others. Recently there has been greater recognition of the former point, the post-napoleonic improvement of world commodity market integration. Jacks (2005) finds decreasing transport costs for wheat before 1860 and Federico (2008) reports declining price dispersion in European wheat prices from the early 1800s on. Kaukiainen (2001) proposes decreasing information costs on why markets could come to function so much better without the widespread use of steam technology and the telegraph. He finds that business letters to and from London traveled on average twice as fast in the middle of the 19th century than at the beginning. Factual transport cost in sail shipping decreased impressively, as Brautaset & Grafe (2005) find. They argue that economies of scale can explain these cost reductions in shipping Broad Trends The country averages are presented in Figures Figure 5.6 contains averages for the major countries such as the U.K., the U.S., France, Austria-Hungary and Germany, while Figure 5.7 features Norway, Sweden, Belgium and Spain. 1 1 The full set of results can be found in the appendix. 122

25 5.6 Results The upper panel shows world market integration, the middle national market integration and the lower one market separation in the sense that it shows the share of price variance, which is explained by neither national nor international comovement. The results are presented in such a fashion that it is possible to get a complete picture of a country s market integration process throughout the 19th century. In the U.K. for example, the upper part of Figure 5.6 shows that it continuously integrated into the world market. At the same time the middle panel shows that conditional on world market integration the importance of national specific shocks decreased continuously. The lower panel shows a very low value of separation. Thus already in the beginning of the 19th century, U.K. markets were either nationally or internationally integrated with local shocks having almost no influence on local prices at all. Three new perspectives on 19th century wheat market integration can be deduced from these results. First, European market integration before 1860 has been neglected in the literature relative to transatlantic market integration. The European market is not the consequence of a return to protectionism under Bismarck & Co in the 1880s but a result of weakening market forces after the end of the Napoleonic Wars. Second, the American Civil War represents a phase of transatlantic market disintegration; therefore postbellum market integration is a candidate for a postwar reconstruction phenomenon. This puts the Grain Invasion into perspective. The third result relates to the relative development of national and international integration, which was different across countries. U.S. Accession to European Integration Levels These results imply that the story to be highlighted in 19th century international wheat market integration is maybe not so much one of America and Europe merging into a common transatlantic market due to reduced transport costs, but rather that of an already integrated European market that opened up to another major supplier. One could also argue that it is likely that the strong increase in comovement between American markets with the world price could have occurred one or two decades earlier but was impeded by the Civil War (see Sharp (2008) 123

26 5.6 Results for a similar argument). The level of integration of the U.S. markets at the end of the century, however, is comparable to that of France s as the top panel of Figure 5.6 shows. While the U.K. s markets are integrated best in the last quarter of the 19th century, major German cities come second, well above the U.S. markets. If the U.S. markets are placed between France and Germany, two protectionist countries in terms of market integration, then a changed perspective relative to the traditional story of protectionist and free market nations after 1880 seems to be justifiable (see O Rourke & Williamson (1999) for the established view). Late National Market Integration and Urban Demand for Wheat The third set of results stems from the ability of the model to truly differentiate between national and international integration. They evolve differently across nations in the 19th century as the middle panel in Figure 5.6 shows. While toward the end of the 19th century we observe decreasing sensitivity to national shocks in the U.K., the U.S. and Germany, in France and especially in Austria- Hungary national shocks explain more than they do in the preceding subperiod (Figure 5.7). A related pattern can be observed in Spain, where national shocks explanatory power decreases but remains on a very high level. On average wheat prices in Spain are explained 50% by national shocks, and about 35% by local shocks (Figure 5.14). While protectionism seems to be the first possible explanation for increasing shares of the national component after 1880, another important factor may have been differences in the timing of industrialization across the nations discussed here. Kopsidis (1998, 2002) argues that integration of agricultural markets may have been spurred by industrialization. The creation of urban demand for agricultural goods led to regional specialization, which made it profitable for farmers to produce for the market and ship their goods to the cities. Even in the absence of technological progress unit transport costs may fall during this process due to economies of scale in market efficiency. 124

27 5.6 Results Country Results The U.K. is certainly the country of which the highest world market integration and lowest separation of markets can be expected. Easy access to most U.K. cities by waterways and early commercialization are some of the often cited reasons for why industrialization first developed in England. Shiue & Keller (2007) for example find that at the beginning of the 1800s the U.K. was much better integrated than Western Europe on average. One still has to bear in mind that the first few decades witnessed the Napoleonic Wars and the Corn Laws that had an insulating effect on the U.K. wheat markets (Sharp, 2006). Thus, one may find some limited world market integration and very low local market separation before Figure 5.8 shows exactly this. Note that there were some differences in the degree to which markets were globally integrated. 76% of wheat prices in Exeter were already determined by the world market, while prices in all other cities stayed below 60%. 1 In the aftermath, all cities converged to the same path of world market integration with shrinking nationally specific price components. Together with Belgium the U.K. emerged clearly as the most integrated wheat market during the second half of the 19th century. In comparison, France s markets experienced a rather marginal development with considerable local heterogeneity. While early in the 1800s cities such as Paris, Lyon and Arras were better integrated into the world market than most U.K. cities, after the 1850s there was rather an anti-global tendency in those cities (Figure 5.9, upper panel). The French average for world integration did not decrease, however, because Toulouse, Mende and Pau followed a path of increasing integration. France ended the early 1900s as an intermediately globally integrated country with the largest national specific component in the sample (Figure 5.6). The figures showing market integration in Belgium speak very much for themselves (Figure 5.10). All three markets behave almost identically. The largest part of price variation is subject to worldwide shocks (upper panel). The only peculiarity is the swing with a short setback to national specific price variation during the third quarter of the 19th century (upper and middle panel). It is 1 In the case of the U.K., France, Spain and Sweden, I did not plot all results to keep the graphs straightforward. All results can be found in the appendix. 125

28 5.6 Results striking that the Belgian markets in the sample were always fully insured against locally specific price shocks, reflected in the extremely low level of separation of all three markets over all periods (lower panel). The disaggregated view allows for detailed insights into the development of German market integration. The upper panel of Figure 5.11 describes the astonishingly strong international integration of the German wheat markets. A global focus arose in the second quarter of the 19th century, too early to be fully explained by railroad or telegraph. It seems plausible to cite the German tariff union as the main driver of market integration. But this would imply strong explanatory power of the national component after 1830, and an increase relative to the first quarter (middle panel). Unfortunately, the small set of markets presented here is not suitable to confirm or refute this hypothesis, because Hamburg never became part of the tariff union, while Königsberg and Berlin have never been divided by tariff boarders. However, it is still interesting to find an already well integrated national market before the tariff union ever existed. Similar results were also found by Fremdling & Hohorst (1979) and Shiue (2005). 1 The evidence rather fits into a picture of increasing Baltic and North Sea trade that included the major Prussian and former Hanseatic cities. The introduction of the railroad was not without effect, though. It became important where waterways were missing: Munich, not connected to the Rhine-Weser-Elbe water network, has not been well integrated into the North-German market, as can be seen in the lower panel. Its separation from the North-German market was attenuated probably with the help of the railroad in last quarter of the 19th century. The United States is the only major nation that developed no strong connections with the world market before Although the U.S. already supplied about half of the U.K. s wheat imports in the 1850s (Sharp, 2006), they were apparently not well integrated in the world market. What matured earlier, however, was a national market. Between the first and the third quarter of the 19th century comovement between the three East Coast cities New York, Philadelphia and Alexandria increased to such an extent that 60-70% of the variance of their prices was explained by the national component (middle panel in Figure 5.12). 1 See Dumke (1991) for details on the German tariff union. 126

Nowcasting Norwegian GDP

Nowcasting Norwegian GDP Nowcasting Norwegian GDP Knut Are Aastveit and Tørres Trovik May 13, 2007 Introduction Motivation The last decades of advances in information technology has made it possible to access a huge amount of

More information

Labor-Supply Shifts and Economic Fluctuations. Technical Appendix

Labor-Supply Shifts and Economic Fluctuations. Technical Appendix Labor-Supply Shifts and Economic Fluctuations Technical Appendix Yongsung Chang Department of Economics University of Pennsylvania Frank Schorfheide Department of Economics University of Pennsylvania January

More information

Identifying the Monetary Policy Shock Christiano et al. (1999)

Identifying the Monetary Policy Shock Christiano et al. (1999) Identifying the Monetary Policy Shock Christiano et al. (1999) The question we are asking is: What are the consequences of a monetary policy shock a shock which is purely related to monetary conditions

More information

Dynamic Factor Models and Factor Augmented Vector Autoregressions. Lawrence J. Christiano

Dynamic Factor Models and Factor Augmented Vector Autoregressions. Lawrence J. Christiano Dynamic Factor Models and Factor Augmented Vector Autoregressions Lawrence J Christiano Dynamic Factor Models and Factor Augmented Vector Autoregressions Problem: the time series dimension of data is relatively

More information

Oil price and macroeconomy in Russia. Abstract

Oil price and macroeconomy in Russia. Abstract Oil price and macroeconomy in Russia Katsuya Ito Fukuoka University Abstract In this note, using the VEC model we attempt to empirically investigate the effects of oil price and monetary shocks on the

More information

1 Bewley Economies with Aggregate Uncertainty

1 Bewley Economies with Aggregate Uncertainty 1 Bewley Economies with Aggregate Uncertainty Sofarwehaveassumedawayaggregatefluctuations (i.e., business cycles) in our description of the incomplete-markets economies with uninsurable idiosyncratic risk

More information

Purchasing power parity over two centuries: strengthening the case for real exchange rate stability A reply to Cuddington and Liang

Purchasing power parity over two centuries: strengthening the case for real exchange rate stability A reply to Cuddington and Liang Journal of International Money and Finance 19 (2000) 759 764 www.elsevier.nl/locate/econbase Purchasing power parity over two centuries: strengthening the case for real exchange rate stability A reply

More information

Endogenous Information Choice

Endogenous Information Choice Endogenous Information Choice Lecture 7 February 11, 2015 An optimizing trader will process those prices of most importance to his decision problem most frequently and carefully, those of less importance

More information

Monetary and Exchange Rate Policy Under Remittance Fluctuations. Technical Appendix and Additional Results

Monetary and Exchange Rate Policy Under Remittance Fluctuations. Technical Appendix and Additional Results Monetary and Exchange Rate Policy Under Remittance Fluctuations Technical Appendix and Additional Results Federico Mandelman February In this appendix, I provide technical details on the Bayesian estimation.

More information

Searching for the Output Gap: Economic Variable or Statistical Illusion? Mark W. Longbrake* J. Huston McCulloch

Searching for the Output Gap: Economic Variable or Statistical Illusion? Mark W. Longbrake* J. Huston McCulloch Draft Draft Searching for the Output Gap: Economic Variable or Statistical Illusion? Mark W. Longbrake* The Ohio State University J. Huston McCulloch The Ohio State University August, 2007 Abstract This

More information

Purchasing power parity: A nonlinear multivariate perspective. Abstract

Purchasing power parity: A nonlinear multivariate perspective. Abstract Purchasing power parity: A nonlinear multivariate perspective Frédérique Bec THEMA, University of Cergy-Pontoise and CREST, France Mélika Ben Salem OEP, Paris-Est University and LEA-INRA (PSE), France

More information

DATABASE AND METHODOLOGY

DATABASE AND METHODOLOGY CHAPTER 3 DATABASE AND METHODOLOGY In the present chapter, sources of database used and methodology applied for the empirical analysis has been presented The whole chapter has been divided into three sections

More information

Econ671 Factor Models: Principal Components

Econ671 Factor Models: Principal Components Econ671 Factor Models: Principal Components Jun YU April 8, 2016 Jun YU () Econ671 Factor Models: Principal Components April 8, 2016 1 / 59 Factor Models: Principal Components Learning Objectives 1. Show

More information

Shortfalls of Panel Unit Root Testing. Jack Strauss Saint Louis University. And. Taner Yigit Bilkent University. Abstract

Shortfalls of Panel Unit Root Testing. Jack Strauss Saint Louis University. And. Taner Yigit Bilkent University. Abstract Shortfalls of Panel Unit Root Testing Jack Strauss Saint Louis University And Taner Yigit Bilkent University Abstract This paper shows that (i) magnitude and variation of contemporaneous correlation are

More information

Identifying Aggregate Liquidity Shocks with Monetary Policy Shocks: An Application using UK Data

Identifying Aggregate Liquidity Shocks with Monetary Policy Shocks: An Application using UK Data Identifying Aggregate Liquidity Shocks with Monetary Policy Shocks: An Application using UK Data Michael Ellington and Costas Milas Financial Services, Liquidity and Economic Activity Bank of England May

More information

Why has globalisation brought such large increases in exports to some countries and not to others?

Why has globalisation brought such large increases in exports to some countries and not to others? by Stephen Redding and Anthony Venables Why has globalisation brought such large increases in exports to some countries and not to others? Stephen Redding and Anthony Venables look at the way internal

More information

University of Pretoria Department of Economics Working Paper Series

University of Pretoria Department of Economics Working Paper Series University of Pretoria Department of Economics Working Paper Series Predicting Stock Returns and Volatility Using Consumption-Aggregate Wealth Ratios: A Nonlinear Approach Stelios Bekiros IPAG Business

More information

Extracting information from noisy time series data

Extracting information from noisy time series data Extracting information from noisy time series data Paul Ormerod (Pormerod@volterra.co.uk) Volterra Consulting Ltd Sheen Elms 135c Sheen Lane London SW14 8AE December 2004 1 Abstract A question which is

More information

Do Markov-Switching Models Capture Nonlinearities in the Data? Tests using Nonparametric Methods

Do Markov-Switching Models Capture Nonlinearities in the Data? Tests using Nonparametric Methods Do Markov-Switching Models Capture Nonlinearities in the Data? Tests using Nonparametric Methods Robert V. Breunig Centre for Economic Policy Research, Research School of Social Sciences and School of

More information

Discussion of Bootstrap prediction intervals for linear, nonlinear, and nonparametric autoregressions, by Li Pan and Dimitris Politis

Discussion of Bootstrap prediction intervals for linear, nonlinear, and nonparametric autoregressions, by Li Pan and Dimitris Politis Discussion of Bootstrap prediction intervals for linear, nonlinear, and nonparametric autoregressions, by Li Pan and Dimitris Politis Sílvia Gonçalves and Benoit Perron Département de sciences économiques,

More information

Euro-indicators Working Group

Euro-indicators Working Group Euro-indicators Working Group Luxembourg, 9 th & 10 th June 2011 Item 9.4 of the Agenda New developments in EuroMIND estimates Rosa Ruggeri Cannata Doc 309/11 What is EuroMIND? EuroMIND is a Monthly INDicator

More information

Explain the impact of location, climate, natural resources, and population distribution on Europe. a. Compare how the location, climate, and natural

Explain the impact of location, climate, natural resources, and population distribution on Europe. a. Compare how the location, climate, and natural SS6G10 Explain the impact of location, climate, natural resources, and population distribution on Europe. a. Compare how the location, climate, and natural resources of Germany, the United Kingdom and

More information

Discussion Papers on Entrepreneurship, Growth and Public Policy

Discussion Papers on Entrepreneurship, Growth and Public Policy Discussion Papers on Entrepreneurship, Growth and Public Policy # 2706 An Empirical Relationship between Entrepreneurship and FDI. A Note. by Miltiades N. Georgiou University of Thessaly Number of Pages:

More information

Econ 423 Lecture Notes: Additional Topics in Time Series 1

Econ 423 Lecture Notes: Additional Topics in Time Series 1 Econ 423 Lecture Notes: Additional Topics in Time Series 1 John C. Chao April 25, 2017 1 These notes are based in large part on Chapter 16 of Stock and Watson (2011). They are for instructional purposes

More information

Financial Factors in Economic Fluctuations. Lawrence Christiano Roberto Motto Massimo Rostagno

Financial Factors in Economic Fluctuations. Lawrence Christiano Roberto Motto Massimo Rostagno Financial Factors in Economic Fluctuations Lawrence Christiano Roberto Motto Massimo Rostagno Background Much progress made on constructing and estimating models that fit quarterly data well (Smets-Wouters,

More information

Warwick Business School Forecasting System. Summary. Ana Galvao, Anthony Garratt and James Mitchell November, 2014

Warwick Business School Forecasting System. Summary. Ana Galvao, Anthony Garratt and James Mitchell November, 2014 Warwick Business School Forecasting System Summary Ana Galvao, Anthony Garratt and James Mitchell November, 21 The main objective of the Warwick Business School Forecasting System is to provide competitive

More information

A Summary of Economic Methodology

A Summary of Economic Methodology A Summary of Economic Methodology I. The Methodology of Theoretical Economics All economic analysis begins with theory, based in part on intuitive insights that naturally spring from certain stylized facts,

More information

Network Connectivity and Systematic Risk

Network Connectivity and Systematic Risk Network Connectivity and Systematic Risk Monica Billio 1 Massimiliano Caporin 2 Roberto Panzica 3 Loriana Pelizzon 1,3 1 University Ca Foscari Venezia (Italy) 2 University of Padova (Italy) 3 Goethe University

More information

Transport as a Large Historical Factor. A Case Study in New Philosophy of History Dan Little University of Michigan-Dearborn

Transport as a Large Historical Factor. A Case Study in New Philosophy of History Dan Little University of Michigan-Dearborn Transport as a Large Historical Factor A Case Study in New Philosophy of History Dan Little University of Michigan-Dearborn New philosophy of history A backwater in philosophy An important set of questions

More information

CROSS-COUNTRY DIFFERENCES IN PRODUCTIVITY: THE ROLE OF ALLOCATION AND SELECTION

CROSS-COUNTRY DIFFERENCES IN PRODUCTIVITY: THE ROLE OF ALLOCATION AND SELECTION ONLINE APPENDIX CROSS-COUNTRY DIFFERENCES IN PRODUCTIVITY: THE ROLE OF ALLOCATION AND SELECTION By ERIC BARTELSMAN, JOHN HALTIWANGER AND STEFANO SCARPETTA This appendix presents a detailed sensitivity

More information

1.2. Structural VARs

1.2. Structural VARs 1. Shocks Nr. 1 1.2. Structural VARs How to identify A 0? A review: Choleski (Cholesky?) decompositions. Short-run restrictions. Inequality restrictions. Long-run restrictions. Then, examples, applications,

More information

Macroeconomics Theory II

Macroeconomics Theory II Macroeconomics Theory II Francesco Franco FEUNL February 2016 Francesco Franco Macroeconomics Theory II 1/23 Housekeeping. Class organization. Website with notes and papers as no "Mas-Collel" in macro

More information

Graduate Macro Theory II: Notes on Quantitative Analysis in DSGE Models

Graduate Macro Theory II: Notes on Quantitative Analysis in DSGE Models Graduate Macro Theory II: Notes on Quantitative Analysis in DSGE Models Eric Sims University of Notre Dame Spring 2011 This note describes very briefly how to conduct quantitative analysis on a linearized

More information

Estimating and Accounting for the Output Gap with Large Bayesian Vector Autoregressions

Estimating and Accounting for the Output Gap with Large Bayesian Vector Autoregressions Estimating and Accounting for the Output Gap with Large Bayesian Vector Autoregressions James Morley 1 Benjamin Wong 2 1 University of Sydney 2 Reserve Bank of New Zealand The view do not necessarily represent

More information

Online Appendix to: Crises and Recoveries in an Empirical Model of. Consumption Disasters

Online Appendix to: Crises and Recoveries in an Empirical Model of. Consumption Disasters Online Appendix to: Crises and Recoveries in an Empirical Model of Consumption Disasters Emi Nakamura Columbia University Robert Barro Harvard University Jón Steinsson Columbia University José Ursúa Harvard

More information

Business Cycle Comovements in Industrial Subsectors

Business Cycle Comovements in Industrial Subsectors Business Cycle Comovements in Industrial Subsectors Michael T. Owyang 1 Daniel Soques 2 1 Federal Reserve Bank of St. Louis 2 University of North Carolina Wilmington The views expressed here are those

More information

Econometría 2: Análisis de series de Tiempo

Econometría 2: Análisis de series de Tiempo Econometría 2: Análisis de series de Tiempo Karoll GOMEZ kgomezp@unal.edu.co http://karollgomez.wordpress.com Segundo semestre 2016 IX. Vector Time Series Models VARMA Models A. 1. Motivation: The vector

More information

A nonparametric test for path dependence in discrete panel data

A nonparametric test for path dependence in discrete panel data A nonparametric test for path dependence in discrete panel data Maximilian Kasy Department of Economics, University of California - Los Angeles, 8283 Bunche Hall, Mail Stop: 147703, Los Angeles, CA 90095,

More information

Global Value Chain Participation and Current Account Imbalances

Global Value Chain Participation and Current Account Imbalances Global Value Chain Participation and Current Account Imbalances Johannes Brumm University of Zurich Georgios Georgiadis European Central Bank Johannes Gräb European Central Bank Fabian Trottner Princeton

More information

Markov Switching Regular Vine Copulas

Markov Switching Regular Vine Copulas Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session CPS057) p.5304 Markov Switching Regular Vine Copulas Stöber, Jakob and Czado, Claudia Lehrstuhl für Mathematische Statistik,

More information

Cross validation of prediction models for seasonal time series by parametric bootstrapping

Cross validation of prediction models for seasonal time series by parametric bootstrapping Cross validation of prediction models for seasonal time series by parametric bootstrapping Robert M. Kunst robert.kunst@univie.ac.at University of Vienna and Institute for Advanced Studies Vienna Prepared

More information

High-dimensional Problems in Finance and Economics. Thomas M. Mertens

High-dimensional Problems in Finance and Economics. Thomas M. Mertens High-dimensional Problems in Finance and Economics Thomas M. Mertens NYU Stern Risk Economics Lab April 17, 2012 1 / 78 Motivation Many problems in finance and economics are high dimensional. Dynamic Optimization:

More information

GREAT BRITAIN: INDUSTRIAL REVOLUTION TO 1851 Student Worksheet

GREAT BRITAIN: INDUSTRIAL REVOLUTION TO 1851 Student Worksheet GREAT BRITAIN: INDUSTRIAL REVOLUTION TO 1851 Student Worksheet Theme of Map This thematic map captures the industrial development of Britain at a time when it had transformed the British economy and society.

More information

A Non-Parametric Approach of Heteroskedasticity Robust Estimation of Vector-Autoregressive (VAR) Models

A Non-Parametric Approach of Heteroskedasticity Robust Estimation of Vector-Autoregressive (VAR) Models Journal of Finance and Investment Analysis, vol.1, no.1, 2012, 55-67 ISSN: 2241-0988 (print version), 2241-0996 (online) International Scientific Press, 2012 A Non-Parametric Approach of Heteroskedasticity

More information

Another Look at the Boom and Bust of Financial Bubbles

Another Look at the Boom and Bust of Financial Bubbles ANNALS OF ECONOMICS AND FINANCE 16-2, 417 423 (2015) Another Look at the Boom and Bust of Financial Bubbles Andrea Beccarini University of Münster, Department of Economics, Am Stadtgraben 9, 48143, Münster,

More information

University of Pretoria Department of Economics Working Paper Series

University of Pretoria Department of Economics Working Paper Series University of Pretoria Department of Economics Working Paper Series On Exchange-Rate Movements and Gold-Price Fluctuations: Evidence for Gold- Producing Countries from a Nonparametric Causality-in-Quantiles

More information

Urban Foundations. Early American Cities. Early American Cities. Early American Cities. Cities in America to 1945

Urban Foundations. Early American Cities. Early American Cities. Early American Cities. Cities in America to 1945 Urban Foundations Cities in America to 1945 Early American Cities The design of early American cities was highly influenced by European traditions The roots of many of these traditions can in turn be traced

More information

TAKEHOME FINAL EXAM e iω e 2iω e iω e 2iω

TAKEHOME FINAL EXAM e iω e 2iω e iω e 2iω ECO 513 Spring 2015 TAKEHOME FINAL EXAM (1) Suppose the univariate stochastic process y is ARMA(2,2) of the following form: y t = 1.6974y t 1.9604y t 2 + ε t 1.6628ε t 1 +.9216ε t 2, (1) where ε is i.i.d.

More information

Cross-sectional space-time modeling using ARNN(p, n) processes

Cross-sectional space-time modeling using ARNN(p, n) processes Cross-sectional space-time modeling using ARNN(p, n) processes W. Polasek K. Kakamu September, 006 Abstract We suggest a new class of cross-sectional space-time models based on local AR models and nearest

More information

Appendix 1. Grain price and tariff data

Appendix 1. Grain price and tariff data Appendix 1. Grain price and tariff data Grain prices The British, French, German and US price data are based on the Board of Trade s (1909) Statistical Tables and Charts Relating to British and Foreign

More information

Signaling Effects of Monetary Policy

Signaling Effects of Monetary Policy Signaling Effects of Monetary Policy Leonardo Melosi London Business School 24 May 2012 Motivation Disperse information about aggregate fundamentals Morris and Shin (2003), Sims (2003), and Woodford (2002)

More information

Using all observations when forecasting under structural breaks

Using all observations when forecasting under structural breaks Using all observations when forecasting under structural breaks Stanislav Anatolyev New Economic School Victor Kitov Moscow State University December 2007 Abstract We extend the idea of the trade-off window

More information

Lecture 2. Business Cycle Measurement. Randall Romero Aguilar, PhD II Semestre 2017 Last updated: August 18, 2017

Lecture 2. Business Cycle Measurement. Randall Romero Aguilar, PhD II Semestre 2017 Last updated: August 18, 2017 Lecture 2 Business Cycle Measurement Randall Romero Aguilar, PhD II Semestre 2017 Last updated: August 18, 2017 Universidad de Costa Rica EC3201 - Teoría Macroeconómica 2 Table of contents 1. Introduction

More information

Fresh perspectives on unobservable variables: Data decomposition of the Kalman smoother

Fresh perspectives on unobservable variables: Data decomposition of the Kalman smoother Fresh perspectives on unobservable variables: Data decomposition of the Kalman smoother AN 2013/09 Nicholas Sander December 2013 Reserve Bank of New Zealand Analytical Note series ISSN 2230 5505 Reserve

More information

Identifying Financial Risk Factors

Identifying Financial Risk Factors Identifying Financial Risk Factors with a Low-Rank Sparse Decomposition Lisa Goldberg Alex Shkolnik Berkeley Columbia Meeting in Engineering and Statistics 24 March 2016 Outline 1 A Brief History of Factor

More information

FaMIDAS: A Mixed Frequency Factor Model with MIDAS structure

FaMIDAS: A Mixed Frequency Factor Model with MIDAS structure FaMIDAS: A Mixed Frequency Factor Model with MIDAS structure Frale C., Monteforte L. Computational and Financial Econometrics Limassol, October 2009 Introduction After the recent financial and economic

More information

Multivariate GARCH models.

Multivariate GARCH models. Multivariate GARCH models. Financial market volatility moves together over time across assets and markets. Recognizing this commonality through a multivariate modeling framework leads to obvious gains

More information

Global Macroeconomic Uncertainty

Global Macroeconomic Uncertainty Global Macroeconomic Uncertainty Tino Berger, University of Cologne, CMR Sibylle Herz, University of Muenster Global Macroeconomic Uncertainty 1 / 1 Goal of the paper 1. Measure global macroeconomic uncertainty

More information

MEASURING THE CORRELATION OF SHOCKS BETWEEN THE UK AND THE CORE OF EUROPE

MEASURING THE CORRELATION OF SHOCKS BETWEEN THE UK AND THE CORE OF EUROPE Measuring the Correlation of Shocks Between the UK and the Core of Europe 2 MEASURING THE CORRELATION OF SHOCKS BETWEEN THE UK AND THE CORE OF EUROPE Abstract S.G.HALL * B. YHAP ** This paper considers

More information

COMMENT ON DEL NEGRO, SCHORFHEIDE, SMETS AND WOUTERS COMMENT ON DEL NEGRO, SCHORFHEIDE, SMETS AND WOUTERS

COMMENT ON DEL NEGRO, SCHORFHEIDE, SMETS AND WOUTERS COMMENT ON DEL NEGRO, SCHORFHEIDE, SMETS AND WOUTERS COMMENT ON DEL NEGRO, SCHORFHEIDE, SMETS AND WOUTERS COMMENT ON DEL NEGRO, SCHORFHEIDE, SMETS AND WOUTERS CHRISTOPHER A. SIMS 1. WHY THIS APPROACH HAS BEEN SUCCESSFUL This paper sets out to blend the advantages

More information

Financial Econometrics

Financial Econometrics Financial Econometrics Multivariate Time Series Analysis: VAR Gerald P. Dwyer Trinity College, Dublin January 2013 GPD (TCD) VAR 01/13 1 / 25 Structural equations Suppose have simultaneous system for supply

More information

Does agglomeration explain regional income inequalities?

Does agglomeration explain regional income inequalities? Does agglomeration explain regional income inequalities? Karen Helene Midelfart Norwegian School of Economics and Business Administration and CEPR August 31, 2004 First draft Abstract This paper seeks

More information

Department of Economics, UCSB UC Santa Barbara

Department of Economics, UCSB UC Santa Barbara Department of Economics, UCSB UC Santa Barbara Title: Past trend versus future expectation: test of exchange rate volatility Author: Sengupta, Jati K., University of California, Santa Barbara Sfeir, Raymond,

More information

On the Distribution Dynamics of Human Development: Evidence from the Metropolitan Regions of Bolivia

On the Distribution Dynamics of Human Development: Evidence from the Metropolitan Regions of Bolivia MPRA Munich Personal RePEc Archive On the Distribution Dynamics of Human Development: Evidence from the Metropolitan Regions of Bolivia Carlos Mendez Kyushu University 24 June 2018 Online at https://mpra.ub.uni-muenchen.de/87571/

More information

An Empirical Analysis of RMB Exchange Rate changes impact on PPI of China

An Empirical Analysis of RMB Exchange Rate changes impact on PPI of China 2nd International Conference on Economics, Management Engineering and Education Technology (ICEMEET 206) An Empirical Analysis of RMB Exchange Rate changes impact on PPI of China Chao Li, a and Yonghua

More information

Inflation Revisited: New Evidence from Modified Unit Root Tests

Inflation Revisited: New Evidence from Modified Unit Root Tests 1 Inflation Revisited: New Evidence from Modified Unit Root Tests Walter Enders and Yu Liu * University of Alabama in Tuscaloosa and University of Texas at El Paso Abstract: We propose a simple modification

More information

Linear Models in Econometrics

Linear Models in Econometrics Linear Models in Econometrics Nicky Grant At the most fundamental level econometrics is the development of statistical techniques suited primarily to answering economic questions and testing economic theories.

More information

4 Shocks and Disturbances to the World Economy

4 Shocks and Disturbances to the World Economy 4 Shocks and Disturbances to the World Economy It is a fact of life that the world economy is subject to continual disturbances or shocks. Some disturbances are large and noticeable, such as the sudden

More information

The regression model with one stochastic regressor (part II)

The regression model with one stochastic regressor (part II) The regression model with one stochastic regressor (part II) 3150/4150 Lecture 7 Ragnar Nymoen 6 Feb 2012 We will finish Lecture topic 4: The regression model with stochastic regressor We will first look

More information

Frederick Wallace Universidad de Quintana Roo. Abstract

Frederick Wallace Universidad de Quintana Roo. Abstract Nonlinear unit root tests of PPP using long-horizon data Frederick Wallace Universidad de Quintana Roo Abstract The Kapetanios, Shin, and Snell (KSS, 2003) test for a nonlinear unit root is used to study

More information

Notes on Random Variables, Expectations, Probability Densities, and Martingales

Notes on Random Variables, Expectations, Probability Densities, and Martingales Eco 315.2 Spring 2006 C.Sims Notes on Random Variables, Expectations, Probability Densities, and Martingales Includes Exercise Due Tuesday, April 4. For many or most of you, parts of these notes will be

More information

Wooldridge, Introductory Econometrics, 4th ed. Chapter 15: Instrumental variables and two stage least squares

Wooldridge, Introductory Econometrics, 4th ed. Chapter 15: Instrumental variables and two stage least squares Wooldridge, Introductory Econometrics, 4th ed. Chapter 15: Instrumental variables and two stage least squares Many economic models involve endogeneity: that is, a theoretical relationship does not fit

More information

Chapter 1. Introduction. 1.1 Background

Chapter 1. Introduction. 1.1 Background Chapter 1 Introduction Science is facts; just as houses are made of stones, so is science made of facts; but a pile of stones is not a house and a collection of facts is not necessarily science. Henri

More information

Will it float? The New Keynesian Phillips curve tested on OECD panel data

Will it float? The New Keynesian Phillips curve tested on OECD panel data Phillips curve Roger Bjørnstad 1 2 1 Research Department Statistics Norway 2 Department of Economics University of Oslo 31 October 2006 Outline Outline Outline Outline Outline The debatable The hybrid

More information

Introduction to Econometrics

Introduction to Econometrics Introduction to Econometrics T H I R D E D I T I O N Global Edition James H. Stock Harvard University Mark W. Watson Princeton University Boston Columbus Indianapolis New York San Francisco Upper Saddle

More information

Cross-Country Differences in Productivity: The Role of Allocation and Selection

Cross-Country Differences in Productivity: The Role of Allocation and Selection Cross-Country Differences in Productivity: The Role of Allocation and Selection Eric Bartelsman, John Haltiwanger & Stefano Scarpetta American Economic Review (2013) Presented by Beatriz González January

More information

Business Cycle Dating Committee of the Centre for Economic Policy Research. 1. The CEPR Business Cycle Dating Committee

Business Cycle Dating Committee of the Centre for Economic Policy Research. 1. The CEPR Business Cycle Dating Committee Business Cycle Dating Committee of the Centre for Economic Policy Research Michael Artis Fabio Canova Jordi Gali Francesco Giavazzi Richard Portes (President, CEPR) Lucrezia Reichlin (Chair) Harald Uhlig

More information

Wage and price setting. Slides for 26. August 2003 lecture

Wage and price setting. Slides for 26. August 2003 lecture 1 B&W s derivation of the Phillips curve Wage and price setting. Slides for 26. August 2003 lecture Ragnar Nymoen University of Oslo, Department of Economics Ch 12.3: The Battle of the mark-ups as a framework

More information

The Dornbusch overshooting model

The Dornbusch overshooting model 4330 Lecture 8 Ragnar Nymoen 12 March 2012 References I Lecture 7: Portfolio model of the FEX market extended by money. Important concepts: monetary policy regimes degree of sterilization Monetary model

More information

Inference in VARs with Conditional Heteroskedasticity of Unknown Form

Inference in VARs with Conditional Heteroskedasticity of Unknown Form Inference in VARs with Conditional Heteroskedasticity of Unknown Form Ralf Brüggemann a Carsten Jentsch b Carsten Trenkler c University of Konstanz University of Mannheim University of Mannheim IAB Nuremberg

More information

The Mixture Approach for Simulating New Families of Bivariate Distributions with Specified Correlations

The Mixture Approach for Simulating New Families of Bivariate Distributions with Specified Correlations The Mixture Approach for Simulating New Families of Bivariate Distributions with Specified Correlations John R. Michael, Significance, Inc. and William R. Schucany, Southern Methodist University The mixture

More information

DISCUSSION PAPERS IN ECONOMICS

DISCUSSION PAPERS IN ECONOMICS STRATHCLYDE DISCUSSION PAPERS IN ECONOMICS UK HOUSE PRICES: CONVERGENCE CLUBS AND SPILLOVERS BY ALBERTO MONTAGNOLI AND JUN NAGAYASU NO 13-22 DEPARTMENT OF ECONOMICS UNIVERSITY OF STRATHCLYDE GLASGOW UK

More information

U n iversity o f H ei delberg

U n iversity o f H ei delberg U n iversity o f H ei delberg Department of Economics Discussion Paper Series No. 585 482482 Global Prediction of Recessions Jonas Dovern and Florian Huber March 2015 Global Prediction of Recessions Jonas

More information

More on Roy Model of Self-Selection

More on Roy Model of Self-Selection V. J. Hotz Rev. May 26, 2007 More on Roy Model of Self-Selection Results drawn on Heckman and Sedlacek JPE, 1985 and Heckman and Honoré, Econometrica, 1986. Two-sector model in which: Agents are income

More information

Parity Reversion of Absolute Purchasing Power Parity Zhi-bai ZHANG 1,a,* and Zhi-cun BIAN 2,b

Parity Reversion of Absolute Purchasing Power Parity Zhi-bai ZHANG 1,a,* and Zhi-cun BIAN 2,b 2016 3 rd International Conference on Economics and Management (ICEM 2016) ISBN: 978-1-60595-368-7 Parity Reversion of Absolute Purchasing Power Parity Zhi-bai ZHANG 1,a,* and Zhi-cun BIAN 2,b 1,2 School

More information

ECO 513 Fall 2009 C. Sims HIDDEN MARKOV CHAIN MODELS

ECO 513 Fall 2009 C. Sims HIDDEN MARKOV CHAIN MODELS ECO 513 Fall 2009 C. Sims HIDDEN MARKOV CHAIN MODELS 1. THE CLASS OF MODELS y t {y s, s < t} p(y t θ t, {y s, s < t}) θ t = θ(s t ) P[S t = i S t 1 = j] = h ij. 2. WHAT S HANDY ABOUT IT Evaluating the

More information

Chapter 4. Explanation of the Model. Satoru Kumagai Inter-disciplinary Studies, IDE-JETRO, Japan

Chapter 4. Explanation of the Model. Satoru Kumagai Inter-disciplinary Studies, IDE-JETRO, Japan Chapter 4 Explanation of the Model Satoru Kumagai Inter-disciplinary Studies, IDE-JETRO, Japan Toshitaka Gokan Inter-disciplinary Studies, IDE-JETRO, Japan Ikumo Isono Bangkok Research Center, IDE-JETRO,

More information

Wooldridge, Introductory Econometrics, 3d ed. Chapter 9: More on specification and data problems

Wooldridge, Introductory Econometrics, 3d ed. Chapter 9: More on specification and data problems Wooldridge, Introductory Econometrics, 3d ed. Chapter 9: More on specification and data problems Functional form misspecification We may have a model that is correctly specified, in terms of including

More information

Finnancial Development and Growth

Finnancial Development and Growth Finnancial Development and Growth Econometrics Prof. Menelaos Karanasos Brunel University December 4, 2012 (Institute Annual historical data for Brazil December 4, 2012 1 / 34 Finnancial Development and

More information

Short T Panels - Review

Short T Panels - Review Short T Panels - Review We have looked at methods for estimating parameters on time-varying explanatory variables consistently in panels with many cross-section observation units but a small number of

More information

Joint Estimation of Risk Preferences and Technology: Further Discussion

Joint Estimation of Risk Preferences and Technology: Further Discussion Joint Estimation of Risk Preferences and Technology: Further Discussion Feng Wu Research Associate Gulf Coast Research and Education Center University of Florida Zhengfei Guan Assistant Professor Gulf

More information

Identifying SVARs with Sign Restrictions and Heteroskedasticity

Identifying SVARs with Sign Restrictions and Heteroskedasticity Identifying SVARs with Sign Restrictions and Heteroskedasticity Srečko Zimic VERY PRELIMINARY AND INCOMPLETE NOT FOR DISTRIBUTION February 13, 217 Abstract This paper introduces a new method to identify

More information

EMERGING MARKETS - Lecture 2: Methodology refresher

EMERGING MARKETS - Lecture 2: Methodology refresher EMERGING MARKETS - Lecture 2: Methodology refresher Maria Perrotta April 4, 2013 SITE http://www.hhs.se/site/pages/default.aspx My contact: maria.perrotta@hhs.se Aim of this class There are many different

More information

A Guide to Modern Econometric:

A Guide to Modern Econometric: A Guide to Modern Econometric: 4th edition Marno Verbeek Rotterdam School of Management, Erasmus University, Rotterdam B 379887 )WILEY A John Wiley & Sons, Ltd., Publication Contents Preface xiii 1 Introduction

More information

DEPARTMENT OF ECONOMICS

DEPARTMENT OF ECONOMICS ISSN 0819-64 ISBN 0 7340 616 1 THE UNIVERSITY OF MELBOURNE DEPARTMENT OF ECONOMICS RESEARCH PAPER NUMBER 959 FEBRUARY 006 TESTING FOR RATE-DEPENDENCE AND ASYMMETRY IN INFLATION UNCERTAINTY: EVIDENCE FROM

More information

Are PPP Tests Erratically Behaved? Some Panel Evidence

Are PPP Tests Erratically Behaved? Some Panel Evidence Are PPP Tests Erratically Behaved? Some Panel Evidence Guglielmo Maria Caporale a, Christoph Hanck b a Brunel University, London b Universität Dortmund October 5, 2006 Abstract This paper examines whether,

More information

State-space Model. Eduardo Rossi University of Pavia. November Rossi State-space Model Fin. Econometrics / 53

State-space Model. Eduardo Rossi University of Pavia. November Rossi State-space Model Fin. Econometrics / 53 State-space Model Eduardo Rossi University of Pavia November 2014 Rossi State-space Model Fin. Econometrics - 2014 1 / 53 Outline 1 Motivation 2 Introduction 3 The Kalman filter 4 Forecast errors 5 State

More information

Time Series Analysis. James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY

Time Series Analysis. James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY Time Series Analysis James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY PREFACE xiii 1 Difference Equations 1.1. First-Order Difference Equations 1 1.2. pth-order Difference Equations 7

More information

ENERGY CONSUMPTION AND ECONOMIC GROWTH IN SWEDEN: A LEVERAGED BOOTSTRAP APPROACH, ( ) HATEMI-J, Abdulnasser * IRANDOUST, Manuchehr

ENERGY CONSUMPTION AND ECONOMIC GROWTH IN SWEDEN: A LEVERAGED BOOTSTRAP APPROACH, ( ) HATEMI-J, Abdulnasser * IRANDOUST, Manuchehr ENERGY CONSUMPTION AND ECONOMIC GROWTH IN SWEDEN: A LEVERAGED BOOTSTRAP APPROACH, (1965-2000) HATEMI-J, Abdulnasser * IRANDOUST, Manuchehr Abstract The causal interaction between energy consumption, real

More information

A Bootstrap Test for Causality with Endogenous Lag Length Choice. - theory and application in finance

A Bootstrap Test for Causality with Endogenous Lag Length Choice. - theory and application in finance CESIS Electronic Working Paper Series Paper No. 223 A Bootstrap Test for Causality with Endogenous Lag Length Choice - theory and application in finance R. Scott Hacker and Abdulnasser Hatemi-J April 200

More information