|
,n∧由使用Gumbel copula和拒绝采样方法进行的校准产生,如表1所示,尺寸d=2、5和25。k 1 2 3 4 5 6 7 8 9 10xk0。000 0.500 0.750 0.875 0.937 0.968 0.984 0.992 0.996 0.998d=2 pk0。100 0.000 0.000 0.000 0.115 0.325 0.206 0.128 0.787 0.048d=5 pk0。100 0.000 0.000 0.000 0.129 0.302 0.202 0.131 0.084 0.053d=25 pk0。100 0.000 0.000 0.000 0.022 0.252 0.216 0.174 0.135 0.102表1:使用Gumbel copula校准的概率权重pk。图2中绘制了权重函数ew(·)=w(·1)和5000个V样本的散点图,其中参考copula为Gumbel,并使用了拒绝采样方法。考虑到本案例研究的设置,很容易检查引理3.2是否满足要求。由于FV的构造,与copula样本相比,更多样本位于上边界或右边界附近。由于目标函数使用S分布函数的估计,我们将样本权重归一化为和1。这进一步减少了Geweke(2005)第4.2.2节或Liu(2008)第2.5.2节中提出的估算误差。为了评估重要性抽样带来的改进,我们对d∈ {2, 5, 25}. 我们使用n=10000的样本量来计算两种算法中每一种算法的重要性抽样器bu,以及所有目标函数的标准蒙特卡罗估计值u。总共500次重复被用来获得这些估计器的经验分布,从而估计它们的方差。虽然样本大小对样本方差的值有影响,但在算法效率的研究中,它不应扮演重要角色。表2、3给出了使用拒绝算法的Gumbel和Clayton情况的结果,表4、5给出了使用直接算法的Gumbel和Clayton情况的结果。
|