Asymptotic distributions for estimated expected functionals of general random elements

Abstract: We consider an estimation problem of expected functionals of a general random element that values in a metric space. If the functional forms an explicit function of some unknown parameters, we can estimate it by plugging-in a suitable estimator in to the function, and we can find the asymptotic distribution. However, if the functional is implicit in the parameters, it causes a problem of specifying asymptotic distribution. This paper gives a general condition to specify the asymptotic distribution even if the functional is implicit in the parameters, and further investigates it in detail when the random elements are semimartingales.

\section{Introduction}

Let $(\Omega,\mathcal{F},{\mathbb{P}})$ be a probability space, and ${\cal X}$ be a metric space with norm $\|\cdot \|$.
Consider a ${\cal X}$-valued random element $X^{\vartheta}$ with an unknown parameter $\vartheta\in \Theta\subset \mathbb{R}^p$, and the distribution of $X^\vartheta$ is $P_{\vartheta}:={\mathbb{P}}\circ (X^\vartheta)^{-1}$.
Suppose that there exists the true value $\vartheta_0\in \Theta$, and we are interested in inference for the following expected functional of $X^\vartheta$:
\begin{align*}
H(\vartheta_0) = \mathbb{E}\left[h(X^{\vartheta_0},\vartheta_0)\right] = \int_{{\cal X}} h(x,\vartheta_0)\,P_{\vartheta_0}(\mathrm{d} x),
\end{align*}
where $h: {\cal X}\times \Theta \to \mathbb{R}$.

Such a expected functional appears in many statistical problem,
where ${\cal X}$ is not only an Euclidean space, but also a functional space.
When $H(\vartheta)$ is written in explicit function of $\vartheta$, we can estimate it as a {\it plug-in estimator} $H(\wh{\vartheta})$ with a suitable estimator of $\vartheta_0$ plugged-in, and it will be easy to evaluate the statistical error by, e.g., the {\it delta method}: for example, if $H$ is differentiable and $\wh{\vartheta}$ is asymptotically normal with asymptotic variance $\sigma^2_0$, then $H(\wh{\vartheta})$ is also asymptotically normal with asymptotic  variance $|\nabla_\vartheta H(\vartheta_0)|^2\sigma^2_0$; see, e.g., Corollary \ref{cor:mc-est}, below.
However, it is not practicable when $H(\vartheta)$ is implicit in $\vartheta$, which is our interest in this paper.

For example, consider a case where $X^\vartheta=(X^\vartheta_t)_{t\in [0,T]}$ is a diffusion process, which is a $C([0,T])$-valued random element. Estimating problem for
\begin{align}
H(\vartheta) = \mathbb{E}\left[e^{-r T}\max (X^{\vartheta}_T  - K,0)\right],\quad r,\, K>0  \label{call-option}
\end{align}
will appear in a typical financial problem for the pricing an {\it European call option} with strike price $K$ and interest rate $r$ when $X^\vartheta$ is a stock price.
This functional $H$ is generally implicit in $\vartheta$ except for some special models of $X^\vartheta$. In such a case, Monte Carlo simulation will be used: based on an estimated value of $\vartheta_0$, say $\wh{\vartheta}$, generate many paths from the estimated distribution $P_{\wh{\vartheta}}$ by simulations, say $X^{\wh{\vartheta}}(i)\ (i=1,2,\dots, B)$, and compute the average $B^{-1}\sum_{i=1}^B e^{-r T}\max (X^{\wh{\vartheta}}_T(i)  - K,0)$, which goes to, not $H(\vartheta_0)$, but $H(\wh{\vartheta})$ as $B\to \infty$, that remains a statistical error. To get information of $H(\vartheta_0)$, we need to know the asymptotic distribution of $H(\wh{\vartheta})$, but such a statistical error seems often ignored in practice since the asymptotic variance of $H(\wh{\vartheta})$ is not clear.

In this paper, we investigate the specification of asymptotic distribution of $H(\wh{\vartheta})$ even in the case where $H(\vartheta)$ is possibly implicit in $\vartheta$ (so the delta method is not available explicitly), when $X^\vartheta$ is a random element that  values in a general metric space $({\cal X},\|\cdot\|)$. More precisely, we will find the asymptotic distribution of
\begin{align}
\gamma_n^{-1}(H(\wh{\vartheta}) - H(\vartheta_0)) \label{eq:target}
\end{align}
when $\gamma_n^{-1}(\wh{\vartheta} - \vartheta_0) \stackrel{d}{\longrightarrow} Z$ as $n\to \infty$ for some random variable $Z$ and norming sequence $\gamma_n\to 0\ (n\to \infty)$.
This problem is a very fundamental problem in statistics, and
it is well known, in the case where $H$ is explicit and differentiable,
that the asymptotic distribution is found by the {\it delta method}:
\[
\gamma_n^{-1}(H(\wh{\vartheta}) - H(\vartheta_0)) \stackrel{d}{\longrightarrow} \nabla_\vartheta H(\vartheta_0) Z.
\]
However, it seems not discussed sufficiently in a statistical context when $H(\vartheta)$ is implicit, and especially when $X^\vartheta$ is a stochastic process.
Formally speaking, we need the following derivative:
\[
\nabla_\vartheta H(\vartheta) = \nabla_\vartheta \mathbb{E}[h(X^\vartheta,\vartheta)] = \mathbb{E}\left[\nabla_xh(X^\vartheta,\vartheta)\nabla_\vartheta X^\vartheta + \dot{h}(X^\vartheta,\vartheta)\right],
\]
where $\dot{h}(x,\vartheta) = \nabla_\vartheta h(x,\vartheta)$, if $\nabla_\vartheta$ and $\mathbb{E}$ are exchangable.
However, the last expectation has a bit problem sisnce the sense of the derivative ``$\nabla_\vartheta X^\vartheta$" is still not clear.

Such a ``derivative" has been considered in sensitivity analysis for expected functionals.
For example, when $X^\vartheta$ is a Euclidian valued random variable with a parameter $\vartheta$ in the distribution, Suri \cite{s83} discusses an expression of a derivative $\nabla_\vartheta X^\vartheta$ based on the distribution function.
Moreover, when $X^\vartheta$ is a stochastic process with the initial value $X^\vartheta_0=\vartheta$,
the map $x\mapsto X^x_t(\omega)$ is called a {\it stochastic flow}, and the continuity and the differentiability with respect to $\vartheta$ can be discussed; see., e.g., Protter \cite{p05}, Chapter V.7 in the case where $X^x$ satisfies a stochastic differential equation.

A similar problem appears in financial computation of {\it Greeks}, which are obtained as derivatives of option prices with respect to some specific parameters since those prices are written by expected functionals of underlying stock prices as in, e.g., \eqref{call-option}.
This problem is recently well studied via {\it Malliavin calculus}, which has a powerful tool such as ``integration-by-parts" formula to compute such derivatives; see, e.g., Davis and Johansson \cite{dj06}, Fourni\'e {\it et al.} \cite{fetal99, fetal01},  Gobet and Kohatsu-Higa \cite{gk03}, among others. see also Kohatsu-Higa and Montero \cite{km04} as a good guidance.  t
On the other hand, Chen and Glasserman \cite{cg07},  Glasserman and Liu \cite{gl10} take the {\it path-wise derivative approach} to compute $\mathbb{E}\left[\nabla_xh(X^\vartheta,\vartheta)\nabla_\vartheta X^\vartheta\right]$ by Monte Carlo simulations; see also Glasserman \cite{gl04}, Chapter 7.

In this paper, we will take a different approach. We consider a differentiability of $X^\vartheta$ with respect to $\vartheta$ in the $L^q$-sense for $q>0$ to evaluate errors in higher order terms; see the condition A4($q$), below. Under this approach,  the regularity conditions for \eqref{eq:target} can be an expectation-based and easy to check,

First, we will discuss general conditions to yield asymptotic distributions in both cases where ${\cal X}=\mathbb{R}^d$ and ${\cal X}$ is a functional space, and then we see the each case in detail.
The former is a standard situation, where the delta method is obtained as a special case; see Corollary \ref{cor:mc-est},  but the latter is performed with a kind of {\it derivative} of $X^\vartheta$ with respect to $\vartheta$; see Theorem \ref{thm:cond}. The case where $X^\vartheta$ are semimartingale that values in $\mathbb{D}$-space is important in applications. In this case, the asymptotic distribution can be described in terms of the {\it derivative process} of $X^\vartheta$ with respect to $\vartheta$ in $L^q$-sense.
The derivative process is essentially the same as the path-wise derivative discussed in  Chen and Glasserman \cite{cg07}, but we will give a different approach and an evaluation for not only a continuous diffusion processes, but also semimartingales with jumps, which is a new contribution on the derivative process because Glasserman and Liu \cite{gl10} just discuss from the simulations point of view.
Our investigation on jump processes indicates that an “error” in \eqref{eq:target} may get worse when $X^\vartheta$ values in $\mathbb{D}$-space than the case where $X^\vartheta$ values in $\mathbb{C}$-space; see Remark \ref{rem:unknown2}.

The paper is organized as follows. In Section \ref{sec:mc}, we shall state fundamental conditions to get the asymptotic distribution in a general formulation, and a special case where ${\cal X} = \mathbb{R}^d$ is described there as a corollary of the general statement as well as the usual delta method.
In later sections, we will consider more specific cases. In Section \ref{sec:C-space}, we consider the case where $X^\vartheta$ is functional valued, and a sufficient condition to ensure that the asymptotic normality of $H(\wh{\vartheta})$ is given in terms of the norm of the functional space ${\cal X}$.
We shall check the condition in each specific form of the functional.
Section \ref{sec:sde} is devoted to the case where $X^\vartheta$ is described by stochastic differential equations with jumps. The situation differs to a large extent when $X^\vartheta$ does not have a jump in the path (${\cal X}$ is a $\mathbb{C}$-space), compared to when $X^\vartheta$ does (${\cal X}$ is a $\mathbb{D}$-space).
The result indicates that we should be careful to use $H(\wh{\vartheta})$ since it may not be asymptotically normal in the case where $X^\vartheta$ is a jump process.

Throughout the paper, we use the following notation.
\begin{itemize}
\item $A\lesssim B$ means that there exists a universal constant $c>0$ such that $A \le c\cdot B$.
\item A $d$-dim Gaussian variable (distribution) with mean $0$ and variance-covariance matrix $\Sigma$ is denoted by $N_d(0,\Sigma)$. We omit the index $d=1$.
\item For a function $f: \mathbb{R}^d\to \mathbb{R}$ and $x=(x_1,\dots,x_d)\in \mathbb{R}^d$,
\[
\nabla_x f = \left(\frac{\partial f}{\partial x_1},\dots, \frac{\partial f}{\partial x_d} \right)^\top,
\]
and $\nabla_x^k = \nabla_x\otimes \nabla_x^{k-1},\ (k=2,3,\dots)$, constitutes a multilinear form.

\item For a function $f:\mathbb{R}^d\times \Theta \to \mathbb{R}$ and an integer $k$,
\begin{align*}
\dot{f}(x,\vartheta) = \nabla_\vartheta f(x,\vartheta); \quad f^{(k)}(x,\vartheta) = \nabla_x^k f(x,\vartheta).
\end{align*}
Note that $\nabla_x^k f$ is a $k$-th order tensor.

\item For a $k$-th order tensor $x=(x_{i_1,i_2,\dots,i_k})_{i_1,\dots,i_k = 1,\dot d} \in \mathbb{R}^d\otimes \dots \otimes \mathbb{R}^d$,
\[
|x| = \sqrt{\sum_{i_1=1}^d \dots\sum_{i_k=1}^d x^2_{i_1,i_2,\dots,i_k}}.
\]
\item For a ${\cal X}$-valued random element $X$, $\|X\|_{L^p}=\left(\mathbb{E}\|X\|^p\right)^{1/p}$ for $p>0$, where
$\|\cdot \|$ a norm on ${\cal X}$, and write $X\in L^p$ if $\|X\|_{L^p} < \infty$.
\end{itemize}

\section{Fundamental conditions for asymptotic distributions}\label{sec:mc}

\subsection{Basic results in general formulation}
Assume that a realization of $X^{\vartheta_0}$ from $P_{\vartheta_0}$, say $X^{\vartheta_0,n}$,  is given, where $n$ is supposed to be a parameter on which the sample size depends. For example, when we observe $n$-samples of i.i.d. variables $\{X_k\}_{k\in \mathbb{N}}$, it can be regarded as $X^{\vartheta_0, n} =(X_1,X_2,\dots,X_n)$, so $n$ represents the number of samples. When $X^{\vartheta_0}$ is a stochastic process $X=(X_t)_{t\ge 0}$, $X^{\vartheta_0,n}$ can be a time-continuous observation in a $[0,n]$-time interval: $X^{\vartheta_0,n}=(X_t)_{t\in [0,n]}$, or it can be discrete samples such as $X^{\vartheta_0,n}=(X_0,X_{t_1},\dots,X_{t_n})$, among others.
We assume that a ``good" estimator of $\vartheta_0$ is given based on the observations $X^{\vartheta_0,n}$, say
\[
\wh{\vartheta}_n:=\wh{\vartheta}(X^{\vartheta_0,n}).
\]
We assume that some estimator of $\vartheta_0$, say $\wh{\vartheta}_n$, is given in a suitable manner.
We shall investigate a fundamental condition under which $H(\wh{\vartheta}_n)$ is the asymptotic distribution is specified.

