Firm Heterogeneity and Aggregate Fluctuations: a Functional VAR model for multidimensional distributions