rstanで最高事後密度区間（HDPI）を出したい

stanは通常、事後期待値が計算されてくる。でも、HDPI最高事後密度区間（HDI最高密度区間）がいいらしいのでrstanでだしてみたい。
HDIntervalパッケージで簡単にできるらしい。詳細は以下の通り。ちなみにfitをextractしとかないと（前回のブログ参照）fitの中はいじれません。