Generalized Fractional Processes with Long Memory and Time Dependent Volatility Revisited

In recent years, fractionally-differenced processes have received a great deal of attention due to their flexibility in financial applications with long-memory. This paper revisits the class of generalized fractionally-differenced processes generated by Gegenbauer polynomials and the ARMA structure (GARMA) with both the long-memory and time-dependent innovation variance. We establish the existence and uniqueness of second-order solutions. We also extend this family with innovations to follow GARCH and stochastic volatility (SV). Under certain regularity conditions, we give asymptotic results for the approximate maximum likelihood estimator for the GARMA-GARCH model. We discuss a Monte Carlo likelihood method for the GARMA-SV model and investigate finite sample properties via Monte Carlo experiments. Finally, we illustrate the usefulness of this approach using monthly inflation rates for France, Japan and the United States.


Introduction
Consider the well-known ARFIMA(p, d, q) model given by: where Y t = (1 − B) d X t , d ∈ (−1, 0.5), { t } is a sequence of uncorrelated (not necessarily independent) random variables, such that Var( t ) = σ 2 , and φ(B) (stationary AR(p)) and θ(B) (invertible MA(q)) polynomials respectively.This standard case of constant variance innovations has been considered in many traditional time series analysis with applications.However, in recent years, there has been a great number of developments based on time-dependent instantaneous innovation variance (or volatility) such that Var( t ) = σ 2 t .In particular, the following cases have been considered: (i) σ 2 t is a deterministic function of t or σ 2 t = f (t), (ii) t follows the family of (G)ARCH process (see, [1,2]), (iii) log(σ 2 t ) is another stochastic process.
These cases (i) to (iii) can be analysed with emphasis on different practical issues.However, in applications, we need additional assumptions on σ 2 t , such as: (a) 0 < m < σ 2 t < M < ∞ to ensure Var(Y t ) is finite, in (i) (b) stationarity or stability of both { 2 t } and log({ 2 t }) in (ii) and (iii) Assumption (a) is imposed as it is natural to set bounds for the deterministic function σ 2 t = f (t).In fact, Assumption (b) is required when σ 2 t is a stochastic process.For case (i), the effect of non-stochastic and time-dependent instantaneous variance when d = 0 was studied by Niemi [3] under the standard AR and MA regularity conditions on the zeros of φ(B) and θ(B), respectively.In his work, Peiris [4] argues that the results of Niemi [3] can be extended to the ARFIMA family when d ∈ (−0.5, 0.5).In (ii), it is known that X t is strictly stationary if t is strictly stationary, and in particular, X t is 2m-th order stationary if t is 2m-th order stationary.For the Integrated GARCH (IGARCH) model, Nelson [5] and Bougerol and Picard [6] have argued that IGARCH is strictly stationary under additional regularity conditions.For case (iii) with stochastic volatility, it is obvious that if log σ 2  t is m-th order stationary, then t is 2m-th order stationary.Therefore, it can be argued that when σ 2  t is a stationary stochastic process, the conditional likelihood estimation carried out by taking σ 2 t = σ 2 cannot be significant since σ 2 t is bounded.This would be useful in the estimation of parameters, especially in both cases (ii) and (iii).See the survey papers of McAleer [7] and Shephard [8] for the various extensions of (ii) and (iii).
An alternative way of modelling time-dependent volatilities has been extensively studied using ARMA models with time-dependent coefficients driven by constant variance innovations.See, for example, [4,[9][10][11][12][13] and the references therein for details.However, this approach is not very attractive in applications, as it involves too many parameters to estimate.
Turning to applications in economic and financial time series, there are many popular directions of modelling and analysis of long-memory.Among others, the analysis of long-memory in inflation has been considered by Backus and Zin [14], Hassler and Wolters [15], Baillie, Chung and Tieslau [16], and Caporale and Gil-Alana [17].In their paper, Delgado and Robinson [18] considered a number of methods for the analysis of long-memory time series using ARFIMA.An alternative and a general approach is to use the ARFIMA family with conditional and stochastic volatility as considered by Baillie et al. [19], Bollerslev and Mikkelsen [20], Ling and Li [21], Breidt et al. [22], Deo and Hurvich [23], and Bos et al. [24].In their recent paper, Bos et al. [24] accommodate the stochastic volatility in ARFIMA modelling.Empirical evidence confirms that such models are very satisfactory in practice.Therefore, the aim of this paper is to extend the ARFIMA models with time-varying volatility to a general flexible class of time series models based on Gegenbauer polynomials together with the ARMA structure.The Gegenbauer ARMA (GARMA) model is a generalization of the ARFIMA model.Clearly, the former encompasses the latter as a special case.We will also extend the class of k-factor Gegenbauer process following Woodward et al. [25], Ferrara and Gueganand [26], and Caporale and Gil-Alana [17], by accommodating time-dependent volatility.
The organization of the paper is as follows.Section 2 reviews the family of GARMA with constant variance (volatility).Section 3 shows the existence and uniqueness of second-order solutions for the GARMA model with time-dependent volatility and develops new classes of GARMA-GARCH and GARMA-SV.Section 4 presents the asymptotic results for the maximum likelihood estimator for the GARMA-GARCH model and reports a Monte Carlo likelihood method for estimating the GARMA-SV model.Section 5 presents an illustrative example via simulation data, while Section 6 demonstrate an empirical example using inflation data in France, Japan and the United States.Section 7 gives concluding remarks.

Basic Results on GARMA with Constant Volatility
In this section, we review the family of GARMA processes with constant volatility.Based on the work of Gray et al. [27] and Chung [28], consider the family of time series generated by: where the polynomials φ(B) and θ(B) are as defined before, |u| ≤ 1, |d| < 1 are real parameters and t is white noise with zero mean and variance σ 2 .
This process in ( 2) is known as Gegenbauer ARMA of order (p, d, q) or GARMA (p, d, q; u) and has the following properties: • The power spectrum is given by: where 2 and i = √ −1.
In their recent papers, a modified class of generalized fractional processes has been studied by Shitan and Peiris [29,30].Now, consider the following special case or the GARMA(0, d, 0; u) process and its properties for later reference.That is, when φ(B) = θ(B) = 1, we have: Under the regularity conditions in R1, there exists the Wold representation to (4) given by: where with ψ 0 = 1 and the Gegenbauer coefficients ψ j have the explicit representation: ) is the Gamma function; see [31] for details).These coefficients ψ j , j ≥ 2 are recursively related by: with initial values ψ 0 = 1 and ψ 1 = 2du.
An Invertible Solution to GARMA(0, d, 0; u) Model Under the MA regularity conditions, there exists an invertible solution to (4) given by: where π j , j ≥ 0 are obtained from (6) replacing d by −d with corresponding initial values.The next section develops the class of GARMA(0, d, 0; u) driven by time-dependent or stochastic innovations for later reference.

