Researcher profile

Nakahiro Yoshida

Nakahiro Yoshida contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
19works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2026arXiv

Statistical inference for highly correlated stationary point processes and noisy bivariate Neyman-Scott processes

Motivated by estimating the lead-lag relationships in high-frequency financial data, we propose noisy bivariate Neyman-Scott point processes with gamma kernels (NBNSP-G). NBNSP-G tolerates noises that are not necessarily Poissonian and has an intuitive interpretation. Our experiments suggest that NBNSP-G can explain the correlation of orders of two stocks well. A composite-type quasi-likelihood is employed to estimate the parameters of the model. However, when one tries to prove consistency and asymptotic normality, NBNSP-G breaks the boundedness assumption on the moment density functions commonly assumed in the literature. Therefore, under more relaxed conditions, we show consistency and asymptotic normality for bivariate point process models, which include NBNSP-G. Our numerical simulations also show that the estimator is indeed likely to converge.

preprint2022arXiv

Order estimate of functionals related to fractional Brownian motion and asymptotic expansion of the quadratic variation of fractional stochastic differential equation

We derive an asymptotic expansion for the quadratic variation of a stochastic process satisfying a stochastic differential equation driven by a fractional Brownian motion, based on the theory of asymptotic expansion of Skorohod integrals converging to a mixed normal limit. In order to apply the general theory, it is necessary to estimate functionals that are a randomly weighted sum of products of multiple integrals of the fractional Brownian motion, in expanding the quadratic variation and identifying the limit random symbols. To overcome the difficulty, we introduce two types of exponents by means of the "weighted graphs" capturing the structure of the sum in the functional, and investigate how the exponents change by the action of the Malliavin derivative and its projection.

preprint2021arXiv

Global jump filters and quasi-likelihood analysis for volatility

We propose a new estimation scheme for estimation of the volatility parameters of a semimartingale with jumps based on a jump-detection filter. Our filter uses all of data to analyze the relative size of increments and to discriminate jumps more precisely. We construct quasi-maximum likelihood estimators and quasi-Bayesian estimators, and show limit theorems for them including $L^p$-estimates of the error and asymptotic mixed normality based on the framework of the quasi-likelihood analysis. The global jump filters do not need a restrictive condition for the distribution of the small jumps. By numerical simulation we show that our "global" method obtains better estimates of the volatility parameter than the previous "local" methods.

preprint2021arXiv

Global jump filters and realized volatility

For a semimartingale with jumps, we propose a new estimation method for integrated volatility, i.e., the quadratic variation of the continuous martingale part, based on the global jump filter proposed by Inatsugu and Yoshida [8]. To decide whether each increment of the process has jumps, the global jump filter adopts the upper $α$-quantile of the absolute increments as the threshold. This jump filter is called global since it uses all the observations to classify one increment. We give a rate of convergence and prove asymptotic mixed normality of the global realized volatility and its variant "Winsorized global volatility". By simulation studies, we show that our estimators outperform previous realized volatility estimators that use a few adjacent increments to mitigate the effects of jumps.

preprint2021arXiv

Simplified quasi-likelihood analysis for a locally asymptotically quadratic random field

The asymptotic decision theory by Le Cam and Hajek has been given a lucid perspective by the Ibragimov-Hasminskii theory on convergence of the likelihood random field. Their scheme has been applied to stochastic processes by Kutoyants, and today this plot is called the IHK program. This scheme ensures that asymptotic properties of an estimator follow directly from the convergence of the random field if a large deviation estimate exists. The quasi-likelihood analysis (QLA) proved a polynomial type large deviation (PLD) inequality to go through a bottleneck of the program. A conclusion of the QLA is that if the quasi-likelihood random field is asymptotically quadratic and if a key index reflecting identifiability the random field has is non-degenerate, then the PLD inequality is always valid, and as a result, the IHK program can run. Many studies already took advantage of the QLA theory. However, not a few of them are using it in an inefficient way yet. The aim of this paper is to provide a reformed and simplified version of the QLA and to improve accessibility to the theory. As an example of the effects of the theory based on the PLD, the user can obtain asymptotic properties of the quasi-Bayesian estimator by only verifying non-degeneracy of the key index.

preprint2020arXiv

Adaptive and non-adaptive estimation for degenerate diffusion processes

