Stochastic Systems: Modeling, Analysis, Synthesis, Control, and their Applications to Engineering
View this Special IssueResearch Article  Open Access
Sumei Zhang, Lihe Wang, "A Fast Fourier Transform Technique for Pricing European Options with Stochastic Volatility and Jump Risk", Mathematical Problems in Engineering, vol. 2012, Article ID 761637, 17 pages, 2012. https://doi.org/10.1155/2012/761637
A Fast Fourier Transform Technique for Pricing European Options with Stochastic Volatility and Jump Risk
Abstract
We consider European options pricing with double jumps and stochastic volatility. We derived closedform solutions for European call options in a double exponential jumpdiffusion model with stochastic volatility (SVDEJD). We developed fast and accurate numerical solutions by using fast Fourier transform (FFT) technique. We compared the density of our model with those of other models, including the BlackScholes model and the double exponential jumpdiffusion model. At last, we analyzed several effects on option prices under the proposed model. Simulations show that the SVDEJD model is suitable for modelling the longtime realmarket changes and stock returns are negatively correlated with volatility. The model and the proposed option pricing method are useful for empirical analysis of asset returns and managing the corporate credit risks.
1. Introduction
The classical BlackScholes (BS) model [1] has long been known to result in systematically biased option valuation. By adding jumps to the archetypal price process with Gaussian innovations Merton [2] is able to partly explain the observed deviations from the benchmark model which are characterized by fat tail and excess kurtosis in the returns distribution. For an overview of “stylized facts” on asset returns see Cont [3]. Statistical properties of implied volatilities are summarized in Cont et al. [4]. In the sequel also other authors develop more realistic models, for example, the pure jump models of Eberlein and Keller [5], Madan et al. [6], and Duffie et al. [7], stochastic volatility models of Steven [8], and stochastic volatility model with normal jumps of Bates [9] and Keppo et al. [10]. The double exponential jumpdiffusion (DEJD) model, recently proposed by Kou [11], generates a highly skewed and leptokurtic distribution and is capable of matching key features of stock and index returns. Moreover, the DEJD model leads to tractable pricing formulas for exotic and path dependent options [12]. Accordingly, the DEJD model has gained wide acceptance. However, the DEJD model cannot capture the volatility clustering effects, which can be captured by stochastic volatility models [13]. Jumpdiffusion models and the stochastic volatility model complement each other: the stochastic volatility model can incorporate dependent structures better, while the DEJD model has better analytical tractability, especially for pathdependent options. Since allowing interest rates to be stochastic does not improve pricing performance any further [14], the model that combines stochastic volatility and double exponential jumpdiffusion (SVDEJD) may be more reasonable.
In the BS setting, the probability measure has a wellknown analytic form [15], but, under stochastic volatility, it can only be obtained numerically. Monte Carlo simulation and the finite difference method are usually used to value the options. But, the two techniques require substantially more computing time and thus are difficult to be applied in real option pricing. Recently, being fast, accurate, and easy to implement, Fourier transforms have been widely used in valuing financial derivatives, for example, Carr and Madan [16] propose Fourier transforms with respect to logstrike price; Duffie et al. [7] offer a comprehensive survey that the Fourier methods are applicable to a wide range of stochastic processes; Carr and Wu [17] apply the transforms to timechanged Lévy processes and the class of generalized affine models. Hurd and Zhou [18] express the spread option payoff in terms of the gamma function and FFT technique. For an overview of option pricing using Fourier transforms, see Schmelzle [19].
The current paper extends the study of option pricing under the DEJD model in three ways. First, we propose a model which combines the double jumps and stochastic volatility. Second, using the martingale method, Fourier transform formula, and FeynmanKac theorem, we obtain a closedform solution for European call options pricing under the proposed model. Third, we obtain fast and accurate numerical solutions for European call options pricing by FFT technique.
The rest of the paper is organized as follows. Section 2 develops the underlying pricing model. Section 3 derives a closedform solution for European call options pricing under the proposed model. Section 4 provides approximation solutions for European call options pricing by FFT technique. Section 5 numerically compares the density of the solutions to the alternative models and analyzes several effects on potion prices. Section 6 concludes. Applied program codes in Matlab package are presented in the appendix.
2. The Model
We consider an arbitragefree, frictionless financial market where only riskless asset and risky asset are traded continuously up to a fixed horizon date . Let be a complete probability space with a filtration satisfying the usual conditions, that is, the filtration is continuous on the right and contains all null sets. Suppose , are both standard Brownian motion, which is adapted, and has correlation with .
Let represent the price for a stock or a stock portfolio. Generally, instantaneous variance of asset returns in financial markets shows randomness; thus, the continuous part of the price process, defined as , is where is riskfree rate and is nonnegative constant, and suppose , which can be set equal to 1 without any loss of generality. The size of the diffusion component is determined by , which represents, absent of any jump occurring, the level of (stochastic) return variance attributable to diffusion variations. For tractability, let follow a squareroot process: where nonnegative constants , and , respectively, reflect the speed of adjustment, the longrun mean, and the variation coefficient of , and suppose .
It has been suggested from extensive empirical studies that markets tend to have both overreaction and underreaction to various good or bad news. One may interpret the jump part of the model as the market response to outside news. Good or bad news arrives according to a Poisson process, and the asset price changes in response according to the jump size distribution. According to Kou [11], the jumps in the logprice are modeled as a sequence of i.i.d. nonnegative random variables that occur at times determined by an independent Poisson process with constant intensity such that has an asymmetric double exponential distribution with the density where denotes the indicator function, so equals 1 if , but 0 otherwise. are upmove jump and downmove jump, respectively. Except for which has correlation with , all sources of randomness, , and , are assumed to be independent.
Because of jumps and stochastic volatility, the riskneutral probability measure is not unique. Following Naik and Lee [20] and Kou [11], by using the rational expectations argument with a HARAtype utility function for the representative agent, one can choose a particular riskneutral measure so that the equilibrium price of an option is given by the expectation under this riskneutral measure of the discounted option payoff. Throughout this paper, we assume that there exists a martingale probability measure being equivalent to . Let be the sum of all the jumps which occur up to and including time , , we have Obviously, is a martingale. Finally, the price process is defined as
Remark 2.1. The model contains most existing models as special cases. For example, we obtain (1) the BS model by setting and ; (2) the SV model by setting ; (3) the DEJD model by setting .
Let , where is standard Brownian motion that is adapted, and independent of , , and random variables . From It’s formula, we have
3. A ClosedForm Solution of European Option Pricing
In this section, we derive closedform solution of a European call option pricing under the SVDEJD model. For a European put option, we can obtain easily corresponding result by the putcall parity [1]. For this purpose, we need the following results.
Lemma 3.1. Supposing the variance process follows (2.2) and are any complex, one has where
Proof. Let . Because of the affine structure of the variance process (2.2), we obtain that has a solution of the following form: From the FeynmanKac formula, is the solution of the following backward Parabolic partial differential equation with the Cauchy problem: Putting (3.3) in (3.4), we have Solving (3.5), we can obtain the result of Lemma 3.1.
Lemma 3.2. Supposing the asset price follows (2.6) and is any complex, one has where
Proof. Let . Because is independent of , and, we have From (2.2) and (2.3), we have Because is a standard Brownian motion, we have Then, Let , and . From Lemma 3.1, we have From (3.8), (3.9), (3.10) and (3.12), we can obtain the required Lemma 3.2.
Lemma 3.3. Suppose is the characteristic function of ; then where
Proof. Let . Because from Lemma 3.2, we can obtain the required Lemma 3.3.
Theorem 3.4. Let denote the log of the strike price , , and the desired value of a maturity call option with strike . Assume that, under , the underlying nondividendpaying stock price and its components are given by (2.1)–(2.5), is the characteristic function of , is the density of ; then the initial call value is written as where represents real part.
Proof. From the riskneutral theory, we have Introducing a change of measure from to by a RadonNikodym derivative, we get With this new measure , the Fourier transform of is defined as Because of the noarbitrage condition, we can obtain From the Fourier transform formula, the probability density for our model is given by Hence, Changing the order of integration, we have From (3.17), (3.20), and (3.23), we can obtain the required Theorem 3.4.
Remark 3.5. In (3.16), tends to not zero as goes to . Hence, is not (absolutely integrable) and a Fourier transform does not exist.
4. Fast Fourier Transform for European Option Pricing
Since the integrand in (3.16) is singular at the required evaluation point , the FFT cannot be applied directly to evaluate the integrals we mentioned above. Therefore, instead of solving for the riskneutral exercise probabilities of finishing intheMoney (ITM), Carr and Madan [16] introduce a new technique with the key idea to calculate the Fourier transform of a modified call option price with respect to the logarithmic strike price. With this specification and a FFT routine, a whole range of option prices can be obtained within a single Fourier inversion. In this section, we develop the numerical solutions of the prices by using the idea of Carr and Madan [16].
4.1. Fourier Transform of ITM and attheMoney (ATM) Option Prices
By introducing an exponential damping factor with , it is possible to make the integrand in (3.16) be square integrable. We modified the pricing function (3.16) by where .
This method is viable when is chosen in a way that the damped option price is well behaved. Damping the option price with makes it integrable for the negative axis . On the other hand, for the option prices increase by the exponential , which influences the integrability for the positive axis. A sufficient condition of to be integrable for both sides (square integrability) is given by being finite, that is, Thus we need , which is equivalent to Therefore, is well behaved when the moments of order of the underlying asset exist and are finite. If not all moments of exist, this will impose an upper bound on . We find that one quarter of this upper bound serves as a good choice for .
Using the trapezoid rule for the integral on the righthand side of (4.1) and setting , an approximation for is The FFT returns values of , and we employ a regular spacing of size so that our values for are This gives us logstrike levels ranging from to , where
In order to apply FFT we define
To obtain an accurate integration with larger values of , we incorporate Simpson’s rule weightings into our summation. From (4.1)–(4.7)and Simpson’s rule weightings, we obtain ATM and ITM call value as where is the Kronecker delta function that is unity for and zero otherwise. The summation in (4.8) is an exact application of the FFT.
4.2. Fourier Transform of outoftheMoney (OTM) Option Prices
In the previous section call values are calculated by an exponential function to obtain square integrable function whose Fourier transform is an analytic function of the characteristic function of the logprice. But, for very short maturities, the call value approaches its non analytic intrinsic value causing the integrand in the Fourier inversion to be high oscillate, and therefore difficult to integrate numerically. We introduce an alternative approach that works with time values only, which is quite similar to the previous approach. But in this case the call price is obtained via the Fourier transform of a modified time value, where the modification involves a hyperbolic sine function instead of an exponential function.
Let denote the time value of an OTM option, that is, for we have the put price for and for we have the call price. Scaling for simplicity, is defined by where is the riskneutral density of the logprice . Let be the Fourier transform of : By considering a damping function , the time value of an option follows a Fourier inversion: where .
The use of the FFT for calculating OTM option prices is similar to (4.8). The only differences are that they replace the multiplication by with a division by and the function call to is replaced by a function call to .
5. Simulation Studies
In this section, to compare across the BS, DEJD, and SVDEJD models, we analyze the probability densities of these models. Then, we analyze mainly the impact of and volatility of volatility on option pricing under the SVDEJD model. For our FFT methods, we used points in our quadrature, implying a logstrike spacing of , which is adequate for practice. For the choice of the dampening coefficient in the transform of the modified call price, we used a value of . For the modified time value, we used . Other parameter values used in the computation are listed in Table 1. (We have used analytic moments to set plausible parameter values for the model. For a formal econometric estimator, one could use these moments to develop a generalized method of moments estimator within the framework of Hall and Inoue [21].)

