Small drift limit theorems for random walks
We show analogs of the classical arcsine theorem for the occupation time of a random walk in $(-\infty,0)$ in the case of a small positive drift. To study the asymptotic behavior of the total time spent in $(-\infty,0)$ we consider parametrized classes of random walks, where the convergence of the parameter to zero implies the convergence of the drift to zero. We begin with shift families, generated by a centered random walk by adding to each step a shift constant $a>0$ and then letting $a$ tend to zero. Then we study families of associated distributions. In all cases we arrive at the same limiting distribution, which is the distribution of the time spent below zero of a standard Brownian motion with drift 1. For shift families this is explained by a functional limit theorem. Using fluctuation-theoretic formulas we derive the generating function of the occupation time in closed form, which provides an alternative approach.In the course also give a new form of the first arcsine law for the Brownian motion with drift.
\section{Introduction}
For the classical symmetric random walk with $\pm 1$ steps it is well known that the three random variables ``time spent on the positive axis", ``position of the first maximum"
and ``last exit from zero'' are identically distributed and (suitably normalized)  asymptotically arcsine-distributed. 
Here the norming factor is the length of the time interval the random walk has been observed, so that the limiting statements refer to ``relative" times. 

Consider now a classical random walk with drift ${\delta}\neq 0$. Clearly the same ``relative" variables can be studied. The asymptotic distribution of the random variable 
``(fraction of) time spent in $(-\infty,{\alpha}]$
has been determined by Tak\'acs \cite{Tak}, by applying a functional limit theorem. 

But if ${\delta}\neq 0$ there is also another, ``absolute'' perspective. If for example ${\delta} > 0$ for a general random walk, 
it is clear that $Z({\delta})=$ ``number of visits in $(-\infty,0)$ ''  is almost surely finite, and that $Z({\delta})\longrightarrow \infty$ in probability as ${\delta}\searrow 0$.
One may ask  if $Z({\delta})$, after multiplication with some deterministic function $a({\delta} )$, has a non-degenerate limit distribution. 
This paper aims to answer these and related questions for random walks in the heavy-traffic regime, i.e., when the drift converges to zero. \color{black}  In all cases the limiting distribution for the occupation time in $(-\infty, 0)$, properly rescaled, turns out to have the density 
\begin{align}
\label{density} 
p(t)=2\,{\varphi(\sqrt{2t})\over \sqrt{2t}}-2\,\Phi(-\sqrt{2t}),\ \ t>0 
\end{align} 
where $\varphi$ and $\Phi$ are the density and the distribution function of $N(0,1)$, respectively. 

The distribution of the occupation time  in $(-\infty,0)$ \color{black} of Brownian motion with  positive \color{black} drift also has density (\ref{density}), and  in Section 2 
we begin with related results for Brownian motion. We  show \color{black} for example  that \color{black} the  distribution of the  time of the last exit from $0$ \color{black}
of Brownian motion with drift
during a finite time interval is composed of the arcsine and a truncated exponential distribution. 
In Section 2 we derive the limiting occupation time distribution for shift families generated from a centered random  walk by adding to each step a shift constant $a>0$ and then letting $a$ tend to zero. The proof that (\ref{density}) gives the asymptotic distribution is based on Donsker's invariance principle. 
In Section 4 we give the key fluctuation-theoretic formulas for the distribution for the occupation time in $(-\infty, 0)$ for general random walks. 


The arcsine law and its ramifications are a classical topic but there are always recent contributions, for example some new explicit distributions \cite{L}, new proofs 
\cite{Hof}, or asymptotic considerations \cite{Mar}. Interesting results on the number of visits to one point by skipfree random walks and related questions can be found in \cite{Ross}. The problem considered in this paper is also connected to the heavy traffic approximation problem in queueing theory, in which the growth of  the all-time maximum of $S_n-na$ (where $S_n$ is the $n$th partial sum of iid random variables with mean zero) is studied as $a\searrow 0$. In the queueing context this is equivalent to the growth of the steady-state waiting time in a $GI/G/1$ system when the traffic load tends to 1. 
This question was first posed by Kingman (see \cite{Kin}) and was investigated by many authors (e.g. \cite{BoxCoh, Kos, Pro, ResSam, Shnwa}). 

\section{Occupation times and last exit from $0$ \color{black} for Brownian motion with drift}
We start by presenting two results on occupation times for Brownian motion  with positive drift ${\delta}>0$ and variance ${\sigma^2}$, one known and one new. 
Let $B_t$ be a standard Brownian motion and $X_t=\sigma B_t+{\delta} t$.



\begin{lem}\label{lem}  
(1) Let $z>0$ and $T_z=\inf\{ t\geq 0\;:\;X_t\geq z\}$ be the first time when $X_t$ reaches level $z$.
Then $T_z$ has Laplace transform
$$\ell_{T_z}(s)=\mathbb{E} e^{-sT_z}=\exp\left(-{ z \over {\sigma^2}}(\sqrt{{\delta}^2+2{\sigma^2} s}-{\delta})\right).$$
(2) Let $V_0=V_0({\delta})=\int_{0}^\infty 1_{(-\infty,0)}(X_t)\,dt$ be the total time that $X_t$ spends below zero.
Then $V_0$ has Laplace transform
$$\ell_{V_0}(s)=\mathbb{E} e^{-sV_0}={ 2{\delta} \over {\delta} +\sqrt{{\delta}^2+2{\sigma^2} s}}.$$
\end{lem}
Proofs for (1) resp. (2) (for ${\sigma^2}=1$) can be found in \cite{KarT} resp. \cite{Imh}. 
Note $({\delta}^2/2{\sigma^2}) V_0$ has the Laplace transform $ 2/(1 +\sqrt{1+ s})$. We call $A$ a generic random variable with this  
Laplace transform. 

The density of $A$ is given by (\ref{density}). 
To see this, note that $1 / \sqrt{1+s}$ is the Laplace transform of the gamma distribution $\Gamma_{1,{1\over 2}}$, which has density 
$$\gamma_{1,{1\over 2}}(t)=1_{(0,\infty)}(t)\,{e^{-t} \over \sqrt{\pi t}}.$$
Therefore $[1-(1 / \sqrt{1+s})]/s$ is the Laplace transform of $1-\Gamma_{1,{1\over 2}}(t)=\int_t^\infty \gamma_{1,{1\over 2}}(x)\,dx$. The equality 
$${1 \over 1+\sqrt{1+s}}={1 \over \sqrt{1+s} }- {1 \over s}\Big(1-{1\over \sqrt{1+s}}\Big)$$
now yields density (\ref{density}). 

