とある病んだ院生の体内時計（サーカディアンリズム）StanのToDoリストが長すぎて若干引き気味のberoberoです。 BUGS/Stan勉強会#2の発表した内容の一つを紹介をしたいと思います。 スライドは埋め込んで、ソースコードのコピペ＆解説をメインにします。 使用したデータは以下。自由に使ってください。 ・元データ: data.txt ・起床時刻だけ抜き出したもの: data_sleepout.txt ・起床時刻から起床時刻の時間: data_out2out.txt ・起床時刻の累積時刻: data_cum.txt まずはじめに起床時間の分布のパラメータを求めたStanコードは以下になります。 7行目・11行目、Stanでvon Mises分布を使うときはsigmaをある程度小さくしましょう。 kickするRコードは以下になります。 9行目: von Mises分布に渡すために0-24を0-2πに変換しています。 次に起床時刻から起
