Data-based Automatic Discretization of Nonparametric Distributions
|
|
- Ralf Newman
- 6 years ago
- Views:
Transcription
1 Data-based Automatic Discretization of Nonparametric Distributions Alexis Akira Toda arxiv: v1 [q-fin.ec] 2 May 2018 May 3, 2018 Abstract Although using non-gaussian distributions in economic models has become increasingly popular, currently there is no systematic way for calibrating a discrete distribution from the data without imposing parametric assumptions. This paper proposes a simple nonparametric calibration method that combines the kernel density estimator and Gaussian quadrature. Applications to an asset pricing model and an optimal portfolio problem suggest that assuming Gaussian instead of nonparametric shocks leads to a 15 25% reduction in the equity premium and overweighting in the stock portfolio because the investor underestimates the probability of crashes. Keywords: calibration, discrete approximation, Gaussian mixture, Gaussian quadrature, kernel density estimator. JEL codes: C63, C65, G11. 1 Introduction This paper studies the following problem, which applied theorists often encounter. A researcher would like to calibrate the parameters of a stochastic model. One of the model inputs is a probability distribution of shocks, which is to be approximated by a discrete distribution. Due to computational considerations, the researcher would like this distribution to have as few support points (nodes) as possible, say five. Given the data of shocks, how should the researcher calibrate the nodes and probabilities of this five-point distribution? While there are many established methods for discretizing processes with Gaussian shocks such as Tauchen(1986), Tauchen and Hussey(1991), and Rouwenhorst (1995), 1 discretizing non-gaussian distributions remains relatively unexplored. However, it has become increasingly common in economics to study models with non-gaussian shocks. For example, the rare disasters model(rietz, 1988; Barro, 2006; Gabaix, 2012) uses rare but large downward jumps to explain asset pricing puzzles. One issue with discretizing non-gaussian distributions is how to calibrate them. If we have a parametric density, it is possible to discretize it using the Gaussian quadrature as in Miller and Rice (1983) or the maximum entropy Department of Economics, University of California San Diego. atoda@ucsd.edu. 1 See Farmer and Toda (2017) and the references therein for a detailed literature review. 1
2 methodasintanaka and Toda(2013,2015)providedthatwecancomputesome moments. However, it is not obvious how to obtain an N-point distribution that approximates the data well without imposing parametric assumptions. Given the data, this paper proposes a new and simple method for automatically calibrating a discrete distribution with a specified number of grid points. The method combines two well-known techniques, the (Gaussian) kernel density estimation and Gaussian quadrature. 2 In the first step, we nonparametrically estimate the probability distribution using the kernel density estimator with the Gaussian (normal) kernel. Because the kernel density estimator can be viewed as a Gaussian mixture, we can compute the moments of order up to 2N efficiently. Using these moments, in the second step we compute the nodes and weights of the N-point Gaussian quadrature. These nodes and weights give us the N-point discretization of the nonparametric distribution. Since this method does not involve optimization (it is a matter of recursively computing moments and solving for the eigenvalues/vectors of a sparse matrix), the implementation is easy and fast. As applications, I discretize the U.S. historical data on aggregate consumption growth and stock returns and solve an asset pricing model and an optimal portfolio problem with constant relative risk aversion utility. I consider two cases in which the investor uses the nonparametric and Gaussian densities. I show that when the investor incorrectly believes that the consumption growth/stock returns distribution is lognormal, the equity premium shrinks by 16 24% and the stock portfolio is overweighted by 16 26% because he underestimates the probability of recessions/crashes. These examples show that the choice of the calibration method may matter quantitatively. 1.1 Related literature The closest paper to mine is Miller and Rice (1983), who use the Gaussian quadrature to discretize distributions. While they consider only the discretization of parametric distributions, my focus is on the discretization of nonparametric distributions estimated from data. Tanaka and Toda (2013) consider the discretization of distributions on preassigned nodes by matching the moments using the maximum entropy principle, and Tanaka and Toda(2015) prove convergence and obtain an error estimate. Farmer and Toda (2017) consider the discretization of general non-gaussian Markov processes by applying the Tanaka-Toda method to conditional distributions. In one of the applications, they discretize a nonparametric density on a preassigned grid by approximating it with a Gaussian mixture. Since computing the kernel density estimator does not require optimizing over parameters (unlike the maximum likelihood estimation of Gaussian mixture parameters or solving the maximum entropy problem), my method is easier and faster to implement, and the grid is chosen endogenously. 2 In this paper the term Gaussian has two different meanings. Gaussian quadrature refers to the quadrature (numerical integration) formula that integrates the maximum number of monomials. Gaussian mixture/distribution refers to the normal mixture/distribution in statistics. The widely used Gaussian quadrature for the Gaussian (normal) distribution is called Gauss-Hermite quadrature. 2
3 2 Method 2.1 Discretizing a density using Gaussian quadrature Suppose for the moment that the nonparametric density f(x) is known. Since stochasticmodelsofteninvolveexpectations,wewouldliketofindnodes{x n } N n=1 and weights {w n } N n=1 such that E[g(X)] = g(x)f(x)dx w n g(x n ), (2.1) where g is a general integrand and X is a random variable with density f(x). The right-hand side of (2.1) defines an N-point quadrature formula. When (2.1) is exact (i.e., becomes =) for all polynomials of degree up to D, we say that the quadrature formula has degree of exactness D. Since the degree of freedom in an N-point quadrature formula is 2N (because there are N nodes and N weights), we cannot expect to integrate more than 2N monomials f(x) = 1,x,...,x 2N 1 exactly. When the quadrature formula (2.1) is exact for these monomials, or equivalently when it has degree of exactness 2N 1, we call the formula Gaussian. The following Golub-Welsch algorithm provides an efficient way to compute the nodes and weights of the Gaussian quadrature. (Appendix B provides more theoretical background.) n=1 Algorithm 1 (Golub and Welsch, 1969). 1. Select a number of quadrature nodes N N. 2. For k = 0,1,...,2N, compute the k-th moment of the density m k = x k f(x)dx. 3. DefinethematrixofmomentsM = (M ij ) 1 i,j N+1 bym ij = m i+j Compute the Cholesky factorization M = R R. Let R = (r ij ) 1 i,j N Define α 1 = r 12 /r 11, α n = rn,n+1 r nn rn 1,n r n 1,n 1 (n = 2,...,N), and β n = rn+1,n+1 r nn (n = 1,...,N 1). Define the N N symmetric tridiagonal matrix α 1 β β 1 α 2 β.. 2. T N =. 0 β 2 α (2.2) βn β N 1 α N 6. Compute the eigenvalues {x n } N n=1 of T N and the correspondingeigenvectors {v n } N n=1. {x n} N n=1 are the nodes of the Gaussian quadrature and the weights {w n } N n=1 are given by w n = m 0 vn1 2 / v n 2 > 0, where v n = (v n1,...,v nn ). 3
4 Once we compute the nodes {x n } N n=1 and weights {w n} N n=1, we can use them as the discrete approximation of the density f. 2.2 Computing moments using kernel density estimator The only inputs to the Golub-Welsch algorithm are the number of nodes N and the moments m k = x k f(x)dx of the density f, where k = 0,...,2N. We can compute these quantities using the kernel density estimator, as follows. Given the data {x i } I i=1, the kernel density estimator (with Gaussian kernel) is defined by f h (x) = 1 Ih I ( ) x xi φ, (2.3) h i=1 where h > 0 is the bandwidth and φ(x) = 1 2π e x2 /2 is the standard normal density. I do not take a stance on the choice of the bandwidth h. When f is itself Gaussian, it has been shown that ( ) 1/5 4 h = σ (2.4) 3I is optimal in the sense of the mean integrated squared error, where σ is the sample standard deviation (Silverman, 1986). Note that the kernel density estimator(2.3) can be viewed as a Gaussian mixture density with I components, where the proportion, mean, and standard deviation of the i-th component are 1/I, x i, and h, respectively. As such, it is straightforward to compute its moments recursively. First, assume that the number of mixtures is 1, so the Gaussian mixture is just N(µ,σ 2 ). Let m k = E[X k ] be the k-th noncentral moment. Clearly m 0 = 1 and m 1 = µ. For k 2, using integration by parts we can compute m k = = =µ x k 1 2πσ 2 e (x µ)2 /2σ 2 dx x k 1 1 [ σ 2( e (x µ)2 /2σ 2) ] +µe (x µ) 2 /2σ 2 dx 2πσ 2 x k 1 1 2πσ 2 e (x µ)2 /2σ 2 dx σ 2 [ [ x k 1 2πσ 2 e (x µ)2 /2σ 2 ] ] (x k ) 1 /2σ 2 dx 2πσ 2 e (x µ)2 =µm k 1 +σ 2 (k 1)m k 2. (2.5) Therefore we can iterate this equation to compute m 2,m 3,...,m 2N. The k-th moment of the kernel density estimator (2.3) is m k = x k fh (x)dx = 1 I I m k (x i,h), (2.6) where m k (µ,σ) is the k-th moment of N(µ,σ 2 ) computed as in (2.5). Since the moments of the Gaussian mixture can be explicitly calculated this way, it is straightforward to compute the nodes and weights of the corresponding Gaussian quadrature using the Golub-Welsch algorithm 1. i=1 4
5 2.3 Summary of algorithm Summarizing the above observations, we obtain the following algorithm for the data-based automatic discretization of nonparametric distributions. Algorithm 2 (Automatic discretization of nonparametric distributions). 1. Given the data {x i } I i=1, choose a bandwidth h (for example, Silverman s rule of thumb (2.4)). 2. Given the number of discrete points N, for k = 0,...,2N compute the k-th moment of the kernel density estimator recursively as in (2.6) and (2.5). 3. Usethese moments {m k } 2N k=0 in the Golub-Welschalgorithm1to compute the nodes { x n } N n=1 and weights{w n} N n=1. The desired discretization assigns probability w n on the point x n. Appendix A shows that the accuracy of the proposed method exceeds that of using a parametric distribution when the latter is misspecified. 3 Applications In this section I illustrate the usefulness of the proposed method using minimal economic examples. In the first application, I solve a consumption-based asset pricing model by calibrating the consumption growth distribution from the data. In the second application, I solve a static optimal portfolio problem by calibrating the stock returns distribution from the data. In each case I compare the nonparametric solution to the case when the distribution is assumed to be lognormal and show that the choice of the discretization method matters quantitatively. 3.1 Asset pricing model Consider a plain-vanilla Lucas (1978) consumption-based asset pricing model, where a representative agent with constant relative risk aversion (CRRA) utility consumes the aggregate endowment. Endowment growth is assumed to be independent and identically distributed (i.i.d.) over time with a general distribution. A stock is modeled as a claim to the aggregate endowment, and the supply of the risk-free asset is zero. Let R > 0 be the gross stock return (which is an i.i.d. random variable) and R f > 0 be the gross risk-free rate. Letting γ > 0 be the relative risk aversion coefficient of the representative agent, it is well known that the log equity premium is given by E[logR] logr f = logm(1)+logm( γ) logm(1 γ) > 0, (3.1) 5
6 where M(s) = E[e sxt ] is the moment generating function of the log endowment growth x t = log(e t /e t 1 ). 3 I obtain the annual data on U.S. real per capita consumption for the period from the spreadsheet of Robert Shiller 4 and apply Algorithm 2 to the log consumption growth logc t /C t 1 to obtain a discrete distribution with nodes { x n } N n=1 and weights {w n} N n=1, where I choose the number of points N = 5 (increasing the number of points further does not change the results). I then compute the moment generating function as M(s) = N n=1 w ne s xn and substitute into (3.1) for s = 1, γ,1 γ to compute the log equity premium for the range γ (0,10]. Probability density Histogram Nonparametric Gaussian Log equity premium (%) Nonparametric Gaussian Log consumption growth (a) Histogram and densities Relative risk aversion (b) Equity premium. 25 Equity premium error (%) Relative risk aversion (c) Relative error. Figure 1: Consumption growth distribution and equity premium. Figure 1a shows the histogram of the log consumption growth distribution as well as the nonparametric kernel density estimator and the Gaussian distribution fitted by maximum likelihood. We can see that the histogram and the nonparametric density have slightly longer tail (expansions and recessions) than the Gaussian distribution. Figure 1b shows the log equity premium (3.1) for the two models. 5 With nonparametric shocks, the equity premium increases by 3 Mylecturenote athttps://sites.google.com/site/aatoda111/file-cabinet/272_l04.pdf contains a proof of (3.1). With i.i.d. endowment growth, it is easy to show that the log equity premium depends only on the risk aversion and it does not matter whether the utility is additive of recursive (Epstein-Zin) To make the results comparable, I discretize the Gaussian density using the Golub-Welsch algorithm with the corresponding Gaussian density, which is equivalent to the well-known 6
7 0.3 percentage points when γ = 10. Figure 1c shows the relative increase in the equity premium when using the nonparametric distribution, which ranges 16 24%. Thus using nonparametric distributions help to explain asset pricing puzzles. 3.2 Optimal portfolio problem Consider a CRRA investor with relative risk aversion γ > 0. Letting R > 0 be the gross stock return, R f > 0 be the gross risk-free rate, and θ be the fraction of wealth (portfolio share) invested in the stock, the investor s optimal portfolio problem is max θ 1 1 γ E[(Rθ +R f(1 θ)) 1 γ ]. (3.2) I obtain the annual data on U.S. nominal stock returns, risk-free rate, and inflation for the period from the spreadsheet of Amit Goyal. 6 For the stock returns I use the CRSP volume-weighted index including dividends. I convert these returns into real log returns and calibrate the log risk free rate as the sample average. The result is R f = I then apply Algorithm 2 to the log excess returns logr logr f to obtain a discrete distribution with nodes { x n } N n=1 and weights {w n} N n=1, where I choose the number of points N = 5 (increasing the number of points further does not change the results). The gross stock return in state n is defined by R n = R f e xn, which occurs with probability w n. Finally, I numerically solve the optimal portfolio problem (3.2). Figure 2 shows the results when we changethe relative risk aversionin the range γ [1,7]. Figure 2a shows the histogram of the log excess returns distribution as well as the nonparametric kernel density estimator and the Gaussian distribution fitted by maximum likelihood. We can see that the histogram and the nonparametric density have a long left tail corresponding to stock market crashes, which the Gaussian distribution misses. Figure 2b shows the optimal portfolio θ for the two models. We can see that when the investor incorrectly believes that the stock returns distribution is lognormal, he overweights the stock portfolio because he underestimates the probability of crashes. Figure 2c shows the percentage of this overweight(portfolioerror)θ G /θ NP 1, wheregand NPstandforGaussian and nonparametric densities. The portfolio error is substantial, in the range of 16 26%. 4 Concluding remarks This paper has proposed a simple, automatic method for discretizing a nonparametric distribution, given the data. Using an asset pricing model and an optimal portfolio problem as a laboratory, I showed that the error from using a parametric distribution (such as the Gaussian distribution) can be substantial. A natural extension is to consider the discretization of Markov processes with nonparametric shocks. For example, one may be tempted to apply the Gauss-Hermite quadrature. 6 The spreadsheet is at Using monthly or quarterly data give qualitatively similar results, though slightly less extreme quantitatively. 7
8 Probability density Histogram Nonparametric Gaussian Optimal portfolio Nonparametric Gaussian Log excess returns (a) Histogram and densities Relative risk aversion (b) Optimal portfolio. 25 Portfolio error (%) Relative risk aversion (c) Portfolio error. Figure 2: Excess returns distribution and numerical solution. kernel density estimation and Gaussian quadrature in the Tauchen and Hussey (1991) method to discretize the AR(1) process x t = ρx t 1 +ε t, where ρ < 1 and the innovations {ε t } t=0 are independent and identically distributed according to some probability density function f. However, it is well known that the Tauchen-Hussey method is not accurate when the persistence ρ is moderately high (Flodén, 2008). Using the AR(1) asset pricing model in the Online Appendix of Farmer and Toda (2017) to evaluate the solution accuracy, I found that the Gaussian quadrature-based methods for discretizing Markov processes is even less accurate when the shock distribution is nonparametric. Therefore for such processes, it is preferable to use the Farmer and Toda (2017) maximum entropy method with an even-spaced grid (see their Section for an example). Fortunately, since the inputs to the Farmer-Toda method are the probability density function and some low order moments (which we can compute exactly as in Section 2.2), the proposed method based on the kernel density estimator is still applicable. Finally, although I proposed my method as a tool for discretization, it can also be used as a quadrature method. For example, Pohl et al. (2018) solve the Bansal and Yaron (2004) long run risks model using the projection method and Gauss-Hermite quadrature, but that is because the model is assumed to have Gaussian shocks. If instead a researcher wishes to use nonparametric shocks, my method can be directly used to construct a quadrature rule from data. 8
9 References Ravi Bansal and Amir Yaron. Risks for the long run: A potential resolution of asset pricing puzzles. Journal of Finance, 59(4): , August doi: /j x. Robert J. Barro. Rare disasters and asset markets in the twentieth century. Quarterly Journal of Economics, 121(3): , doi: /qjec Leland E. Farmer and Alexis Akira Toda. Discretizing nonlinear, non-gaussian Markov processes with exact conditional moments. Quantitative Economics, 8(2): , July doi: /qe737. Martin Flodén. A note on the accuracy of Markov-chain approximations to highly persistent AR(1) processes. Economics Letters, 99(3): , June doi: /j.econlet Xavier Gabaix. Variable rare disasters: An exactly solved framework for ten puzzles in macro-finance. Quarterly Journal of Economics, 127(2): , doi: /qje/qjs001. Gene H. Golub and John H. Welsch. Calculation of Gauss quadrature rules. Mathematics of Computation, 23(106): , may doi: /s Robert E. Lucas, Jr. Asset prices in an exchange economy. Econometrica, 46 (6): , November doi: / Allen C. Miller, III and Thomas R. Rice. Discrete approximations of probability distributions. Management Science, 29(3): , March doi: /mnsc Walter Pohl, Karl Schmedders, and Ole Wilms. Higher-order effects in asset pricing models with long-run risks. Journal of Finance, doi: /jofi Thomas A. Rietz. The equity risk premium: A solution. Journal of Monetary Economics, 22(1): , July doi: / (88) K. Geert Rouwenhorst. Asset pricing implications of equilibrium business cycle models. In Thomas F. Cooley, editor, Frontiers of Business Cycle Research, chapter 10, pages Princeton University Press, Bernard W. Silverman. Density Estimation for Statistics and Data Analysis, volume 26 of Monographs on Statistics and Applied Probability. Chapman & Hall/CRC, New York, Ken ichiro Tanaka and Alexis Akira Toda. Discrete approximations of continuous distributions by maximum entropy. Economics Letters, 118(3): , March doi: /j.econlet Ken ichiro Tanaka and Alexis Akira Toda. Discretizing distributions with exact moments: Error estimate and convergence analysis. SIAM Journal on Numerical Analysis, 53(5): , doi: /
10 George Tauchen. Finite state Markov-chain approximations to univariate and vector autoregressions. Economics Letters, 20(2): , doi: / (86) George Tauchen and Robert Hussey. Quadrature-based methods for obtaining approximate solutions to nonlinear asset pricing models. Econometrica, 59 (2): , March doi: / A Accuracy As in any numerical method, evaluating the accuracy is very important. In this section I evaluate the accuracy of the proposed method using the optimal portfolio problem in Section 3.2 as a laboratory. I design the numerical experiment as follows. First I fit a Gaussian mixture distribution with two components to the annual log excess returns data. The proportion, mean, and standard deviation of each mixture components are p = (p j ) = (0.1392,0.8608), µ = (µ j ) = ( ,0.1064), and σ = (σ j ) = (0.2164, ), respectively. Next, I assume that the true excess returns distribution is this Gaussian mixture and solve the optimal portfolio problem for relative risk aversion γ {2, 4, 6} using the Gaussian quadrature (Golub-Welsch algorithm 1) for Gaussian mixtures with 11 points. Finally, I generate random numbers from this Gaussian mixture with various sample sizes, discretize these distributions with various methods, and compute the optimal portfolio. I repeat this procedure with M = 1,000 Monte Carlo replications and compute the relative bias and mean absolute error (MAE) Bias = 1 M MAE = 1 M M m=1 ( θm /θ 1), (A.1a) M θ m /θ 1, (A.1b) m=1 where θ m is the optimal portfolio from simulation m and θ is the theoretical optimal portfolio. For the sample size I consider T = 100,1,000,10,000, and for the number of quadrature nodes I consider N = 3,5,7,9. For the discretization method I consider three cases. The first is the kernel density estimator-gaussian quadrature method (Algorithm 2), which I refer to as KDE-GQ. The second is the Gauss-Hermite quadrature, where the mean and standard deviation are estimated by maximum likelihood. This is the most natural method if the returns distribution is lognormal. The third is the maximum entropy method proposed by Tanaka and Toda (2013, 2015) and Farmer and Toda (2017) where the kernel density estimatoris fed into, which I referto as KDE-ME.For this method one needs to assign the grid and the number of moments to match. Following Corollary 3.5 of Farmer and Toda (2017), I use an even-spaced grid centered at the sample mean that spans 2(N 1) times the sample standard deviation at both sides, where N is the number of grid points. I match 4 moments whenever possible for N 5, and otherwise I match 2 moments (mean and variance). For more details on the exact algorithm, please refer to Tanaka and Toda (2013, 2015) and Sections 2 and 3.2 of Farmer and Toda (2017). 10
11 Tables 1 and 2 show the relative bias and mean absolute error of the optimal portfolio, respectively. As expected, the optimal portfolio computed using Gauss-Hermite is biased upwards because it uses the Gaussian distribution, which underestimates the probability of crashes. KDE-GQ and KDE-ME perform similarly for N 5 grid points, although KDE-GQ has slightly better small sample properties (T = 100). For N = 3 grid points, in which case it is impossible to match 4 moments with KDE-ME, the proposed KDE-GQ method performs better. Finally, increasing N beyond 5 does not improve the bias or the mean absolute error, which suggests that using a five-point distribution is enough (at least for solving this portfolio problem). Table 1: Relative bias of the optimal portfolio. Method KDE-GQ Gauss-Hermite KDE-ME T N γ = γ = γ = ,000 10, Note: the table reports the relative bias of the optimal portfolio defined by (A.1a). For discretization methods, KDE-GQ uses Algorithm 2, Gauss-Hermite uses the Gauss- Hermite quadrature (with mean and standard deviation estimated by maximum likelihood), and KDE-ME uses the maximum entropy method with the kernel density estimator. T is the sample size in each simulation. N is the number of nodes in the quadrature formula. γ is the relative risk aversion. All results are based on 1,000 Monte Carlo replications. B Gaussian quadrature (not for publication) In this appendix we prove some properties of the Gaussian quadrature. For notational simplicity let us omit a,b (so means b a ) and assume that w(x)x n dx exists for all n 0. For functions f,g, define the inner product (f,g) by (f,g) = b a w(x)f(x)g(x) dx. (B.1) As usual, define the norm of f by f = (f,f). The first step is to construct orthogonal polynomials {p n (x)} N n=0 corresponding to the inner product (B.1). Definition 1 (Orthogonalpolynomial). The polynomials {p n (x)} N n=0 are called orthogonal if (i) degp n = n and the leading coefficient of p n is 1, and (ii) for all m n, we have (p m,p n ) = 0. 11
12 Table 2: Relative mean absolute error of the optimal portfolio. Method KDE-GQ Gauss-Hermite KDE-ME T N γ = γ = γ = ,000 10, Note: the table reports the relative mean absolute error of the optimal portfolio defined by (A.1b). See Table 1 for the definition of variables. Some authors require that the polynomials are orthonormal, so (p n,p n ) = 1. In this paper we normalize the polynomials by requiring that the leading coefficient is 1, which is useful for computation. The following three-term recurrence relation (TTRR) shows the existence of orthogonal polynomials and provides an explicit algorithm for computing them. Proposition2(Three-termrecurrencerelation,TTRR). Let p 0 (x) = 1, p 1 (x) = x (xp0,p0), and for n 1 define p 0 2 p n+1 (x) = ( x (xp n,p n ) p n 2 ) p n (x) p n 2 Then p n (x) is the degree n orthogonal polynomial. p n 1 2p n 1(x). (B.2) Proof. Let us show by induction on n that (i) p n is an degree n polynomial with leading coefficient 1, and (ii) (p n,p m ) = 0 for all m < n. The claim is trivial for n = 0. For n = 1, by construction p 1 is a degree 1 polynomial with leading coefficient 1, and since p 0 (x) = 1, we obtain (p 1,p 0 ) = (( x (xp 0,p 0 ) p 0 2 )p 0,p 0 ) = (xp 0,p 0 ) (xp 0,p 0 ) = 0. Suppose the claim holds up to n. Then for n + 1, by (B.2) the leading coefficient of p n+1 is the same as that of xp n, which is 1. If m = n, then (( ) ) (p n+1,p n ) = x (xp n,p n ) p n 2 p n p n 2 p n 1 2p n 1,p n = (xp n,p n ) (xp n,p n ) p n 2 p n 1 2(p n 1,p n ) = 0. 12
13 If m = n 1, then (p n+1,p n 1 ) = (( x (xp n,p n ) p n 2 ) ) p n p n 2 p n 1 2p n 1,p n 1 = (xp n,p n 1 ) (xp n,p n ) p n 2 (p n,p n 1 ) p n 2 = (p n,xp n 1 ) p n 2. Since the leading coefficients of p n,p n 1 are 1, we can write xp n 1 (x) = p n (x)+ q(x), where q(x) is a polynomial of degree at most n 1. Clearly q can be expressed as a linear combination of p 0,p 1,...,p n 1, so (p n,q) = 0. Therefore (p n+1,p n 1 ) = (p n,p n +q) p n 2 = p n 2 +(p n,q) p n 2 = 0. Finally, if m < n 1, then (( (p n+1,p m ) = x (xp n,p n ) p n 2 ) ) p n p n 2 p n 1 2p n 1,p m = (xp n,p m ) (xp n,p n ) p n 2 (p n,p m ) p 2 n p n 1 2(p n 1,p m ) = (p n,xp m ) = 0 because xp m is a polynomial of degree 1+m < n. The following lemma shows that an degree n orthogonal polynomial has exactly n real roots (so they are all simple). Lemma 3. p n (x) has exactly n real roots on (a,b). Proof. By the fundamental theorem of algebra, p n (x) has exactly n roots in C. Suppose on the contrary that p n (x) has less than n real roots on (a,b). Let x 1,...,x k (k < n) those roots at which p n (x) changes its sign. Let q(x) = (x x 1 ) (x x k ). Since p n (x)q(x) > 0 (or < 0) almost everywhere on (a,b), we have (p n,q) = w(x)p n (x)q(x)dx 0. On the other hand, since degq = k < n, we have (p n,q) = 0, which is a contradiction. The following theorem shows that using the N roots of the degree N orthogonal polynomial p N (x) as quadrature nodes and choosing specific weights, we can integrate all polynomials of degree up to 2N 1 exactly. Thus Gaussian quadrature always exists. Theorem 4 (Gaussian quadrature). Let a < x 1 < < x N < b be the N roots of the degree N orthogonal polynomial p N and define w n = w(x)l n (x)dx 13
14 for n = 1,...,N, where L n (x) = m n x x m x n x m is the degree N 1 polynomial that takes value 1 at x n and 0 at x m (m {1,...,N}\n). Then w(x)p(x) dx = w n p(x n ) n=1 for all polynomials p(x) of degree up to 2N 1. Proof. Since degp 2N 1 and degp N = N, we can write p(x) = p N (x)q(x)+r(x), (B.3) where degq,degr N 1. Since q can be expressed as a linear combination of orthogonal polynomials of degree up to N 1, we have (p N,q) = 0. Hence w(x)p(x)dx = (p N,q)+ w(x)r(x) dx = w(x)r(x) dx. On the other hand, since {x n } N n=1 are roots of p N, we have for all n, so in particular p(x n ) = p N (x n )q(x n )+r(x n ) = r(x n ) w n p(x n ) = n=1 w n r(x n ). Therefore it suffices to show (B.3) for polynomials r of degree up to N 1. Let us show that r(x) = r(x n )L n (x) n=1 identically. To see this, let r be the right-hand side. Since L n (x m ) = δ mn (Kronecker s delta), we have r(x m ) = r(x n )L n (x m ) = n=1 n=1 δ mn r(x n ) = r(x m ), so r and r agree on N distinct points {x n } N n=1. Since each L n(x) is a degree N 1 polynomial, we have deg r N 1. Therefore it must be r = r. Since r can be represented as a linear combination of L n s, it suffices to show (B.3) for all L n s. But since by definition the claim is true. w(x)l n (x)dx = w n = n=1 w m δ mn = m=1 w m L n (x m ), m=1 14
15 In practice, how can we compute the nodes {x n } N n=1 and weights {w n} N n=1 of the N-point Gaussian quadrature? The solution is given by the following Golub-Welsch algorithm. Theorem 5 (Golub and Welsch, 1969). For each n 1, define α n,β n by α n = (xp n 1,p n 1 ) p n 1 2, β n = p n p n 1 > 0. Define the N N symmetric tridiagonal matrix T N as in (2.2). Then the Gaussian quadrature nodes {x n } N n=1 are eigenvalues of T N. Letting v n = (v n1,...,v nn ) be an eigenvector of T N corresponding to eigenvalue x n, then the weights {w n } N n=1 are given by w n = v2 n1 v n 2 w(x)dx > 0. (B.4) Proof. By (B.2) and the definition of α n,β n, for all n 0 we have p n+1 (x) = (x α n+1 )p n (x) β 2 np n 1 (x). Note that this is true for n = 0 by defining p 1 (x) = 0 and β 0 = 0. For each n, let p n(x) = p n (x)/ p n be the normalized orthogonal polynomial. Then the above equation becomes p n+1 p n+1 (x) = p n (x α n+1 )p n (x) p n 1 β 2 n p n 1 (x). Dividing both sides by p n > 0, using the definition of β n,β n+1, and rearranging terms, we obtain β n p n 1 (x)+α n+1p n (x)+β n+1p n+1 (x) = xp n (x). In particular, setting x = x k (where x k is a root of p N ), we obtain β n p n 1(x k )+α n+1 p n(x k )+β n+1 p n+1(x k ) = x k p n(x k ). for all n and k = 1,...,N. Since β 0 = 0 by definition and p N (x k) = 0 (since x k is a root of p N and hence p N = p N/ p N ), letting P(x) = (p 0(x),...,p N 1 (x)) and collecting the above equation into a vector, we obtain T N P(x k ) = x k P(x k ) for k = 1,...,N. Define the N N matrix P by P = (P(x 1 ),...,P(x N )). Then T N P = diag(x 1,...,x N )P, so x 1,...,x N are eigenvalues of T N provided that P is invertible. Now since {p n} N 1 n=0 are normalized and Gaussian quadrature integrates all polynomials of degree up to 2N 1 exactly, we have δ mn = (p m,p n ) = N w(x)p m (x)p n (x)dx = w k p m (x k)p n (x k) for m,n N 1. Letting W = diag(w 1,...,w N ), this equation becomes PWP = I. Therefore P,W are invertible and x 1,...,x N are eigenvalues of T N. Solving for W and taking the inverse, we obtain k=1 W 1 = P P 1 N 1 = p k w (x n) 2 > 0 n 15 k=0
16 for all n. To show (B.4), let v n be an eigenvector of T N corresponding to eigenvalue x n. Then v n = cp(x n ) for some constant c 0. Taking the norm, we obtain N 1 v n 2 = c 2 P(x n ) 2 = c 2 k=0 p k(x n ) 2 = c2 w n w n = c2 v n 2. Comparing the first element of v n = cp(x n ), noting that p 0 (x) = 1 and hence p 0 = p 0/ p 0 = 1/ p 0, we obtain c 2 = vn1 2 p 0 2 = vn1 2 w(x)p 0 (x) 2 dx = vn1 2 w(x) dx, which implies (B.4). 16
Higher-Order Dynamics in Asset-Pricing Models with Recursive Preferences
Higher-Order Dynamics in Asset-Pricing Models with Recursive Preferences Walt Pohl Karl Schmedders Ole Wilms Dept. of Business Administration, University of Zurich Becker Friedman Institute Computational
More informationMarkov-Chain Approximations for Life-Cycle Models
Markov-Chain Approximations for Life-Cycle Models Giulio Fella Giovanni Gallipoli Jutong Pan May 3, 2017 Abstract Non-stationary income processes are standard in quantitative life-cycle models, prompted
More informationUniversity of California Berkeley
Working Paper #2018-02 Infinite Horizon CCAPM with Stochastic Taxation and Monetary Policy Revised from the Center for Risk Management Research Working Paper 2018-01 Konstantin Magin, University of California,
More informationDiscretizing Markov Processes with Exact Conditional Moments
Discretizing Markov Processes with Exact Conditional Moments Leland E. Farmer Alexis Akira Toda This version: January 17, 2016 Abstract Approximating stochastic processes by finite-state Markov chains
More informationDiscrete State Space Methods for Dynamic Economies
Discrete State Space Methods for Dynamic Economies A Brief Introduction Craig Burnside Duke University September 2006 Craig Burnside (Duke University) Discrete State Space Methods September 2006 1 / 42
More informationIdentifying 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 informationHigh-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 informationOrthogonal Polynomials and Gaussian Quadrature
Orthogonal Polynomials and Gaussian Quadrature 1. Orthogonal polynomials Given a bounded, nonnegative, nondecreasing function w(x) on an interval, I of the real line, we consider the Hilbert space L 2
More informationThe Hansen Singleton analysis
The Hansen Singleton analysis November 15, 2018 1 Estimation of macroeconomic rational expectations model. The Hansen Singleton (1982) paper. We start by looking at the application of GMM that first showed
More informationEquilibrium Valuation with Growth Rate Uncertainty
Equilibrium Valuation with Growth Rate Uncertainty Lars Peter Hansen University of Chicago Toulouse p. 1/26 Portfolio dividends relative to consumption. 2 2.5 Port. 1 Port. 2 Port. 3 3 3.5 4 1950 1955
More informationUncertainty Per Krusell & D. Krueger Lecture Notes Chapter 6
1 Uncertainty Per Krusell & D. Krueger Lecture Notes Chapter 6 1 A Two-Period Example Suppose the economy lasts only two periods, t =0, 1. The uncertainty arises in the income (wage) of period 1. Not that
More informationLecture 3: Huggett s 1993 model. 1 Solving the savings problem for an individual consumer
UNIVERSITY OF WESTERN ONTARIO LONDON ONTARIO Paul Klein Office: SSC 4044 Extension: 85484 Email: pklein2@uwo.ca URL: http://paulklein.ca/newsite/teaching/619.php Economics 9619 Computational methods in
More informationLecture 6: Recursive Preferences
Lecture 6: Recursive Preferences Simon Gilchrist Boston Univerity and NBER EC 745 Fall, 2013 Basics Epstein and Zin (1989 JPE, 1991 Ecta) following work by Kreps and Porteus introduced a class of preferences
More informationDynamic Programming. Macro 3: Lecture 2. Mark Huggett 2. 2 Georgetown. September, 2016
Macro 3: Lecture 2 Mark Huggett 2 2 Georgetown September, 2016 Three Maps: Review of Lecture 1 X = R 1 + and X grid = {x 1,..., x n } X where x i+1 > x i 1. T (v)(x) = max u(f (x) y) + βv(y) s.t. y Γ 1
More informationMultivariate 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 informationAsset pricing in DSGE models comparison of different approximation methods
Asset pricing in DSGE models comparison of different approximation methods 1 Introduction Jan Acedański 1 Abstract. There are many numerical methods suitable for approximating solutions of DSGE models.
More information1. Using the model and notations covered in class, the expected returns are:
Econ 510a second half Yale University Fall 2006 Prof. Tony Smith HOMEWORK #5 This homework assignment is due at 5PM on Friday, December 8 in Marnix Amand s mailbox. Solution 1. a In the Mehra-Prescott
More informationA NEW NONLINEAR FILTER
COMMUNICATIONS IN INFORMATION AND SYSTEMS c 006 International Press Vol 6, No 3, pp 03-0, 006 004 A NEW NONLINEAR FILTER ROBERT J ELLIOTT AND SIMON HAYKIN Abstract A discrete time filter is constructed
More informationSpeculation and the Bond Market: An Empirical No-arbitrage Framework
Online Appendix to the paper Speculation and the Bond Market: An Empirical No-arbitrage Framework October 5, 2015 Part I: Maturity specific shocks in affine and equilibrium models This Appendix present
More informationSolving non-linear systems of equations
Solving non-linear systems of equations Felix Kubler 1 1 DBF, University of Zurich and Swiss Finance Institute October 7, 2017 Felix Kubler Comp.Econ. Gerzensee, Ch2 October 7, 2017 1 / 38 The problem
More informationDiscretizing Stochastic Processes with Exact Conditional Moments
Discretizing Stochastic Processes with Exact Conditional Moments Leland E. Farmer Alexis Akira Toda This version: August 25, 2015 Abstract Approximating stochastic processes by finite-state Markov chains
More informationA Robust Approach to Estimating Production Functions: Replication of the ACF procedure
A Robust Approach to Estimating Production Functions: Replication of the ACF procedure Kyoo il Kim Michigan State University Yao Luo University of Toronto Yingjun Su IESR, Jinan University August 2018
More informationAssessing Markov Chain Approximations: A Minimal Econometric Approach
Assessing Markov Chain Approximations: A Minimal Econometric Approach Pedro Silos Dept. of Economics University of Iowa October 18, 2004 Abstract Markov chain approximations of continuous state-space processes
More informationDo 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 informationSources of Entropy in Representative Agent Models
Sources of Entropy in Representative Agent Models David Backus (NYU), Mikhail Chernov (LSE), and Stanley Zin (NYU) Conference in Honor of T. Sargent and C. Sims May 2012 This version: May 5, 2012 1 / 21
More informationA nonparametric method of multi-step ahead forecasting in diffusion processes
A nonparametric method of multi-step ahead forecasting in diffusion processes Mariko Yamamura a, Isao Shoji b a School of Pharmacy, Kitasato University, Minato-ku, Tokyo, 108-8641, Japan. b Graduate School
More informationTeruko Takada Department of Economics, University of Illinois. Abstract
Nonparametric density estimation: A comparative study Teruko Takada Department of Economics, University of Illinois Abstract Motivated by finance applications, the objective of this paper is to assess
More informationDepartment 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 informationUniversity 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 informationOptimal Investment Strategies: A Constrained Optimization Approach
Optimal Investment Strategies: A Constrained Optimization Approach Janet L Waldrop Mississippi State University jlc3@ramsstateedu Faculty Advisor: Michael Pearson Pearson@mathmsstateedu Contents Introduction
More informationDiscrete State-Space Methods for the Study of Dynamic Economies
Discrete State-Space Methods for the Study of Dynamic Economies Craig Burnside Duke University September 25, 2006 Updated Version 1 Introduction In this chapter I discuss a number of numerical methods
More informationComprehensive Exam. Macro Spring 2014 Retake. August 22, 2014
Comprehensive Exam Macro Spring 2014 Retake August 22, 2014 You have a total of 180 minutes to complete the exam. If a question seems ambiguous, state why, sharpen it up and answer the sharpened-up question.
More informationIntroduction to Algorithmic Trading Strategies Lecture 10
Introduction to Algorithmic Trading Strategies Lecture 10 Risk Management Haksun Li haksun.li@numericalmethod.com www.numericalmethod.com Outline Value at Risk (VaR) Extreme Value Theory (EVT) References
More informationFinite State Markov-chain Approximations to Highly. Persistent Processes
Finite State Markov-chain Approximations to Highly Persistent Processes Karen A. Kopecky y Richard M. H. Suen z This Version: November 2009 Abstract The Rouwenhorst method of approximating stationary AR(1)
More informationWeek 1 Quantitative Analysis of Financial Markets Distributions A
Week 1 Quantitative Analysis of Financial Markets Distributions A Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036 October
More informationClass notes: Approximation
Class notes: Approximation Introduction Vector spaces, linear independence, subspace The goal of Numerical Analysis is to compute approximations We want to approximate eg numbers in R or C vectors in R
More informationAsset Pricing. Question: What is the equilibrium price of a stock? Defn: a stock is a claim to a future stream of dividends. # X E β t U(c t ) t=0
Asset Pricing 1 Lucas (1978, Econometrica) Question: What is the equilibrium price of a stock? Defn: a stock is a claim to a future stream of dividends. 1.1 Environment Tastes: " # X E β t U( ) t=0 Technology:
More informationA Test of Cointegration Rank Based Title Component Analysis.
A Test of Cointegration Rank Based Title Component Analysis Author(s) Chigira, Hiroaki Citation Issue 2006-01 Date Type Technical Report Text Version publisher URL http://hdl.handle.net/10086/13683 Right
More informationArea I: Contract Theory Question (Econ 206)
Theory Field Exam Summer 2011 Instructions You must complete two of the four areas (the areas being (I) contract theory, (II) game theory A, (III) game theory B, and (IV) psychology & economics). Be sure
More informationIntroduction to Orthogonal Polynomials: Definition and basic properties
Introduction to Orthogonal Polynomials: Definition and basic properties Prof. Dr. Mama Foupouagnigni African Institute for Mathematical Sciences, Limbe, Cameroon and Department of Mathematics, Higher Teachers
More informationMonte Carlo Methods. Leon Gu CSD, CMU
Monte Carlo Methods Leon Gu CSD, CMU Approximate Inference EM: y-observed variables; x-hidden variables; θ-parameters; E-step: q(x) = p(x y, θ t 1 ) M-step: θ t = arg max E q(x) [log p(y, x θ)] θ Monte
More informationMultivariate Distributions
IEOR E4602: Quantitative Risk Management Spring 2016 c 2016 by Martin Haugh Multivariate Distributions We will study multivariate distributions in these notes, focusing 1 in particular on multivariate
More informationThe Discretization Filter: A Simple Way to Estimate Nonlinear State Space Models
The Discretization Filter: A Simple Way to Estimate Nonlinear State Space Models Leland E. Farmer Department of Economics University of Virginia This Version: August 7, 217 First Version: August 214 Link
More informationIn manycomputationaleconomicapplications, one must compute thede nite n
Chapter 6 Numerical Integration In manycomputationaleconomicapplications, one must compute thede nite n integral of a real-valued function f de ned on some interval I of
More informationMethod of Finite Elements I
Method of Finite Elements I PhD Candidate - Charilaos Mylonas HIL H33.1 and Boundary Conditions, 26 March, 2018 Institute of Structural Engineering Method of Finite Elements I 1 Outline 1 2 Penalty method
More informationErrata for Campbell, Financial Decisions and Markets, 01/02/2019.
Errata for Campbell, Financial Decisions and Markets, 01/02/2019. Page xi, section title for Section 11.4.3 should be Endogenous Margin Requirements. Page 20, equation 1.49), expectations operator E should
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More informationPolitical Cycles and Stock Returns. Pietro Veronesi
Political Cycles and Stock Returns Ľuboš Pástor and Pietro Veronesi University of Chicago, National Bank of Slovakia, NBER, CEPR University of Chicago, NBER, CEPR Average Excess Stock Market Returns 30
More information1 The Basic RBC Model
IHS 2016, Macroeconomics III Michael Reiter Ch. 1: Notes on RBC Model 1 1 The Basic RBC Model 1.1 Description of Model Variables y z k L c I w r output level of technology (exogenous) capital at end of
More informationJoint 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 informationDiscretizing Stochastic Processes with Exact Conditional Moments
Discretizing Stochastic Processes with Exact Conditional Moments Leland E. Farmer Alexis Akira Toda This version: March 25, 2015 Abstract Approximating stochastic processes by finite-state Markov chains
More information1 Description of variables
1 Description of variables We have three possible instruments/state variables: dividend yield d t+1, default spread y t+1, and realized market volatility v t+1 d t is the continuously compounded 12 month
More informationNonparametric Methods
Nonparametric Methods Michael R. Roberts Department of Finance The Wharton School University of Pennsylvania July 28, 2009 Michael R. Roberts Nonparametric Methods 1/42 Overview Great for data analysis
More informationMacroeconomic Topics Homework 1
March 25, 2004 Kjetil Storesletten. Macroeconomic Topics Homework 1 Due: April 23 1 Theory 1.1 Aggregation Consider an economy consisting of a continuum of agents of measure 1 who solve max P t=0 βt c
More informationECONOMICS 7200 MODERN TIME SERIES ANALYSIS Econometric Theory and Applications
ECONOMICS 7200 MODERN TIME SERIES ANALYSIS Econometric Theory and Applications Yongmiao Hong Department of Economics & Department of Statistical Sciences Cornell University Spring 2019 Time and uncertainty
More informationPrice and Wealth Dynamics in a Speculative Market with an Arbitrary Number of Generic Technical Traders
Price and Wealth Dynamics in a Speculative Market with an Arbitrary Number of Generic Technical Traders Mikhail Anufriev Giulio Bottazzi S.Anna School for Advanced Studies, Pisa, Italy February 20, 2005
More informationAmbiguous Business Cycles: Online Appendix
Ambiguous Business Cycles: Online Appendix By Cosmin Ilut and Martin Schneider This paper studies a New Keynesian business cycle model with agents who are averse to ambiguity (Knightian uncertainty). Shocks
More informationNumerical integration and importance sampling
2 Numerical integration and importance sampling 2.1 Quadrature Consider the numerical evaluation of the integral I(a, b) = b a dx f(x) Rectangle rule: on small interval, construct interpolating function
More informationLarge-scale eigenvalue problems
ELE 538B: Mathematics of High-Dimensional Data Large-scale eigenvalue problems Yuxin Chen Princeton University, Fall 208 Outline Power method Lanczos algorithm Eigenvalue problems 4-2 Eigendecomposition
More informationComputing Gauss-Jacobi quadrature rules
.... Computing Gauss-Jacobi quadrature rules Alex Townsend Joint work with Nick Hale NA Internal Seminar Trinity Term 2012 . Gauss Jacobi Quadrature An n-point Gauss Jacobi quadrature rule: 1 1 w(x)f (x)dx
More informationDynamic Risk Measures and Nonlinear Expectations with Markov Chain noise
Dynamic Risk Measures and Nonlinear Expectations with Markov Chain noise Robert J. Elliott 1 Samuel N. Cohen 2 1 Department of Commerce, University of South Australia 2 Mathematical Insitute, University
More informationGeneralized Impulse Response Analysis: General or Extreme?
MPRA Munich Personal RePEc Archive Generalized Impulse Response Analysis: General or Extreme? Kim Hyeongwoo Auburn University April 2009 Online at http://mpra.ub.uni-muenchen.de/17014/ MPRA Paper No. 17014,
More informationConcavity of the Consumption Function with Recursive Preferences
Concavity of the Consumption Function with Recursive Preferences Semyon Malamud May 21, 2014 Abstract Carroll and Kimball 1996) show that the consumption function for an agent with time-separable, isoelastic
More informationProblem 1A. Find the volume of the solid given by x 2 + z 2 1, y 2 + z 2 1. (Hint: 1. Solution: The volume is 1. Problem 2A.
Problem 1A Find the volume of the solid given by x 2 + z 2 1, y 2 + z 2 1 (Hint: 1 1 (something)dz) Solution: The volume is 1 1 4xydz where x = y = 1 z 2 This integral has value 16/3 Problem 2A Let f(x)
More informationA Novel Nonparametric Density Estimator
A Novel Nonparametric Density Estimator Z. I. Botev The University of Queensland Australia Abstract We present a novel nonparametric density estimator and a new data-driven bandwidth selection method with
More informationNonparametric Econometrics
Applied Microeconometrics with Stata Nonparametric Econometrics Spring Term 2011 1 / 37 Contents Introduction The histogram estimator The kernel density estimator Nonparametric regression estimators Semi-
More informationThe Laplacian PDF Distance: A Cost Function for Clustering in a Kernel Feature Space
The Laplacian PDF Distance: A Cost Function for Clustering in a Kernel Feature Space Robert Jenssen, Deniz Erdogmus 2, Jose Principe 2, Torbjørn Eltoft Department of Physics, University of Tromsø, Norway
More informationNonparametric Identification of a Binary Random Factor in Cross Section Data - Supplemental Appendix
Nonparametric Identification of a Binary Random Factor in Cross Section Data - Supplemental Appendix Yingying Dong and Arthur Lewbel California State University Fullerton and Boston College July 2010 Abstract
More informationON MEASURING HEDGE FUND RISK
ON MEASURING HEDGE FUND RISK Alexander Cherny Raphael Douady Stanislav Molchanov Department of Probability Theory Faculty of Mathematics and Mechanics Moscow State University 119992 Moscow Russia E-mail:
More informationBehavioral Economics. Lecture 10
14.127 Behavioral Economics. Lecture 10 Xavier Gabaix April 15, 2004 1 Hyperbolic discounting Luttmer and Mariotti (JPE 2003) hyperbolics does not make much difference/improvement over exponential discounting.
More informationProblem Set 4. Graduate Macro II, Spring 2011 The University of Notre Dame Professor Sims
Problem Set 4 Graduate Macro II, Spring 2011 The University of Notre Dame Professor Sims Instructions: You may consult with other members of the class, but please make sure to turn in your own work. Where
More informationMonte Carlo Solution of Integral Equations
1 Monte Carlo Solution of Integral Equations of the second kind Adam M. Johansen a.m.johansen@warwick.ac.uk www2.warwick.ac.uk/fac/sci/statistics/staff/academic/johansen/talks/ March 7th, 2011 MIRaW Monte
More informationChapter 6. Maximum Likelihood Analysis of Dynamic Stochastic General Equilibrium (DSGE) Models
Chapter 6. Maximum Likelihood Analysis of Dynamic Stochastic General Equilibrium (DSGE) Models Fall 22 Contents Introduction 2. An illustrative example........................... 2.2 Discussion...................................
More informationMonte Carlo Integration I [RC] Chapter 3
Aula 3. Monte Carlo Integration I 0 Monte Carlo Integration I [RC] Chapter 3 Anatoli Iambartsev IME-USP Aula 3. Monte Carlo Integration I 1 There is no exact definition of the Monte Carlo methods. In the
More informationDynamic System Identification using HDMR-Bayesian Technique
Dynamic System Identification using HDMR-Bayesian Technique *Shereena O A 1) and Dr. B N Rao 2) 1), 2) Department of Civil Engineering, IIT Madras, Chennai 600036, Tamil Nadu, India 1) ce14d020@smail.iitm.ac.in
More information1 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 informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond January 18, 2017 Contents 1 Batch and Recursive Estimation 2 Towards Bayesian Filtering 3 Kalman Filter and Bayesian Filtering and Smoothing
More informationAPPENDIX B GRAM-SCHMIDT PROCEDURE OF ORTHOGONALIZATION. Let V be a finite dimensional inner product space spanned by basis vector functions
301 APPENDIX B GRAM-SCHMIDT PROCEDURE OF ORTHOGONALIZATION Let V be a finite dimensional inner product space spanned by basis vector functions {w 1, w 2,, w n }. According to the Gram-Schmidt Process an
More informationIndeterminacy and Sunspots in Macroeconomics
Indeterminacy and Sunspots in Macroeconomics Friday September 8 th : Lecture 10 Gerzensee, September 2017 Roger E. A. Farmer Warwick University and NIESR Topics for Lecture 10 Tying together the pieces
More informationMA Advanced Macroeconomics: Solving Models with Rational Expectations
MA Advanced Macroeconomics: Solving Models with Rational Expectations Karl Whelan School of Economics, UCD February 6, 2009 Karl Whelan (UCD) Models with Rational Expectations February 6, 2009 1 / 32 Moving
More informationCurrency Risk Factors in a Recursive Multi-Country Economy
Currency Risk Factors in a Recursive Multi-Country Economy Ric Colacito Max Croce F. Gavazzoni Rob Ready 1 / 28 Motivation The literature has identified factor structures in currency returns Interest Rates
More informationAre Probabilities Used in Markets? 1
Journal of Economic Theory 91, 8690 (2000) doi:10.1006jeth.1999.2590, available online at http:www.idealibrary.com on NOTES, COMMENTS, AND LETTERS TO THE EDITOR Are Probabilities Used in Markets? 1 Larry
More informationOn orthogonal polynomials for certain non-definite linear functionals
On orthogonal polynomials for certain non-definite linear functionals Sven Ehrich a a GSF Research Center, Institute of Biomathematics and Biometry, Ingolstädter Landstr. 1, 85764 Neuherberg, Germany Abstract
More informationStochastic processes. MAS275 Probability Modelling. Introduction and Markov chains. Continuous time. Markov property
Chapter 1: and Markov chains Stochastic processes We study stochastic processes, which are families of random variables describing the evolution of a quantity with time. In some situations, we can treat
More informationVolume 30, Issue 1. Measuring the Intertemporal Elasticity of Substitution for Consumption: Some Evidence from Japan
Volume 30, Issue 1 Measuring the Intertemporal Elasticity of Substitution for Consumption: Some Evidence from Japan Akihiko Noda Graduate School of Business and Commerce, Keio University Shunsuke Sugiyama
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationMultivariate Normal-Laplace Distribution and Processes
CHAPTER 4 Multivariate Normal-Laplace Distribution and Processes The normal-laplace distribution, which results from the convolution of independent normal and Laplace random variables is introduced by
More informationWealth, Information Acquisition and Portfolio Choice: A Correction
Wealth, Information Acquisition and Portfolio Choice: A Correction Joel Peress INSEAD There is an error in our 2004 paper Wealth, Information Acquisition and Portfolio Choice. This note shows how to correct
More informationA note on multivariate Gauss-Hermite quadrature
A note on multivariate Gauss-Hermite quadrature Peter Jäckel 1th May 5 1 Introduction Gaussian quadratures are an ingenious way to approximate the integral of an unknown function f(x) over a specified
More information1: PROBABILITY REVIEW
1: PROBABILITY REVIEW Marek Rutkowski School of Mathematics and Statistics University of Sydney Semester 2, 2016 M. Rutkowski (USydney) Slides 1: Probability Review 1 / 56 Outline We will review the following
More informationDiscussion Papers in Economics
Discussion Papers in Economics No. 10/11 A General Equilibrium Corporate Finance Theorem for Incomplete Markets: A Special Case By Pascal Stiefenhofer, University of York Department of Economics and Related
More informationInflation 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 informationOrthogonal Polynomials, Quadratures & Sparse-Grid Methods for Probability Integrals
1/31 Orthogonal Polynomials, Quadratures & Sparse-Grid Methods for Probability Integrals Dr. Abebe Geletu May, 2010 Technische Universität Ilmenau, Institut für Automatisierungs- und Systemtechnik Fachgebiet
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationIdentifying 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 informationAdvanced SDPs Lecture 6: March 16, 2017
Advanced SDPs Lecture 6: March 16, 2017 Lecturers: Nikhil Bansal and Daniel Dadush Scribe: Daniel Dadush 6.1 Notation Let N = {0, 1,... } denote the set of non-negative integers. For α N n, define the
More informationSpatially Smoothed Kernel Density Estimation via Generalized Empirical Likelihood
Spatially Smoothed Kernel Density Estimation via Generalized Empirical Likelihood Kuangyu Wen & Ximing Wu Texas A&M University Info-Metrics Institute Conference: Recent Innovations in Info-Metrics October
More informationA Polynomial Chaos Approach to Robust Multiobjective Optimization
A Polynomial Chaos Approach to Robust Multiobjective Optimization Silvia Poles 1, Alberto Lovison 2 1 EnginSoft S.p.A., Optimization Consulting Via Giambellino, 7 35129 Padova, Italy s.poles@enginsoft.it
More informationA Moment-Matching Method for Approximating Vector Autoregressive Processes by Finite-State Markov Chains
A Moment-Matching Method for Approximating Vector Autoregressive Processes by Finite-State Markov Chains Nikolay Gospodinov Concordia University and CIREQ Damba Lkhagvasuren Concordia University and CIREQ
More information