We discuss parametric estimation of a degenerate diffusion system from time-discrete observations. The first component of the degenerate diffusion system has a parameter $θ_1$ in a non-degenerate diffusion coefficient and a parameter $θ_2$ in the drift term. The second component has a drift term parameterized by $θ_3$ and no diffusion term. Asymptotic normality is proved in three different situations for an adaptive estimator for $θ_3$ with some initial estimators for ($θ_1$ , $θ_2$), an adaptive one-step estimator for ($θ_1$ , $θ_2$ , $θ_3$) with some initial estimators for them, and a joint quasi-maximum likelihood estimator for ($θ_1$ , $θ_2$ , $θ_3$) without any initial estimator. Our estimators incorporate information of the increments of both components. Thanks to this construction, the asymptotic variance of the estimators for $θ_1$ is smaller than the standard one based only on the first component. The convergence of the estimators for $θ_3$ is much faster than the other parameters. The resulting asymptotic variance is smaller than that of an estimator only using the increments of the second component.

preprint2020arXiv

Asymptotic expansion of a variation with anticipative weights

Asymptotic expansion of a variation with anticipative weights is derived by the theory of asymptotic expansion for Skorohod integrals having a mixed normal limit. The expansion formula is expressed with the quasi-torsion, quasi-tangent and other random symbols. To specify these random symbols, it is necessary to classify the level of the effect of each term appearing in the stochastic expansion of the variable in question. To solve this problem, we consider a class ${\cal L}$ of certain sequences $({\cal I}_n)_{n\in{\mathbb N}}$ of Wiener functionals and we give a systematic way of estimation of the order of $({\cal I}_n)_{n\in{\mathbb N}}$. Based on this method, we introduce a notion of exponent of the sequence $({\cal I}_n)_{n\in{\mathbb N}}$, and investigate the stability and contraction effect of the operators $D_{u_n}$ and $D$ on ${\cal L}$, where $u_n$ is the integrand of a Skorohod integral. After constructed these machineries, we derive asymptotic expansion of the variation having anticipative weights. An application to robust volatility estimation is mentioned.

preprint2020arXiv

Marked point processes and intensity ratios for limit order book modeling

This paper extends the analysis of Muni Toke and Yoshida (2020) to the case of marked point processes. We consider multiple marked point processes with intensities defined by three multiplicative components, namely a common baseline intensity, a state-dependent component specific to each process, and a state-dependent component specific to each mark within each process. We show that for specific mark distributions, this model is a combination of the ratio models defined in Muni Toke and Yoshida (2020). We prove convergence results for the quasi-maximum and quasi-Bayesian likelihood estimators of this model and provide numerical illustrations of the asymptotic variances. We use these ratio processes in order to model transactions occuring in a limit order book. Model flexibility allows us to investigate both state-dependency (emphasizing the role of imbalance and spread as significant signals) and clustering. Calibration, model selection and prediction results are reported for high-frequency trading data on multiple stocks traded on Euronext Paris. We show that the marked ratio model outperforms other intensity-based methods (such as "pure" Hawkes-based methods) in predicting the sign and aggressiveness of market orders on financial markets.

preprint2020arXiv

Rate of Estimation for the Stationary Distribution of Stochastic Damping Hamiltonian Systems with Continuous Observations

We study the problem of the non-parametric estimation for the density $π$ of the stationary distribution of a stochastic two-dimensional damping Hamiltonian system $(Z_t)_{t\in[0,T]}=(X_t,Y_t)_{t \in [0,T]}$. From the continuous observation of the sampling path on $[0,T]$, we study the rate of estimation for $π(x_0,y_0)$ as $T \to \infty$. We show that kernel based estimators can achieve the rate $T^{-v}$ for some explicit exponent $v \in (0,1/2)$. One finding is that the rate of estimation depends on the smoothness of $π$ and is completely different with the rate appearing in the standard i.i.d.\ setting or in the case of two-dimensional non degenerate diffusion processes. Especially, this rate depends also on $y_0$. Moreover, we obtain a minimax lower bound on the $L^2$-risk for pointwise estimation, with the same rate $T^{-v}$, up to $\log(T)$ terms.

preprint2013arXiv

Edgeworth expansion for functionals of continuous diffusion processes

This paper presents new results on the Edgeworth expansion for high frequency functionals of continuous diffusion processes. We derive asymptotic expansions for weighted functionals of the Brownian motion and apply them to provide the second order Edgeworth expansion for power variation of diffusion processes. Our methodology relies on martingale embedding, Malliavin calculus and stable central limit theorems for semimartingales. Finally, we demonstrate the density expansion for studentized statistics of power variations.

