![起床時刻の分布は？
von Mises分布 （フォン?ミーゼス分布）
??? ???? ? ? ?
? ?; ?, ? =
2??0 ?
ただし， ?? ?
0 ≤ ? < 2?
?>0
0 ≤ ? < 2?
は第1種変形Bessel関数
[wikipediaより]](https://image.slidesharecdn.com/random-131221232219-phpapp01/85/-9-320.jpg)
More Related Content
More from . . (11)
とある病んだ院生の体内时计（サーカディアンリズム）
- 4. day wday 1 2 3 4 4 5 6 … 591 591 592 592 593 sleep.in 6 0 1 2 2 3 4 … 1 1 2 2 3 9 8.67 3 6 17.33 10.67 2 … 0 22 0 21.33 0 sleep.out 17 18.17 11.87 14.75 22.08 13.75 4.67 … 7.17 24 6.83 24 6.5 注: ここは一続き 曜日 0: 日 1: 月 … 6: 土 ? 2003.3.17～2004.10.27 ? 593日間のデータ, 674回寝てる
- 9. 起床時刻の分布は？ von Mises分布 （フォン?ミーゼス分布） ??? ???? ? ? ? ? ?; ?, ? = 2??0 ? ただし， ?? ? 0 ≤ ? < 2? ?>0 0 ≤ ? < 2? は第1種変形Bessel関数 [wikipediaより]
- 10. 厂迟补苍コード sigmaの上限が大きいと 実行時にエラーが出ます. Error in function boost::math::cyl_bessel_i<d>(d,d): numeric overflow
- 12. 結果 0 18 + 6 12 ? ?????? = 11.03 ? ?????? = 3.28 plotには {circular}パッケージを使いました. 95%??: 10.55 ?? 11.52 95%??: 2.83 ?? 3.71
- 15. 厂迟补苍コード
- 18. 厂迟补苍コード 変化点を表すパラメーター こんな感じで入れる. cf. The BUGS Book 11.7.1
- 19. 結果 ? ? ?????? = 22.5 ? ?????? ?????? = 20.5 ? ?????? = 21.4 = 21.4 日 時刻 収束はしたけどもっと範囲を狭めて iterを多くした方がいいかも.
- 22. 厂迟补苍コード 範囲が各々異なる潜在変数zの作り方 cf. Stan v2.0.1のマニュアル 49.2-49.4 ARモデル
- 25. まとめ ? Stanで周期データを分析してみた. von Mises分布もある！ ? 前処理につぐ前処理が意外と近道. ? 規則正しい生活をしましょう！
- 26. Web上の情報（予定含む） ? 生態学会大会 (2014) の自由集会 – 角度や時間など循環する変数の統計モデリング – http://hosho.ees.hokudai.ac.jp/~kubo/ce/EcoSj2014.html ? 石原氏による「角度データの統計処理基礎? 」 – http://www.qbio.jp/wiki/%E7%AC%AC%E5%9B%9B%E5%9B%9E%E5%B9%B4%E4%BC% 9A_%E3%83%81%E3%83%A5%E3%83%BC%E3%83%88%E3%83%AA%E3% 82%A2%E3%83%AB_%E7%B5%B1%E8%A8%88 – http://qbio.jp/images/5/53/%E8%A7%92%E5%BA%A6%E7%B5%B1%E8%A8%8 8%E9%85%8D%E5%B8%83_qbio4th.pdf ? 睡眠関係のナショジオの連載記事 – http://nationalgeographic.jp/nng/article/20121127/332007/