5.1. Probability Densities under Alternative Models
We compare the probability densities of the SVDEJD model, the BS model, and the DEJD model to verify the rationality of our model. Suppose is the characteristic function of and the probability density of our model. From FFT algorithm, can be approximated by The density has the mean and variance given by
Figure 1 shows the figures of the probability density , compared with the normal density with the same mean and variance given by (5.2). The first figure compares the overall shapes of the densities of the SVDEJD model and the BS model, the second one details the shapes around the peak areas, and the last one shows the right tail. From Figure 1, we can see that the leptokurtic and skewness feature of the density of our model is quite evident. Moreover, additional numerical plots suggest that the feature of skewness becomes more significant if increases, which is impossible for the DEJD model.
We also compare the shortterm and longterm densities of the SVDEJD model, the BS model, and the DEJD one. Figure 2 shows their densities under months and years. From Figure 2, we can see that the SVDEJD model and the DEJD model generate virtually identical densities for shortterm options, with a slight departure occurring between the two densities in the upper tail. This means that differential pricing performance between the SVDEJD model and the DEJD model is unlikely to occur when they are applied to price shortterm OTM puts and that only when they are applied to deep ITM puts (and deep OTM calls) can differences be observed between these models. Yet, compared to the BS model density, the densities of the two models are distinctly different: they all have leptokurtic and skewness feature. Therefore, the two models can potentially correct the BS model’s tendency to underprice deep OTM puts and overprice deep OTM calls. The longterm density curves in Figure 2 still show significantly different pricing structures between the BS and its two alternatives. But, more importantly, the densities of the SVDEJD model and the DEJD model also exhibit different shapes now. The SVDEJD density has higher peak and assigns more weight to both the entire lower tail and the far upper tail, but less weight to those payoffs than the DEJD.
Our simulation studies have demonstrated that the SVDEJD model has better performance than the DEJD one on pricing longterm options, while both the DEJD model and the SVDEJD model have better performance than the BS model.
5.2. Effects of the Main Parameter on Option Values
In Table 2, we use the SVDEJD model to examine the effects of volatility of volatility , the correlation coefficient , exercise price , and maturity time on option values. We analyze the prices of threemonth call options and twoyear call options. To examine the effect of the negative correlation coefficient, we have calculated the model with , , and . The prices for threemonth call options associated with volatility of volatility and are relatively close. With , the largest price difference is 0.0149; with , the largest price difference is only 0.0008; with , the largest price difference is 0.0213. However, the difference is significantly larger when longertime horizons such as twoyear call options are valued. With , the largest price difference is an increase on 0.0149 to 0.0545 and the effect of volatility of volatility is an increase for ITM calls and a decrease for OTM calls. With , the largest price difference is increase of 0.0008 to 0.01 and the effect of volatility of volatility is a small decreas that is negligible for option values. With , the largest price difference is an increase of 0.0213 to 0.0633 and the effect of volatility of volatility is a decrease for ITM and ATM calls and an increase for OTM calls. The correlation parameter has several effects depending on the relation between the strike price and the current stock price. A negative tends to produce higher values for ITM calls and lower values for OTM money calls.

