#########################################
#  第13章 時変分散と時変係数ARモデル(TSSSパッケージ)
#########################################
#
# 地震波データの入力
#
data(MYE1F) # an earthquake wave data 
#
# 時変分散モデルの推定とデータの等分散化
# トレンド次数=2，システムノイズ分散の初期値と探索幅
#
tvvar(MYE1F, 2, 6.6e-06, 1.0e-06)


#　地震データの入力
#
data(MYE1F) # an earthquake wave data 
#
#　時変係数ARモデルの推定
#　AR次数=4，AR係数トレンド次数=2，
#
data(MYE1F) # an earthquake wave data 
z <- tvar(MYE1F, 4, 2, 20,tau20=6.6e-06, delta=1.0e-06) 
z$tau2 
z$sigma2 
z$lkhood 
z$aic

#
#　n=630と1026の2か所で構造変化があったと仮定
#
z <- tvar(MYE1F, 4, 2, 20, c(630,1026), 
　　　6.6e-06, 1.0e-06) 
z$tau2 
z$sigma2 
z$lkhood 
z$aic


#
#　時変係数ARモデルから時変スペクトルを計算
#　時変スペクトルを3次元表示
#
spec <- tvspc(z$arcoef, z$sigma2,, 20) 
persp(spec$x, spec$y, spec$z, expand=0.5, theta=20, 
　　col="lightblue", ticktype="detailed", xlab="f", ylab="n", 
　　zlab="log p(f)")