The effect of intraday periodicity on realized volatility measures

We focus on estimating daily integrated volatility (IV) by realized measures based on intraday returns following a discrete-time stochastic model with a pronounced intraday periodicity (IP). We demonstrate that neglecting the IP-impact on realized estimators may lead to invalid statistical inference concerning IV for a common finite number of intraday returns. For a given IP functional form, we analytically derive robust IP-correction factors for realized measures of IV as well as their asymptotic distributions. We show both in Monte Carlo simulations and empirically that the proposed bias corrections are the robust way to account for IP by computing realized estimators.


Introduction
For the majority of financial markets, the pervasive intraday periodicity (IP), which often takes a U -or mirrored J -form during the daily trading time, is a well documented empirical feature of intraday absolute returns (Wood et al. 1985;Harris 1986). As it appears to be highly correlated with intraday variation of trading volume, Admati and Pfleiderer (1988) propose to explain daily IP-shape by strategic interaction of traders around market openings and closures, whereas the periodicity at weekly or monthly horizons could be attributed to the impact of slowly varying macroeconomic fundamentals (Andersen and Bollerslev 1998b;Andersen et al. 2001Andersen et al. , 2003. The availability of high-frequency data allows for the construction of precise estimators of daily integrated volatility (I V ) for risky asset returns. The realized volatility (RV ) defined as a sum of squared intraday returns is known to be a consistent estimator of daily I V in absence of jumps. Other realized measures such as the bipower variation (BV ) should be used for I V estimation in presence of finite-activity jumps during the day (cf. Aït-Sahalia and Jacod 2014). Barndorff-Nielsen and Shephard (2004) derive the asymptotic properties of these quantities for the number of intraday returns M → ∞ under mild assumptions on the corresponding pricing process. In practice, however, the number of intraday returns available for computing realized estimators often remains limited due to insufficient liquidity and/or irregular trading activity. Even for highly liquid stocks a practitioner could prefer to rely on simple realized estimators based e.g. on 5 min returns primarily in order to escape from adverse effects of market microstructure noise (MMN) which are particularly pronounced at (ultra) high sampling frequencies (cf. Aït-Sahalia and Jacod 2014).
In this paper we analyze and quantify the impact of IP on the finite M properties of RV and BV estimators by providing corresponding formal statements given the IP functional form. To the best of our knowledge, this research agenda has not been explored yet, although there is a vast amount of literature devoted to modeling and estimating IP (cf. Engle et al. 1990;Andersen et al. 2019;Christensen et al. 2018). Thus, our investigation provides useful insights for exploring differences between the asymptotic theory and the practical finite sample performance of realized measures based on intraday data.
To model intraday returns, we presume a discrete time stochastic specification as in Andersen and Bollerslev (1997), where the variance of intraday returns is written as a product of the deterministic periodic and stochastic volatility (SV) components. The IP is assumed to be constant for all days, whereas the SV part is slowly changing over time. Our framework is motivated by the empirical evidence that the IP captures a vast part of intraday volatility variation whereas the SV impact is of a smaller order (cf. Christensen et al. 2018). We show that for the commonly available finite number of intraday returns M neglecting the impact of IP would lead to non-valid statistical inference concerning daily I V . For a given IP, we compute the first and the second moments of RV and BV , moreover, we establish the asymptotic bivariate distribution of these measures as M → ∞. We also quantify the impact of IP on realized tri-power (T P) and quad-power (Q P) estimators of daily integrated quarticity (I Q) required for statistical inference about I V .
Our major finding is that for the commonly available finite number of intraday returns the impact of IP should be explicitly addressed when making statistical inference concerning I V . While the RV estimator of I V is unaffected by IP, BV has a finite sample bias negligible only for large sample sizes which are not always available in practice. Moreover, by estimating I Q one should account -at least for small M values -for scaling factors, which depend on the functional form of the IP. We derive the explicit expressions for these IP-correction factors and provide the asymptotic distribution for their estimators.
Our theoretical results are illustrated in a Monte Carlo study, where we investigate the impact of IP on various realized measures. We find that the IP-correction procedures proposed in this paper are helpful against the adverse impact of IP on realized measures in finite samples. Moreover, our IP-corrected estimates are advantageous compared to immediate removing of the estimated IP due to their robustness with respect to IP misspecifications at days with unusual pattern of intraday volatility (cf. Gabrys et al. 2013;Kokoszka and Reimherr 2013). In such situations our approach is preferable in terms of relative bias and mean squared error (MSE) compared to the standard procedure with immediate scaling of intraperiod returns by the estimated IP profile as in Boudt et al. (2011), for example. In the empirical application we estimate the IP for the daily volatility of the Dow Jones Industrial Average Index with and without IP bias corrections.
The remaining part of the paper is organized as follows. In Sect. 2, we introduce the model for intraday returns and discuss the realized estimators of daily I V and I Q.
The theoretical results are derived in Sect. 3 where we establish both finite sample and asymptotic stochastic properties of commonly applied realized estimators for a given IP form. Moreover, we provide expressions for IP correction factors and derive the asymptotic distributions of their estimators. Our approach is illustrated in Sect. 4 by means of a simulation study and in Sect. 5 by an empirical application. Sect. 6 concludes, whereas the proofs are placed in the Appendix.

Measuring daily volatility based on intraday information
Before we introduce our model in Sect. 2.2, we provide definitions of the objects which are of importance for our analysis. For this purpose we start from a general jump-diffusion model for log price increments in order to define the objects of our interest, namely the daily I V and daily I Q. Next, we present the corresponding realized estimators which are based on M intraday returns. These realized measures are consistent estimators of I V for M → ∞, however, in practice we often face M ≤ 10 2 primarily because ultra high frequency returns are contaminated by MMN. Thus, it is of importance to explore the finite sample stochastic properties of realized estimators. For this purpose we then consider a discrete time model for intraday returns with an explicit functional specification of intraday periodicity (IP) and study the impact of IP on the realized measures of I V and I Q for finite M.

Model for intraday returns and realized measures
Assume that log-prices of risky assets p(t) = ln P t follow a continuous time process with (possible) additive jump components. We consider a day t as the period of interest with the daily return r t = p(t) − p(t −1) and focus on the integrated volatility (I V ), which is defined for day t as where σ (u) is a spot volatility. In order to make statistical inference about I V measures one also needs statements concerning the daily integrated quarticity (I Q) defined by The availability of intraday returns allows to construct precise realized estimators (Andersen and Bollerslev 1998a) for the daily I V which are of immense practical importance for estimation and inferences concerning daily volatility. Assume that M equally spaced intraday returns are available for day t. We denoted them by r t,m = Then the daily return r t is the sum of intraday returns with r t = M m=1 r t,m . The most popular I V estimator is the realized volatility (RV ) measure given as Barndorff- Nielsen and Shephard (2004) show that RV t is a consistent estimator of I V t without jumps at day t, i.e. RV t p −→ I V t as M → ∞. Although the estimator RV t possesses a set of appealing properties, it is not appropriate in the presence a non-zero jump component.
The bipower variation (BV ) proposed by Barndorff-Nielsen and Shephard (2004) is a jump-robust estimator of I V . It is consistent even in presence of jumps, i.e. BV t p −→ I V t as M → ∞. However, RV has a smaller variance than BV if there is no jumps, so the common practice is first to test for a jump during day t. Then, in case of a significantly large positive distance between RV and BV indicating jumps, one should apply BV ; otherwise RV is to use.
Intraday returns are also suitable for the purpose of estimating the unknown daily I Q required for computing variances of RV and BV measures. The realized quarticity (R Q) is a consistent estimator of I Q in the case of no jumps with R Q t p −→ I Q t as M → ∞. However, as R Q t measure is not robust (cf. Andersen et al. 2014), Barndorff-Nielsen and Shephard (2004) suggest to use the realized tri-power (T P) and quad-power (Q P) measures defined by |r t,m−3 ||r t,m−2 ||r t,m−1 ||r t,m |, Although both RV and BV have appealing stochastic properties as M → ∞, their practical implementation is often based on (say) 5 min intraday returns which makes M = 78 intraday observations for a usual 6.5-hours trading day, because of MMN which hinders the use of ultra high frequency data for construction of realized estimators (McAleer and Medeiros 2008). To overcome these problems, the recent research has been focused on making realized estimators more robust to these features. However, in many situations the common practice still remains to sample returns at a lower frequency, i.e. to consider 5-, 10-, or even 15-min intraday returns (Andersen et al. 2011). We follow this strand of literature and concentrate on profound understanding of stochastic properties of realized estimators for comparatively small values M.
Although we focus in our analysis primarily on the classical RV , BV , Q P, and T P measures, we also provide a discussion of recently proposed further realized measures in Sect. 5.2 and in Dette et al. (2022).

Discrete time model for intraday returns
The IP in absolute intraday returns is one of the most important stylized facts characterizing high frequency data. In order to investigate the IP impact on realized measures for a fixed number of intraday returns M, next we introduce a discrete time model in (4) which is central for study. There is a substantial scope of recent literature concerning discrete-time modeling of intraday returns whereas the IP is assumed to be a multiplicative scaling component (Boudt et al. 2011;Engle and Sokalska 2012;Bekierman and Gribisch 2021). Following Andersen and Bollerslev (1997), we focus on a discrete stochastic model for intraday return without jumps, which is given as where s t,m > 0 is the deterministic IP component and γ 2 t,m > 0 is the stochastic part. Note that there is no leverage effect in (4), as it is mostly of importance for daily returns but much less pronounced for high-frequency intraday returns, see e.g. Bollerslev et al. (2006).
For our theoretical derivations we presume that the stochastic part remains constant within day t (Andersen and Bollerslev 1998b;Hecq et al. 2012), i.e. γ t,m = σ t for all m = 1, . . . , M, but may change from one day to another. This assumption is justified by the empirical evidence that IP commonly accounts for a vast part of intraday volatility variation (cf. Christensen et al. 2018). In Sect. 4 we relax this assumption in the Monte Carlo simulation study by considering the intraday SV which follows a diffusion process as e.g. in Goncalves and Meddahi (2009). Based on the results of our Monte Carlo simulations, we conclude that our major findings also hold in the SV setting.
In line with the current literature (Hecq et al. 2012), we set the IP as constant at different days so that we further skip the time index with s t,m = s m . Moreover, the periodic component is standardized such that it sums up to M over the day with M m=1 s 2 m = M. Of course, a special case s m = 1 for all m = 1, . . . , M corresponds to no IP. Putting all together, we separate the intraday periodic component s m which is solely responsible for intraday heteroskedasticity, and interday stochastic component σ t which could change from one day to another by writing Although the model in (4) and (5) is fairly simple, it allows a detailed analysis of the IP impact on popular realized measures of the objects of our interest, which are the I V for day t as well as the I Q written by (Andersen et al. 2014) Note that as in general it holds that M m=1 s 4 m ≥ M, the I Q is directly influenced by IP. We aim to investigate the impact of IP on realized measures of I V and I Q.

The impact of intraday periodicity on RV and BV
To gain results on the IP impact on realized measures, re-write the normalized IP {s m } M m=1 as where g : [0, 1] → R + with R + := (0, +∞) is a given non-normalized function. This normalization is very common in IP-literature (cf. Andersen and Bollerslev 1997, p. 153). The functional form of g(·) could be quite flexible and is subject to very general regularity conditions specified in the following propositions.

Bias and variance of realized estimators
For the discrete model of intraday returns (4)-(5) and the IP from {g(m/M)} M m=1 , we derive expectation, bias and variance of RV and BV estimators of daily I V in the next proposition. (8) for some function g :

BV t is an asymptotically unbiased estimator of I V . (B) The (co)variances of RV t and BV t are given as
and Cov(RV t , (1)). The asymptotic scaling factor ξ is defined by (10) there is no IP.
The property R M ∈ (0, 1) follows from the proof of Proposition 1 in the Appendix. More precisely, by the Cauchy-Schwarz inequality, we have Thus, in the case of IP, RV is an unbiased estimator for I V but BV has a finite M bias which should be corrected in applications. Since the expectation of BV is given by we suggest the following bias-corrected measure Hence, it holds that BV t = (M − 1) which should be replaced by its empirical counterpartζ M based on IP estimatesŝ m in practice. By the same principle, we define the IP factor ξ M in R Q for finite M values as with lim M→∞ ξ M = ξ as in (10) of Proposition 1. Note that in case of no IP it holds that ξ M = ξ M,R Q = 1. Of course, in applications we replace ξ M,R Q by its estimator ξ M =ξ M,R Q which is discussed below in Sect. 3.3. Hence, we show analytically that although for BV it holds that lim M→∞ ζ M = 1, its counterpart ξ for the realized measures of I Q is still present and could be (depending on the data) rather substantial.
In the context of the IP-bias corrections for further realized measures, we investigate by means of numerical analysis several popular MMN-robust realized estimators such as min RV or med RV in the follow-up paper of Dette et al. (2022).
In the next proposition we provide the expectations of realized estimators R Q, T P and Q P serving as the measures for I Q under the model defined in (4) and (5) with the function g(·).
Proposition 2 Assume that the IP is given by (8), then the expectations of R Q t , Q P t and T P t are given as E[ is square integrable, where the asymptotic scaling factor ξ which comprises the impact of IP is given by (10).
Proposition 2 suggests the following factorization of Q P and T P measures of I Q in finite samples: whereby ξ M,Q P and ξ M,T P are the IP scaling factors for Q P and T P, respectively. This factorization appears to be useful both in the simulations in Sect. 4 and in the empirical study in Sect. 5.

Asymptotic distribution in case of intraday periodicity
Next, we provide the corresponding bivariate limit distribution for RV and BV as M → ∞ for our discrete time model of intraday returns with IP. (4) and (5) without jumps and assume that the IP component is given by (8) with a continuously differentiable function g : [0, 1] → R. Then, as M → ∞,

Theorem 1 Consider model
The integrated quarticity I Q t = ξ · σ 4 t can be consistently estimated by R Q t , Q P t , or T P t .
Thus, a pronounced IP with ξ > 1 causes more variability of I V estimators compared to the case of no IP where ξ = 1. The asymptotic (1 − α)-confidence interval for daily I V and I Q based on RV and R Q measures is given according to Theorem 1 as where z α/2 is the (α/2)-quantile of the standard normal distribution. The asymptotic confidence intervals based on BV , T P or Q P measures are constructed similarly. The results in Proposition 2 and Theorem 1 are useful for statistical inference on I V . As R Q is not robust, one could estimate I Q in presence of jumps by either T P or Q P directly as in (2) or (3), i.e. without estimating ξ separately. However, as we show in the Monte Carlo simulation, the approximation T P t ≈ I Q t is only precise for fairly large M. For this reason, for finite M we recommend to use the IP-scaled estimators Q P t or T P t as well as the estimated scaling factorξ M,R Q from Eq. (15) below for the construction of confidence intervals as

Estimation of IP correction factors
and use the statisticŝ as the estimators of ξ M and ζ M , respectively. If the variance exists, it follows from (5) (14). As a consequence, we expect that the estimators in (14) and (15) are consistent for s m , m = 1, . . . , M, as well as for ξ M and ζ M , respectively. Note that as long as the ergodicity condition (16) is met, the assumption of constant intraday volatility as in Sect. 2 could be relaxed for the following analysis in Sect. 3.3, e.g. by allowing intraday stochastic volatility.
In order to make the intuitive arguments above more precise we investigate in the following the asymptotic distribution of the statisticsζ M andξ M for finite M and T → ∞. Note that under common assumptions, such as mixing conditions (see, for example Dehling et al. 1986, among many others) or physical dependence conditions   (16) and (17) hold, then, as T → ∞, that

Proposition 3 Assume that
where the asymptotic variances are given by where we put s 0 = s M+1 = 0.
When estimating IP elements s m , we require for consistency that the number of days T → ∞, whereas the number of intraday observations M is here fixed. Hence, we derive the distributions of the correction factorsζ M andξ M for fixed M and T → ∞ in Proposition 3. Differently, inference about the asymptotic correction factor ξ in Theorem 1 would require that M → ∞ because thenξ M is a consistent estimator of ξ . These results allow to make statistical inferences and conduct tests for IP correction factors ξ M and ζ M with T → ∞. Of course, an extension of our theoretical findings for IP estimated from a finite sample T is also of interest. However, this is also a quite challenging task, for which Christensen et al. (2018) present with some theoretical considerations in their Proposition 3.1. We investigate this issue in the Monte Carlo simulations in Sect. 4.2.
In Proposition 3 we provide the explicit asymptotic results for IP scaling factors given the SD estimator. In general, our results for the SD estimator could be extended for any consistent estimator of IP. In this paper both in our simulations and empirical studies we use a more robust WSD estimator of Boudt et al. (2011) which is described in the Appendix. However, it is much more difficult to get analytical results such as in Proposition 3 for this more complicated WSD estimator which is based on order statistics. For this reason we recommend to make statistical inferences for the WSD approach by using bootstrap procedures, as e.g. in Goncalves and Meddahi (2009) or in Dette et al. (2022).

Simulation study
We illustrate our theoretical findings by means of an extensive Monte Carlo simulation study which is structured as follows. First, we introduce the U -shaped IP functional form and discuss the parameter choice for the model in (4)-(5). We generate intraday returns for estimation of IP shape and construction of various realized measures. In Sect. 4.1 we study the impact of IP on various realized estimators. In Sect. 4.2 we investigate our IP-corrections in terms of MSE whereby we compare our approach with those of Boudt et al. (2011).
We consider M = 26, M = 78, or M = 390 intraday returns, which roughly correspond to sampling at 15 min, 5 min, or 1 min for 6.5-hour trading days, respectively. For constant intraday volatility, we fix I V = σ 2 = 1 and I Q = ξ · σ 4 = ξ . We generate M intraday returns for each of T = 10 4 days with r t,m ∼ N (0, γ 2 t,m s 2 m ) where γ 2 and s 2 m are the respective SV and IP components. As a baseline, we set γ 2 t,m = 1/M which corresponds to I V = 1 and I Q = ξ . Later we also consider a SV model where we assume that γ 2 t,m is governed by the process where t = 1/M and u t,m ∼ iid N (0, 1), with γ 2 1,1 = 1. This model is a discretized version of the GARCH (1,1) diffusion used by Andersen and Bollerslev (1998a), Goncalves and Meddahi (2009) with parameters implying an autoregressive persistence; in the time series context it is related to state-space models for realized volatilities (cf. Golosnoy et al. 2021). In presence of intraday SV there is no exact analytical expression for the impact of IP on realized measures. However, since the SV in (24) follows a highly persistent process, the empirical contribution of γ 2 t,m to intraday heteroskedasticity is of a smaller order compared to the impact of IP (see e.g., Christensen et al. 2018;Bekierman and Gribisch 2021). Moreover, from the empirical perspective, one could precisely estimate the SV components γ t,m ex post, see Bekierman and Gribisch (2016). Then one could calculate the product componentŝ s t,m =ŝ mγt,m and compute the IP correction factors for this day t based on the obtained estimatesŝ t,m .
The IP components follow a quadratic convex U -shaped given by  (25) is motivated by our empirical findings, see Fig. 5. Note that our theoretical results are applicable for any functional form of IP including asymmetric mirrored Jshape specifications, as e.g. a more flexible asymmetric U -shaped IP specification as in Hasbrouck (1999) and Andersen et al. (2012). Because of M m=1 s 2 m = M, it holds for (25) that c 2 = 12 (1 − c 1 ) / (M 2 + 2). We select c 1 ∈ {0.01, 0.11, . . . , 0.91, 1} with the most IP curvature for c 1 → 0 + whereas c 1 = 1 is the no-IP case. Next we focus on the value c 1 = 0.71 corresponding to the evidence from U.S. stock market. Note that the IP form could be very distinct at 'special' days characterized by specific announcements and/or unexpected events where the IP curvature could be much more (or less) pronounced.
The asymptotic scaling factor ξ defined in (10) is plotted as a function of c 1 in Fig. 1 where we observe that ξ is substantially larger than one even for c 1 close to one. As the IP is unknown in practice, we construct estimatorsĝ(·) andŝ 2 (·) by applying the W S D estimator of Boudt et al. (2011) which does not require a-priori specification of the IP functional form; the implementation details for the W S D estimator are provided in the Appendix. Then we compute the average estimates of the scaling factorsξ M,R Q , ξ M,T P ,ξ M,Q P which are plotted in Fig. 1. We observe that the factorξ M,R Q appears to be very close to ξ M even for a fairly small value M = 26.

The impact of IP on realized measures
As RV measure is not affected by IP, we study the IP impact on other measures, such as BV for I V , and T P, Q P for I Q. After generating IID intraday returns as specified above, we calculate BV t , T P t , and Q P t as well as the IP-corrected estimators denoted by BV t , etc. for each day t = 1, . . . , T . Additionally, we compute the jump-robust medRV and minRV measures of Andersen et al. (2012). Then we build time averages, e.g. BV = (1/T ) · T t=1 BV t , for all measures. These averages of BV , BV , medRV , and minRV are shown in Fig. 2; whereas of T P, Q P, T P, and Q P in Fig. 3 for different M and c 1 values. All measures should be equal to one for no IP with c 1 = 1. In Fig. 2 the IP-bias in BV is quite pronounced for M = 26 and M = 78 for large and medium curvatures and is still visible even for M = 390. Remarkably, the biases in medRV and minRV are even stronger than in BV . As expected, the bias-corrected mean of BV is close to the true I V for all M, so the suggested correction functions properly.
The averages of the original T P, Q P and scaled T P, Q P are reported in Fig. 3. The original measures are downward biased for finite M = 26, 78, compared to almost unbiased measures for M = 390 which is close to the asymptotic value ξ , see Fig. 1. Remarkably, in case of 5 minute returns with M = 78 the bias is still quite substantial for empirical relevant values of IP curvature parameter c 1 ∈ [0.6, 0.8]. As expected, Q P is more biased than TP for finite M due to longer lags involved in its computing. Our scaled T P and Q P measures are equal to the value σ 4 = 1 for all considered values of c 1 and M which is an appealing property. Summarizing, the IP has a substantial influence on BV , minRV , medRV , so that IP-corrections are needed to get valid statistical inference on I V .
In order to shed light on the finite sample validity of the result in Proposition 3, we provide the exemplary QQ-plots in Fig. 4 for the statisticsζ M andξ M for the case of M = 78 and c 1 = 0.6 which are based on the W S D estimator of IP (cf. Boudt et al. 2011). Being standardized properly, both statistics seem to approach normality with the increase of estimation period T , in particular, even T = 250 the QQ plots show a decent fit.

The comparison of IP-bias corrected estimators
In Proposition 1 we show analytically that the original BV is downward biased so that the true level of risk measured by I V is underestimated. This is a rather undesired scenario from the risk management point of view making an MSE comparison of biased and unbiased measures not reasonable because the MSE is symmetric for upward and downward biases. For this reason, we provide an MSE comparison only for IP-bias corrected estimators.
In particular, we contrast the relative bias and MSE of our corrected estimator BV t in (11) with those of Boudt et al. (2011) where the estimated IP componentŝ m is immediately removed from intraday returns by computing r * t,m = r t,m /ŝ m . Then the BCL-estimator BV * t is given as This immediate removing of IP appears to be a common approach in the current literature (cf. Golosnoy et al. 2012;Bekierman and Gribisch 2021;Christensen et al. 2018), whereby the IP estimatesŝ m are based on historical data. To contrast the relative biases and MSEs of BV * t and BV t , we generate intraday returns for T = 250 pre-sample days with the IP parameter value c h 1 = 0.5 as in (4) and (25) and use them to get IP estimatesŝ m . We denote by c h 1 the 'historical' value of c 1 assumed to be constant during the pre-sample period. Then, we focus on the next day's IP which functional form is described by the 'current value' of c 1 denoted by c c 1 . That means, for this next (single) day of our interest, intraday returns are generated with either unchanged current IP parameter c c 1 = 0.5 = c h 1 or changed parameter c c 1 = c h 1 . The latter case could occur at some 'special' days characterized by announcements, unexpected events etc. Note that the values c c 1 = 0.1 (extremely pronounced IP) and c c 1 = 0.9 (almost no IP) are not empirically relevant but considered for the illustration purposes only. We calculate BV * and BV for this day of interest using the historical IP estimatesŝ m , so there is an IP misspecification when c c 1 = c h 1 . We repeat the procedure (generating T = 250 pre-sample days and an additional 'day of interest') 10 4 times and put the computed relative biases and MSEs in Table 1, where we show results for constant intraday volatility in Block A and intraday SV is Block B.
Hence, we study both the effect of estimation risk in case of unchanged IP and the effect of a change in the IP form. The latter is of much practical importance, as there are many empirical confirmations for time variability of IP even after excluding days with important macroeconomic announcements, see Andersen et al. (2001), Hecq et al. (2012), or Andersen et al. (2019. We observe in Table 1 that in case of no IP change with c h 1 = c c 1 , the realized measure BV * of Boudt et al. (2011) is preferable in terms of relative bias and MSE. These findings are in line with the recent theoretical results of Ghysels et al. (2021) who show under similar assumptions that the most efficient quarticity estimators (in terms of MSE) can be obtained by an immediate adjustment of intraday returns for the IP as e.g. by Boudt et al. (2011) which is also confirmed by our evidence. However, even a comparatively small change in IP, e.g. from c h 1 = 0.5 to c c 1 = 0.7, leads to a substantial increase in the relative bias and MSE of BV * . This evidence remains also for intraday SV in Block B. Hence, our IP-correction approach is more robust compared to the procedure of Boudt et al. (2011) which is rather sensible to even small changes in the IP form.

Empirical study
In our application we work with intraday returns for the Dow Jones Industrial Average Index with the focus on measuring daily I V . Our dataset consists of intraday observations from January 1996 to December 2010 with non-regular trading time days skipped. Days with non-regular trading times are those where the trading time have been substantially shorter than the common 6.5 hours because of these or that reasons. We consider 15 min, 10 min, 5 min and 2 min intraday returns; 5 min returns is the most popular choice in practice. To avoid the opening bias effects, we skip the first daily observation which is the common practice to escape from the impact of noisy overnight quotes. The final sample consists of 3329 days with M = 24, 37, 76 or 193 observations for 15, 10, 5 or 2 min frequency, respectively. We estimate the intraday IP s m , and calculate the IP correction factors ζ and ξ . Then we present descriptive statistics for both uncorrected and IP-corrected realized measures. Finally, we discuss the impact of IP bias in some further realized measures which are proposed in the literature.

Estimation of intraday pattern and descriptive statistics
We estimate the IP with the non-parametric W S D estimator of Boudt et al. (2011) and show them in Fig. 5 with componentsŝ m normalized such that M m=1ŝ 2 m = M. The IP pattern has a convex U -shape for all considered sampling frequencies. It is high during morning and afternoon hours and low during the lunch break. In numerical terms, it is about twice as high during the peak in the morning compared to the trough in the middle of the day. Hence, in our empirical application for Dow Jones index the IP estimate approximately corresponds to c 1 ≈ 0.7 in our simulation study. However, the IP curvature could be more pronounced for individual stocks. For example, in   Table 2 with 95% confidence intervals obtained by the bootstrap procedure with 10 3 − 1 replications outlined in Dette et al. (2022). The estimated factorξ M,R Q is larger than one and numerically similar for all sampling frequencies which corresponds to the evidence of similar patterns in Fig. 5. We also estimate finite sample corrections for BV , T P, and Q P. The correction estimatê ζ M for BV gets closer to one with increasing sampling frequencies. The same holds for other finite M scaling factors so that BV is biased downward. For example, we getξ M,Q P ≈1.23 for 15 min frequency but onlyξ M,Q P ≈1.03 for 5 min frequency returns. This evidence indicates that the IP bias problem should not be very acute for the U.S. stock market during the considered period of time.
Next, we provide descriptive statistics of realized estimators. In Table 3 we report the full sample averages of both uncorrected and IP-corrected realized measures of  daily I V and I Q. The average RV is stable for all sampling frequencies, whereas the average BV is biased downwards compared to RV even after the IP correction, however, the corrected BV gets numerically closer to RV . The average values for T P and Q P are almost the same for 15 and 10 min frequencies, but are much higher for 5 min and 2 min returns. Note that there are reported difficulties in estimating I Q based on frequencies of 5 min or higher (Andersen et al. 2014). Another important quantity is the average relative component (RV − BV )/RV which, as expected, gets substantially smaller with the IP-bias corrected BV both for all sampling frequencies.

Discussion and further extensions
Our study is primarily focused on the popular realized measures which are widely used in practice. Recently, there have been much developments in estimation of I V and I Q based on high frequency observations. In particular, we would like to mention the robust threshold power variations (Corsi et al. 2010), the minRV and medRV estimators (Andersen et al. 2012), and the nearest neighborhood generalizations (Andersen et al. 2014). Moreover, in order to mitigate MMN-problems which arise by using ultra high frequency returns, several noise robust methods have been developed, as e.g. the pre-averaged BV (Podolskij and Vetter 2009), pre-averaged threshold RV (cf. Aït-Sahalia and Jacod 2014), or pre-averaged threshold BV (Christensen et al. 2014). The IP bias problem is also present for some of these more advanced realized measures as it could be shown by a simple Monte Carlo simulation exercise. To illustrate this point, we plot the IP-biases for minRV and medRV measures of Andersen et al. (2012) in Fig. 2 and observe that they are even larger than for the original BV estimator. The same argument applies to the minR Q and medR Q measures of I Q which are elaborated by Andersen et al. (2014). Of course, the IP-correction factors should be newly derived for each particular measure, however, this challenging task is beyond the scope of our paper and is left for future research, as e.g. the task of deriving IP correction factors for realized portfolio weights (Golosnoy et al. 2019, 2020, or Golosnoy and Gribisch 2022. The impact of IP on these additional estimators is investigated in the follow-up paper of Dette et al. (2022).
(BV ) which are commonly applied for estimation of daily integrated volatility (I V ), or tri-power (T P) and quad-power (Q P) variations which are serving as measures for daily integrated quarticity (I Q).
In this paper we investigate the impact of intraday periodicity (IP) on the finite sample properties of these realized measures. For our analysis we assume a discrete time model for intraday returns on risky assets and postulate a multiplicative deterministic IP component which is often of U -shape empirically. For a number of intraday returns M → ∞ the impact of IP is asymptotically negligible, however, we show that the IP-impact should be taken into account for a practically relevant situation with finite M. In particular, we prove that finite sample corrections of BV as well as of T P and Q P measures are necessary to obtain valid statistical inferences concerning daily I V . We derive analytically the factors for IP-correction and analyse their stochastic properties. Our results are illustrated by means of a Monte Carlo simulation study for both constant and stochastic intraday volatility models. Finally, we evaluate IP correction factors empirically for daily I V of the Dow Jones Industrial Average Index.
where S M = (s 2 1 , . . . , s 2 M ) and we have used the fact that M m=1 S D 2 m,T = σ 2 T → σ 2 . Now a componentwise Taylor expansion gives  (17) and (16) we have