|
这两个区域之间的转换由标准的马尔可夫转换模型描述,转换矩阵为Q。观测方程[rE]t=gsi+Lt1- Lt(gsi)- rt+TT~ N(0,σ)), (33)其中t∈ Z+和si=si(t)∈ {s,s}表示时间t的状态。状态方程slt+1=Lt+(Tt- Lt)gsi- rtLt+a(1)- Lt)1- Tt+k(1)- Lt)Ttt、 (34)Tt+1=Tt+kTt(Tt- Lt)(1)- (Tt)t、 (35)gsi=(cin s,cin s,(马尔可夫转换模型)。(36)先验σ~ IG(10-2, 10-2) ,c,c~ U(-0.25, 0.25).(IG表示逆伽马分布,U表示均匀分布)。转移率矩阵:Q=-λ λu -u, λ, u ~ U(01100)。初始条件:L~ U(0.2,0.3),T~ U(0.3,0.4)。参数值:a=0.05,k=0.05,t=0.1。我们取样200次,选择100个老化样品(即丢弃前100个样品)。选择tc,c,λ,u。逆伽马分布通常用作噪声项方差的先验。仅对and施加限制,以确保所有变量保持有限且T>L。控制时间尺度的参数(a、k、,t) 被选为常数,所有变量都是有限的。时间步长t=0.1(年)相当于大约一个月,与我们每月提供一个数据点(每年12个数据点)的事实一致。只要KandareMain为正值且不会导致变量不确定,我们本可以选择估算模型参数A和Kas。然而,估算和K(例如~ U(0,10)和k~ U(0,1))将导致杠杆和信任的动态变得有些多余,因为它可以有效地减少用于仅拟合等式(33)和等式(36)的等式系统,而RTS作为数据和常数提供。为了防止杠杆和信任很快达到一个稳定的状态,我们给ANDK赋值。
|