高松市の最大降水量って増えてる？

『Stan と R でベイズ統計モデリング』読書会@関西学院大学Osaka.Stan #6でのLT資料です。一般化極値分布をStanで定義して，状態空間モデルを使って最大降水量のモデリングをしてみました。

Transcript

1. ߴদࢢͷ࠷େ߱ਫྔͬͯ
   ૿͑ͯΔʁ ಙԬɹେ
   ʢNUPLVPLB!UBLBNBUTVVBDKQʣ
   ʰ4UBO
   ͱ
   3
   ͰϕΠζ౷ܭϞσϦϯάʱಡॻձ
   
   !ؔ੢ֶӃେֶ
   0TBLB4UBO
   
   
   ೥݄೔
   
2. ࠓ೔ͷ͓͠ͳ͕͖ w ࣗݾ঺հ
   w ࠷େ߱ਫྔͷ࣌ܥྻมԽͷਪఆ
 ʢ̍֊ࠩ෼Ͱ͔͠΍ͬͯ·ͤΜʣ
   w ΋͏গ͠΍Γ͔ͨͬͨ͜ͱ 2
3. ࣗݾ঺հ w ໊લɿಙԬɹେʢͱ͓͔͘ɹ·͞Δʣ
   w ॴଐɿߴদେֶൃୡՊֶ෦
   w ઐ໳ɿڭҭ৺ཧֶɼ઀ۙճආಈػ͚ͮɼଞऀࢤ޲త ಈػ͚ͮ
   w ͦͷଞؔ৺ɿ৺ཧ౷ܭɼϕΠζ౷ܭɼֶڃن໛
   w 5XJUUFSɿɹɹɹ!U@NBDZB
   w ;ͭ͏ͷڭҭ৺ཧֶऀ
   w ϓϩάϥϛϯάॳ৺ऀʂ 3
4. Ϟνϕʔγϣϯ w ঢ়ଶۭؒϞσϧΛ࢖ͬͯΈ͍ͨ
   w ҰൠԽۃ஋෼෍Λ࢖ͬͯΈ͍ͨ
   w XFCʹ͋ΔσʔλΛ࢖ͬͯΈ͍ͨ
   w ࡢ೥ͷ4UBO
   "EWFOU
   $BMFOEFS
   ͕ෆ׬શ೩ম 4
5. ࠷େ߱ਫྔ w ؾ৅ிͰ೥͝ͱͷ࠷େ߱ਫྔΛެදͯ͘͠Ε͍ͯΔ 5
6. ࠷େ߱ਫྔͷ෼෍ •ਖ਼ن෼෍ͯ͠ͳ͍ؾ͕͢Δ •࠷େ஋ͳͷͰҰൠԽۃ஋෼෍ ʹै͏ؾ͕͢Δ •ҰൠԽۃ஋෼෍(GEV)͸3ύϥ
 gevmxt ( mean, sigma, xi ) •GEV͸Stanʹ࣮૷͞Εͯͳ͍ 6
7. ͜ΕΛ࣮૷ •OsakaStan#5ͷLTʹͯ͜ͷ͋ͨΓ͸࣮૷ࡁΈ 7
8. ͜Μͳײ͡ 8
9. ࠷େ߱ਫྔͷ࣌ܥྻมԽ w ී௨ʹϓϩοτ͢Δͱ͜Μͳײ͡
   w ͔͔ͤͬͩ͘Βঢ়ଶۭؒϞσϧΛ࢖͍͍ͨ 9
10. HFWNYU Ͱঢ়ଶۭؒϞσϧ 10
11. ਪఆ݁Ռ w ऩଋͯ͠ΔͬΆ͍͚ͲO@F⒎গͳ͍Α͏ͳɾɾɾ 11 <தུ>
12. ൺ΂Δͱ͜͏ 12 ϩʔσʔλΛ ϓϩοτ ঢ়ଶۭؒϞσϧͷ ਪఆ݁ՌΛϓϩοτ
13. ਖ਼ن෼෍ͱ΋ൺ΂ͯΈΔ 13 ਖ਼ن෼෍ ϩʔσʔλ ҰൠԽۃ஋෼෍
14. ೥ޙΛ༧ଌͯ͠Έͨ w ͱΓ͋͑ͣಉ͡৚݅ͷαϯϓϦϯάΛ͢Δͱऩଋͤͣɾɾɾ 14
15. ͜͜·Ͱͷ·ͱΊ w ߴদࢢͷ࠷େ߱ਫྔ͕૿͑ͯΔ͔͸Α͘Θ͔Βͳ͍
    w ঢ়ଶۭؒϞσϧΛ࢖ͬͨํ͕มಈ͕׈Β͔ʹͳΔ
    w ҰൠԽۃ஋෼෍Λ࠷େ஋ʹ࢖༻ͯ͠Έ͕ͨɼਖ਼ن෼ ෍ͱൺֱͯ͠ɼฏۉ஋ͷมಈ͸ΑΓ׈Β͔ʹͳͬͨ ͕ɼ۠ؒਪఆͷ෯͸ͻΖ͕ͬͨʢࣄલ෼෍Λݕ౼͢ Δඞཁ͸͋Δ͔΋ʣ 15
16. ΋͏গ͠΍Γ͔ͨͬͨ͜ͱ w ҰൠԽۃ஋෼෍Ͱਪఆͨ͠ͱ͖ʹ͸ɼਪఆ஋͕ෛʹ ͳ͍ͬͯΔͷͰθϩ੾அ෼෍Λ૊Έ߹ΘͤΔ
    w ೥ޙ͘Β͍·Ͱͷ༧ଌ۠ؒͷਪఆʢճࠩ෼ͳͲ ͍Ζ͍Ζͱࢼͯ͠Έ͍ͨʣ
    w εΫϨΠϐϯάͯ͠શࠃͷσʔλΛҾͬு͖ͬͯͯ ϚϧνϨϕϧͷঢ়ଶۭؒϞσϧ 16
17. &OKPZ 17