GARMA with Time-Dependent Innovations
Suppose that Var( t ) in ( 4) is time-dependent.Consider the class of regular (in both AR and MA) GARMA(0, d, 0; u) driven by time-dependent innovations satisfying: where σ t is the time-dependent volatility.Below, we establish the existence and uniqueness of second-order solutions to (4) with innovations in (8) under certain additional regularity conditions.

Unique Stable Solutions
We use the following general approach: Let (Ω, A, P) be a probability space, and let L r 0 (Ω, A, P) be the space of all real-valued random variables on (Ω, A, P) with finite r-th order moments.Suppose that {ζ t } is sequence of random variables in L r 0 (Ω, A, P) and M s (ζ) is the closed linear subspace of L r 0 (Ω, A, P) spanned by the elements ζ t , t ≤ s.Let L r 0 (Ω, A, P) be the closed linear subspace spanned by all of the elements ζ t ; t ≥ 1 by M(ζ).
In the case of r = 2, let ξ and ζ be any two random elements in L 2 0 (Ω, A, P) such that the inner product and the norm satisfy respectively.It is easy to verify that L 2 0 (Ω, A, P) is a Hilbert space.Intuitively, if the volatility process is stationary, it guarantees the existence of the second moment of t , which enables us to establish the following two lemmas.Lemma 1.Under the AR regularity conditions R1, if the innovation process is stationary, then the solution in (5) belongs to L 2 0 (.), i.e., X t ∈ L 2 0 (.).
Proof.Since the innovation process t is stationary, var( t ) = σ 2 t is bounded.Hence the solution Lemma 2. Under the MA regularity conditions R2, if the innovation process is stationary, then an invertible solution to (7) belongs to L 2 0 (.) or t ∈ L 2 0 (.).
The proof is similar to that of Lemma 1.
Lemma 3.Under the both AR and MA regularity conditions and stationarity of the innovation process, one has M t (X) = M t ( ).
The proof to this follows from Lemmas 1 and 2. Next, we consider two special cases useful in applications.