For $z\geq 0$ let $V_z = 
\int_{0}^\infty 1_{(-\infty,z)}(X_t)\,dt$ the total time the process spends below $z$.
Then the obvious decomposition (obtained by conditioning on $T_z$ ) $V_z=T_z+V_0^\prime$ (where $V_0^\prime $ is independent of  $T_z$ and dsitributed as $V_0$) yields

\begin{lem}\label{lem2} $V_z$ has Laplace transform
$$\ell_{V_z}(s)=\mathbb{E}(e^{-sV_z})=\ell_{T_z}(s)\,\ell_{V_0}(s).$$
\end{lem}
The density and distribution function are given in \cite{Imh}.\\

We focus in the sequel on the time spent on  the negative axis, but it is also of interest to  look at the other classical arcsine variable, i.e., the 
time of the last exit from $0$. Here we determine its distribution. \color{black}   
 Let ${\delta} \in\mathbb{R}\setminus\{0\}, {\sigma^2} =1$, so that $X_t=B_t+{\delta} t$, and consider  $W=\sup\{t\in [0,1]\,:\,X_t=0\}$, the last time
$X_t$ visits $0$ in $[0,1]$. 

Recall that for ${\delta}=0$, i.e., for the standard Brownian motion, the standard arcsine distribution (which has density $1_{(0,1)}(t)(1 / \sqrt{t(1-t)})$ 
and distribution function $(2 / \pi)\arcsin(\sqrt{t})$ on $[0,1]$) is the distribution of the last exit time from zero in the interval $[0,1]$. 

\color{black} The distribution of $W$ turns out to have a nice representation in terms of the standard arcsine distribution
and a truncated exponential distribution.  
As this result seems new, we provide a proof. 
 
\begin{thm}\label{thmle} $W\stackrel{d}{=}C\cdot \min \{1,D_{\delta}\}$ where $C$ and $D_{\delta}$ are independent, $C$ is arcsine-distributed, and $D_{\delta}$ is $\exp({\delta}^2/2)$-distributed. 
The moments of $W$ are given by 
$$\mathbb{E}W^k=  {2k \choose k}{1 \over 2^{2k}}\,\int_0^1 ky^{k-1} e^{-{\delta}^2y/2}\,dy, \  \ k\ge 1.$$
\end{thm} 


\begin{proof} We use a random walk approximation in the style of Tak\'acs \cite{Tak}.
Let $Y_1,Y_2,\ldots$ be iid with $$\mathbb{P}(Y_i=1)=p={1 \over 2}+{{\delta} \over 2\sqrt{n}}, \ \ \mathbb{P}(Y_i=-1)=q=1-p$$ 
($p$ and $q$ depend on $n$, but this is suppressed in the notation) and partial sums
$S_0=0,\,S_k=\sum_{i=1}^k Y_i$. 

It it easy to see that the  processes $X^{(n)}$ defined by 
$$X^{(n)}(t)={1\over \sqrt{n}} S_{\lfloor nt \rfloor}, \ \ 0\le t \le 1 $$ 
converge in distribution to $X=(X_t)_{t\in [ 0,1]}$ in $D[0,1]$. 

Furthermore, the last-exit time from $0$ is continuous in the Skorohod topology on $D[0,1]$ 
on a set of $P_X$-measure 1, and
 $$T_n=\sup\{t\in[0,1]\;:\;X^{(n)}(t)=0\} ={1\over n} \max\{0\leq k\leq n\;:\; S_k=0\} =:{M_n / n}\color{black}$$
 Then it suffices to show that ${M_N/ N}\longrightarrow C\cdot \min\{1,D_{\delta}\}$ as $N\longrightarrow \infty$. 
 
 Since ${1/ \sqrt{1-4pqz^2}}$ and $(\sqrt{1-4pqz^2})/(1-z)$ are the generating functions for the sequences of probabilities $\mathbb{P}(S_n=0)$ and $\mathbb{P}(S_1\not = 0,\ldots,S_n\not = 0)$, \color{black}
 respectively, 
the generating function of $M_N$ is 
\begin{align*}
\mathbb{E}t^{M_N}  &\ =\sum_{k=0}^N t^{k}\, \mathbb{P}(S_k=0, S_{k+1} \not=0,\ldots,S_{N}\not=0)\\
                   & =\sum_{k=0}^N t^{k}\, \mathbb{P}(S_k=0)\,\mathbb{P}(S_{1} \neq 0,\ldots,S_{N-k}\neq 0)\\ 
                   &=[z^{N}] {1\over \sqrt{1-4pqt^2z^2}}{\sqrt{1-4pqz^2} \over 1-z}\\
                   &=[z^{N}] {1\over \sqrt{1-4pqt^2z^2}}{\sqrt{1-4pqz^2} \over 1-z^2}(1+z). 
\end{align*}
(Here and in the following $[z^{N}] f(z)$ denotes the coefficient of $[z^{N}]$ in the Taylor expansion of the function $f(z)$ around zero.) 
Thus the generating functions for $N=2n+1$ and $N=2n$ are identical and it is enough to consider even $N$.
Let $N=2n$ be even (and  $n>{\delta}^2$) and $U_n={M_N / 2}$. 
Then the generating function of $U_n$ is
\begin{align*}
\mathbb{E}t^{U_n}&=[z^{2n}] {1\over \sqrt{1-4pqtz^2}}{\sqrt{1-4pqz^2} \over 1-z^2}\\
                   &=[z^{n}] {1\over \sqrt{1-4pqtz}}{\sqrt{1-4pqz} \over 1-z}
\end{align*}
so that the $k-$th factorial moment $u_{k,n}=\mathbb{E}\left(U_n(U_n-1)\cdots(U_n-k+1)\right)$ of $U_n$ is
given by 
\begin{align*}
u_{k,n}&= k! (-1)^k {-\frac{1}{2} \choose k}(4pq)^k\,[z^{n-k}]{1\over ({1-4pqz})^k\,(1-z)}\\
       &= k (-1)^k {-\frac{1}{2} \choose k}(4pq)^k\,[z^{n-k}]{1\over (1-z)}\int _{0}^{\infty} x^{k-1} e^{-(1-4pqz)x}\,dx\\
       &= (-1)^k {-\frac{1}{2} \choose k}(4pq)^k\,\int _{0}^{\infty} kx^{k-1} e^{-x}\left(\sum_{j=0}^{n-k}\frac{(4pqx)^j}{j!}\right)\,dx. 