We have also compared the model with the BS model, which can be interpreted as a firstorder approximation with no jumps, and . A common practice is to set the implied volatility in the BS model so that the model matches the price for the option with a strike price closest to the current stock price. For some comparisons not reported here, we have set the implied volatility in the BS model so that it matches the price generated by the stochastic volatility model for an ATM option. The BS implied volatility is very close to the expected volatility under the riskneutral distribution when shortterm options are valued. When longerterm options are used, there is a significant difference between the BS implied volatility and the expected volatility. As an approximation, the BS model tends to undervalue ITM calls and overvalue OTM calls.
6. Conclusion
The SVDEJD model incorporates several important features of stock returns. We derive a closedform solution for European call options in the model by using the martingale method, Fourier inversion transform formula, and FeynmanKac theorem. Using FFT, we obtain fast and accurate numerical solution to European option under the model. The comparison of densities of the alternative models shows that the SVDEJD model has better pricing performance on longtime options. An analysis of the model reveals that volatility of volatility and the correlation coefficient have significant impact on option values, especially longtime option, stock returns are negatively correlated with volatility, and these negative correlations are important for option valuation.
Appendix
A.1. Matlab Codes for ITM and ATM Options Pricing by FFT
function CV =inSVDexpJ(ata1, ata2, lamta, sigma, thetav, alphav,rho, sigmav, r, p, s0, v0, strike, T)x0 = log(s0)alpha = 2.55N = 4096c = 600eta = c/Nb =pi/etau = [0:N1]*etalamda = 2*b/Nposition = (log(strike) + b)/lamda + 1v = u  (alpha+1)*ik=p*ata1/(ata11)+(1p)*ata2/(ata2+1)1l=p*ata1./(ata1i*v)+(1p)*ata2./(ata2+i*v)m=sqrt((alphavi*v*rho*sigma*sigmav).^{∧}2+i*v.*(1i*v)*(sigma * sigmav)^{∧}2)n=2*m+(alphavmi*v*rho*sigma*sigmav).*(1exp(m*T))A=(2*m./n).^{∧}(2*thetav/sigmav^{∧}2)B=i*v*x0+(thetav*(alphavm)*T)/sigmav^{∧}2(i*v*rho*sigma*thetav*T)/sigmav…+lamta*T*(li*v*k1)+i*v*r*TC=(i*v.*(i*v1)*sigma^{∧}2.*(1exp(m*T)))./ncharFunc=A.*exp(C*v0+B)ModifiedCharFunc = charFunc*exp(r*T)./(alpha^{∧}2+alpha  u.^{∧}2 + i*(2*alpha +1)*u)SimpsonW = 1/3*(3 + (−1).^{∧}[1:N]  [1, zeros(1,N1)])FftFunc = exp(i*b*u).*ModifiedCharFunc*eta.*SimpsonWpayoff = real(fft(FftFunc))CallValueM = (exp(log(strike)*alpha))’*payoff/piformat shortCV= CallValueM(round(position)).
A.2. Matlab Codes for OTM Options Pricing by FFT
function CV=outSVDexpJ(ata1,ata2,lamta,sigma,thetav,alphav, rho,sigmav,r,p,s0,v0,strike,T)x0 = log(s0)alpha = 1.55N=4096c = 600eta = c/Nb =pi/etau = [0:N1]*etalamda = 2*b/Nposition = (log(strike) + b)/lamda + 1w1 = ui*alphaw2 = u+i*alphav1 = ui*alpha iv2 = u+i*alpha ik=p*ata1/(ata11)+(1p)*ata2/(ata2+1)1l1=p*ata1./(ata1i*v1)+(1p)*ata2./(ata2+i*v1)m1=sqrt((alphavi*v1*rho*sigma*sigmav).^{∧}2+i*v1.*(1i*v1)*(sigma*sigmav)^{∧}2)n1=2*m1+(alphavm1i*v1*rho*sigma*sigmav).*(1exp(m1*T))A1=(2*m1./n1).^{∧}(2*thetav/sigmav^{∧}2)B1=i*v1*x0+(thetav*(alphavm1)*T)/sigmav^{∧}2(i*v1*rho*sigma*thetav*T)/sigmav…+lamta*T*(l1i*v1*k1)+i*v1*r*TC1=(i*v1.*(i*v11)*sigma^{∧}2.*(1exp(m1*T)))./n1charFunc1=A1.*exp(C1*v0+B1)ModifiedCharFunc1 = exp(r*T)*(1./(1+i*w1)…exp(r*T)./(i*w1)  charFunc1./(w1.^{∧}2  i*w1))l2=p*ata1./(ata1i*v2)+(1p)*ata2./(ata2+i*v2)m2=sqrt((alphavi*v2*rho*sigma*sigmav).^{∧}2+i*v2.*(1i*v2)*(sigma*sigmav)^{∧}2)n2=2*m2+(alphavm2i*v2*rho*sigma*sigmav).*(1exp(m2*T))A2=(2*m2./n2).^{∧}(2*thetav/sigmav^{∧}2)B2=i*v2*x0+(thetav*(alphavm2)*T)/sigmav^{∧}2(i*v2*rho*sigma*thetav*T)/sigmav…+lamta*T*(l2i*v2*k1)+i*v2*r*TC2=(i*v2.*(i*v21)*sigma^{∧}2.*(1exp(m2*T)))./n2charFunc2=A2.*exp(C2*v0+B2)ModifiedCharFunc2 = exp(r*T)*(1./(1+i*w2) exp(r*T)./(i*w2)… charFunc2./(w2.^{∧}2  i*w2))ModifiedCharFuncCombo = (ModifiedCharFunc1  ModifiedCharFunc2)/2SimpsonW = 1/3*(3 + (−1).^{∧}[1:N]  [1, zeros(1,N1)])FftFunc = exp(i*b*u).*ModifiedCharFuncCombo*eta.*SimpsonWpayoff = real(fft(FftFunc))CallValueM = payoff/pi/sinh(alpha*log(strike))format shortCV= CallValueM(round(position)).
Acknowledgments
This work is supported by the National Natural Science Foundation of China (Grant no. 11171266) and the Science Plan Foundation of the Education Bureau of Shanxi Province (nos. 11JK0491 and 11JK0493). The authors are grateful to the reviewers as well as the editor for their valuable comments and suggestions, which led to a substantial improvement of the paper.
References
 F. Black and M. S. Scholes, “The pricing of options and corporate liabilities,” Journal of Political Economy, vol. 81, no. 3, pp. 637–654, 1973. View at: Google Scholar
 R. C. Merton, “Option pricing when underlying stock returns are discontinuous,” Journal of Financial Economics, vol. 3, no. 12, pp. 125–144, 1976. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 R. Cont, “Empirical properties of asset returns: stylized facts and statistical issues,” Quantitative Finance, vol. 1, no. 2, pp. 223–236, 2001. View at: Publisher Site  Google Scholar
 R. Cont, J. da Fonseca, and V. Durrleman, “Stochastic models of implied volatility surfaces,” Economic Notes, vol. 31, no. 2, pp. 361–377, 2002. View at: Publisher Site  Google Scholar
 E. Eberlein and U. Keller, “Hyperbolic distributions in finance,” Bernoulli, vol. 1, no. 3, pp. 281–299, 1995. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 D. B. Madan, P. P. Carr, and E. C. Chang, “The variance gamma process and option pricing,” European Finance Review, vol. 2, no. 1, pp. 79–105, 1998. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 D. Duffie, J. Pan, and K. Singleton, “Transform analysis and asset pricing for affine jumpdiffusions,” Econometrica, vol. 68, no. 6, pp. 1343–1376, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 L. H. Steven, “A closedform solution for options with stochastic volatility with applications to bond and currency options,” Review of Financial Studies, vol. 6, no. 2, pp. 327–343, 1993. View at: Publisher Site  Google Scholar
 D. S. Bates, “Jumps and stochastic volatility: exchange rate processes implicit in deutsche mark options,” Review of Financial Studies, vol. 9, no. 1, pp. 69–107, 1996. View at: Publisher Site  Google Scholar
 J. Keppo, X. Meng, S. Shive, and M. Sullivan, “Modelling and hedging options under stochastic pricing parameters,” in Proceedings of the Industrial and Operations Engineering, University of Michigan at Ann Arbor, 2003. View at: Publisher Site  Google Scholar
 S. G. Kou, “A jumpdiffusion model for option pricing,” Management Science, vol. 48, no. 8, pp. 1086–1101, 2002. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 S. G. Kou and H. Wang, “Option pricing under a double exponential jump diffusion model,” Management Science, vol. 50, no. 9, pp. 1178–1192, 2004. View at: Publisher Site  Google Scholar
 J. R. Birge and V. Linetsky, Handbooks in Operations Research and Management Science, vol. 15, NorthHolland, 2008.
 G. Bakshi, C. Cao, and Z. Chen, “Pricing and hedging longterm options,” Journal of Econometrics, vol. 94, no. 12, pp. 277–318, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 T. W. Epps, Pricing Derivative Securities, World Scientific, Singapore, 2nd edition, 2007.
 P. P. Carr and D. B. Madan, “Option valuation using the fast Fourier transform,” Journal of Computational Finance, vol. 2, no. 4, pp. 61–73, 1999. View at: Google Scholar
 P. P. Carr and L. Wu, “Timechanged Lévy processes and option pricing,” Journal of Financial Economics, vol. 71, no. 1, pp. 113–141, 2004. View at: Publisher Site  Google Scholar
 T. R. Hurd and Z. Zhou, “A Fourier transform method for spread option pricing,” SIAM Journal on Financial Mathematics, vol. 1, pp. 142–157, 2010. View at: Publisher Site  Google Scholar
 M. Schmelzle, “Option pricing formulae using Fourier transform: theory and application,” 2010, http://pfadintegral.com. View at: Google Scholar
 V. Naik and M. Lee, “General equilibrium pricing of options on the market portfolio with discontinuous returns,” Review of Financial Studies, vol. 3, no. 4, pp. 493–521, 1990. View at: Publisher Site  Google Scholar
 A. R. Hall and A. Inoue, “The large sample behaviour of the generalized method of moments estimator in misspecified models,” Journal of Econometrics, vol. 114, no. 2, pp. 361–394, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH
Copyright
Copyright © 2012 Sumei Zhang and Lihe Wang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.