Two Special Cases
Consider two popular cases where the innovations t follow GARCH or SV processes.
3.2.1.GARMA(0, d, 0; u)-GARCH(r, s) Suppose that σ 2 t in (8) follows a GARCH(r, s) process, such that where r and s are positive integers (both not simultaneously zero), α 0 > 0, α i ≥ 0 and β j ≥ 0. It is well-known that an equivalent representation of the GARCH(r, s) process is: where , and the disturbances z t and ξ t are mutually independent for all t.Further assume that the roots of ρ(z) = 0 and ν(z) = 0 lie outside the unit circle.Note that σ 2 ξ measures the conditional volatility of the log-volatility.

we have the following lemma:
Lemma 5.Under the AR regularity condition on ρ(θ), the log-volatility process in (10) has uniquely determined L 2 0 (.) solution, such that is the mean of the log-volatility process.
It is clear from (11) that the log-volatility process {h t } converges in the mean square to κ and: In their paper, Chesney and Scott [32] have considered the case where h t follows an AR(1) process, such that where κ * is a positive constant, and |ρ| < 1.
Remark.Clearly, {w t } in ( 12) is not a martingale difference series, and hence, it is not useful in applications.However, this can be written as an ARMA(1,1) in the form: This result can be extended to any general ARMA structure for the log-volatility process.
Lemma 7. Suppose that the SV process follows an ARMA(l, m) model as in (10).Then, the corresponding {ln( 2 t )} process satisfies and ARMA(k, k) in the form: Proof.The proof follows from extending the approach in Lemma 6.
Section 4 discusses the estimation of parameters.

Estimation of Parameters
This section discusses the estimation for the GARMA (p,d,q; u) model with time-dependent volatility.We divide the section into to two parts, namely, (i) the GARMA-GARCH and (ii) the GARMA-SV.
For the ARFIMA (p, d, q)-GARCH (r, s) model, Breidt et al. [22] suggested the above approach.However, Ling and Li [21] established the consistency and asymptotic normality of the corresponding ML estimator and showed that the information matrix is block-diagonal.Combining the results of Ling and Li [21] and Chung [28,33], we can obtain the asymptotic results for the corresponding estimator of the GARMA-GARCH model.Proposition 1.Let û and ( γ, δ) be approximate ML estimators of u and (γ, δ) based on a sample {X t } n t=1 from a GARMA-GARCH model under the conditions in Lemma 4.Then, û is asymptotically independent of ( γ, δ) and: and ( W1 (t), W2 (t)) and (W 1 (t), W 2 (t)) are two independent Brownian motions with mean zero and covariance where Y 1 is a random variable defined as As discussed in Chung [28,33], the convergence rates of û are faster than those of the remaining parameters.The off-diagonal blocks in the information matrix (with respect to the parameter u and the remaining parameters) approach zero.Hence, the distribution of û is asymptotically independent of the remaining parameters.
Below, we report the asymptotic result of the remaining parameters.
Proposition 2. Based on the sample {X t } n t=1 from the GARMA-GARCH model and under the conditions in Lemma 4, we have where u = 1 and Proof.See Appendix A.3.
For the case of u to be one, we can use the asymptotic result of Ling and Li [21].Following Gray et al. [27] and Chung [28], in practice, we use the grid search procedure for different value of u over the range [−1, 1] to minimize the likelihood function.For selecting the order of the GARMA-GARCH model, we can use an information criterion, such as AIC and BIC.Furthermore, we can use the conventional t test for the parameters except for u.For testing the null hypothesis regarding u, we can use the approach of Chung [33] to obtain percentiles via simulations.