\end{align*}
Now denote by $Poiss (\lambda)$ a random variable having  the Poisson  distribution with parameter $\lambda$. As $4pq=1-({{\delta}^2}/{2n})$, we obtain 
\begin{align*}
&\int _{0}^{\infty} k x^{k-1} e^{-x}\left(\sum_{j=0}^{n-k}\frac{(4pqx)^j}{j!}\right)\,dx\\
&=\int _{0}^{\infty} k x^{k-1} e^{-x(1-4pq)}\,\mathbb{P}\left(Poiss(4pqx)\leq n-k \right)\,dx\\
&=n^k\,\int _{0}^{\infty} k y^{k-1} e^{-{\delta}^2 y/2}\,\mathbb{P}\left(Poiss((n-\frac{{\delta}^2}{2})y)\leq n-k\right)\,dy. 
\end{align*}
By the central limit theorem,  
$$\mathbb{P}\left(Poiss((n-\frac{{\delta}^2}{2})y)\leq n-k\right)\longrightarrow\left\{
\begin{array}{cll}  &1  &\mbox{ for } 0\leq y<1\\ &\frac{1}{2} &\mbox{ for } y=1 \\ &0 &\mbox{ for } y>1
\end{array}\right.$$
so that for every $k$ we have 
$$\frac{u_{n,k}}{n^k}\longrightarrow (-1)^k {-\frac{1}{2} \choose k} \int_{0}^1 ky^{k-1}e^{-y{\delta}^2/2}\,dy.$$ 
 Hence ${\mathbb{E}T_N^k}/{N^k}$
 tends to the same limit. This shows the second assertion. Finally,  $$\mathbb{E} C^k = (-1)^k {-\frac{1}{2} \choose k}={2k \choose k} {1 \over 2^{2k}}$$  and 
integration by parts shows that $\int_{0}^1 ky^{k-1}e^{-y{\delta}^2/2}\,dy=\mathbb{E} \min\{1,D_{\delta}^k \}$.
Thus all moments of ${T_N / N}$ converge to the corresponding moments of $C\cdot\min\{1,D_{\delta} \}$. Since the distribution of 
$C\cdot\min\{1,D_{\delta} \}$ is clearly determined by its moments the first assertion follows.
\end{proof} 

\begin{rem} 
As an immediate consequence of the scaling properties of Brownian motion we see that  
the distribution of 
$$W_T=\sup\{ t\leq T\,:\, \sigma B_t+{\delta} t=0\}$$ is the same as that of $C\cdot\min\{T,D_{{\delta}/\sigma} \}$. 
The time of the last zero of $\;\sigma B_t+{\delta} t\;$ in the interval
$[0,\infty)$ is thus distributed as $C\cdot D_{{\delta}/\sigma} $, which is the gamma distribution with parameters ${{\delta}^2 / 2\sigma^2}$ and ${1}/{2}$.
\hfill\vrule height5pt width5pt depth0pt 
\end{rem}

\begin{rem} Clearly $V_0$ (the occupation time on the negative axis) is stochastically smaller than $W_\infty$ (the last exit time from zero), and the results above quantify this
precisely. We find e.g. that $$\mathbb{E}(V_0)=\dfrac {{\sigma^2}}{2{\delta}^2}=\dfrac{1}{2}\,\mathbb{E}(W_\infty).$$  
\hfill\vrule height5pt width5pt depth0pt 
\end{rem}
\color{black}

\begin{rem} Last-exit times of Brownian motion from moving boundaries have been studied intensively,
and  more complicated  expressions for the density of the last-exit time from a linear boundary were derived in \cite{Sal} and \cite{Imh2}. The representation in (\ref{thmle})  appears to be new, as it
is not mentioned in the encyclopedic monograph \cite{BorS}. For the density of the sojourn time found by Tak\'acs by a random walk limit two ``purely Brownian" explanations have 
been given in \cite{DonY}. It is natural to ask for such an explanation for the representation in (\ref{thmle}).
  \hfill\vrule height5pt width5pt depth0pt  
\end{rem}


\section{Limit of occupation times for shifted random walks}


In this section we consider a shifted random walk. Specifically, let $(X_{{\delta},1},X_{{\delta},2},\ldots)$ be a  parametrized sequence
of iid random variables with $\mathbb{E}(X_{{\delta},i})=0$, ${\mathrm{Var}}(X_{{\delta},i})={\sigma^2}({\delta}) \in (0, \infty)$. Let ${\delta}>0$ and 
$$ Y_i^{\delta} =X_{{\delta},i}+{\delta}, \ \ 
  S_{{\delta},n}=\sum_{i=1}^n X_{{\delta},i}, \ \ 
 S_n^{\delta}=\sum_{i=1}^n Y_i^{\delta}.$$ 
We are interested in the occupation time 
  $$Z_0^{\delta}=\sum_{i=1}^\infty 1_{(-\infty,0)}(S_n^{\delta}).$$ 
Throughout this section we assume that ${\sigma^2}({\delta})\longrightarrow {\sigma^2}>0$ as ${\delta}\longrightarrow 0$ and that 
 the following Lindeberg-type condition holds:  for every ${\varepsilon}>0$,
\begin{equation}\label{lind}
\lim_{{\delta}\longrightarrow 0} \int_{|{\delta} X_{{\delta},1}|>{\varepsilon}} X_{{\delta},1}^2\,d\mathbb{P}=0.
\end{equation} 
These conditions are chosen such that for the triangular array with the variables 
$$Z_{{\delta},k}=\dfrac{{\delta}}{\sigma({\delta})}X_{{\delta},k}, \ \ k=1,\ldots,\lfloor\dfrac{1}{{\delta}^2}\rfloor $$ 
the central limit theorem holds: indeed, ${\mathrm{Var}}(Z_{{\delta},1})={\delta}^2$ and the Lindeberg condition for this triangular array reads as 
$$\lim_{{\delta}\longrightarrow 0} \frac{1}{{\delta}^2}\int_{|Z_{{\delta},1}|>{\varepsilon}{\delta}^2\lfloor\frac{1}{{\delta}^2}\rfloor} Z_{{\delta},1}^2\,d\mathbb{P}=
\lim_{{\delta}\longrightarrow 0} \frac{1}{{\sigma^2}({\delta}))}\int_{|{\delta} X_{{\delta},1}|>{\varepsilon} \sigma({\delta}) {\delta}^2\lfloor\frac{1}{{\delta}^2}\rfloor} X_{{\delta},1}^2\,d\mathbb{P}=
0 \ \  \text{for every} \  {\varepsilon}>0,   
$$ 
which is clearly true under the conditions above. 

