风险度量VaR (Value at risk) 是指在正常的市场环境下, 金融资产或证劵组合在一定的持有期内和一定的置信水平下, 预期的最大损失称为风险价值或在险价值, 其定义一般表述为设$\{Y_t\}_{t=0}^n$是$n$个时间段资产价格序列, $X_t=-\log({Y_t}/{Y_{t-1}})$是第$t$个时间段对数收益.假设$\{X_t\}_{t=1}^n$是相依严平稳序列, 其边缘分布函数为$F(x)=P(X_t\leq x),\ x\in R.$对任意给定$p\in (0,1),$ $F(x)$的$p$-分位数定义为$\xi_p=F^{-1}(p)=\inf\{u: F(u)\geq p\},$其中函数$F^{-1}(t)$, $0﹤t﹤1$是$F(\cdot)$的广义逆函数.易知$\xi_p$满足$F(\xi_p-)\leq p\leq F(\xi_p)$.置信水平为$1-p$的风险度量VaR值定义为$ v_p=-\xi_p=-\inf\{u: F(u)\geq p\}.$
在许多金融模型中, 风险度量VaR方法能够较准确地度量由不同风险来源及其相互作用而产生的潜在损失, 而且与其它风险度量方法相比, VaR方法的计算也相对简便.自1990年以来, 国内外学者对风险度量VaR方法进行了一定程度的研究.如: Dowd [1]给出的风险度量VaR样本分位数估计为$Q_{n,p}=-X_{([np]+1)},$其中$X_{(r)}$为样本$X_1, \cdots, X_n,$的第$r$个次序统计量. Koji et al [2]提出了用样本分位数$-X_{([np])}$估计VaR.谢佳利和杨善朝等[3]采用加权样本分位数估计量估计VaR, 该方法的估计精度相对文献[1, 2]有较好的改进.
记$F_n(x)=\frac{1}{n}\sum\limits_{i=1}^{n}I(X_i\leq x)$为总体$X$的经验分布函数, 其样本分位数记为$Z_{n,p}$.显然$Z_{n,p}=-Q_{n,p}$.由于风险度量VaR与分位数只相差一个负号, 所以分位数的各种非参数估计自然也是VaR的非参数估计.而对于样本分位数估计的研究, 自Bahadur[4]给出了独立样本下分位数估计的Bahadur表示之后, 许多学者对混合和负相协 (NA) 以及负象限相依 (NOD) 样本下研究了分位数估计的Bahadur表示及渐近性质对此进行了研究, 也取得了不少成果, 这里就不一一列出.
定义1.1 称随机变量$\{X_1, \cdots, X_n, n\geq 2\}$是正相协 (PA) 的, 设$A_1$与$A_2$为$\{1,\cdots ,n\}$任何两个不相交的非空子集, 如果$ {\rm Cov}(f_1(X_i,i\in A_1),f_2(X_j,j\in A_2))\geq 0$成立, 其中$f_1$和$f_2$是任何两个使得协方差存在且对每个变元均非降 (或对每个变元均非升) 的函数.称随机变量序列$\{X_i,i\geq1\}$是PA序列, 如果对任何$n\geq 2, X_1,\cdots ,X_n$是PA的.
PA序列的一个基本性质: PA序列经过单调函数变换得到的序列仍为PA序列.
对于NA序列的Rosenthal型矩不等式和Bernstein型概率不等式的研究已经有很好的成果, 这为研究负相协下统计大样本性质提供了很好的研究工具.然而, PA序列受到其协方差系数影响, 比较理想的Rosenthal型矩不等式和Bernstein型概率不等式目前比较少见, 这给研究PA下估计量的统计大样本性质带来一定程度的不便.
由于PA序列包括正相关的正态随机变量族,而且PA序列在与实际应用有关的模型中(如可靠性理论, 渗透理论及多元统计分析等) 有广泛的应用, 所以受到了学者们的广泛关注, 具体文献在此就不一一列举.但是PA样本下风险度量VaR分位数估计的大样本性质及其Bahadur表示的研究尚未见文献涉及.因此本文将采用文献[5]中PA序列指数不等式以及PA序列协方差性质, 讨论PA序列下风险度量VaR分位数估计的强相合性和渐近正态性及其Bahadur表示.所得结果推广了相关文献的结果.
下面我们给出PA序列指数不等式以及一些辅助结论.文中$C, C_1, C_2, \cdots$表示不依赖于$n$的正常数, 且在不同的地方取值可能相同也可能不同.
引理2.1[5] (i) 设$X_1,X_2,\cdots $是零均值正相协序列, $\max\limits_{1\leq i\leq n}|X_i|\leq c_n <\infty {\rm a.s.}\ (n=1,2,\cdots)$;
(ii) 设$u(n)=\sup\limits_{i\geq 1}\sum\limits_{j:j-i\geq n}{\rm cov}(X_i,X_j)$, 满足$\sum\limits_{i=1}^\infty u^{1/2}(2^i)<\infty$;
(iii) 设$\{p_n:n\geq 1\}$是一串正整数序列, 满足$p_n\leq \frac{n}{2}$成立.又设$\theta >0$, 满足$0<\theta p_nc_n\leq 1$.则存在一个不依赖于$n$的正常数$C_1$, 使得对任意的$\varepsilon >0$, 有
引理2.2[6] 设$F(x)$是右连续分布函数, 则广义逆函数$F^{-1}(t)$, 在$0<t<1$非降且左连续, 并且满足
(i) $F^{-1}(F(x))\leq x,-\infty <x<+\infty;$
(ii) $F(F^{-1}(t))\geq t,0<t<1;$
(iii) $F(x)\geq t\Longleftrightarrow x\geq F^{-1}(t).$
引理2.3[7] 令$p\in (0,1),\ Z_{p,n}=F_n^{-1}(p)=\inf \{x:F_n(x)\geq p\}$.假设$P(X_i=X_j)=0,i\neq j$, 那么$p<F_n(Z_{p,n})<p+1/n\ \ {\rm a.s.}$.
引理2.4 (i) 如果$\{X_n,n\geq 1\}$具有相同分布函数$F(x)$和有界密度函数$f(x)$的严平稳PA序列, $F(x)$在$\xi_p$的邻域$N_p$内连续可导 ($p\in(0,1)$), 且$0<d=\sup\{f(x):x\in N_p\}<\infty$;
(ii) 设$v(n)=\sup\limits_{i\geq 1}\sum\limits_{j:j-i\geq n}{\rm cov}^{1/3}(X_i,X_j)$满足$v(n)=O(e^{-2\sqrt{n}})$;
(iii) 对任意满足$d_n\rightarrow 0, n^\frac{1}{4}d_n/\log n \rightarrow \infty,n\rightarrow \infty$的正实数序列$\{d_n\}_{n\geq 1}$, 记$D_n=[\xi_p-d_n, \xi_p+d_n]$.则
注2.1 如果$v(n)=O(e^{-2\sqrt{n}})$, 则必有
事实上, 由于$\sum\limits_{n=1}^\infty v^{1/2}(2^n)\leq C\sum\limits_{i=1}^\infty e^{-2^{n/2}}$, 易知当$n>N$有$2^{n/2}>n$.故$\sum\limits_{n=N}^\infty e^{-2^{n/2}} <\infty, $从而$\sum\limits_{n=1}^\infty v^{1/2}(2^n)<\infty$.
引理2.4的证明 设$t_n=d_nn^{-\frac{1}{4}}, S_{r,n}=\xi_p+rt_n,$
由于
当$x\in [\xi_p+rt_n,\xi_p+(r+1)t_n]$时, 由$F_n(x)$是非降函数以及微分中值定理得
和
根据 (2.1) 和 (2.2) 式得
而
令
由PA序列性质知$\eta_{n1},\cdots,\eta_{nn}$仍是平稳PA随机变量, 且$E\eta_{ni}=0,\ |\eta_{ni}|\leq 2\stackrel {\Delta}{=}c_n,\ i=1,\cdots,n$, 故引理2.1中的条件 (i) 满足.又令$w(n)=\sup\limits_{i\geq 1}\sum\limits_{j:j-i\geq n} {\rm Cov}(\eta_{ni},\eta_{nj})$, 根据引理2.4条件 (i) 及Roussas [8]中的引理2.6 (也可见文献[9]的 (4) 式) 知, 存在某正常数$M$有
再根据引理2.4的条件 (ii) 以及注2.1得$w(n)\leq v(n)<\infty$, 且$\sum\limits _{i=1}^\infty w^{1/2}(2^i)<\infty,$从而引理2.1条件 (ii) 满足.又令$\theta =\frac{1}{\sqrt{n}}>0,\ p_n\leq \frac{\sqrt{n}}{2}\leq n/2$, 则$0<\theta p_nc_n\leq 1 $, 故引理2.1条件 (iii) 满足.由此, 根据引理2.1可得
类似于$I_1$的计算方法, 得
故由 (2.4), (2.7) 和 (2.8) 式得
从而$\sum\limits_{n=1}^\infty P(\max\limits_{-m_n\leq r\leq m_n}|\Delta_{r,n}|>t_n)<\infty.$故由Borel-Cantelli引理及 (2.3) 式得
引理2.4得证.
下面给出本文的主要结果.
定理3.1 (1) 设$\{X_n,n\geq 1\}$是具有相同分布函数$F(x)$和有界密度函数$f(x)$的严平稳PA序列, $F(x)$在$\xi_p$处可导且$F'(\xi_p)=f(\xi_p)>0$;
(2) 设$v(n)=\sup\limits_{i\geq 1}\sum\limits_{j:j-i\geq n}{\rm Cov}^{1/3}(X_i,X_j)$满足$v(n)=O(e^{-2\sqrt{n}})$;
(3) 设$\{d_n\}_{n\geq 1}$为满足$d_n\rightarrow 0,n^\frac{1}{4}d_n/\log n\rightarrow \infty,\ n\rightarrow \infty$的任意正实数序列.则
定理3.2 设$\{X_n,n\geq 1\}$是具有相同分布函数$F(x)$和有界密度函数$f(x)$的严平稳PA序列, $F(x)$在$\xi_p$的邻域$N_p$内连续可导 ($p\in(0,1)$), 且$0<d=\sup\{f(x):x\in N_p\}<\infty$.如果定理3.1的条件 (2)-(3) 满足, 则
定理3.3 在定理3.2的条件下, 如果$\forall x,y\in R$, $j \geq 1$, 有
其中$F_{1,j+1}(x,y)$为$(X_1, X_j)$的联合分布函数.令$\sigma_0^2=F(x)[1-F(x)]+\sum\limits_{j=1}^{\infty}E(Z_{n1}Z_{nj+1})$, 记$\stackrel {d}{\to} $为依分布收敛.则
定理3.1的证明 对于任意的$\varepsilon >0$有
由引理2.2可得
其中$w_i=I(X_i>\xi_p+\varepsilon n^{-\frac{1}{4}}d_n),\ \delta_1=F(\xi_p+\varepsilon n^{-\frac{1}{4}}d_n)-p$.同理可得
其中$v_i=I(X_i\leq \xi_p-\varepsilon n^{-\frac{1}{4}}d_n), \ \delta_2=p-F(\xi_p-\varepsilon n^{-\frac{1}{4}}d_n).$
类似于 (2.5) 式关于$\{\eta_{ni}, i\geq 1\}$的性质讨论, 知$\{w_i-Ew_i\}_{1\leq i\leq n}$和$\{v_i-Ev_i\}_{1\leq i\leq n}$都为0均值平稳PA序列, 且满足$|w_i-Ew_i|\leq 2\stackrel {\Delta}{=}c_n$, $|v_i-Ev_i|\leq 2\stackrel {\Delta}{=}c_n.$令$\theta =1/\sqrt{n}>0$, $p_n\leq \sqrt{n}/2$, 有$\theta p_nc_n\leq 1.$从而由引理2.1及 (4.1)-(4.3) 式可得
其中$c,c_1$为常数.又因为$F(x)$在$\xi_p$点处连续, $F'(\xi_p)>0$, 此处$\xi_p$是不等式$F(x-)\leq p\leq F(x)$的唯一解且$F(\xi_p)=p$, 由Taylor展开得
从而$\min(\delta_1,\delta_2)=f(\xi_p)\cdot \varepsilon n^{-\frac{1}{4}}d_n$, ($n\rightarrow \infty$), 代入 (4.4) 式得
于是
故由Borel-cantelli引理得$v_p-Q_{n,p}=o(n^{-\frac{1}{4}}d_n) \ \ {\rm a.s.}. $定理得证.
定理3.2的证明 由定理3.1知
由引理2.4得
又由引理2.3得
令$\theta_n$是介于$Q_{p,n}$与$v_p$之间的随机变量, 结合 (4.5)-(4.7) 式, 利用Taylor展开得
由此可得
定理3.2证明完毕.
定理3.3的证明 由定理3.2得$v_p-Q_{n,p}=(p-F_n(\xi_p))/f(\xi_p)+O(d_nn^{-\frac{1}{4}})\ \ {\rm a.s.}. $故要证明定理3.3成立, 只需证明$F_n(\xi_p)-p$具有渐近正态性.记$Z_{ni}=I_{(X_i\leq \xi_p)}-EI_{(X_i\leq \xi_p)}$, 则
令$r_{n2}<r_{n1}$为正整数, $r_{n2}/r_{n1}\to 0$, 且当$n\to \infty$时, $r_{n1}/n\to 0$成立.又取$k=[n/(r_{n1}+r_{n2}]$, 满足$k(r_{n1}+r_{n2})/n\to 1$, $k r_{n2}/n\to 0$, 且$r_{n1}^2/n\to 0$.利用Bernstein分块方法, $S_n$可分解为
其中
$k_m=(m-1)(r_{n1}+r_{n2})+1,\ \ l_m=(m-1)(r_{n1}+r_{n2})+r_{n1}+1,\ \ m=1,\cdots, k. $则在定理的条件下, 由文献Li and Yang[10]的引理3.4, 得到
再由文献Li et al.[11]的引理A.3及文献Li and Yang[10]的引理4.2和4.3得
结合 (4.9)-(4.11) 式可得
由此知$F_n(\xi_p)-p \stackrel {d}{\to} N (0,\sigma_0^2)$.这样定理3.3证明完毕.