|
,N,(4)一阶条件(4)是通过对对数似然函数方程(3)中ai的微分得到的。非线性方程组H(a)=0,可使用标准数值算法求解。此后,a的计算解(即最大似然估计)由a表示*≡ (a)*, . . . , 一*N) 。通过逐段重复这个过程,我们得到了活动向量SNA的估计值*, 一*, . . . , 一*tmaxo公司。3.3关系借出的统计检验我们提出了两种统计检验;一个用于识别参与关系借贷的银行对,另一个用于检测依赖关系的银行。在samemanner中,当我们估计银行的活动水平时(第3.2节),我们将每日数据集拆分为TMAX聚合时段,并逐时段实施测试。3.3.1基于边缘的关系借贷测试如果i银行对贸易伙伴没有偏好,因此按照能力模型中的建议以随机方式找到合作伙伴,那么i银行和j银行在给定时期内的双边交易数量mij应服从以下二项分布:g(mij | a*i、 a*j)=τmiju(a*i、 a*j) mij(1- u(a*i、 a*j) )τ-mij, i、 j=1,N、 (5)相比之下,如果i银行与j银行有较强的(即非随机)合作关系,则MIJ的分布将偏离二项分布。设McIj表示第c百分位(0≤ c≤ 100)g(mij | a*i、 a*j) (即c/100=G(mcij | a*i、 a*j) ,其中G是G(mij | a)的累积分布函数(CDF)*i、 a*j) 。如果c值接近100时,mij>MCIJ,则交易的经验数量太大,无法用随机机会来解释,这表明我们使用基于改进牛顿法(称为信赖域折线法)的Matlab函数fsolve解决了问题。
|