We use similar ideas as Prohorov \cite{Pro}, who proved the following:   

\begin{thm} (Prohorov) In the situation above let
$M^{\delta} =\min\{ S_n^{\delta}\;:\;n\geq 0\}$. Then 
$$\mathbb{P}({\delta} M^{\delta}>x)\longrightarrow e^{-2x/{\sigma^2}} \ \mbox{ for all } x>0.$$ 
\end{thm} 
In \cite{Pro} 
 the maximum in the case of negative drift was considered instead of $M^{\delta}$. The result had been proved earlier 
by Kingman under the assumption of the existence of an exponential moment.

The following lemma will be needed to obtain tightness bounds.
\begin{lem}
 In the situation above let $z\geq 0$ and let ${\delta}_k>0$ be a sequence of positive numbers satisfying $\sup_{k\geq 1} {\sigma^2}({\delta}_k)< \infty$.
 Then for every ${\varepsilon}>0$ we can find a $T$ such that for all $k$
$$\mathbb{P}(\sup_{n\geq T/{\delta}_k^2} (|S_{{\delta}_k,n}|-n{\delta}_k)\geq -\frac{z}{{\delta}_k})<{\varepsilon}.$$ 
\end{lem}
 \begin{proof} First consider a sequence $S_n$ of partial sums
of an arbitrary iid sequence $(X_i)$ with $\mathbb{E}(X_1)=0$ and ${\mathrm{Var}}(X_1)={\sigma^2}$. Let $a,b>0,Na>b$ and consider the event $E_N =\{\sup_{n\geq N} (|S_n|-na)\geq -b\}$. 
Clearly \begin{align*}E_N&=\bigcup_{j=0}^\infty\left\{\max_{2^jN\leq n< 2^{j+1}N} (|S_n|-na)\geq -b\right\}\\
                         &\subseteq \bigcup_{j=0}^\infty\left\{\max_{2^jN\leq n< 2^{j+1}N}|S_n|\geq 2^jNa- b \right\}\\
                         &\subseteq \bigcup_{j=0}^\infty\left\{\max_{n\leq 2^{j+1}N}|S_n|\geq 2^jNa- b \right\}.\\
\end{align*}
By Kolmogorov's inequality,  
$$\mathbb{P}(\max_{n\leq 2^{j+1} N}|S_n|\geq 2^jNa-b)\leq \frac{2^{j+1} N{\sigma^2}}{(2^jNa-b)^2}. $$
Now set $N=\frac{T}{{\delta}^2}, a={\delta},  b=\frac{z}{{\delta}}, X_i=X_{{\delta},i}$. It follows that 
$$\mathbb{P}(\sup_{n\geq T/{\delta}^2} (|S_{{\delta},n}|-n{\delta})\geq -\frac{z}{{\delta}})\leq\sum_{j=0}^{\infty}\frac{2^{j+1} T {\sigma^2}({\delta})}{(2^jT-z)^2}.$$
The bound on the right side depends on ${\delta}$ only via ${\sigma^2}({\delta})$ and can clearly be made arbitrarily small (under the assumptions above).
\end{proof}

\begin{cor}\label{cor1}
 In the situation above let $z\geq 0$ and let ${\delta}_k>0$ be a sequence of positive numbers satisfying $\sup_{k\geq 1} {\sigma^2}({\delta}_k)<\infty$. 
Then for every ${\varepsilon}>0$ one can find a $T$ such that for all $k$
$$\mathbb{P}(\min_{n\geq T/{\delta}_k^2} {\delta}_k(S_{{\delta}_k,n}+n{\delta}_k)\leq z)<{\varepsilon}. $$ 
\end{cor}

\begin{thm}\label{thmconv}
 $${{\delta}^2 \over 2{\sigma^2}({\delta})\color{black}}Z_0^{\delta} \longrightarrow A \ \mbox{ in distribution} \mbox{ as }{\delta}\searrow 0. $$
\end{thm}
\begin{proof}
By the remark following \ref{lem} it suffices to show that ${\delta}^2 Z_0^{\delta}\longrightarrow V_0$ in distribution, where $V_0$ is the distribution of the time
the process $W_t=\sigma B_t+t$ spends below zero. 

Let $T>0$ and consider the sequence of processes $$U^{\delta}(t)={\delta}\,\sum_{i=1}^{\lfloor t/{\delta}^2\rfloor}Y_i^{\delta}, \;\;\;0\leq t \leq T.$$
By Donsker's limit theorem (in the version for triangular arrays, see e.g. \cite{Bil}, p.147), 
the sequence $U^{\delta}\longrightarrow \sigma B + id$ in distribution in $D[0,T]$, where $\sigma B + id$
denotes the Brownian motion with variance ${\sigma^2}$ and drift $1$, i.e., with coordinate variables $\sigma B_t +t$. For any bounded Borel function
$v$ on $[0,T]$ the functional 
$x\mapsto \int_{0}^T v(x_t)  \,dt$ on $D[0,T]$ is Skorohod-measurable and continuous except on a set of $B$-measure $0$ (see e.g. \cite{Bil}, p. 247).
Thus, 
\begin{align*}
{\delta}^2 \mbox{card}(\{n\; :\;  S_n^{\delta} < 0, 1\leq n \leq T/{\delta}^2\})&=\int_0^{{\delta}^2\lfloor T/{\delta}^2\rfloor} 1_{(-\infty,0)}(U^{\delta}(t))\,dt\\
 & \longrightarrow \int_0^T 1_{(-\infty,0)}(X_t)\,dt \mbox{ as }{\delta}\searrow 0\,
\end{align*} 
in distribution and we will be done if we can justify the interchange of the limits $T\longrightarrow \infty$ and $ {\delta} \searrow 0$.
Let ${\delta}_k > 0$ be a sequence decreasing to zero and let ${\varepsilon}>0$. By corollary \ref{cor1} 
we can find an $N$ such that $\mathbb{P}(\min_{n\geq N/{\delta}_k^2} S_n^{{\delta}_k}\leq 0)<{\varepsilon}$ for all $k$.
Thus,  
\begin{equation}\label{lim}\lim_{T\longrightarrow \infty} \sup_{k\ge 1}\mathbb{P}(\min_{n\geq T/{\delta}_k^2} S_n^{{\delta}_k}\leq 0)=0
\end{equation} 
and the assertion follows since, by the monotone convergence theorem,  
$$ \lim_{T\longrightarrow \infty}\int_0^T 1_{(-\infty,0)}(X_t)\,dt = \int_0^\infty 1_{(-\infty,0)}(X_t)\,dt.$$ 
\end{proof}