preprint2013arXiv

Martingale Expansion in Mixed Normal Limit

The quasi-likelihood estimator and the Bayesian type estimator of the volatility parameter are in general asymptotically mixed normal. In case the limit is normal, the asymptotic expansion was derived in Yoshida (1997) as an application of the martingale expansion. The expansion for the asymptotically mixed normal distribution is then indispensable to develop the higher-order approximation and inference for the volatility. The classical approaches in limit theorems, where the limit is a process with independent increments or a simple mixture, do not work. We present asymptotic expansion of a martingale with asymptotically mixed normal distribution. The expansion formula is expressed by the adjoint of a random symbol with coefficients described by the Malliavin calculus, differently from the standard invariance principle. Applications to a quadratic form of a diffusion process ("realized volatility") are discussed.

preprint2012arXiv

Asymptotic expansion for the quadratic form of the diffusion process

In [8], asymptotic expansion of the martingale with mixed normal limit was provided. The expansion formula is expressed by the adjoint of a random symbol with coefficients described by the Malliavin calculus, differently from the standard invariance principle. As an application, an asymptotic expansion for a quadratic form of a diffusion process was derived in the same paper. This article gives some details of the derivation, after a short review of the martingale expansion in mixed normal limit.

preprint2012arXiv

Nondegeneracy of Random Field and Estimation of Diffusion

We construct a quasi likelihood analysis for diffusions under the high-frequency sampling over a finite time interval. For this, we prove a polynomial type large deviation inequality for the quasi likelihood random field. Then it becomes crucial to prove nondegeneracy of a key index chi_0. By nature of the sampling setting, chi_0 is random. This makes it difficult to apply a naive sufficient condition, and requires a new machinery. In order to establish a quasi likelihood analysis, we need quantitative estimate of the nondegeneracy of chi_0. The existence of a nondegenerate local section of a certain tensor bundle associated with the statistical random field solves this problem.

preprint2012arXiv

Quasi-Likelihood Analysis for Stochastic Regression Models with Nonsynchronous Observations

We consider nonsynchronous sampling of parameterized stochastic regression models, which contain stochastic differential equations. Constructing a quasi-likelihood function, we prove that the quasi-maximum likelihood estimator and the Bayes type estimator are consistent and asymptotically mixed normal when the sampling frequency of the nonsynchronous data becomes large.

preprint2010arXiv

Second-order asymptotic expansion for a non-synchronous covariation estimator

In this paper, we consider the problem of estimating the covariation of two diffusion processes when observations are subject to non-synchronicity. Building on recent papers \cite{Hay-Yos03, Hay-Yos04}, we derive second-order asymptotic expansions for the distribution of the Hayashi-Yoshida estimator in a fairly general setup including random sampling schemes and non-anticipative random drifts. The key steps leading to our results are a second-order decomposition of the estimator's distribution in the Gaussian set-up, a stochastic decomposition of the estimator itself and an accurate evaluation of the Malliavin covariance. To give a concrete example, we compute the constants involved in the resulting expansions for the particular case of sampling scheme generated by two independent Poisson processes.

preprint2007arXiv

Moment estimation for ergodic diffusion processes

We investigate the moment estimation for an ergodic diffusion process with unknown trend coefficient. We consider nonparametric and parametric estimation. In each case, we present a lower bound for the risk and then construct an asymptotically efficient estimator of the moment type functional or of a parameter which has a one-to-one correspondence to such a functional. Next, we clarify a higher order property of the moment type estimator by the Edgeworth expansion of the distribution function.

preprint2007arXiv

Parametric estimation for partially hidden diffusion processes sampled at discrete times

For a one dimensional diffusion process $X=\{X(t) ; 0\leq t \leq T \}$, we suppose that $X(t)$ is hidden if it is below some fixed and known threshold $τ$, but otherwise it is visible. This means a partially hidden diffusion process. The problem treated in this paper is to estimate finite dimensional parameter in both drift and diffusion coefficients under a partially hidden diffusion process obtained by a discrete sampling scheme. It is assumed that the sampling occurs at regularly spaced time intervals of length $h_n$ such that $n h_n=T$. The asymptotic is when $h_n\to0$, $T\to\infty$ and $n h_n^2\to 0$ as $n\to\infty$. Consistency and asymptotic normality for estimators of parameters in both drift and diffusion coefficients are proved.