#########################################
#  第8章 局所定常ARモデルの推定
#########################################
# テストデータの読み込み
  data(MYE1F)

# LSAR1　任意個の区間への自動分割
#
  lsar1(MYE1F, 10, 100)

#LSAR2　変化時点の精密な推定
#
# P波到着時刻の推定
#　
  lsar2(MYE1F, 10, c(400,800), c(600,700))
#
# S波到着時刻の推定
#　
  lsar2(MYE1F, 10, c(800,1200), c(1000,1100))


#
# 変化点の事後確率の計算
#
# P波の到着時刻分布
#
x <- lsar2(MYE1F, 10, c(400,800), c(600,700))
AICP <- x$aic
post <- exp( -(AICP-min(AICP))/2 )
post <- post/sum(post)
plot(post,type="l",col="red",lwd=2)


#
# S波の到着時刻分布
#
x <- lsar2(MYE1F, 10, c(800,1200), c(1000,1100))
AICP <- x$aic
post <- exp( -(AICP-min(AICP))/2 )
post <- post/sum(post)
plot(post,type="l",col="red",lwd=2)