\begin{rem}
A related discussion can be found in \cite{Shnwa}. In that paper, Shneer and Wachtel derived an extension of Kolmogorov's inequality and treated the maximum of random walks with negative drift and step size distributions attracted to a stable law of index ${\alpha}\in(1,2]$.
In the case of finite variance (${\alpha}=2$) they already remarked that their results (including in particular the crucial relation \eqref{lim}) remain valid 
  if the conditions assumed above hold.
\hfill\vrule height5pt width5pt depth0pt 
\end{rem} 

\begin{rem} 
Assume that the $X_i$ are independent with $\mathbb{E}(X_i)=0$ and variances ${\mathrm{Var}}(X_i)=\sigma_i^2$
and satisfy Lindeberg's condition. Let $s_i^2=\sum_{k=1}^i \sigma_k^2$. Then the step processes
$X_n(t)$ which jump to the value ${S_i / s_n}$ at time ${s_i^2 / s_n^2}$ converge weakly to
a standard Brownian in $D[0,1]$ (by Prohorov's extension of Donsker's theorem). One may thus expect that they exhibit a similar limiting behavior.
\hfill\vrule height5pt width5pt depth0pt 
\end{rem} 

 Finally, replacing $0$ by $z/{\delta}$ and repeating the steps in the proof of \ref{thmconv} yields 

\begin{thm}\label{thmz}
In the situation above let $z>0$ and $Z_z^{\delta}=\sum_{n=1}^\infty 1_{(-\infty,z)}(S^{\delta}_n)$ . Then 
${\delta}^2 Z^{\delta}_{z/{\delta}} \longrightarrow V_z$ in distribution, where the Laplace transform of $V_z$ is given in Lemma \ref{lem2}  
with ${\delta}=1$. 
\end{thm} 

 If here $z$ depends on ${\delta}$ such that  ${\delta} z({\delta})\longrightarrow 0$ as ${\delta}\longrightarrow 0$ we find \color{black}

\begin{prp} In the situation above let $(z({\delta}))$ a sequence of positive numbers with $z({\delta})=\mathrm{o}(1/{\delta})$ and $\sup_{\delta} z({\delta})<\infty$. Then 
	$${\delta}^2 Z^{\delta}_{z({\delta})}\longrightarrow V_0=2{\sigma^2} A \mbox{  as  } {\delta} \longrightarrow 0.$$
\end{prp}
\begin{proof} Clearly 
$V_0$ is stochastically smaller than any distributional limit of ${\delta}^2 Z^{\delta}_{z({\delta})}$ 
(because $Z_0^{\delta}$ is stochastically smaller than $Z_y^{\delta}$ for $y\geq 0$),
furthermore $V_y=T_y+V_0$ is stochastically smaller than $V_z$ for $y\leq z$. 
Let ${\varepsilon}>0$ and $C=\sup_{\delta} z({\delta})$, then $C<\infty$ and ${\delta}^2 Z^{\delta}_{{\varepsilon} C/{\delta}}\longrightarrow V_{{\varepsilon} C}$ in distribution as ${\delta} \longrightarrow 0$ (by theorem \ref{thmz}).
Since $Z^{\delta}_{z({\delta})}=Z^{\delta}_{{\varepsilon} z({\delta})/{\varepsilon}}$ is stochastically smaller than $Z^{\delta}_{{\varepsilon} C/{\delta}}$ for ${\delta}\leq {\varepsilon}$, 
any distributional limit of ${\delta}^2 Z^{\delta}_{z({\delta})}$ is stochastically smaller than $V_{{\varepsilon} C}$. Thus the distributional limit
exists and equals  $V_0$.
\end{proof}


We close this section with an application of Theorem \ref{thmconv} 
in a frequently encountered situation.

\begin{exl}\label{exl1} (Expectation shift in exponential families.)\\
Let $U$ be a non-constant real random variable such that the moment generating function $$m(s)=\mathbb{E} e^{sU}$$ is finite in an open interval $I$ around  $0$, and $E(U)=m^\prime(0)=0$,
$\mathrm{Var}(Y)={\sigma^2}$. 
 
For $p\in I$ let $U_p$ have the ``associated'' distribution with moment generating function $m_p(s)={ m(p+s) \over m(p)}$, clearly $U_p$ has expectation 
$\mathbb{E}(U_p)={ m^\prime(p) \over m(p)}$ and variance ${\sigma^2}(p)={ m^{\prime\prime}(p)m(p)- (m^\prime(p))^2 \over m(p)^2}$.\\ Let $Z_0(p)$ denote the random variable
``time spent in $(-\infty,0)$'' by the random walk generated by iid variables with distribution $U_p$. Then 
$${(\mathbb{E}(U_p))^2 \over 2 {\sigma^2}(p)}\,Z_0(p)\longrightarrow  A  \mbox{ in distribution for } p\searrow 0 \;\;\;.$$
\begin{proof}
It is well known that $s\mapsto \log m(s)$ is strictly convex on $I$, thus $p\mapsto {m^\prime(p) \over m(p)}=\mathbb{E}(U_p)$ is strictly increasing. Thus we may parameterize the
distributions by ${\delta}(p)= \mathbb{E}(U_p)$. We have  ${\delta}(p)\searrow 0$ for $p \searrow 0$ and 
${\sigma^2}(p)\longrightarrow {\sigma^2}$ as $p\searrow 0$. Let $X_{{\delta}(p)}=U_p-\mathbb{E}(U_p)$ and $Y^{{\delta}(p)}=X_{{\delta}(p)}+{\delta}(p)=U_p$.
Then the Lindeberg condition \eqref{lind} is satisfied, since by Chebyshev's inequality   
$$\int_{|{\delta}(p) X_{{\delta}(p)}|>{\varepsilon}} X_{{\delta}(p)}^2\,d\mathbb{P}\leq {{\delta}^2(p){\sigma^2}(p)\over {\varepsilon}^2}$$
and the claim follows from Theorem \ref{thmconv}. 
\end{proof}



\end{exl} 



\color{black}

\section{The fluctuation theoretic approach}

The topics investigated here belong to the fluctuation theory of random
walks. We recall some basic facts, which will be used in 
the sequel and can e.g. be found in Section XII.7 of \cite{Fel2}. 