Estimation of GARMA-SV Model
Suppose that X = {X 1 , • • • , X n } are generated from the GARMA (p,d,q; u)-SV model ( 2), ( 8) with the SV structure where . Since {h t } is a latent process, the evaluation of the likelihood function requires integrating it with respect to (h 1 , . . ., h n ).
As mentioned earlier, the evaluation of the likelihood function involves high-dimensional integration, which is difficult to calculate.Nevertheless, among others, Danielsson and Richard [34], Shephard and Pitt [35], Durbin and Koopman [36,37] and Liesenfeld and Richard [38] suggested evaluating high-dimensional integrals using simulation methods and then maximizing the corresponding likelihood function.In this case, we use the Monte Carlo likelihood (MCL) estimator of Durbin and Koopman [37].In their recent paper, Bos et al. [24] extended the MCL estimator for the estimation of parameters of ARFIMA-SV parameters.This approach creates a set of realized values for h = (h 1 , h 2 , • • • , h n ) by 'importance sampling'.Conditional on h and using the prediction error decomposition, it is clear that the density is given by: log(X|h , d, φ, θ, κ.ρ, σ where a t is the one-step ahead prediction error, f t is its variance and h is from importance sampling.From (19), we evaluate the simulated likelihood function based on the true density and the importance density using the results of Durbin and Koopman [37].We also extend the work of Bos et al. [24], by replacing the autocovariance functions of the ARFIMA by those of the GARMA and use the results of McElroy and Holan [39] to obtain the ML estimator.As a practical issue, we use the grid search procedure for different values of u over the range [−1, 1] for minimizing the likelihood function.
Noting that ∂h t /∂u = 0, ∂h t /∂γ = 0, we can consider that the information matrix of (u, γ , δ 2 ) has a block diagonal structure similar to the GARMA-GARCH case, where δ 2 = (κ, σ ξ , ρ) .If the MCL approximates the true likelihood accurately, we can use the conventional t test for the parameters except for u.For testing the null hypothesis regarding u, we can use the approach of Chung [33] to obtain percentiles via simulations, based on Proposition 1.For this purpose, we need to replace K in Proposition 1 with K * = E[exp(−2h t )], by noting ∂h t /∂u = 0 unlike the GARMA-GARCH process.
In the next section, we will investigate the finite sample properties of the MCL estimator.In selecting the order of the GARMA-SV model, we can use the MCL to calculate the information criterion, such as AIC and BIC.
It is worth mentioning a semi-parametric estimation procedure for the long-memory parameter under heteroskedasticity.See, for example, [40].For the GARMA-SV model, we extend this in two directions: one is to to replace the conditional heteroskedasticity by SV, and the other is to extend the ARFIMA to the GARMA.For the latter case regarding the GARMA process, Hidalgo and Soulier [41] developed a log-periodogram regression estimator extending the work of Robinson and Henry [40].
Now, look at a simulation study in order to illustrate certain properties of this GARMA-SV class.

Simulation Results
In this section, we show two kinds of simulation results for the GARMA-SV model.One is an illustrative example to show the pattern of a GARMA-SV process, and the other is the Monte Carlo results for the finite sample properties of the ML estimator for the GARMA-SV.
Take the time series of 300 observations (discarding the first 200) for further analysis.It is clear that the Gegenbauer frequency is ω G = 0.2532π = 0.7954.
For the simulated GARMA(0, d, 0; u), Figure 1 shows the time series plot, the sample periodogram and the autocorrelation functions (ACF) of series.In addition, 95% confidence intervals for the squared series are presented.The periodogram indicates that the peak is at ω = 0.2733π, which is close to the actual value.The ACFs of the raw series decay hyperbolically and periodically, confirming the features of the Gegenbauer process.On the other hand, the ACFs of the squared series are insignificant for higher lags.Figure 2 presents the plots for a simulated GARMA(0, d, 0; u)-SV process.The time series plot indicates the pattern of the time-dependent volatility.Figure 2b shows three major peaks at frequencies of 0.0467π, 0.1600π and 0.2733π.Clearly, the last peak is close to the true Gegenbauer frequency.Figure 2c shows that the ACF of the raw series decays hyperbolically and periodically (similar to the result in Figure 1c of the Gegenbauer process).The ACF in Figure 2d shows a geometric decay, which is different from Figure 1d, which shows constant volatility.As shown in Figure 1b, the periodogram of the GARMA process has one major peak at the Gegenbauer frequency, which is further from the origin.These plots show that there is a possibility of having multiple peaks under the time-dependent volatility model.The next section gives an empirical study based on a monthly consumer price index as an example.

Monte Carlo Experiments
In order to examine the finite sample properties of the MCL estimator for the GARMA-SV model, we conduct a simulation study.For the ARFIMA-SV model, Bos et al. [24] conducted Monte Carlo experiments for the MCL estimator for the ARFIMA (1, d, 0)-SV model, and found (i) downward bias in d, (ii) upward bias in φ, (iii) downward bias in ρ and (iv) upward bias in σ ξ , with the sample size T = 500.We consider the GARMA (1, d, 0; u)-SV model in order to check the effects of the additional parameter u.As explained in the previous section, we use the grid-search for estimating u over the interval [−1, 1].
Table 1 shows the sample means, standard deviations and root mean squared errors (RMSE) of the MCL estimators.For the parameter u, the bias and RMSE are smaller than those for d and φ 1 , implying a faster convergence rate for the MCL estimator of u, as expected from Theorem 1 of Chung [33], for the simple GARMA process and from Proposition 1 of the current paper for the GARMA-GARCH model.Figure 3 shows the estimated density of the û, indicating that 94% of the estimates are located between (0.65,0.75) and that 2.8% of them are greater than 0.95.For the remaining parameters, Table 1 agrees with the simulation results of Bos et al. [24], except for φ 1 .The MCL estimator for φ 1 gives a downward bias, which may be explained by the fact that not only (d, φ 1 ), but also u, affect the autocorrelation structure of the GARMA part.

Applications of GARMA-SV Models
We use the monthly consumer price index (CPI) P t in three countries: France, Japan and the United States (U.S.) to illustrate the modelling described in this paper.The sample period starts from 1960M11 to 2015M11 in all three countries.The data source is the IMF's International Financial Statistics published on the IMF web page.The CPIs are normalized at 100 for the year 2010.>From the price index series, we calculate the monthly inflation rates, Π t = 100 × log(P t /P t−1 ) giving T = 660 observations.By removing the seasonal effects, regress the inflation rates on a series of monthly dummies D t as Π t = D t δ + e t .Define X t = Π + êt , where Π is the average of inflation rates and êt is the residuals through the method of OLS.Table 2 shows the corresponding descriptive statistics.For seasonally-adjusted inflation rates, Figure 4 shows the estimated periodograms.As in Figure 4, the highest peak is at ω = 0.0033π for all three countries.Apart from the pole, there are at least two small peaks.The second highest peak is at ω = 0.6697π for France and Japan, while the peak in the second mass is at ω = 0.1606π for the U.S.Although it is difficult to distinguish, the third highest peak for France is at ω = 0.9970π.Based on the ACF and spectral properties, we examine the reasons for these multiple peaks.The first candidate is time-dependent volatility as discussed in the previous section.The second one is the general GARMA process with multi-factors, and the third is a nested model of the multi-factor GARMA and the SV.In this case, we estimate the three-factor GARMA-SV model given by: based on ( 8) and ( 18) and using the MCL estimation technique.By the sample periodogram in Figure 4, we expect that one of the three Gegenbauer frequencies is zero, corresponding to the highest peak at the point close to the origin.We do not report the higher orders of (p, q), as the estimates were insignificant (using the mean subtracted series following the recommendation of Chung [28]).The estimates of d i are located in (0.05, 0.2) and are significant at five percent, except for d 2 for the U.S. series.For all three countries, the null hypothesis of u 1 = 1 cannot be rejected.A positive value of u i indicates that the Gegenbauer frequency is close to the origin, while a negative value implies that it is close to π.Hence, the Gegenbauer frequency for U.S. data is close to the origin, while the estimates of (u 2 , u 3 ) of the other two countries are close to π.It is interesting to note that the estimates of ρ are positive and significant, indicating the appropriateness of accommodating the time-dependent volatility structure in this multi-factor GARMA model.The empirical results support this three-factor GARMA-SV model for France and Japan, while the U.S. series favours the two-factor GARMA-SV model.Since the estimates of ρ are close to one for Japan and the U.S. series, we consider a GARMA model for the SV process, {h t }.Using the MCL estimates from Table 3, we obtain the residuals to satisfy Now, we consider the interpretation of the parameter u i .First of all, the behaviour of a periodic long-memory process is different from a periodic short memory process (for example, a seasonal ARMA).As shown by Chung [28,33], the j-th autocovariance function for the GARMA process is approximated by Kj 2d−1 cos(ω g j), where K is a constant.Hence, the operator (1 − u i B + B 2 ) d i produces a periodic long-memory with cycles of every 2π/ arccos(ω g,i ).For France, the values for i = 2, 3 are 2.986 and 2.018, respectively.These results indicate that there are periodic long-range dependences with respect to every two and three months.The values of Japan are 2.986 and 3.976, implying the periodic long-memory with respect to every three and four months.For the U.S., the estimate for u 3 is 5.946, producing the periodic long-memory with respect to every six months.As u 1 is equal to one in three countries, there is no periodic long-memory for the first factor.Figure 5 shows the periodograms for the log of the squared residuals, log ˆ 2 t , which can be considered as the proxy of log-volatility.Figure 5 shows that there is a distinguished peak at zero frequency for the case of Japan and the U.S., implying that a short memory model model is adequate for France, while long-memory models are appropriate for the remaining two countries.Furthermore, Japan also has five other periodic peaks, implying a multiple periodic long-memory.For estimating the GARMA-SV model when h t follows a GARMA model, the MCL technique requires a smoother simulation to be applied to the second GARMA process, which is an extension of the work of de Jong and Shephard [43].Another task is to reduce the computational time under the multiple grid-search method for finding optimal values of u i 's in mean and volatility.We need to wait for further research for these problems.
In this section, we found that in addition to non-periodical long-range dependence, the empirical results indicate the existence of the periodic long memory under the time-dependent volatility.

Conclusions
This paper considers the class of GARMA-GARCH and GARMA-SV models in detail.We have established the existence and uniqueness of second order solutions and investigated the asymptotic results of the approximated ML estimator for the GARMA-GARCH model.We have explained the Monte Carlo likelihood (MCL) technique for estimating the GARMA-SV model.Using a set of simulated data, we have presented an illustrative example with possible multiple peaks in the GARMA-SV model.We also conducted a Monte Carlo experiment to investigate finite sample properties of the MCL estimator.In the empirical analysis, we have estimated a three-factor GARMA-SV model for monthly inflation rates of France, Japan and the United States and found that the data favour multi-factor GARMA-SV models showing that periodic long-memory exists.
There are many possible extensions of this new approach.In a future paper, we will develop models allowing both the mean and variance to follow Gegenbauer processes, as discussed in the previous section.For this purpose, we may extend the MCL estimation method by considering the simulation smoother for the GARMA process.By extending the work of Robinson and Henry [40], we plan to develop an estimation method for the Gegenbauer frequency and the long-memory parameter in the presence of long-memory in volatility.
where  [16], Ling and Li [21] and Weiss [44].This will not affect the asymptotic efficiency and other asymptotic properties.
We obtain the second derivatives, by directly differentiating ∂l t /∂γ, ∂l t /∂u and ∂l t /∂δ.For using later, we only give: Let F t be the past information of t available up to time t.Then, we obtain the conditional expectation of the second derivatives: Below, we will deal with u and (γ , δ ) separately, since their speeds of convergence are different, as shown by Chung [33].
Noting that we apply the proof of Theorem 3.1 of Ling and Li (1997) for (γ , δ ) to obtain the asymptotic properties of the information matrix: as n → ∞, and Ω γ and Ω δ are positive definite matrices, where Ω γ and Ω δ are defined in (17).We can show that the information matrix is a block diagonal, in the following way.We can express each element of ∂σ 2 t /∂δ as the infinite sum 2 t−i (i = 1, 2, . . ., ), while each element of ∂σ 2 t /∂γ has a representation of the infinite sum of t−i t−i−j (j = 1, 2, , . ..).Thus, we have E[(∂σ  [21], we can also show that −n −1 ∑ n t=1 (∂ 2 l t /∂δ∂u) −→ a.s.O.The block diagonal structure of the information matrix with respect to (u, γ ) and δ implies that the distribution of ( û, γ ) is asymptotically independent of that of δ.Hence, we focus on the first derivatives and the information matrix regarding (u, γ ) for deriving the asymptotic properties of û, in the next subsection.
For deriving the asymptotic properties of the information matrix regarding û, we present the results of Chung [33].
and the results indicate that when |u| < 1 and when u = ±1: As discussed in [21], we also show that every element of I γu = E [(∂ t /∂γ) (∂ t /∂u)] is finite.Using a similar approach as in Theorem 3.1 of [21], we have E[∂ 2 L(λ)/∂γ∂u] = O(1).By the above results for I u , we show that E[∂ 2 L(λ)/∂u 2 ] does not exist, indicating that the usual asymptotic theory based on O p (n −1/2 ) convergence will not work as discussed in [33].
Appendix A.2. Proof of Proposition 1 Since the information matrix is block diagonal, we focus on the parts related to (u, γ ) in order to prove the proposition.
We consider the following Taylor series expansion of the first conditions for the maximization of the approximate likelihood function around the true value of λ , where are the rates suggested by Chung [33].
then we obtain since û is asymptotically independent of γ.In the remainder of the proof, we derive the limiting distributions of the two series in (A6) and then that of c n ( û − u) in (A7).For this purpose, we use the theory of nonstationary ARMA process with GARCH errors, derived by Ling and Li [45].By the conditions of Lemma 4, we can write Assuming pre-sample values, 0 = 1 = 2 = • • • = 0 and combining with (A1), we have Then, from (A2), we have Below, we apply the theorems of [45] to derive their limiting distributions.For this purpose, we separately consider three cases: |u| < 1, u = 1 and u = −1.Applying (A5) and the proof of Theorem 3.2 of [21], we can show that the conditions provided by Basawa et al. [46] are satisfied.The only difference is on the third derivative of t with respect to d: which has the finite second moment as:

Figure 3 .
Figure 3.Estimated density of the MCL estimator of u.

Figure 4 .
Figure 4. Sample periodograms of inflation rates.The inflation rates are seasonally adjusted via monthly dummy variables.

Table 1 .
Finite sample performance of the Monte Carlo likelihood (MCL) estimator for T = 500.

Table 2 .
Descriptive statistics of inflation rates.
Note: Π t denotes the original inflation rate, while X t denotes the seasonally-adjusted series.

Table 3 .
QML estimates of the multi-factor GARMA-SV model.