We make the following conditions.
\begin{as} \label{as:h-deriv}
For any $\vartheta'\in \Theta$,  $\displaystyle \nabla_\vartheta \mathbb{E}[h(X^{\vartheta'},\vartheta)] = \mathbb{E}[\dot{h}(X^{\vartheta'},\vartheta)]$.
\end{as}

\begin{as}\label{as:h-cont}
The function $\vartheta \mapsto \mathbb{E}[\dot{h}(X^\vartheta,\vartheta)]$ is continuous.
\end{as}

\begin{as}\label{as:Z}
There exists a diagonal matrix  $\Gamma_n=\mathrm{diag}(\gamma_n^{(1)},\dots,\gamma_n^{(p)})$ with $\gamma_n^{(k)} > 0$ and
$\gamma_{n*}:=\max_{1\le k\le p} \gamma_n^{(k)} \downarrow 0\ (n\to \infty)$ such that the estimator $\wh{\vartheta}_n$ satisfies
\[
\Gamma_n^{-1}(\wh{\vartheta}_n - \vartheta_0) \stackrel{d}{\longrightarrow} Z;\qquad \gamma_{n*}^{-1}(\wh{\vartheta}_n - \vartheta_0) \stackrel{d}{\longrightarrow} Z^*,
\]
as $n\to \infty$, for $p$-dim random variables $Z$ and $Z^*$.
\end{as}

\begin{flushleft}
{\bf A4($\mb{q}$).}\ \ There exists  a ${\cal X}^p$-valued random element $Y^\vartheta$ such that
$Y^\vartheta \in L^q$ for $q>0$ and
\begin{align*}
\|X^{\vartheta + u} - X^\vartheta - u^\top Y^\vartheta\|_{L^q} =o(|u|),\quad |u|\to 0,
\end{align*}
uniformly in $\vartheta\in \Theta$.
\end{flushleft}

\begin{remark}
As for condition A\ref{as:Z}, it usually holds that $\gamma_n^{(k)}=1/\sqrt{n}$ for all $k$ in i.i.d.-cases, but there are some examples where the rates of convergence are different among parameters, e.g., for a sequence such that $T_n/n\to 0$ as $n\to \infty$ and constants $\sigma_1^2,\sigma_2^2\ne 0$,
\[
\mathrm{diag}(\sqrt{T_n},\sqrt{n})(\wh{\vartheta}_n^{(1)} - \vartheta_0^{(1)}, \wh{\vartheta}_n^{(2)} - \vartheta_0^{(2)})^\top \stackrel{d}{\longrightarrow} N_2\left(0,\mathrm{diag}(\sigma_1^2,\sigma_2^2)\right)=Z.
\]
In such a case, $\gamma_{n*}=1/\sqrt{T_n}$ and we have
\[
\sqrt{T_n}(\wh{\vartheta}_n^{(1)} - \vartheta_0^{(1)}, \wh{\vartheta}_n^{(2)} - \vartheta_0^{(2)})^\top \stackrel{d}{\longrightarrow} (N(0,\sigma_1^2),0)^\top = Z^*,
\]
which is a degenerate random variable; see also Examples \ref{ex:levy} and \ref{ex:ou2}.
\end{remark}

\begin{remark}
In condition A4($q$), the random element $Y^\vartheta$ is interpreted as the first derivative of $X^\vartheta$ with respect to $\vartheta$ in the sense of $L^q$.
\end{remark}

Although the following seems to be a simple result, we shall claim it here since it is a basis of the discussion below.
\begin{thm}\label{thm:mc-est}
Suppose that A\ref{as:h-deriv} -- A\ref{as:Z} hold true, and that there exists a constant vector $C_{\vartheta_0} \in \mathbb{R}^p$ such that
for $\gamma_{n*}:=\max_{1\le k\le p}\gamma_n^{(k)}$,
\begin{align}
\gamma_{n*}^{-1} \mathbb{E}\left[h(X^{\vartheta},\vartheta_0) - h(X^{\vartheta_0},\vartheta_0)\right]\Big|_{\vartheta = \wh{\vartheta}_n} = C_{\vartheta_0}^\top \gamma_{n*}^{-1}(\wh{\vartheta}_n - \vartheta_0) + o_p(1),\label{term2}
\end{align}
as $n\to \infty$. Then, it holds that
\begin{align*}
\gamma_{n*}^{-1} [H(\wh{\vartheta}_n) - H(\vartheta_0)] \stackrel{d}{\longrightarrow} (\mathbb{E}[\dot{h}(X^{\vartheta_0},\vartheta_0)] + C_{\vartheta_0})^\top \,Z^*, \quad n\to \infty.
\end{align*}
\end{thm}

\begin{proof}
Let $X_*^{\vartheta_0}\sim P_{\vartheta_0}$, which is independent of the data $X^{\vartheta_0,n}$. Then we have that
\begin{align*}
H(&\wh{\vartheta}_n) - H(\vartheta_0) \notag \\
&= \mathbb{E}\left[h(X_*^{\wh{\vartheta}_n},\wh{\vartheta}_n) - h(X_*^{\vartheta_0},\vartheta_0) \big| X^{\vartheta_0,n}\right]  \notag \\
&= \mathbb{E}\left[h(X_*^{\wh{\vartheta}_n},\wh{\vartheta}_n) - h(X_*^{\wh{\vartheta}_n},\vartheta_0) \big| X^{\vartheta_0,n}\right]
+ \mathbb{E}\left[h(X_*^{\wh{\vartheta}_n},\vartheta_0) - h(X_*^{\vartheta_0},\vartheta_0) \big| X^{\vartheta_0,n}\right] \notag \\
&= \mathbb{E}\left[\dot{h}(X_*^{\vartheta},\vartheta_0 + \eta_n(\wh{\vartheta} - \vartheta_0)\right]\Big|_{\vartheta = \wh{\vartheta}_n} (\wh{\vartheta}_n - \vartheta_0)
+ \mathbb{E}\left[h(X^{\vartheta},\vartheta_0) - h(X^{\vartheta_0},\vartheta_0)\right]\Big|_{\vartheta = \wh{\vartheta}_n},
\end{align*}
where $\eta_n$ is a random variable values in $[0,1]$.
We use the mean value theorem in the last equality.
Then, under A\ref{as:h-cont}, the continuous mapping theorem yields the result.
\end{proof}

This theorem immediately leads us a version of the {\it delta method} when ${\cal X} = \mathbb{R}^d$.
\begin{cor}\label{cor:mc-est}
Consider the case where ${\cal X} = \mathbb{R}^d$ and $X^\vartheta$ be a random variable with probability density $f:{\cal X}\times \Theta \to \mathbb{R}$ in Theorem \ref{thm:mc-est}. Suppose that $f$ is twice differentiable with respect to $\vartheta\in \Theta$ with $\int_{\cal X} h(x,\vartheta_0)\dot{f}(x,\vartheta_0)\,\mathrm{d} x < \infty$. Moreover, suppose A\ref{as:Z} holds, and that it holds for the second derivative of $f$ in $\vartheta$, say $\ddot{f}$, such that
\begin{align}
\sup_{\vartheta\in \Theta} \left|\int_{\cal X} h(x,\vartheta_0)\ddot{f}(x,\vartheta)\,\mathrm{d} x\right| < \infty. \label{f-dd}
\end{align}
Then $C_\vartheta$ in \eqref{term2} is given by $C_\vartheta=\int_\mathbb{R} h(x,\vartheta)\dot{f}(x,\vartheta)\,\mathrm{d} x$, and it follows that
\[
\gamma_{n*}^{-1} [H(\wh{\vartheta}_n) - H(\vartheta_0)] \stackrel{d}{\longrightarrow}  \nabla_\vartheta H(\vartheta_0)^\top \,Z^*, \quad n\to \infty,
\]
where
\[
\nabla_\vartheta H(\vartheta_0) = \int_{\cal X} \left[\dot{h}(x,\vartheta_0)f(x,\vartheta_0)+ h(x,\vartheta_0)\dot{f}(x,\vartheta_0)\right]\,\mathrm{d} x.
\]
\end{cor}

\begin{proof}
For $\vartheta \in \Theta$ and $u\in \mathbb{R}^p$ with $\vartheta + u\in \Theta$, it follows from Taylor's formula that
\begin{align*}
\mathbb{E}\left[h(X^{\vartheta + u},\vartheta) - h(X^\vartheta,\vartheta)\right] &= \int_{\cal X} h(x,\vartheta)\left[f(x,\vartheta + u) - f(x,\vartheta)\right]\,\mathrm{d} x \\
&= \int_{\cal X} h(x,\vartheta)\left[u^\top \dot{f}(x,\vartheta) + u^\top \ddot{f}(x,\vartheta_u) u\right]\,\mathrm{d} x
\end{align*}
where $\vartheta^u:= \vartheta + \eta_u u$ for some $\eta_u\in [0,1]$.
Hence, when $\vartheta=\vartheta_0$ and $u=\wh{\vartheta}_n - \vartheta_0$ and both sides are multiplied by $\gamma_{n*}^{-1}$, we have that
\begin{align*}
\gamma_{n*}^{-1}\mathbb{E}\left[h(X^{\wh{\vartheta}_n},\vartheta_0) - h(X^{\vartheta_0},\vartheta_0)\right]
&= \left(\int_{\cal X} h(x,\vartheta_0)\dot{f}(x,\vartheta_0)\,\mathrm{d} x\right)^\top \cdot \gamma_{n*}^{-1}(\wh{\vartheta}_n - \vartheta_0) \\
&\quad + O_p(|\wh{\vartheta}_n - \vartheta_0|),\quad n\to \infty.
\end{align*}
Therefore, $C_{\vartheta_0} = \int_{\cal X} h(x,\vartheta_0)\dot{f}(x,\vartheta_0)\,\mathrm{d} x$, and Theorem \ref{thm:mc-est} and the condition A\ref{as:h-deriv} yield that
\begin{align*}
\gamma_{n*}^{-1} [H(\wh{\vartheta}_n) - H(\vartheta_0)] &\stackrel{d}{\longrightarrow}  \left( \int_{\cal X} \left[\dot{h}(x,\vartheta_0)f(x,\vartheta_0)+ h(x,\vartheta_0)\dot{f}(x,\vartheta_0)\right]\,\mathrm{d} x\right)^\top Z^*\\
&= \nabla_\vartheta H(\vartheta_0)^\top \,Z^*,\quad n\to \infty.
\end{align*}
\end{proof}

When ${\cal X}$ is not Euclidean, but some functional spaces, the following theorem will be useful to specify the value of $C_{\vartheta_0}$ in Theorem \ref{thm:mc-est}.

\begin{thm}\label{thm:cond}
Suppose that assumptions A\ref{as:Z} and A4($q$) hold for a constant $q>1$, and that there exists a $\mathbb{R}^p$-valued random variable $G_{\vartheta_0} \in L^1$ such that for each $u\in \mathbb{R}$ with $\vartheta_0 + u\in \Theta$,
\begin{align}
\left|\mathbb{E}[h(X^{\vartheta_0 + u},\vartheta_0) - h(X^{\vartheta_0},\vartheta_0)  -  u^\top G_{\vartheta_0}]\right| \lesssim \|X^{\vartheta_0 + u} - X^{\vartheta_0} - u^\top Y^{\vartheta_0}\|_{L^q}  +  r_u,  \label{mvt}
\end{align}
where $r_u=o(|u|)$ as $|u|\to 0$. Then the equality \eqref{term2} holds true with $C_{\vartheta_0}=\mathbb{E}[G_{\vartheta_0}]$.
\end{thm}

\begin{proof}
The assumption A4($q$) with $q>1$ implies that
\begin{align*}
\left|\mathbb{E}[h(X^{\vartheta_0 + u},\vartheta_0) - h(X^{\vartheta_0},\vartheta_0)  -  u^\top G_{\vartheta_0}]\right| = o(|u|), \quad |u|\to 0.
\end{align*}
Then, it follows that
\[
\mathbb{E}[h(X^{\vartheta_0 + u},\vartheta_0) - h(X^{\vartheta_0},\vartheta_0) ] =  \mathbb{E}[G_{\vartheta_0}]^\top u + o(|u|), \quad |u|\to 0.
\]
When $u = \wh{\vartheta}_n - \vartheta_0$ and both sides are multiplied by $\gamma_{n*}^{-1}$,
we obtain
\begin{align*}
\gamma_{n*}^{-1}\mathbb{E}\left[h(X^{\vartheta},\vartheta_0) - h(X^{\vartheta},\vartheta_0)\right] \Big|_{\vartheta=\wh{\vartheta}_n}
&= \mathbb{E}[G_\vartheta]^\top  \gamma_{n*}^{-1}(\wh{\vartheta}_n - \vartheta_0) + o_p(|\gamma_{n*}^{-1}(\wh{\vartheta}_n - \vartheta_0)|).
\end{align*}
The last term converges to zero in probability under A\ref{as:Z}. This ends the proof.
\end{proof}

\begin{example}
Consider a random variable $X^\vartheta$ values on ${\cal X}=\mathbb{R}$ with distribution function $F_\vartheta:\mathbb{R} \to [0,1]$ and a parameter $\vartheta\in \Theta \subset \mathbb{R}$. Suppose that a positive density $\nabla_x F_\vartheta(x) = f(x,\vartheta)$ exists, and that $\nabla_\vartheta^2 F_\vartheta^{-1}$ is bounded for simplicity:
\begin{align}
D:=\sup_{x\in [0,1],\vartheta\in \Theta}|\nabla_\vartheta^2 F_\vartheta^{-1}(x)| < \infty. \label{bdd}
\end{align}
We shall consider a ``derivative $\nabla_\vartheta X^\vartheta$" after the idea by Suri \cite{s83}: we may set
\[
X_\vartheta = F^{-1}_\vartheta(U);\quad \nabla_\vartheta X^\vartheta :=\nabla_\vartheta F_\vartheta^{-1}(U)
\]
where $U$ is a uniform random variable on $[0,1]$ independent of $\vartheta$.
This leads us to $F_\vartheta(F^{-1}_\vartheta(U)) = U$. Differentiating the both sides, we have that $\nabla_\vartheta F(X^\vartheta) + f(X^\vartheta,\vartheta) \nabla_\vartheta X^\vartheta = 0$, and that
\[
\nabla_\vartheta X^\vartheta = -\frac{\nabla_\vartheta F_\vartheta(X^\vartheta)}{f(X^\vartheta, \vartheta)}.
\]
Under the assumption \eqref{bdd}, we can easily see by Taylor's formula that, for any $q>0$,
\[
\|X^{\vartheta + u} - X^\vartheta - u^\top\nabla_\vartheta X^\vartheta\|^q \le \frac{D}{2}|u|^2 = o(u),\quad |u|\to 0,
\]
which yields the condition A($q$) with $Y^\vartheta = \nabla_\vartheta X^\vartheta$.

Now we suppose that $\nabla_x h(X^{\vartheta_0},\vartheta_0) \in L^r$ for some $r>1$ with $1/q + 1/r =1$. Then we can see that
\[
G_{\vartheta_0} = \nabla_x h(X^{\vartheta_0},\vartheta_0)\nabla_\vartheta X^{\vartheta_0}.
\]
Indeed, we see by Taylor's formula that
\begin{align*}
\left| \mathbb{E} \left[h(X^{\vartheta_0 + u},\vartheta_0) - h(X^{\vartheta_0},\vartheta_0) - u G_{\vartheta_0} \right] \right|
&= \left| \mathbb{E}\left[\nabla_x h(X^{\vartheta_0},\vartheta_0) (X^{\vartheta_0+u} - X^{\vartheta_0}) - uG_{\vartheta_0}\right] \right| \\
&\le \mathbb{E}\left[|\nabla_x h(X^{\vartheta_0},\vartheta_0)| \cdot |X^{\vartheta_0+u} - X^{\vartheta_0} - Y^{\vartheta_0}|\right] \\
&\le \|\nabla_x h(X^{\vartheta_0},\vartheta_0)\|_{L^r} \|X^{\vartheta_0+u} - X^{\vartheta_0} - Y^{\vartheta_0}\|_{L^q},
\end{align*}
under the assumption A\ref{as:h-deriv}, which yields the inequality \eqref{mvt}.
Therefore it follows by the integration-by-parts that
\begin{align*}
C_{\vartheta_0} &= \mathbb{E}[G_{\vartheta_0}] = - \int_\mathbb{R} h(x,\vartheta_0) \frac{\nabla_\vartheta F_{\vartheta_0}(x)}{f(x,\vartheta_0)}\,f(x,\vartheta_0)\,\mathrm{d} x
 = \int_\mathbb{R} h(x,\vartheta_0)\dot{f}(x,\vartheta_0)\,\mathrm{d} x,
\end{align*}
which coincides with the expression of $C_{\vartheta_0}$ in Corollary \ref{cor:mc-est}.
\end{example}

\if
\begin{example}
Let $X^\vartheta$ be a real valued random variable with $\Theta \subset \mathbb{R}$, and assume that we have a set of i.i.d. samples
$X^{\vartheta_0,n} = (X_1,\dots,X_n)$ with probability density $f(x,\vartheta)$ with Fisher information
$I_\vartheta := - \int_\mathbb{R} \frac{\dot{f}^2(x,\vartheta)}{f(x,\vartheta)}\,\mathrm{d} x < \infty$ for any $\vartheta\in \Theta$.
Suppose the regularities are such that the maximum likelihood estimator $\wh{\vartheta}_n$ is asymptotically normal:
\[\UTF{00DF}%
\sqrt{n}(\wh{\vartheta}_n - \vartheta_0) \stackrel{d}{\longrightarrow} N(0,I_{\vartheta_0}^{-1}),\quad n\to \infty.
\]
Moreover, suppose that $h:\mathbb{R}\times \Theta \to \mathbb{R}$ is a function such that all of the conditions in Theorem \ref{thm:mc-est} are satisfied, and that $C_\vartheta=\int_\mathbb{R} h(x,\vartheta)\dot{f}(x,\vartheta)\,\mathrm{d} x$ is continuous in $\vartheta$.
Then, as $n\to \infty$,
\[
\sqrt{n}\left(\wh{H}^*(\wh{\vartheta}_n) - H(\vartheta_0)\right) \stackrel{d}{\longrightarrow} N(0,C_{\vartheta_0}^2 I_{\vartheta_0}^{-1}).
\]
Therefore, an $\alpha$-confidence interval for $H(\vartheta_0)$ via an estimator $H(\wh{\vartheta}_n)$ is given by
\[
\left[H(\wh{\vartheta}_n) - \frac{z_{\alpha/2}}{\sqrt{n}} C_{\wh{\vartheta}_n} I_{\wh{\vartheta}_n}^{-1/2},H(\wh{\vartheta}_n) + \frac{z_{\alpha/2}}{\sqrt{n}} C_{\wh{\vartheta}_n} I_{\wh{\vartheta}_n}^{-1/2}\right],
\]
where $z_\alpha$ is the upper $\alpha$ percentile of $N(0,1)$.
\end{example}

\begin{example}
Let us consider a more concrete example than the previous.
Let $X^{\vartheta} = X_{\alpha,\beta}$ be a random variable that follows a gamma distribution with probability density
\[
\gamma(x;\alpha,\beta) = \frac{\beta^\alpha}{\Gamma(\alpha)}x^{\alpha-1}e^{-\beta x}\mathbf{1}_{\{x>0\}},\quad \vartheta = (\alpha,\beta)\in \Theta.
\]
Assuming that $\alpha>1$, we shall estimate $H(\vartheta) = \mathbb{E}[h(X)]$ for a smooth function $h\in C^\infty(\mathbb{R})$ by Monte Carlo simulation based on $n$-i.i.d. samples: $X_1,\dots,X_n$.
We suppose that for the $n$-th order Taylor expansion of $h$ around $x=0$,
say $h_n(x) := \sum_{k=1}^n \frac{h^{(k)(0)}}{k!}x^k$, satisfies
\[
\sup_{x\in \mathbb{R}} |h_n(x) - h(x)| \to 0\quad n\to \infty.
\]
The maximum likelihood estimator $\wh{\vartheta}_n$ is given by the solution to
$\sum_{i=1}^n \nabla_\vartheta \log \gamma(X_i;\wh{\vartheta}_n) = 0$, and it satisfies
\[
\sqrt{n}(\wh{\vartheta}_n - \vartheta) \stackrel{d}{\longrightarrow} N_2(0,I^{-1}(\vartheta)), \quad n\to \infty,
\]
where $I(\vartheta)$ is the Fisher information matrix, given by
\[
I(\vartheta) = \begin{pmatrix}\psi'(\alpha) & 1/\beta \\ 1/\beta & \alpha/\beta^2 \end{pmatrix},
\]
and $\psi$ is the digamma function with derivative $\psi'$.
By direct calculation, we have $\dot{\gamma}(x;\vartheta) = (\nabla_\alpha \gamma,\nabla_\beta \gamma)$:
\begin{align*}
\nabla_\alpha\gamma(x;\vartheta) &= \left(\frac{\log \beta}{\Gamma(\alpha)} - \psi(\alpha)\right)\gamma(x;\alpha,\beta) + \beta \gamma(x;\alpha-1,\beta) \\
\nabla_\beta\gamma(x;\vartheta) &= \frac{1}{\beta}\gamma(x;\alpha,\beta) - \frac{\alpha}{\beta}\gamma(x;\alpha+1,\beta).
\end{align*}
Therefore, it follows from the dominated convergence theorem that
\begin{align*}
C_\vartheta &= \int_0^\infty h(x)\dot{\gamma}(x;\vartheta)\,\mathrm{d} x  = \sum_{k=1}^\infty \frac{h^{(k)}(0)}{k!} \int_0^\infty x^k \dot{\gamma}(x;\vartheta)\,\mathrm{d} x\\
&= \sum_{k=1}^\infty \frac{h^{(k)}(0)}{k!}
\begin{pmatrix}
 \displaystyle \left(\frac{\log \beta}{\Gamma(\alpha)} - \psi(\alpha)\right)\mathbb{E}[X_{\alpha,\beta}^k] + \beta \mathbb{E}[X_{\alpha-1,\beta}^k]) \\
 \displaystyle \frac{1}{\beta}\mathbb{E}[X_{\alpha,\beta}^k]- \frac{\alpha}{\beta}\mathbb{E}[X_{\alpha+1,\beta}^k],
\end{pmatrix}
\end{align*}
which is computable by the formula
\[
\mathbb{E}[X_{\alpha,\beta}^k] = \frac{\alpha(\alpha+1)\dots(\alpha + k-1) }{\beta^k}.
\]
If $\Theta$ is bounded, then all the conditions in Corollary \ref{cor:mc-est} are satisfied, and hence
we have the asymptotic distribution of the estimator $H(\wh{\vartheta}_n)$:
\[
\sqrt{n}(H(\wh{\vartheta}_n) - H(\vartheta)) \stackrel{d}{\longrightarrow} N\left(0, C_\vartheta^\top I^{-1}(\vartheta)C_\vartheta\right),\quad n\to \infty.
\]
\end{example}
\f

\section{Expected functionals for stochastic processes}\label{sec:C-space}
In this section, we consider the case where ${\cal X}$ is a  functional space on a compact set $K\subset \mathbb{R}$,  e.g., $\mathbb{C}(K)$, $\mathbb{D}(K)$, with the sup norm
\[
\|x\| = \sup_{t\in K}|x_t|,\quad x=(x_t)_{t\in K}\in {\cal X}.
\]
Without loss of generality, we assume that $K=[0,1]$ for notational simplicity, so we consider the case where $X^\vartheta$ is a continuous time stochastic process on $[0,1]$.

\subsection{Functionals of expected integrals}
In this section, we are interested in the expected integral-type functionals
\[
H(\vartheta) = \mathbb{E}\left[\int_0^1 V_{\vartheta}(X_t^{\vartheta},t)\,\mathrm{d} t\right],
\]
for a function $V:\mathbb{R}^d\times [0,1]\to \mathbb{R}$. This is the case where $H(\vartheta) = \mathbb{E}[h(X^\vartheta,\vartheta)]$ with
\[
h(x,\vartheta) = \int_0^1V_\vartheta(x_t,t)\,\mathrm{d} t,\quad x\in {\cal X}
\]
The marginal distribution of a stochastic process $X^\vartheta$ is generally not explicit and the expectation
$\mathbb{E}[V_\vartheta(X_t^\vartheta,t)]$ is not clear. In such a case, Theorem \ref{thm:cond} can be useful to the analysis if the assumption A4($q$) can be confirmed.

\begin{example}\label{ex:diff}
Suppose that $X^\vartheta$ satisfies the following 1-dim stochastic differential equation:
\[
X^\vartheta_t = x(\vartheta) + \int_0^t a(X^\vartheta_s,\vartheta)\,\mathrm{d} s + \int_0^t b(X^\vartheta_s,\vartheta)\,\mathrm{d} W_s,
\]
where $W$ is a Wiener process and $a,b$ are functions with some ``good" regularities and $\vartheta\in \mathbb{R}^p$ is the unknown parameter. According to Section \ref{sec:sde}, under some regularities,
the {\it derivative process} $Y^\vartheta = (Y^\vartheta_t)_{t\in [0,1]}$ is given as follows.
\[
Y_t^\vartheta = \dot{x}(\vartheta) + \int_0^t A(X_s,Y_s,\vartheta)\,\mathrm{d} s + \int_0^t B(X_s,Y_s,\vartheta)\,\mathrm{d} W_s,
\]
where $\dot{x}(\vartheta) =\nabla_\vartheta x(\vartheta)$ and $A, B$ are $\mathbb{R}$-valued functions on $\mathbb{R}\times \mathbb{R}^p\times \mathbb{R}$, which are of the form
\[
A(x,y,\vartheta)= \nabla_x a(x,\vartheta) y + \dot{a}(x,\vartheta);\quad B(x,y,\vartheta)= \nabla_xb(x,\vartheta)y + \dot{b}(x,\vartheta).
\]
This $Y^\vartheta$ can satisfy
\[
\mathbb{E}\|X^{\vartheta+u} - X^\vartheta - u^\top Y^\vartheta \|^p \lesssim |u|^{2p}.
\]
for each $u\in \mathbb{R}^p$ and any $p\ge 2$, which implies A4($p$).
\end{example}

\begin{thm}\label{thm:int-V}
Suppose that there exists an integer $n\ge 1$ and $\vartheta\in \Theta$ such that
$V^{(n)}_\vartheta(x,t):=\nabla_x^n V_\vartheta(x,t)$ is Lipschitz continuous with respect to $x$, uniformly in $t\in [0,1]$:
\[
\sup_{t\in [0,1]}|V^{(n)}_\vartheta(x,t) - V^{(n)}_\vartheta(y,t)|\lesssim |x-y|, \quad x,y\in \mathbb{R}.
\]
Moreover, suppose that A4($q$) holds for some $q\ge 2n$, and that
\[
\sup_{t\in [0,1]}|V^{(k)}_\vartheta(X^{\vartheta}_t,t)| \in L^r,\quad k=1,\dots,n,
\]
for some $r>1$ with $1/r + 1/q = 1$.
Then, condition \eqref{mvt} holds with
\[
G_{\vartheta}=\int_0^1 V^{(1)}_\vartheta(X^\vartheta_t,t)Y_t^{\vartheta}\,\mathrm{d} t.
\]
\end{thm}

\begin{proof}
We shall check condition \eqref{mvt} in Theorem \ref{thm:cond}.
In the proof, for notational simplicity we consider only the case where $d=1$.
The general case can be shown in a similar manner.

Let
\[
R_t^{\vartheta,u}:= X_t^{\vartheta + u} - X_t^\vartheta - Y_t^\vartheta, \quad u\in \mathbb{R}^p.
\]
We note that $\|R^{\vartheta,u}\|_{L^q} \lesssim |u|^q$ by A4($q$).
It follows from Taylor's formula that
\begin{align*}
&\mathbb{E}\left[h(X^{\vartheta + u},\vartheta) - h(X^\vartheta,\vartheta) - u^\top \int_0^1 V^{(1)}_\vartheta(X^\vartheta_t,t)Y_t^{\vartheta}\,\mathrm{d} t \right] \\
&= E\left[\int_0^1 \left\{V_\vartheta(X_t^{\vartheta+u},t) - V_\vartheta(X_t^\vartheta,t) - u^\top V^{(1)}_\vartheta(X^\vartheta_t,t)Y_t^{\vartheta}\right\}\,\mathrm{d} t \right] \\
&= \mathbb{E}\left[\int_0^1 V^{(1)}_\vartheta(X_t^\vartheta,t) R_t^{\vartheta,u}\,\mathrm{d} t\right]
+\sum_{k=2}^{n-1}\frac{1}{k!}\mathbb{E}\left[\int_0^1 V^{(k)}_\vartheta(X^\vartheta,t)(R_t^{\vartheta,u} + u^\top Y_t^\vartheta)^k\,\mathrm{d} t\right] \\
&\qquad + \frac{1}{n!}\mathbb{E}\left[\int_0^1 V^{(n)}_\vartheta(\wt{X}_t^{\vartheta,u},t)(R_t^{\vartheta,u} + u^\top Y_t^\vartheta)^n\,\mathrm{d} t\right],
\end{align*}
where $\wt{X}_t^{u} = X^\vartheta + \eta_\vartheta^u (X^{\vartheta + u} - X^\vartheta)$ for some random number $\eta_\vartheta^u \in [0,1]$.

Firstly, it follows from H\"older's inequality that for $q,r>1$ with $1/q + 1/r=1$,
\[
\left|\mathbb{E}\left[\int_0^1 V^{(1)}_\vartheta(X_t^\vartheta,t) R_t^{\vartheta,u}\,\mathrm{d} t\right] \right|
\lesssim \left\|\sup_{t\in [0,1]}V ^{(1)}_\vartheta(X_t^\vartheta,t)\right\|_{L^r} \|R^{\vartheta,u}\|_{L^q}
\]
Secondly, noticing that
\begin{align*}
\|\wt{X}^{\vartheta,u} - X^\vartheta\|_{L^2} \le \|R^{\vartheta,u}\|_{L^2} + |u| \|Y^\vartheta\|_{L^2} = O(|u|),\quad |u|\to 0,
\end{align*}
we see that
\begin{align*}
&\left|\mathbb{E}\left[\int_0^1 V^{(n)}_\vartheta(\wt{X}_t^{u},t)(R_t^{\vartheta,u} + u^\top Y_t^\vartheta)^n\,\mathrm{d} t\right] - \mathbb{E}\left[\int_0^1 V^{(n)}_\vartheta(X_t^\vartheta,t)(R_t^{\vartheta,u} + u^\top Y_t^\vartheta)^n\,\mathrm{d} t\right]\right| \\
&\lesssim 2^{n-1}\int_0^1 \mathbb{E}\left[ |\wt{X}_t^{u} - X_t^\vartheta|\left\{|R_t^{\vartheta,u}|^n + |u|^n |Y_t^\vartheta|^n\right\}\right] \,\mathrm{d} t\\
&\lesssim \|\wt{X}^{\vartheta,u} - X^\vartheta\|_{L^2}\left\{ \| R_\vartheta^u\|_{L^{2n}}^{n} + |u|^n\| Y^\vartheta\|_{L^{2n}}^{n}\right\}
=o(|u|^n),\qquad \mbox{as $|u|\to 0$},
\end{align*}
Finally, from the Schwartz inequality, it is easy to see that for each $k=2,\dots,n$,
\begin{align*}
&\left|\mathbb{E}\left[\int_0^1 V^{(k)}_\vartheta(X_t^\vartheta,t)(R_t^{\vartheta,u} + u^\top Y_t^\vartheta)^k\right]\right| \\
&\le \int_0^1 \mathbb{E}\left[|V^{(k)}_\vartheta(X_t^\vartheta,t)|\cdot \|R^{\vartheta,u} + u^\top Y^\vartheta\|^k\right]\,\mathrm{d} t \\
&\le 2^{k-1}\int_0^1 \|V^{(k)}_\vartheta(X_t^\vartheta,t)\|_{L^s}\,\mathrm{d} t \cdot \left\{ \|R^{\vartheta,u}\|_{L^q}^k + |u|^k\|Y^\vartheta\|_{L^q}^k \right\} \\
&= o(|u|^k),
\end{align*}
where $s>1$ with $1/s + k/q = 1$. Note that such an $s>1$ exists under our assumption
since $(1 - k/q)^{-1} \ge n/(n-1)> 1$ when $q\ge 2n$. As a result, we have that
\begin{align*}
\left|\mathbb{E}\left[h(X^{\vartheta + u},\vartheta) - h(X^\vartheta,\vartheta) - u^\top \int_0^1 V^{(1)}_\vartheta(X^\vartheta_t,t)Y_t^{\vartheta}\,\mathrm{d} t \right] \right|
\lesssim \|R^{\vartheta,u}\|_{L^p} + o(|u|^2),
\end{align*}
which implies condition \eqref{mvt} in Theorem \ref{thm:cond} with
$G_\vartheta = \int_0^1 V^{(1)}_\vartheta(X^\vartheta_t,t)Y_t^{\vartheta}\,\mathrm{d} t$.  Therefore, the proof is completed.
\end{proof}

\begin{remark}
If the function $V_\vartheta$ is a ``good" function such that a ``lower" derivative is Lipschitz continuous,
then Theorem \ref{thm:int-V} requires only a ``small" $q\ge 2$ for A4$(q)$ to hold true.
The more “violent” the function $V$ is, the stronger the integrability condition becomes.
\end{remark}

\begin{example}
Consider a 1-dim (ergodic) diffusion process $X^\vartheta=(X_t)_{t\ge 0}$:  for a constant $x>0$,
\[
X^{\vartheta_0}_t = x + \int_0^t a(X^{\vartheta_0}_s,\vartheta_0)\,\mathrm{d} s + \int_0^t b(X^{\vartheta_0}_t)\,\mathrm{d} W_t,
\]
where $\vartheta_0\in \mathbb{R}$ is unknown, and consider the estimation of
\[
H(\vartheta_0) =  \int_0^T e^{-r t}U(X^{\vartheta_0}_t)\,\mathrm{d} t,
\]
for a constant  $r>0$ and a function $U\in C(\mathbb{R})$,
which is the case where $V_\vartheta(x,t) = e^{-r t}U(x)\mathbf{1}_{[0,T)}(t)$. See also Example \ref{ex:ou2} for practical applications of this example.

Assume that we have continuous data $\{X_t\}_{t\in [0,T]}$, and consider the {\it long term asymptotics}: $T\to \infty$. Then, under some regularities, the maximum likelihood estimator of $\vartheta$, say $\wh{\vartheta}_T$, satisfies
\[
\sqrt{T}(\wh{\vartheta} - \vartheta_0) \stackrel{d}{\longrightarrow} N(0,I^{-1}(\vartheta_0)), \quad T\to \infty,
\]
where $I(\vartheta) = \int_{\mathbb{R}} \frac{a^2(x,\vartheta)}{b^2(x)}\,\pi(\mathrm{d} x)$ for a stationary distribution $\pi$,
and it can be estimated by, e.g.,
\[
\wh{I}_T(\vartheta) = \frac{1}{T}\int_0^T \frac{a^2(X_t,\vartheta)}{b^2(X_t)}\,\mathrm{d} t \stackrel{p}{\longrightarrow} I(\vartheta), \quad T\to \infty,
\]
uniformly in $\vartheta\in \Theta$ (see, e.g., Kutoyants \cite{k04}).
Therefore, considering the derivative process $Y^\vartheta$ given in Example \ref{ex:diff}, we have that
\[
\sqrt{T}(H(\wh{\vartheta}_T) - H(\vartheta_0)) \stackrel{d}{\longrightarrow} N(0, C_{\vartheta_0}^2 I(\vartheta_0)^{-1}), \quad T\to \infty,
\]
where
\[
C_\vartheta = \mathbb{E}\left[\int_0^T e^{-r t}\nabla_x U(X_t^\vartheta)Y^\vartheta_t\,\mathrm{d} t\right]
\]
Therefore an $\alpha$-confidence interval for $H(\vartheta_0)$ is given by
\[
\left[H(\wh{\vartheta}_T) - \frac{z_{\alpha/2}}{\sqrt{T}} C_{\wh{\vartheta}_T}\wh{I}_T(\wh{\vartheta}_T)^{-1/2},
H(\wh{\vartheta}_T) + \frac{z_{\alpha/2}}{\sqrt{T}} C_{\wh{\vartheta}_T}\wh{I}_T(\wh{\vartheta}_T)^{-1/2} \right].
\]
In practice, $H(\wh{\vartheta}_T)$ and $C_{\wh{\vartheta}_T}$ will be computed by Monte Carlo simulations by a suitable discretization if needed. Of course, the same argument is possible in the case where $X^\vartheta$ is discretely observed; cf. Example \ref{ex:ou2}.
\end{example}

\subsection{Functionals of integrated professes}
Let us consider the following quantity: for a function $\varphi_\vartheta:\mathbb{R}\to \mathbb{R}$ and $T\in (0,1]$,
\[
H(\vartheta) = \mathbb{E}\left[\varphi_{\vartheta}\left(\frac{1}{T}\int_0^T X_t^{\vartheta} \,\mathrm{d} t\right)\right].
\]
We use the following notation for simplicity:
\[
X_* = \frac{1}{T}\int_0^T X_t\,\mathrm{d} t,
\]
for a process $X=(X_t)_{t\in [0,1]}$. Then, we have the following theorem.
\begin{thm}\label{thm:asian}
Suppose that there exists an integer $n\ge 1$ and $\vartheta\in \Theta$ such that
$\varphi_{\vartheta}^{(n)}(x)$ is Lipschitz continuous:
\[
|\varphi_{\vartheta}^{(n)}(x) - \varphi_{\vartheta}^{(n)}(y)|\lesssim |x-y|, \quad x,y\in \mathbb{R}.
\]
Moreover, suppose that A4($q$) holds for some $q\ge 2n$, and that
\[
\varphi_{\vartheta}^{(k)}(X^{\vartheta}_*) \in L^r,\quad k=1,\dots,n,
\]
for the constant $r>1$ with $1/r + 1/q = 1$.
Then, condition \eqref{mvt} holds with
\[
G_\vartheta = \varphi^{(1)}_\vartheta (X^\vartheta_*)Y^\vartheta_*.
\]
\end{thm}

\begin{proof}
It follows from Jensen's inequality that
\begin{align*}
|X^{\vartheta+u}_* - X^\vartheta_* - u^\top Y^\vartheta_*| &\le \frac{1}{T}\int_0^T |X^{\vartheta+u}_t - X^\vartheta_t - u^\top Y^\vartheta_t|\,\mathrm{d} t  \le \|X^{\vartheta + u} - X^\vartheta - u^\top Y^\vartheta\|,
\end{align*}
with probability one.
Hence, $Y^\vartheta_* = \frac{1}{T}\int_0^T Y^\vartheta_t\,\mathrm{d} t$ is the derivative of $X^\vartheta$ w.r.t. $\vartheta$.

We can take the same argument as in Theorem \ref{thm:int-V}: we use Taylor's formula and  H\"older's inequality to obtain
\begin{align*}
&\left|\mathbb{E}[h(X^{\vartheta+u}) - h(X^\vartheta) -  \varphi^{(1)}(X^\vartheta_*)u^\top Y^\vartheta_*]\right|  \\
&\le \mathbb{E}\left| \varphi^{(1)}(X^{\vartheta}_*)(X^{\vartheta+u}_* - X^\vartheta_* - u^\top Y^\vartheta_*)\right|
+ \sum_{k=1}^{n-1}\frac{1}{k!}\mathbb{E}\left| \varphi^{(k)}(X^{\vartheta}_*)(X^{\vartheta+u}_* - X^\vartheta_*)^k \right| \\
&\quad + \frac{1}{n!} \mathbb{E}\left| \varphi^{(k)}(\wt{X}^{\vartheta,u}_*)(X^{\vartheta+u}_* - X^\vartheta_*)^n \right|.
\end{align*}
Then, the same argument as in the proof of Theorem \ref{thm:int-V}
enables us to check condition \eqref{mvt} in Theorem \ref{thm:cond}.

\end{proof}

\begin{example}\label{ex:asian}
When $X^\vartheta$ is a stock price, the price of an {\it Asian call option} for $X^\vartheta$ with maturity $T$ and strike price $K$ is given by
\[
{\cal C}_{T,K} = \mathbb{E}\left[\max \left\{ \frac{1}{T}\int_0^T X_t^{\vartheta_0}\,\mathrm{d} t - K,0\right\} \right]
\]
where $\delta>0$ is an interest rate and $\mathbb{E}$ is usually taken as an expectation with respect to the risk-neutral probability. This is approximated as
\[
H_{\epsilon}(\vartheta_0) := \mathbb{E}\left[\varphi_{\epsilon}\left(  \frac{1}{T}\int_0^T X_t^{\vartheta_0}\,\mathrm{d} t \right)\right]
\]
by a function $\varphi_{\epsilon}(x) \in C^\infty(\mathbb{R})$ such that
\[
\sup_{x}|\varphi_{\epsilon}(x) - \max\{x-K,0\}| \to 0,\quad {\epsilon}\to 0.
\]
For example, we can take a function $\varphi_{\epsilon}(x) = 2^{-1}(\sqrt{(x-K)^2 + {\epsilon}^2} + x - K)$.
Then, it follows by the dominated convergence theorem that $H_{\epsilon}(\vartheta) \to {\cal C}_{T,K}$ as ${\epsilon}\to 0$ if $X^\vartheta\in L^1$.

Assume that a suitable estimator of $\vartheta_0\in \mathbb{R}^p$ is obtained, e.g.,
\[
\sqrt{T} (\wh{\vartheta}_T - \vartheta_0) \stackrel{d}{\longrightarrow} N(0,\Sigma),\quad T\to \infty,
\]
for a positive-definite matrix $\Sigma\in \mathbb{R}^p\otimes \mathbb{R}^p$.
Then, we can apply Theorem \ref{thm:asian} to $H_{\epsilon}(\vartheta)$, and letting $T\to \infty$ as well as ${\epsilon}\to 0$,  we have
\[
\sqrt{T}\left(H_{\epsilon}(\wh{\vartheta}_T) - H_{\epsilon}(\vartheta_0)\right) \stackrel{d}{\longrightarrow}  N(0,C_{\vartheta_0}^\top \Sigma \,C_{\vartheta_0}),
\]
where
\[
C_{\vartheta} = \lim_{{\epsilon}\to 0}\mathbb{E}\left[\frac{Y^\vartheta_*}{2}\left\{\frac{(X^\vartheta_* - K)}{\sqrt{(X^\vartheta_* - K)^2 + {\epsilon}^2}} + 1\right\}\right]  = \mathbb{E}\left[\frac{Y^\vartheta_*}{2}\left\{\mathrm{sgn}(X^\vartheta_* - K) + 1\right\} \right],
\]
with $\mathrm{sgn}(z) = \mathbf{1}_{\{z>0\}} - \mathbf{1}_{\{z<0\}}$.
Note that this quantity would be computed by Monte Carlo simulation in practice with $\vartheta_0$ replaced by $\wh{\vartheta}_T$, or some estimators based on discrete samples of $X^\vartheta$ in practice.
We will discuss when the condition A4$(q)$ holds when $X^\vartheta$ is a semimartingale with jumps in Section \ref{sec:sde}.
\end{example}

\begin{remark}\label{rem:general}
According to the proof of Theorem \ref{thm:asian}, we can consider more general functionals for $X^\vartheta_*$ under some smoothness conditions for $\varphi_\vartheta$.
That is, suppose that there exists an $\mathbb{R}^p$-valued random variable $\wt{Y}^\vartheta$ such that the following inequality holds:
\begin{align}
|X^{\vartheta+u}_* - X^\vartheta_* - u^\top \wt{Y}^\vartheta|\lesssim \|X^{\vartheta+u} - X^\vartheta - u^\top Y^\vartheta\| + |u|^{1+\delta}\quad a.s., \label{eq:X-star}
\end{align}
for $\delta>0$, and the derivative is $Y^\vartheta$. Then, the same proof as that of Theorem \ref{thm:asian} works with
\begin{align*}
G_\vartheta = \varphi^{(1)}_\vartheta(X^\vartheta_*) \wt{Y}^\vartheta.
\end{align*}
For example, let
\[
X^\vartheta_* = \int_0^T U(X^\vartheta_t)\,\mathrm{d} t
\]
for $T>0$ and $U\in C^2(\mathbb{R})$ be a function with bounded derivatives. Then we find that
\[
\wt{Y}^\vartheta = \int_0^T U^{(1)}(X^\vartheta_t)Y^\vartheta_t\,\mathrm{d} t,
\]
since it follows that
\begin{align*}
|X^{\vartheta+u}_* - X^\vartheta_*  - u^\top \wt{Y}^\vartheta| &\le \int_0^T |U(X^{\vartheta + u}_t) - U(X^\vartheta_t) - u^\top U^{(1)}(X^\vartheta_t)Y^\vartheta_t|\,\mathrm{d} t \\
&\lesssim \int_0^T |U^{(1)}(X^\vartheta_t)(X^{\vartheta+u}_t - X^\vartheta_t - u^\top Y^\vartheta_t)\,\mathrm{d} t + |u|^2 \\
&\lesssim \|X^{\vartheta+u} - X^\vartheta - u^\top Y^\vartheta\| + |u|^2.
\end{align*}
This argument can include Theorem \ref{thm:int-V}.
\end{remark}

\begin{remark}
You might also be interested in the case where $X^\vartheta_*$ is an extreme-type functional such as
$X_*^\vartheta = \inf_{s \le t} X^\vartheta_s$, which is important when, e.g., $\varphi(x) = \mathbf{1}_{\{x < 0\}}$, the function $H(\vartheta) =\mathbb{E}[\varphi(X_*^\vartheta)]$ stands for the hitting time distribution:
\[
H(\vartheta) = {\mathbb{P}}(\tau^\vartheta \le t),\quad \tau^\vartheta = \inf\{t>0\,|\,X^\vartheta_t < 0\},
\]
or we can approximate $\varphi$ with a bounded smooth function such as, e.g.,  $\varphi_{\epsilon}(x) = [1 + e^{-x/{\epsilon}}]^{-1}\ \to \varphi(x)\ ({\epsilon}\to 0)$,
among others.

When $X^\vartheta$ is a continuous diffusion process, Gobet and Kohatsu \cite{gk03} obtain a derivative of ${\mathbb{P}}(\tau^\vartheta \le t)$ via Malliavin Calculus.
However, in our approach, it is not so easy to find a suitable random variable $\wt{Y}^\vartheta$ satisfying the inequality \eqref{eq:X-star}, except for a trivial case where the derivative process $Y^\vartheta$ is a constant.
One might expect that $\wt{Y}^\vartheta = \sup_{s\le t}Y_s^\vartheta$ in general, but it fails.
This important  case is an open problem.
\end{remark}

\section{Expected functionals of semimartingales}\label{sec:sde}

\subsection{Stochastic differential equations with jumps}

On a stochastic basis $(\Omega,\mathcal{F},\mathbb{F},{\mathbb{P}})$ with a filtration $\mathbb{F}=(\mathcal{F}_t)_{t\ge 0}$,
consider a $1$-dim stochastic process $X=(X_t)_{t\in [0,T]}$ that satisfies the following stochastic differential equation (SDE) with a  multidimensional parameter $\vartheta\in \Theta\subset \mathbb{R}^p$:
\begin{align}
X^\vartheta_t = x(\vartheta) + \int_0^t a(X^\vartheta_s,\vartheta)\,\mathrm{d} s + \int_0^t b(X^\vartheta_s,\vartheta)\,\mathrm{d} W_s
+\int_0^t \int_{E} c(X^\vartheta_{s-},z,\vartheta)\,\wt{N}(\mathrm{d} t,\mathrm{d} z),  \label{sde}
\end{align}
where $E=\mathbb{R} \setminus \{0\}$; $x:\Theta \to \mathbb{R}$; $a:\mathbb{R}\times \Theta\to \mathbb{R}$, $b:\mathbb{R} \times \Theta \to \mathbb{R}\otimes \mathbb{R}$ and $c:\mathbb{R}\times E\times \Theta\to \mathbb{R}$; $W$ is a $\mathbb{F}$-Wiener process.
Moreover, $\wt{N}(\mathrm{d} t,\mathrm{d} z):=N(\mathrm{d} t,\mathrm{d} z) - \nu(z)\,\mathrm{d} z\mathrm{d} t$, which is the {\it compensated Poisson random measure}, where $N$ is a Poisson random measure associated with a $\mathbb{F}$-L\'evy process, say $Z=(Z_t)_{t\ge 0}$ with the L\'evy density $\nu$:
\[
N(A\times (0,t]) = \sum_{s \le t}\mathbf{1}_{\{\Delta Z_s\in A\}},\quad A\subset E,
\]
and $\mathbb{E}[N(\mathrm{d} t,\mathrm{d} z)] = \nu(z)\,\mathrm{d} z\mathrm{d} t$.

In what follows, we assume that $\nu$ is {\it essentially} known:
some cases can be rewritten into a model for a known $\nu$ even if $\nu$ has some unknown parameters (see Remark \ref{rem:unknown} below). However, if it is not the case, the situation may be totally different from ours, and the argument in this section would no longer work; see Remark \ref{rem:unknown2}.

\begin{remark}\label{rem:unknown}
Some cases where the L\'evy density $\nu$ depends on an unknown parameter, say $\nu_\vartheta$, can be rewritten into the form of \eqref{sde} with a known L\'evy process by changing the coefficients $a$ and $c$, suitably.
For example, consider the following SDE:
\begin{align}
\mathrm{d} X_t = a(X_t)\,\mathrm{d} t + b(X_t)\,\mathrm{d} W_t +  \int_E c(X_{t-},z)\,N_\vartheta(\mathrm{d} t,\mathrm{d} z), \label{sde2}
\end{align}
where $N_\vartheta$ is the Poisson random measure associated with a compound Poisson process of the form
$Z^\vartheta_t = \sum_{i=1}^{N_t} U^\vartheta_i$ such that $N$ is a Poisson process with intensity $\lambda_0$, and the $U_i^\vartheta$'s are i.i.d. sequences with probability density $f_\vartheta$ with $\mathbb{E}[U^\vartheta_i] = \eta$ and $Var(U^\vartheta_i) = \zeta^2$.
Suppose that $\lambda_0$ is known, but $\vartheta=(\eta,\zeta)$ is unknown.
In this case, we can rewrite $Z^\vartheta\,(=Z^{(\eta,\zeta)})$ as
\[
Z^{(\eta,\zeta)}_t = \sum_{i=1}^{N_t} (\zeta U_i^{(0,1)} + \eta) = \int_0^t \int_E (\zeta z + \eta)\,N_{(0,1)}(\mathrm{d} s,\mathrm{d} z),
\]
where $N_{(0,1)}$ is the Poisson random measure associated with $Z^{(0,1)}$.
Then, the SDE \eqref{sde2} is written as
\begin{align*}
\mathrm{d} X_t &= a(X_t)\,\mathrm{d} t + b(X_t)\,\mathrm{d} W_t +  \int_E c(X_{t-},\zeta z+\eta)\,N_{(0,1)}(\mathrm{d} t,\mathrm{d} z) \\
&= \left[a(X_t) +  \lambda_0\int_E c(X_t,\zeta z+\eta)f_{(0,1)}(z)\,\mathrm{d} z\right]\,\mathrm{d} t+ b(X_t)\,\mathrm{d} W_t  \\
&\qquad +  \int_E c(X_{t-}, \zeta z+\eta)\,\wt{N}_{(0,1)}(\mathrm{d} t,\mathrm{d} z),
\end{align*}
where the L\'evy density $\lambda_0 f_{(0,1)}(z)$ is known. See also Example \ref{ex:ou}.
\end{remark}

The semimartingale $X^\vartheta$ in \eqref{sde} is a ${\cal X}=\mathbb{D}([0,T])$-valued random element.
In what follows, we consider a metric space $({\cal X},\|\cdot\|)$  with the sup norm:
\[
\|X^\vartheta\|=\|X^\vartheta\|_T := \sup_{t\in [0,T]}|X_t^\vartheta|.
\]

We make some assumptions.

\begin{asb}\label{asb:lg}
For each $x, z\in \mathbb{R}$,
\[
|a(x,\vartheta)| + |b(x,\vartheta)|  \lesssim 1 + |x|;\quad |c(x,z,\vartheta)|\lesssim |z|(1 + |x|),
\]
uniformly in $\vartheta\in \Theta$.
\end{asb}

\begin{asb}\label{asb:smooth}
The functions $a,b$ and $c$ are twice differentiable in $x$, and
the derivatives $\nabla_x^k a$ and $\nabla_x^k b$ $(k=1,2)$ are uniformly bounded.
Moreover, $|\nabla_x^k c(x,z,\vartheta)| \lesssim |z|$.
\end{asb}

\begin{asb}\label{asb:dot}
The functions $a,b$ and $c$ are differentiable in $\vartheta$.
It follows that
\[
|\dot{a}(x,\vartheta)| + |\dot{b}(x,\vartheta)|\lesssim 1 + |x|; \quad |\dot{c}(x,z,\vartheta)|\lesssim  |z|(1 + |x|),
\]
uniformly in $\vartheta\in \Theta$.
\end{asb}

\begin{asb}\label{asb:nu-moment}
For any $p>0$, $\displaystyle \int_{|z|>1} z^p\,\nu(z)\,\mathrm{d} z < \infty.$
\end{asb}

\begin{asb}\label{asb:moment}
For any $p>0$ and $T>0$, $\|X\|_T^p < \infty$.
\end{asb}

\subsection{Derivative processes}

Let $Y^\vartheta =(Y^\vartheta_t)_{t\ge 0}$ be a $p$-dim stochastic process satisfying the following SDE:  $Y^\vartheta_0=\dot{x}(\vartheta)$,
\begin{align}
\mathrm{d} Y^\vartheta_t &= A(X^\vartheta_t,Y^\vartheta_t,\vartheta)\,\mathrm{d} t + B(X^\vartheta_t,Y^\vartheta_t,\vartheta)\,\mathrm{d} W_t + \int_{E} C(X^\vartheta_{t-},Y^\vartheta_{t-},z,\vartheta)\,\wt{N}(\mathrm{d} t,\mathrm{d} z), \label{Y-sde}
\end{align}
for each $\vartheta\in \Theta$, where
\begin{align*}
A(x,y,\vartheta) &= \nabla_x a(x,\vartheta)y + \dot{a}(x,\vartheta); \\
B(x,y,\vartheta) &= \nabla_x b(x,\vartheta)y + \dot{b}(x,\vartheta); \\
C(x,y,z,\vartheta) &= \nabla_x c(x,z,\vartheta)y + \dot{c}(x,z,\vartheta).
\end{align*}

In this section, we will show that the above $Y^\vartheta=(Y^\vartheta_t)_{t\ge 0}$ can be
the {\it derivative process} of $X^\vartheta$ with respect to $\vartheta$ in the sense of $L^q$.
For that purpose, we shall give some preliminary lemmas.

\begin{lemma}\label{lem:M}
Let $g:\mathbb{R}\times \mathbb{R}\to \mathbb{R}$ be of polynomial growth.
Then ,under B\ref{asb:moment}, it holds for $p=2^m\,(m\in \mathbb{N})$ that
\[
\mathbb{E}\left\| \int_0^t \int_E g(X_{s-},z)\,\wt{N}(\mathrm{d} s,\mathrm{d} z) \right\|_T^p \lesssim \mathbb{E}\left[\int_0^T\int_E |g(X_{s-},z)|^p\,\nu(z)\,\mathrm{d} z\mathrm{d} s\right]
\]
\end{lemma}

\begin{proof}
See Shimizu and Yoshida \cite{sy06}, Lemma 4.1.
\end{proof}

\begin{lemma}\label{lem:DX}
Suppose that assumptions B\ref{asb:lg} -- B\ref{asb:moment} hold, and that $\dot{x}(\vartheta)$ is uniformly bounded on $\Theta$.
Then, it follows for any $T>0$, $p\ge 2$ and $u\in \mathbb{R}^p$ with $\vartheta+u\in \Theta$ that
\[
\mathbb{E}\|X^{\vartheta+u} - X^\vartheta\|_T^p \lesssim |u|^p.
\]
\end{lemma}

\begin{proof}
It follows from Jensen's inequality that
\begin{align}
|X^{\vartheta+u}_t - X^\vartheta_t|^p
&\lesssim |x(\vartheta+u) - x(\vartheta)|^p + t^{p-1}\int_0^t |\wt{A}_s(u,\vartheta)|^p\,\mathrm{d} s + \left|\int_0^t \wt{B}_t(u,\vartheta)\,\mathrm{d} W_s\right|^p \notag \\
&\quad + \left|\int_0^t \int_{E}\wt{C}_s(u,z,\vartheta)\,\wt{N}(\mathrm{d} s,\mathrm{d} z)\right|^p, \label{DX}
\end{align}
with
\begin{align*}
\wt{A}_t(u,\vartheta)&:= a(X^{\vartheta+u}_t,\vartheta+u) - a(X^{\vartheta}_t,\vartheta); \\
\wt{B}_t(u,\vartheta)&:= b(X^{\vartheta+u}_t,\vartheta+u) - b(X^{\vartheta}_t,\vartheta); \\
\wt{C}_t(u,z,\vartheta)&:= c(X^{\vartheta+u}_t,z,\vartheta+u) - c(X^{\vartheta}_t,z,\vartheta).
\end{align*}
Then, since it holds that $|x(\vartheta+u) - x(\vartheta)|\lesssim |u|$ from the mean value theorem,
Lemma \ref{lem:M} and Burkholder-Davis-Gundy's inequality yield  that
\begin{align*}
\mathbb{E}\left\|X^{\vartheta+u} - X^\vartheta \right\|_T^p &\lesssim |u|^{p} + \int_0^T \mathbb{E}\left[ |\wt{A}_s(u,\vartheta)|^p + |\wt{B}_s(u,\vartheta)|^p\right]\,\mathrm{d} s +  \mathbb{E}\left[\int_0^T \int_{E}|\wt{C}_s(u,z,\vartheta)|^p\,\nu(z)\,\mathrm{d} z\mathrm{d} s\right]
\end{align*}
It follows from the mean value theorem and assumptions B\ref{asb:lg} -- B\ref{asb:dot} that
\begin{align*}
|\wt{A}_t(u,\vartheta)|^p &= |\nabla_xa(X^*,\vartheta^*)(X_t^{\vartheta+u} - X_t^\vartheta) + \dot{a}(X^*,\vartheta^*)^\top u|^p \\
&\lesssim |X_t^{\vartheta+u} - X_t^\vartheta|^p + \left(1 + \|X\|_T^p\right)|u|^p
\end{align*}
Hence, it follows from B\ref{asb:moment} that
\[
\mathbb{E}|\wt{A}_t(u,\vartheta)|^p \lesssim \mathbb{E}|X_t^{\vartheta+u} - X_t^\vartheta|^p + |u|^p
\]

Similarly, we also have that
\begin{align*}
\mathbb{E}|\wt{B}_t(u,\vartheta)|^p &\lesssim \mathbb{E}|X^{\vartheta+u}_t - X^\vartheta_t|^p +  |u|^p;\\
\mathbb{E}|\wt{C}_t(u,z,\vartheta)|^p &\lesssim|z|^p\left(\mathbb{E}|X^{\vartheta+u}_t - X^\vartheta_t|^p +  |u|^p\right).
\end{align*}
Hence, assumption B\ref{asb:nu-moment} yields that
\[
\mathbb{E}\left\|X^{\vartheta+u} - X^\vartheta \right\|_T^p \lesssim |u|^{p} + \int_0^T \mathbb{E}\left\|X^{\vartheta+u} - X^\vartheta \right\|_t^p \,\mathrm{d} t.
\]
Finally, Gronwall's inequality completes the proof.
\end{proof}

The next theorem is the consequence of this section.

\begin{thm}\label{thm:deriv-X}
Suppose that assumptions B\ref{asb:lg} -- B\ref{asb:moment} hold.
Moreover, suppose that the initial value $x(\vartheta)=X^\vartheta_0$ is twice differentiable with respect to the bounded derivatives, and that
the solution $Y^\vartheta$ to \eqref{Y-sde} satisfies $\|Y^\vartheta\|_T < \infty$ for any $T>0$.
Then, for any $p\ge 2$, there exists a positive constant $C_p$ depending on $p$ such that
\[
\mathbb{E}\left\|X^{\vartheta+u} - X^\vartheta - u^\top Y^\vartheta \right\|_T^p \le C_p|u|^{2p},\quad h \in \mathbb{R}^p.
\]
\end{thm}

\begin{proof}
First, we shall consider the case where $p=2^m\ (m\in \mathbb{N})$.
Applying Jensen's inequality to the $\mathrm{d} t$-integral part, we see that
\begin{align}
|X^{\vartheta+u}_t - X^\vartheta_t - u^\top Y^\vartheta_t|^p
&\lesssim |x(\vartheta+u) - x(\vartheta) - u^\top \dot{x}(\vartheta)|^p + t^{p-1}\int_0^t |\wt{A}_s(u,\vartheta)|^p\,\mathrm{d} s \notag \\
&\quad + \left|\int_0^t \wt{B}_t(u,\vartheta)\,\mathrm{d} W_s\right|^p + \left|\int_0^t \int_{E}\wt{C}_{s-}(u,z,\vartheta)\,\wt{N}(\mathrm{d} s,\mathrm{d} z)\right|^p, \label{DX-Y}
\end{align}
where
\begin{align*}
\wt{A}_t(u,\vartheta)&:= a(X^{\vartheta+u}_t,\vartheta+u) - a(X^{\vartheta}_t,\vartheta) - u^\top [\nabla_x a(X^\vartheta_t,\vartheta)Y^\vartheta_t+ \dot{a}(X^\vartheta_t,\vartheta)]; \\
\wt{B}_t(u,\vartheta)&:= b(X^{\vartheta+u}_t,\vartheta+u) - b(X^{\vartheta}_t,\vartheta) - u^\top [\nabla_x b(X^\vartheta_t,\vartheta)Y^\vartheta_t+ \dot{b}(X^\vartheta_t,\vartheta)]; \\
\wt{C}_t(u,z,\vartheta)&:= c(X^{\vartheta+u}_t,z,\vartheta+u) - c(X^{\vartheta}_t,z,\vartheta) - u^\top [\nabla_x c(X^\vartheta_t,z,\vartheta)Y^\vartheta_t+ \dot{c}(X^\vartheta_t,z,\vartheta)].
\end{align*}

Take $\sup_{t\in [0,T]}$ and the expectation $\mathbb{E}$ on both sides to obtain that
\begin{align*}
\mathbb{E}\|X^{\vartheta+u} - X^\vartheta - u^\top Y^\vartheta\|_T^p
&\lesssim |u|^{2p} + \int_0^T \mathbb{E}|\wt{A}_t(u,\vartheta)|^p\,\mathrm{d} t + \mathbb{E}\left\|\int_0^\cdot \wt{B}_s(u,\vartheta)\,\mathrm{d} W_s\right\|_T^p\\
&\quad + \mathbb{E}\left\|\int_0^\cdot \int_{E}\wt{C}_{s-}(u,z,\vartheta)\,\wt{N}(\mathrm{d} s,\mathrm{d} z)\right\|_T^p.
\end{align*}
Using Burkholder-Davis-Gundy's inequality and Lemma \ref{lem:M}, we have that
\begin{align*}
\mathbb{E}\|X^{\vartheta+u} - X^\vartheta - u^\top Y^\vartheta\|_T^p
&\lesssim |u|^{2p} + \int_0^T \mathbb{E}|\wt{A}_s(u,\vartheta)|^p\,\mathrm{d} s + \mathbb{E}\left|\int_0^T \wt{B}_s^2(u,\vartheta)\,\mathrm{d} s\right|^{p/2} \\
&\quad +  \mathbb{E}\left[\int_0^T \int_{E}\wt{C}_s^p(u,z,\vartheta)\,\nu(z)\,\mathrm{d} z\mathrm{d} s\right] \\
&\lesssim |u|^{2p} + \int_0^T \mathbb{E}\left[ |\wt{A}_s(u,\vartheta)|^p + |\wt{B}_s(u,\vartheta)|^p\right]\,\mathrm{d} s \\
&\quad +  \mathbb{E}\left[\int_0^T \int_{E}|\wt{C}_s(u,z,\vartheta)|^p\,\nu(z)\,\mathrm{d} z\mathrm{d} s\right]
\end{align*}

According to assumptions B\ref{asb:lg}, B\ref{asb:smooth}, and Taylor's formula, we have, e.g.,
\begin{align*}
\wt{A}_t(u,\vartheta) &= \nabla_x a(X^{\vartheta}_t)(X^{\vartheta+u}_t - X^\vartheta_t) + u^\top \dot{a}(X^\vartheta,\vartheta)   + \frac{1}{2}[(X^{\vartheta+u}-X^\vartheta)\nabla_x + u^\top \nabla_\vartheta]^2a(X_t^*,\vartheta^*) \\
&\quad - u^\top [\nabla_x a(X^\vartheta_t,\vartheta)Y^\vartheta_t+ \dot{a}(X^\vartheta,\vartheta)],
\end{align*}
where $X^*$ is a random variable between $X_t^{\vartheta+u}$ and $X^\vartheta$, $\vartheta^* \in [\vartheta,\vartheta+u]$.
Since the second derivatives are bounded, and from B\ref{asb:dot}, we have that
\begin{align*}
|\wt{A}_t(u,\vartheta)|^p &\lesssim|X^{\vartheta+u}_t - X^\vartheta_t - u^\top Y^\vartheta_t|^p + |u|^{2p} + |X^{\vartheta+u}_t-X^\vartheta_t|^{2p}
+ |u|^p|X^{\vartheta+u}_t-X^\vartheta_t|^p
\end{align*}
Similarly, we also have that
\begin{align*}
|\wt{B}_t(u,\vartheta)|^p &\lesssim|X^{\vartheta+u}_t - X^\vartheta_t - u^\top Y^\vartheta_t|^p + |u|^{2p} + |X^{\vartheta+u}_t-X^\vartheta_t|^{2p}
+ |u|^p|X^{\vartheta+u}_t-X^\vartheta_t|^p;\\
|\wt{C}_t(u,z,\vartheta)|^p &\lesssim|z|^p\left(|X^{\vartheta+u}_t - X^\vartheta_t - u^\top Y^\vartheta_t|^p + |u|^{2p} + |X^{\vartheta+u}_t-X^\vartheta_t|^{2p} + |u|^p|X^{\vartheta+u}_t-X^\vartheta_t|^p\right).
\end{align*}

Hence, under B\ref{asb:nu-moment}, it follows from Lemma \ref{lem:DX} that
\[
\mathbb{E}\|X^{\vartheta+u} - X^\vartheta - u^\top Y^\vartheta\|_T^p
\lesssim |u|^{2p} + \int_0^T \mathbb{E}\|X^{\vartheta+u} - X^\vartheta - u^\top Y^\vartheta\|_t^p \,\mathrm{d} t,
\]
and Gronwall's inequality yields the consequence.

For any $p\ge 2$, we write the binomial expansion of $p$ as $p= \sum_{k=1}^m \delta_k 2^k$, where $m$ is an integer and $\delta_k = 0$ or $1$.
Note that we have already proved the consequence for $p$ with $m=1$ and $\delta_1=0,1$.
Next, we assume that the consequence also  holds true for some $m$ and any $\delta_k\ (k=1,2,\dots,m)$.
Then, the Cauchy-Schwartz inequality yields that for $q = \sum_{k=2}^{m} 2^{k}\delta_{k-1}$,
\begin{align*}
&\mathbb{E}\|X^{\vartheta+u} - X^\vartheta - u^\top Y^\vartheta \|_T^p \\
&= \mathbb{E}\left[  \|X^{\vartheta+u} - X^\vartheta - u^\top Y^\vartheta \|_T^{2^m\delta_m}\prod_{k=1}^{m-1} \|X^{\vartheta+u} - X^\vartheta - u^\top Y^\vartheta \|_T^{\delta_k2^k}\right] \\
&\le \sqrt{\mathbb{E}\left[ \left\|X^{\vartheta+u} - X^\vartheta - u^\top Y^\vartheta \right\|_T^{2^{m+1}\delta_m}\right] }
\sqrt{\mathbb{E}\left[  \left\|X^{\vartheta+u} - X^\vartheta - u^\top Y^\vartheta \right\|_T^{\sum_{k=2}^{m} 2^{k}\delta_{k-1}}\right] } \\
&\le \sqrt{C_{2^m\delta_m}|u|^{2\cdot 2^{m+1} \delta_m}} \sqrt{C_q|u|^{2q}} \\
&\lesssim |u|^{2^{m+1}\delta_m + q} = |u|^{2p}.
\end{align*}
This completes the proof.
\end{proof}

\begin{remark}\label{rem:unknown2}
If the random measure $N$ essentially includes unknown parameters, then the derivative process in the sense of $L^q$ cannot exist. To see this, consider a simple case where $X^\vartheta$ is a Poisson process with (unknown) intensity $\vartheta$: $X^\vartheta\sim Po(\vartheta t)$, which is not the case described in Remark \ref{rem:unknown}. In this case, we cannot compute the expectation $\mathbb{E}\|X^{\vartheta+u} - X^\vartheta \|_T^p$ since we do not know the joint distribution of $(X^{\vartheta+u},X^\vartheta)$.
 This consideration indicates that we should be careful when we compute expected functionals of $X^\vartheta$ by Monte Carlo simulation when it has an unknown jump part.
\end{remark}

\begin{example}[L\'evy processes]\label{ex:levy}
Consider a 1-dim L\'evy process $X^\vartheta$ starting at $x>0$:
\[
X^\vartheta_t = x + \mu t + \sigma W_t + \eta S_t,
\]
where $S$ is a known {\it L\'evy process} with $\mathbb{E}[S_1]=1$ and $\eta\ne 0$.
We set $\vartheta=(\mu,\sigma,\eta) \in \Theta\subset \mathbb{R}^3$.
Then, this is the case of \eqref{sde} with
\[
a(x,\vartheta) = \mu + \eta,\quad b(x,\vartheta) = \sigma,\quad c(x,z,\vartheta)=\eta z,\quad X_0=x,
\]
Hence, the derivative process $Y^\vartheta$ is a 3-dim  L\'evy process of the form
\[
Y^\vartheta_t =\left(t, W_t, S_t\right)
\]
\end{example}

\begin{example}\label{ex:ou}
Consider an  O-U process $X=(X_t)_{t\ge 0}$ written as
\begin{align}
\mathrm{d} X^\vartheta_t = -\mu X^\vartheta_t\,\mathrm{d} t + \sigma \,\mathrm{d} W_t + \mathrm{d} Z^\eta_t, \quad X_0 = x\ (\mbox{const.}) \label{ou}
\end{align}
where $\vartheta = (\mu,\sigma,\eta)$, $W$ is a Wiener process, and $Z^\eta$ is a compound Poisson process with
known intensity, and  the mean of the jumps  is $\eta$
Then, the SDE \eqref{ou} is rewritten as
\[
X^\vartheta_t = x +  \int_0^t (-\mu X^\vartheta_s + \eta)\,\mathrm{d} s + \sigma W_t + \int_0^t \int_E (z+\eta)\wt{N}(\mathrm{d} t,\mathrm{d} z),
\]
where $\wt{N}$ is the compensated Poisson random measure associated with $Z^0\,(\eta =0)$ (see Remark \ref{rem:unknown}).

Then, the derivative process $Y^\vartheta=(Y^1_t,Y^2_t,Y^3_t)_{t\ge 0}$ satisfies the following SDE:
\begin{align*}
Y^1_t &= \int_0^t (-\mu Y_s^1 - X^\vartheta_s)\,\mathrm{d} s; \quad
Y^2_t = - \mu \int_0^t Y^2_s\,\mathrm{d} s + W_t; \\
Y^3_t &=  \int_0^t (1 - \mu Y^3_s)\,\mathrm{d} s + Z_t^0,
\end{align*}
since $\int_E z\nu_0(z)\,\mathrm{d} z = 0$.
The equation for $Y^1$ is an ordinary differential equation for almost all $\omega \in \Omega$, and
the equations for $Y^2$ and $Y^3$ are O-U type SDEs.
Therefore, we can solve these equations explicitly, as follows:
\begin{align*}
Y^1_t &= -\int_0^t X^\vartheta_s e^{-\mu(t-s)}\,\mathrm{d} s,\quad
Y^2_t = \int_0^t e^{-\mu(t-s)}\,\mathrm{d} W_s,\\
Y^3_t &= \frac{1}{\mu} (1-e^{-\mu t}) + \int_0^t e^{-\mu (t-s)}\,\mathrm{d} Z^0_s,
\end{align*}
and
\begin{align}
X^\vartheta_t = xe^{-\mu t} + \int_0^t e^{-\mu(t-s)}[\sigma\,\mathrm{d} W_s + \mathrm{d} Z^\eta_s]. \label{eq:ou}
\end{align}
\end{example}

\subsection{Expected functionals for semimartingales}

For each $\vartheta\in \Theta$, let $\varphi_\vartheta:\mathbb{R}\to \mathbb{R}$ and
\[
H(\vartheta) = \mathbb{E}\left[\varphi_\vartheta(X^{\vartheta}_*)\right],
\]
where $X^\vartheta_*$ is a $\mathbb{R}$-valued random functional of $X^\vartheta$ such that the inequality
\begin{align}
|X^{\vartheta+u}_* - X^\vartheta_* - u^\top \wt{Y}^\vartheta|\lesssim \|X^{\vartheta+u} - X^\vartheta - u^\top Y^\vartheta\| + |u|^{1+\delta}\quad a.s.,  \label{Y-tilde}
\end{align}
holds true for some $\wt{Y}^\vartheta$ and $\delta>0$; see Remark \ref{rem:general} for some examples.
Summing up our results in Sections \ref{sec:mc}, \ref{sec:C-space} and \ref{sec:sde} with Remark \ref{rem:general}, we can immediately obtain the following result.
\begin{thm}\label{thm:final}
Suppose that the same assumptions as in Theorem \ref{thm:deriv-X} hold.
Moreover, suppose that there exists an integer $n\ge 1$  such that
$\varphi_{\vartheta_0}^{(n)}(x)$ is Lipschitz continuous:
\[
|\varphi_{\vartheta_0}^{(n)}(x) - \varphi_{\vartheta_0}^{(n)}(y)|\lesssim |x-y|, \quad x,y\in \mathbb{R},
\]
and that for some constant $r>2$,
\[
\varphi_{\vartheta_0}^{(k)}(X^{\vartheta_0}_*) \in L^r,\quad k=1,\dots, n.
\]
Furthermore, assume that we have an estimator of $\vartheta_0$ based on some observations depending on a parameter $n$, say $\wh{\vartheta}_n$, such that assumption A\ref{as:Z} holds true.
Then the asymptotic distribution of $H(\wh{\vartheta}_n)$ is specified:
\[
\gamma_{n*}^{-1}(H(\wh{\vartheta}_n) - H(\vartheta_0)) \stackrel{d}{\longrightarrow} \left(\mathbb{E}\left[\dot{\varphi}_{\vartheta_0}(X^{\vartheta_0})\right] + C_{\vartheta_0}\right)^\top Z^*,\quad n\to \infty,
\]
and the deterministic vector $C_\vartheta$ is given by
\[
C_\vartheta = \mathbb{E}\left[\varphi^{(1)}(X^\vartheta_*)\wt{Y}^\vartheta\right],
\]
where $\wt{Y}^\vartheta$ is given in \eqref{Y-tilde}.
\end{thm}

\begin{example}[Ornstein-Uhlenbeck type processes]\label{ex:ou2}
This is a continuation of the previous Example \ref{ex:ou}.
Let us consider the same SDE as \eqref{ou}, and  consider the {\it expected discounted functional} for a constant $\delta>0$,
\[
H(\vartheta) = \mathbb{E}\left[\int_0^T e^{-\delta t} V(X_t) \,\mathrm{d} t\Big| X_0=x\right],
\]
which is an important quantity in insurance and finance because such a functional can represent an option price when $X$ is a stock price (see, e.g., Karatzas and Shereve \cite{ks91}), or it can represent some aggregated “costs” or “risks” in insurance businesses when $X$ is an asset process of the company; see, e.g., Feng and Shimizu \cite{fs13}. The constant $\delta>0$ is interpreted as an interest rate.

Here, we shall consider a simple case where $V(x) = x$:
\[
H(\vartheta) = \int_0^t e^{-\delta t}\mathbb{E}[X_t]\,\mathrm{d} t,
\]
Noticing that from expression  \eqref{eq:ou},
\[
\mathbb{E}[X_t] = xe^{-\mu t} + \mathbb{E}\left[\int_0^t e^{-\mu(t-s)} \,\mathrm{d} Z^\eta_s\right] = xe^{-\mu t} + \frac{\eta}{\mu}(1-e^{-\mu t}),
\]
we can compute $H(\vartheta)$ explicitly as
\[
H(\vartheta) = \frac{x}{\mu + \delta}(1-e^{-(\mu + \delta)T})
+ \frac{\eta}{\mu}\left[\frac{1}{\delta} (1-e^{-\delta T}) - \frac{1}{\mu + \delta}(1-e^{-\mu + \delta)T})\right].
\]

Suppose that $Z$ is a compound Poisson process, and that we have a set of discrete samples
$(X_{t_1},X_{t_2},\dots,X_{t_n})$ with $t_k = kh_n$ for $h_n>0$, and
assume some asymptotic conditions on $n$ and $h_n$, e.g., $h_n\to 0$ and $nh_n^2\to 0$.
Although we omit the details of the regularity conditions here, we can construct an asymptotic normal (efficient) estimator of $\vartheta = (\mu,\sigma,\eta)$, say $\wh{\vartheta}_n$, such that
\[
\Gamma_n^{-1} (\wh{\vartheta}_n - \vartheta) \stackrel{d}{\longrightarrow} N_3(0,\Sigma), \quad n\to \infty
\]
with $\Gamma_n = {\rm diag} (1/\sqrt{nh_n}, 1/\sqrt{n}, 1/\sqrt{nh_n})$ and a diagonal matrix $\Sigma=\mathrm{diag}(\Sigma_1,\Sigma_2,\Sigma_3)$ (see, e.g., Shimizu and Yoshida \cite{sy06}).
In this case, we have  $\gamma_{n*} = 1/\sqrt{nh_n}$, and Theorem \ref{thm:int-V} says that
\[
\sqrt{nh_n}[H(\wh{\vartheta}_n) - H(\vartheta_0)] \stackrel{d}{\longrightarrow} N\left( 0, C_{\vartheta_0}^\top \mathrm{diag}(\Sigma_1,0,\Sigma_3) C_{\vartheta_0}\right)
\]
where
\[
C_\vartheta = \left(\int_0^T e^{-\delta t}\mathbb{E}[Y^\vartheta_t]\,\mathrm{d} t\right) =:(C_\vartheta^1,C_\vartheta^2,C_\vartheta^3)^\top;
\]
with $C_\vartheta^2=0$ and
\begin{align*}
C_\vartheta^1 &= \frac{\eta - \mu x}{\mu(\delta +\mu)^2}\left[1 - (\mu + \delta)e^{-(\mu + \delta)T - e^{-(\mu+\delta)T}}\right] + \frac{\eta}{\delta \mu^2}(1-e^{-\delta T})  \\ &\qquad  + \frac{\eta}{\mu^2(\mu + \delta)}(1-e^{-(\mu + \delta)T}); \\
C_\vartheta^3 &= \frac{1}{\mu}\left[\frac{1}{\delta}(1-e^{-\delta T}) - \frac{1}{\mu + \delta}(1-e^{-(\mu + \delta)T})\right].
\end{align*}
\end{example}

\section{Numerical experiments}
In this section, we shall illustrate the result of Theorem \ref{thm:final} by an example described in Introduction. That is, evaluating the statistical error of Monte Carlo estimation of a European call option based on a diffusion process of the underlying asset.
In the experiments, we assume that the asset process is observed discretely in time with the {\it small noise asymptotics}, which philosophically corresponds to a kind of long term observations; see, e.g., Shimizu \cite{s17}, Section 2.3.
We first compute estimators of unknown parameters in the process from discrete samples, and compute the price of the European call option by Monte Carlo simulations based on the estimated process, and evaluate the statistical error.

In Section \ref{sec:sim-setting}, we will describe a general framework of the simulations as well as how to construct estimators of unknown parameters in the process. In Section \ref{sec:BS}, we will particularly consider the Black-Scholes model, and investigate the asymptotic distribution of the estimated European call price by the Monte Carlo method.

We used the {\tt yuima} package in {\tt R} for simulating paths and discrete sampling from diffusion processes; see Brouste {\it et al.} \cite{yuima} for details.

\subsection{European call-type functionals under the small noise diffusions}\label{sec:sim-setting}

We assume that the stochastic process $X^\vartheta=(X^\vartheta_t)_{t \ge 0}$ satisfies the following stochastic differential equation: for $\vartheta=(\mu,\sigma)$ and ${\epsilon}>0$,
\begin{align}
\mathrm{d} X_t^{\vartheta_0,{\epsilon}} = b(X^{\vartheta_0,{\epsilon}}_t, \mu_0)\,\mathrm{d} t + {\epsilon}\cdot a(X^{\vartheta_0,{\epsilon}}_t,\sigma_0) \,\mathrm{d} W_t,\quad
X^{\vartheta_0,{\epsilon}}_0=x,
\label{small-diff}
\end{align}
and suppose that we observe $X^{\vartheta_0,{\epsilon}}$ discretely at time points $t_k=k/n\, (k=0,1,\dots, n)$ in $[0,1]$-interval, and write the samples $X_{t_k}\,(k=0,1,\dots, n)$ and denote by $\Delta_kX = X_{t_k} - X_{t_{k-1}}$.
Our purpose is to estimate the following expected functional with parameters $r, K, T>0$:
\begin{align}
{\cal C}_{T,K}^\vartheta :=\mathbb{E}\left[e^{-r T}\max \{X^\vartheta_T - K, 0\}\right]. \label{E-call}
\end{align}
Under the {\it small noise asymptotics}:
\[
({\epsilon}\sqrt{n})^{-1} = O(1), \quad {\epsilon} \to 0, \quad n\to \infty,
\]
we can see that a minimum contrast estimator
\begin{align}
\wh{\vartheta}_n &= \arg\min_\vartheta M_n(\vartheta), \label{M-estimator}
\end{align}
with the contrast function
\begin{align*}
M_n(\vartheta) &=\sum_{k=1}^n  \left[\frac{n}{{\epsilon}^2}\frac{\left(\Delta_kX - \frac{1}{n} b(X_{t_{k-1}},\mu)\right)^2}{a^2(X_{t_{k-1}},\sigma)} + \log a^2(X_{t_{k-1}},\sigma)\right],
\end{align*}
is asymptotically normal:
\begin{align}
\left({\epsilon}^{-1} (\wh{\mu}_n - \mu_0),\sqrt{n}(\wh{\sigma}_n - \sigma_0)\right) \stackrel{d}{\longrightarrow} N_2\left(0, I^{-1}_{\vartheta_0}\right),\quad n\to \infty,
\label{asymptotics}
\end{align}
where
\[
I_\vartheta = \mathrm{diag}\left( \int_0^1 \left(\frac{\nabla_\mu b(X^{\vartheta,0}_s,\mu)}{a(X^{\vartheta,0}_s,\sigma)}\right)^2\mathrm{d} s, \
\frac{1}{2}\int_0^1 \left(\frac{\nabla_\sigma a^2(X^{\vartheta,0}_s,\sigma)}{a^2(X^{\vartheta,0}_s,\sigma)}\right)^2\,\mathrm{d} s\right);
\]
see S{\o}rensen and Uchida \cite{su03}. There is an another type of estimators as in Uchida \cite{u08}.
See also Long {\it et al.} \cite{letal13} and Shimizu \cite{s17} if $X$ is a jump process.

We would like to numerically demonstrate the result of Theorem \ref{thm:final} with
\[
\varphi_\vartheta(x) = e^{-rT}\max\{x - K,0\},
\]
and discretely observed diffusions as in \eqref{small-diff}, but $\varphi_\vartheta$ is not differentiable at $x=K$.
In order to apply the theorem, we can use the approximation as in Example \ref{ex:asian} by
\[
\varphi^\delta(x) = \frac{e^{-rT}}{2} \left(\sqrt{(x-K)^2 + \delta^2} + x - K\right).
\]
Setting $H_\delta(\vartheta) = \mathbb{E}[\varphi^\delta(X^\vartheta_T)]$, we have by Theorem \ref{thm:final} that
\begin{align}
{\epsilon}^{-1} ({\cal C}^{\wh{\vartheta}_n}_{T,K} - {\cal C}^{\vartheta_0}_{T,K}) \stackrel{d}{\longrightarrow} N(0,C_{\vartheta_0}^\top I^{-1}_{\vartheta_0} \,C_{\vartheta_0}), \label{normal}
\end{align}
under the asymptotics \eqref{asymptotics} as well as $\delta\to 0$,
where
\[
C_{\vartheta} = \mathbb{E}\left[\frac{Y^{\vartheta,0}_T}{2}\left\{\mathrm{sgn}(X^{\vartheta,0}_T - K) + 1\right\} \right],
\]
with $\mathrm{sgn}(z) = \mathbf{1}_{\{z>0\}} - \mathbf{1}_{\{z<0\}}$ and the derivative process $Y^{\vartheta,{\epsilon}}$.
We will compute $C_{\vartheta_0}$ by Monte Carlo simulations later.

In the next section, we shall numerically illustrate \eqref{normal} by a more concrete model.
\begin{remark}\label{rem:P*}
We claim that our setting above is a bit different from the financial practical problem since we do not care about the expectation in \eqref{E-call}, which should be the one with respect to a {\it risk neutral probability} although we should consider the parameter estimation under the physical probability measure ${\mathbb{P}}$ in financial problem. However, we shall try it just to numerically confirm our theoretical results.
\end{remark}

\subsection{The Black-Scholes model with small noise}\label{sec:BS}
For simulations, we shall consider the standard Black-Scholes assumption:
\begin{align}
X_t^{\vartheta_0,{\epsilon}} = x_0 + \int_0^t \mu_0 X_s^{\vartheta_0,{\epsilon}}\,\mathrm{d} s + {\epsilon} \int_0^t \sigma_0  X_s^{\vartheta_0,{\epsilon}}\,\mathrm{d} W_s. \label{BSmod}
\end{align}
We assume that $X^{\vartheta_0,{\epsilon}}$ is observed at $t_k = k/n\ (k=0,1,2,\dots)$, and set the true parameters
\begin{align}
\vartheta_0=(\mu_0,\sigma_0) = (0.2,1.0),\quad x_0=1.0,\quad {\epsilon} = 1/\sqrt{n} \label{truth}
\end{align}
for generating discrete samples from $X^{\vartheta_0,{\epsilon}}$. Then the minimum contrast estimator \eqref{M-estimator} is given in explicit form:
\begin{align}
\wh{\vartheta}_n = (\wh{\mu}_n,\wh{\sigma}^2_n) = \left(\sum_{k=1}^n \frac{\Delta_kX}{X_{t_{k-1}}},\
{\epsilon}^{-2}\sum_{k=1}^n \frac{(\Delta_k X - \frac{\wh{\mu}_n}{n}X_{t_{k-1}})^2}{X^2_{t_{k-1}}}\right) \label{estimators}
\end{align}
and the Fisher information matrix becomes
\[
I_{\vartheta_0} = \mathrm{diag}\left(\sigma_0^{-4}, 2\sigma_0^{-2}\right) = \mathrm{diag}(1,2).
\]

The derivative process $Y^\vartheta_t = (Y^1_t,Y^2_t)$ is given by
\begin{align*}
\mathrm{d} Y^{1,{\epsilon}}_t &= (\mu Y^{1,{\epsilon}}_t + X_t^{\vartheta,{\epsilon}})\,\mathrm{d} t + {\epsilon}\cdot \sigma Y^{1,{\epsilon}}_t\,\mathrm{d} W_t, \\
\mathrm{d} Y^{2,{\epsilon}}_t &= \mu Y^{2,{\epsilon}}_t\,\mathrm{d} t + {\epsilon} (\sigma Y^{2,{\epsilon}}_t + X_t^{\vartheta,{\epsilon}})\,\mathrm{d} W_t
\end{align*}
with $Y^{\vartheta,{\epsilon}}_0 = (0,0)$. A sample path of $(X^{\vartheta_0,{\epsilon}}, Y^{1,{\epsilon}}, Y^{2,{\epsilon}})$ is given in Figure \ref{fig:path}.

\begin{figure}
\begin{center}
\includegraphics[height=8cm, width=10cm]{paths.eps}
\caption{Sample path of $(X^{\vartheta_0,{\epsilon}}, Y^{1,{\epsilon}}, Y^{2,{\epsilon}})$ with $(\mu_0, \sigma_0, x_0)=(0.5,1.0,1.0)$ with ${\epsilon} =1$.}
\label{fig:path}
\end{center}
\end{figure}

We can find the explicit formula for ${\cal C}_{T,K}^{\vartheta_0}$, \eqref{E-call}
by using the well-known formula of the European call option price:
\[
{\cal C}_{T,K}^{\vartheta_0} = e^{-(r-\mu_0)T}\left[x \Phi(d_1) - K e^{-\mu_0 T}\Phi(d_2)\right],
\]
where $\Phi(x) = \frac{1}{\sqrt{2\pi}}\int_{-\infty}^x e^{-z^2/2}\,\mathrm{d} z$ and
\[
d_1 = \frac{\log (x/K) + ( r + {\epsilon}^2\sigma_0^2/2)T}{{\epsilon}\sigma\sqrt{T}};\quad
d_2 = d_1- {\epsilon}\sigma\sqrt{T}.
\]
One may notice that the above formula is a bit different from the usual Black-Scholes formula.
As is pointed out in Remark \ref{rem:P*},  we must be careful that the usual formula is under the {\it risk neutral probability}, under which the drift parameter $\mu_0$ corresponds to the interest rate.
Since we are now ignoring the risk neutral transform, we have to make a minor modification to the formula.

On the other hand, computing ${\cal C}_{T,K}^{\vartheta_0}$ by Monte Carlo simulations,
we see from Figure \ref{fig:C0} showing relative errors for Monte Carlo estimators, say ${\cal C}_{T,K}^*$, that 10,000 samples seem to be enough to compute ${\cal C}_{T,K}^{\vartheta_0}$.
Therefore, we also take 10,000 samples when we compute an estimator ${\cal C}_{T,K}^{\wh{\vartheta}_n}$ by Monte Carlo simulations, below.

\begin{figure}
\begin{center}
\includegraphics[height=8cm]{Rerror.eps}
\caption{The relative error $({\cal C}_{T,K}^* -{\cal C}_{T,K}^{\vartheta_0})/{\cal C}_{T,K}^{\vartheta_0}$, where ${\cal C}_{T,K}^*$ is a Monte Carlo estimator. The horizontal axis is the number of iteration for Monte Carlo simulations. }
\label{fig:C0}
\end{center}
\end{figure}

Numerical experiments are done by the following steps:
\begin{itemize}
\item[(1)] Generate a path of \eqref{BSmod}, and get discrete samples.
\item[(2)] Compute $\wh{\vartheta}_n$ in \eqref{estimators}, and generate paths of $X^{\wh{\vartheta}_n,{\epsilon}}$.
\item[(3)] Based on paths of (2), compute ${\cal C}_{T,K}^{\wh{\vartheta}_n}$ by a Monte Carlo method:
\[
{\cal C}_{T,K}^{\wh{\vartheta}_n} \approx \frac{1}{B}\sum_{k=1}^B e^{-r T}\max\{X_T^{\wh{\vartheta}_n,{\epsilon}}(k) - K,0\},
\]
where $X_T^{\wh{\vartheta}_n,{\epsilon}}(k)$ is the value of $X^{\wh{\vartheta}_n,{\epsilon}}_T$ starting from the initial value $X^{\wh{\vartheta}_n,{\epsilon}}_0=1.0$ for the $k$-th sample path.
\item[(4)] Compare $N(0,1)$ and the (estimated) distribution of
\begin{align}
\wh{Z}_n:={\epsilon}^{-1} ({\cal C}^{\wh{\vartheta}_n}_{T,K} - {\cal C}^{\vartheta_0}_{T,K})/\sqrt{C_{\vartheta_0}^\top I^{-1}_{\vartheta_0} \,C_{\vartheta_0}}
\label{normed.estimator}
\end{align}
for different sample sizes: $n=50, 100,$ and $300$.
\end{itemize}
In our experiments, we put
\[
T = 1.0,\quad K=0.75, \quad r = 0.05,\quad x=1.0.
\]
Then,  we had
\[
C_{\vartheta_0}=\begin{pmatrix}1.64937\\ 0.00585\end{pmatrix},
\]
by Monte Carlo simulations, and that the asymptotic variance in \eqref{normed.estimator} is
\[
C_{\vartheta_0}^\top I^{-1}_{\vartheta_0} \,C_{\vartheta_0} = 1.649396.
\]
We iterate the steps (1)--(4) 300 times, and show the histograms of $\wh{Z}_n$ and their estimated densities (by the kernel method, which were done by {\tt density()} in {\tt R}) as well as the normal QQ-plots in Figures \ref{fig:n=50} and \ref{fig:n=500}. Then, from their graphs, we can confirm that the asymptotic normality holds true even in the case where sample size $n$ is relatively small.

\begin{center}
\begin{figure}[htbp]
\begin{tabular}{c}
\begin{minipage}{0.5\hsize}
\includegraphics[width=7cm, height=5.5cm]{n=50.eps}
\end{minipage}
\begin{minipage}{0.5\hsize}
\includegraphics[width=7cm, height=5.5cm]{qq50.eps}
\end{minipage}
\end{tabular}
\caption{As $n=50$: Histogram of $\wh{Z}_n$ and its estimated density (blue line) as well as the standard normal density (red line) (left);
Normal QQ-plot (right) .}\label{fig:n=50}
\label{fig:theoretical}
\end{figure}
\end{center}

\begin{center}
\begin{figure}[htbp]
\begin{tabular}{c}
\begin{minipage}{0.5\hsize}
\includegraphics[width=7cm, height=5.5cm]{n=500.eps}
\end{minipage}
\begin{minipage}{0.5\hsize}
\includegraphics[width=7cm, height=5.5cm]{qq500.eps}
\end{minipage}
\end{tabular}
\caption{As $n=500$: Histogram of $\wh{Z}_n$ and its estimated density (blue line) as well as the standard normal density (red line) (left);
Normal QQ-plot (right).}\label{fig:n=500}
\label{fig:theoretical}
\end{figure}
\end{center}

\ \vspace{3mm}\\
\begin{flushleft}
{\bf\large Acknowledgements.}
The author expresses the sincere thanks to anonymous referees for detailed suggestions and proposals that makes the paper improve extensively.
This research was partially supported by JSPS KAKENHI Grant-in-Aid for Scientific Research (A) \#17H01100; (B) \#18H00836 and JST CREST \#PMJCR14D7, Japan.
\end{flushleft}