We consider a random walk $(S_n)_{n\ge 1}$, i.e., a sequence of partial sums of iid random variables 
and let $R=\inf\{ n\geq 1\;:\; S_n<0\}$ and  $W=\inf\{n \geq 1\;:\; S_n\geq 0\}$ be the lengths of the 
first strictly descending and weakly ascending ladder epochs of the random walk, 
respectively.We denote by $r(z)$ and $a(z)$
denote the corresponding  probability generating functions  and set $\mu = \mathbb{E} W $.   
The occupation time of interest is $Z_0=\sum_{n=1}^\infty 1_{(-\infty, 0)}(S_n)$. 
\begin{thm}\label{thmSpA}(Sparre Andersen)
For $|z|<1$
\begin{align*}
{1 \over 1-r(z)}&=\exp \left\{\sum_{n=1}^\infty {z^n \over n} \mathbb{P}(S_n<0)\right\} \mbox{  }\\ 
{1 \over 1-a(z)}&=\exp \left\{\sum_{n=1}^\infty {z^n \over n} \mathbb{P}(S_n \ge 0)\right\} \mbox{  }\end{align*}
\end{thm}

An immediate consequence is the factorization theorem. 
\begin{thm}(Duality) 
For $|z|<1$
\begin{align*}
(1-r(z))(1-a(z))=1-z. 
\end{align*}
\end{thm}

It follows from the factorization theorem is that $W$ ($R$) has a finite expected value if and only if $R$ ($W$) is defective, and that
the relations $\mathbb{E}(R)\mathbb{P}(W=\infty)=1$ and $\mathbb{E}(W)\mathbb{P}(R=\infty)=1$ hold. 


At the combinatorial heart of fluctuation theory is the ``Sparre Andersen transformation'' (made explicit by Feller and 
refined by Bizley and Joseph) given in Lemma 3 of XII.8 of \cite{Fel2}:

\begin{lem} Let $x_1,\ldots,x_n$ be real numbers with exactly $k\geq 0$ negative partial sums
$s_{i_1},\ldots,s_{i_k}$, where $i_1>\ldots> i_k$. Write down $x_{i_1},...,x_{i_k}$ followed
by the remaining $x_i$ in their original order. (If $k=0$, the sequence remains unchanged).  The transformation thus defined is invertible, and
the first (absolute) minimum of the partial sums of the new arrangement  
occurs at the $k$-th place.
\end{lem}

Clearly this extends to {\it infinite sequences} with exactly $k$ negative partial sums: just apply the bijection above 
to an initial section large enough to contain all the negative partial sums, and leave the rest unchanged.


The following formulas express 
the generating function of $Z_0$ in terms of $r(z)$ or of $a(z)$, respectively. 


\begin{thm
\begin{align}\label{asc}
\mathbb{E} z^{Z_0}={1-r(1) \over 1-r(z)}  = {1 \over \mu}{ 1-a(z) \over 1-z}
 =\exp\left\{-\sum_{k=1}^\infty (1-z^k)\frac{\mathbb{P}(S_k<0)}{k} \right\}. 
\end{align} 
\end{thm}



\begin{proof}
According to Lemma 2.3, each sequence $x_1,x_2\ldots$ with exactly $k$ negative partial sums there corresponds (by a finite reordering) a unique sequence with first (absolute) minimum  at the $k$th
place. The partial sums $s_0=0,s_1,s_2,\ldots$ of the rearranged sequence consist of a first part $s_0,s_1,\ldots,s_k$ 
and a second part $s_{k+1},s_{k+2},\ldots $ such that  the partials sums satisfy $s_i>s_k$ for $i\leq k$ and $s_i-s_k\geq 0$ for $i>k$. 
For a random walk the joint distribution of the $X_i$ is invariant under finite permutations, and the two parts are independent. The first part has probability $$\mathbb{P}(0>S_k,S_1>S_k,\ldots,S_{k-1}>S_k)=\mathbb{P}(S_1<0,\ldots,S_k<0)$$ 
(by reversing the order of the variables), the second part has  probability 
$$\mathbb{P}(S_{k+1}-S_k \geq 0,S_{k+2}-S_k\geq 0,\ldots)=\mathbb{P}(S_1\geq 0,S_2\geq 0,\ldots)=1-r(1).$$
This yields the first equation of (\ref{asc}). 
The second   
one follows immediately from the factorization identity $(1-a(z))(1-r(z))=1-z$ (recall Theorem 4.2) and the third one from Sparre Andersen's theorem. 
\end{proof}

In some cases $r(z)$ can be computed in closed form, and the asymptotics of $Z_0 $ can be obtained from an explicit formula. 
An example is the normal random walk.
Let the iid steps $X_i$ be $N({\delta}, {\sigma^2})$-distributed. Here we only assume  
that ${\delta} \neq 0$, i.e., we consider the cases of positive and negative ${\delta}$ simultaneously and let $d:=|{\delta}|, q:=\frac{{\delta}^2}{2{\sigma^2}}$. 
\begin{exl} For the normal random walk we have\\
(a) $ r(z)=1-(1-z)^{{1 \over 2}}\,\exp \left(\mbox{ sign}({\delta})  {d^2 \over \pi{\sigma^2}} {\displaystyle \int_{0}^{1}\int_0^{\infty} \,{ e^{-d^2(y^2+x^2) / {2{\sigma^2}}} \over 1-ze^{-d^2(y^2+x^2) / {2{\sigma^2}}}}}\,dy\,dx \right). $ 

(b) $qZ_0\longrightarrow A$ in distribution as ${\delta}^2/{\sigma^2}\searrow 0,{\delta} \searrow 0$.  

(c) $r(e^{-qs})^{1/\sqrt{q}}\longrightarrow e^{-(\sqrt{1+s} -1)}$ as $q \searrow 0,{\delta} \nearrow 0$.
\end{exl}
Note that here ${\sigma^2}$ may vary with ${\delta}$, it is only essential that ${\delta}/\sigma \longrightarrow 0$.
\begin{proof} 
Directly from Sparre Andersen's theorem we find that 
\begin{align*}
\log\left({1 \over 1-r(z)}\right)&=\sum_{n=1}^\infty {z^n \over n} \mathbb{P}(S_n<0) 
                                    =\sum_{n=1}^\infty {z^n \over n} \int_{-\infty}^{-n{\delta}} {1 \over \sqrt{2n\pi{\sigma^2}}}\, e^{-x^2 / {2n{\sigma^2}}}\,dx
\\ &                                    =\sum_{n=1}^\infty {z^n \over n}\left({1 \over 2} -   \mbox{sign}({\delta})\, \int_{0}^{n d} {1 \over \sqrt{2n\pi{\sigma^2}}}\, e^{-x^2 / {2n{\sigma^2}}}\,dx\right).
\end{align*} 
Hence, 
\begin{align}\label{normal}
 1-r(z)=(1-z)^{{1 \over 2}}\,\exp(\mbox{ sign}({\delta}) G(z)),
 \end{align} 
where 
$$G(z)=\sum_{n=1}^\infty {z^n \over n}\, \int_{0}^{nd} {1 \over \sqrt{2n\pi{\sigma^2}}}\, e^{{-x^2 / 2n{\sigma^2}}}\,dx.$$
We have 
\begin{align*}
\int_{0}^{nd} {1 \over \sqrt{2n\pi{\sigma^2}}}\, e^{-x^2 / {2n{\sigma^2}}}\,dx&= \int_{0}^{d} \sqrt{ {n \over 2\pi{\sigma^2}}}\, e^{-ny^2 / {2{\sigma^2}}}\,dy
 =
{n \over \pi{\sigma^2}} \int_{0}^{d}\int_0^{\infty} \, e^{-n(y^2+x^2) /{2{\sigma^2}}}\,dy\,dx
\\ 
&= {nd^2 \over \pi{\sigma^2}} \int_{0}^{1}\int_0^{\infty} \, e^{-nd^2(y^2+x^2) / {2{\sigma^2}}}\,dy\,dx 
\end{align*}
and therefore  
$$G(z)={d^2 \over \pi{\sigma^2}} \int_{0}^{1}\int_0^{\infty} \,{ e^{-d^2(y^2+x^2) / {2{\sigma^2}}} \over 1-ze^{-d^2(y^2+x^2) / {2{\sigma^2}}}}\,dy\,dx, $$
proving (a). Note that $G(z)$ depends only on the ratio $q={d^2/ 2{\sigma^2}}$. Fix $s>0$. Setting $z=e^{-qs}$ we obtain for $q\searrow 0$ (by dominated convergence):
\begin{align*} 
G(e^{-qs})&={2\over \pi} \int_{0}^{1}\int_0^{\infty} \,{ q e^{-q(y^2+x^2) } \over 1-e^{-q(s+y^2+x^2) }}\,dy\,dx \\ & 
 \longrightarrow 
{2\over \pi} \int_{0}^{1}\int_0^{\infty} \,{ 1 \over s+y^2+x^2}\,dy\,dx 
=\log \left({1+\sqrt{1+s}\over \sqrt{s}}\right). 
\end{align*}
From this (b) and (c) follow easily. 
\end{proof} 
It is of methodological interest to have also a purely fluctuation-theoretic proof of Theorem 3.4., i.e., a proof which does not rely on the ``functional limit theorem'' approach used above. \color{black}  The reviewer suggested the following alternative derivation of \ref{thmconv} 
based on Theorem 4.4.   Assume the conditions introduced in Section 3.  
\begin{thm} = {\bf Theorem 3.4}\color{black}  
 $${{\delta}^2 \over 2{\sigma^2}({\delta})\color{black}}Z_0^{\delta} \longrightarrow A \ \mbox{ in distribution} \mbox{ as }{\delta}\searrow 0. $$
\end{thm}
\begin{proof}
In principle, we follow the line of argument used for a similar proof in \cite{Shnwa}. Let ${\varepsilon}>0$ and split the series
in the exponent of the right-hand side of \eqref{asc} into three parts:
$$\sum_{k=1}^\infty =\sum_{k=1}^{{\varepsilon}/{\delta}^2}\; + \; \sum_{{\varepsilon}/{\delta}^2}^{T/{\delta}^2}\; + \;\sum_{T/{\delta}^2}^\infty=\sum\nolimits_1 +\sum\nolimits_2 +\sum\nolimits_3\;.$$
Let $s>0$ and set $$z=e^{-s{\delta}^2/2{\sigma^2}({\delta})}.$$ 
We consider the different sums separately, starting with  $\sum\nolimits_1$:
$$\sum_{k=0}^{{\varepsilon}/{\delta}^2} (1-z^k) \frac{\mathbb{P}(S^{\delta}_k<0)}{k}\leq \frac{s{\delta}^2}{2{\sigma^2}({\delta})} \sum_{k=0}^{{\varepsilon}/{\delta}^2} \mathbb{P}(S^{\delta}_k<0)\leq \frac{s{\varepsilon}}{2{\sigma^2}({\delta})}.$$
Furthermore, $\mathbb{P}(S_k^{\delta}<0)= \mathbb{P}(\sum_{j=1}^k X_{{\delta},j}<-k{\delta})\leq {\sigma^2}({\delta})/(k{\delta}^2)$ by Chebyshev's inequality.
 Therefore we  obtain for ${\varepsilon}>{\delta}^2$ 
$$\sum_{k\geq {\varepsilon}/{\delta}^2}\frac{(1-z^k)}{k}\mathbb{P}(S^{\delta}_k<0) \leq \frac{{\sigma^2}({\delta})}{{\delta}^2}\sum_{k\geq {\varepsilon}/{\delta}^2}\frac{1}{k^2} \leq\frac{{\sigma^2}({\delta})}{{\delta}^2} \int_{{\varepsilon}/{\delta}^2}^\infty \frac{1}{(x-1)^2}\,dx=\frac{{\sigma^2}({\delta})}{{\varepsilon}-{\delta}^2}.$$
Since $\sigma({\delta})\longrightarrow {\sigma^2}\in(0,\infty)$ as ${\delta}\longrightarrow 0$, there is a ${\delta}_0$ such that $2{\delta}_0^2<{\varepsilon}$ and ${\sigma^2}({\delta})$ is bounded for ${\delta}\leq {\delta}_0$. 
Without loss of generality assume in the sequel ${\delta}\leq{\delta}_0$. 
Then $\sum\nolimits_3$ can be made arbitrarily small by a suitable choice of $T$, and $\sum\nolimits_2\leq 2C/{\varepsilon}$ for a suitable constant $C$. 

For $\sum\nolimits_2$ we use the asymptotic normality of ${\delta} S^{\delta}_{t/{\delta}^2}$ (which is implied by the Lindeberg condition, see  the beginning of section 3)\color{black}:
$$\mathbb{P}({\delta} S_k^d<0) \longrightarrow \mathbb{P}(N(t,{\sigma^2} t)<0)=\Phi(-\sqrt{\frac{t}{{\sigma^2}}})\; \mbox{ as }\; {\delta} \longrightarrow 0,\,k{\delta}^2\longrightarrow t$$
(uniformly for $t\in [{\varepsilon},T]$), and by the dominated convergence we conclude that 
 $$\sum\nolimits_2 \longrightarrow \int_{\varepsilon}^T \frac{1-e^{-t/2{\sigma^2}}}{t}\, \Phi(-\sqrt{\frac{t}{{\sigma^2}}})\,dt.$$
 Letting ${\varepsilon}\longrightarrow 0, T\longrightarrow \infty$ we finally arrive at
 \begin{align} 
\label{domi} 
\mathbb{E} e^{-s\frac{{\delta}^2}{2{\sigma^2}}}\longrightarrow \exp\{-\int_0^\infty \frac{1-e^{-su}}{u}\,\Phi(-\sqrt{2u})\,du\}.
\end{align} 
Evaluating the integral finishes the proof. Avoiding the calculation, it suffices to notice that the right side of (\ref{domi})  
is independent of the underlying 
distribution of the random walk so that one can look at the example of the normal random walk computed above, which leads to the 
conclusion that the right side of (\ref{domi})is equal to $2/(1+\sqrt{1+s})$. 
\end{proof}  

The advantage of this proof is that it generalizes to the ${\alpha}$-stable case ($1<{\alpha}<2$) essentially unchanged - the main difficulties (the corresponding  estimates
for these cases) can be overcome using inequality $(6)$ in \cite{Shnwa}. 

 We close this section with a few remarks on the simple random walk \color{black}
taking step +1 with probability $p>1/2$ and step -1 with  probability $q=1-p$. It is well-known that in this example 
$$r(z)={1 -\sqrt{1-4pqz^2} \over 2pz},$$ 
so that a quick calculation shows that 
$$ \mathbb{E} z^{Z_0}={1-r(1) \over 1-r(z)} \color{black}
= \dfrac{(p-q)(1+\sqrt{ 1-4pqz^2})}{p(1-2z^2 +\sqrt{ 1-4pqz^2})}$$  

and  
$2(p-\frac{1}{2})^2Z_0 \longrightarrow A$ in distribution as $p\searrow 1/2$. 



\begin{rem}  
 Let $T_0(p)=\sup\{n\ge 0\,:\,S_n^{(0)}=0\}$ the time of the last return to the origin. 
In the symmetric case $p=1/2$  the walk is persistent and  $T_0(1/2)=\infty$ almost surely.  In the transient case $p>1/2$,  
 $T_0(p)$ has generating function 
$$h(z)={ p-q \over \sqrt{1-4pqz^2}}.$$
A short computation yields that ${1 \over 2}(p-q)^2T_0(p)$ converges in distribution as $p\searrow 1/2$, the limiting distribution 	
having Laplace transform $\dfrac{1}{\sqrt{ 1+s}}$, i.e., being the $\Gamma_{1,{1 \over 2}}$ distribution with density $\gamma_{1,{1 \over 2}}(t)$ as above.
\hfill\vrule height5pt width5pt depth0pt 
\end{rem}
  \begin{rem}     
Let $N_0(p)$ denote the number of zeros of the random walk. Then
$$\mathbb{P}(N_0(p)=r,\,T_0(p)=2n)={r \over n-r}{2n- r \choose n} 2\,(pq)^n$$
and $(\delta N_0(p),{1 \over 2}\delta^2 T_0(p))$ converges weakly to the distribution with density
$$f(y,t)=1_{(0,\infty)}(y)\,1_{(0,\infty)}(t)\,{y \over 2t}{1 \over \sqrt{2\pi t}}\, e^{-({y^2 / 4t})\,-t}\;\;.$$
In particular,  $\delta N_0(p)$ is asymptotically $\exp(1)$. 
 For the symmetric random walk let $N_{0, 2n}$ denote the number of zeros up to time $2n$. A classical theorem of Chung-Hunt \cite{ChH} states that 
$\sqrt{{2 / n}} N_{0, 2n} $  is asymptotically distributed as $
|N(0,1)|$. All these results show that deviations from the symmetric random walk become clearly visible after $n\approx \delta^{-2}$ steps. 
 While characteristics like the positive sojourn time and the last exit time from zero are in both cases of approximately the same size their distributions differ. 
For the last exit time from zero a precise description is given in Theorem 2.3.  \color{black}. 
\hfill\vrule height5pt width5pt depth0pt 
\end{rem}  

Apparently the distribution of $A$ occurs naturally as a limit of occupation times for random walks with drift.
It is well-known (see e.g. Section XIV.3 in \cite{Fel2}) that the deeper reason for the frequent occurrence of the (generalized) arcsine distributions 
lies in their intimate connection to distribution functions with regularly varying tails. The same explanation applies here.
 In the case of drift zero the distribution functions of the ladder epochs are attracted to the standard positive stable distribution
of index 1/2 and the positive (negative) sojourn times are asymptotically arcsine-distributed. In the cases with small drift (and finite variance) the ladder epochs are
attracted to an associated distribution of this stable distribution, and therefore the positive (negative) sojourn times have asymptotically the distribution of $A$. 

{\bf Acknowledgement}. We would like to thank the referee for valuable remarks and in particular for suggesting the alternative proof
of Theorem \ref{thmconv} given in Section 4.
\color{black}





2016-05-31T02:06:15.000
1601
1601.00609
en
https://arxiv.org/abs/1601.00609
We show analogs of the classical arcsine theorem for the occupation time of a random walk in $(-\infty,0)$ in the case of a small positive drift. To study the asymptotic behavior of the total time spent in $(-\infty,0)$ we consider parametrized classes of random walks, where the convergence of the parameter to zero implies the convergence of the drift to zero. We begin with shift families, generated by a centered random walk by adding to each step a shift constant $a>0$ and then letting $a$ tend to zero. Then we study families of associated distributions. In all cases we arrive at the same limiting distribution, which is the distribution of the time spent below zero of a standard Brownian motion with drift 1. For shift families this is explained by a functional limit theorem. Using fluctuation-theoretic formulas we derive the generating function of the occupation time in closed form, which provides an alternative approach.In the course also give a new form of the first arcsine law for the Brownian motion with drift.
Probability (math.PR)
Small drift limit theorems for random walks
Qwen/Qwen-72B
1. YES
2. YES