Upgrade to Pro
— share decks privately, control downloads, hide ads and more …
Speaker Deck
Features
Speaker Deck
PRO
Sign in
Sign up for free
Search
Search
ぼくの町の不安定
Search
tanaka_marimo
December 01, 2018
0
730
ぼくの町の不安定
ぼくの町の不安定の話です。
tanaka_marimo
December 01, 2018
Tweet
Share
More Decks by tanaka_marimo
See All by tanaka_marimo
ベイズABテスト
tanaka_marimo
0
300
状態空間モデルによる予測と補間
tanaka_marimo
0
190
時変係数モデル
tanaka_marimo
0
360
brmsの使い方
tanaka_marimo
0
470
不連続回帰デザイン
tanaka_marimo
0
500
真夜中には何が起こっているのか
tanaka_marimo
0
970
佐渡島を見積もる
tanaka_marimo
0
960
Featured
See All Featured
How GitHub Uses GitHub to Build GitHub
holman
467
290k
5 minutes of I Can Smell Your CMS
philhawksworth
199
19k
Imperfection Machines: The Place of Print at Facebook
scottboms
257
12k
Let's Do A Bunch of Simple Stuff to Make Websites Faster
chriscoyier
501
140k
Fireside Chat
paigeccino
19
2.6k
Done Done
chrislema
178
15k
What's new in Ruby 2.0
geeforr
335
31k
Designing for humans not robots
tammielis
247
25k
Writing Fast Ruby
sferik
619
59k
GitHub's CSS Performance
jonrohan
1023
450k
Build The Right Thing And Hit Your Dates
maggiecrowley
23
1.9k
Typedesign – Prime Four
hannesfritz
36
2k
Transcript
Japan.R 2018 ぼくの町の不安定 2018年12⽉1⽇(⼩雪) Japan.R 2018:LT
Japan.R 2018 誰?
Japan.R 2018 名乗るほどの者では ございません。
Japan.R 2018
Japan.R 2018
Japan.R 2018
Japan.R 2018 別に
Japan.R 2018 別に
Japan.R 2018 不安定
Japan.R 2018 調べてみた n 2001年〜2018年の6~9⽉(18年) n 関東平野32地点の午後の不安定降⽔を推定 n よく晴れた暑い⽇を抽出 n
関東平野で午前に降⽔がない⽇ n 東京・横浜・熊⾕のどこかで最⾼気温30℃以上の⽇ n 台⾵の影響を受けない、且つ、梅⾬期間以外の⽇ n つくば上空の気温・湿度から不安定指数を計算
Japan.R 2018 地上データ(気象庁HP) http://www.data.jma.go.jp/gmd/risk/obsdl/index.php
Japan.R 2018 梅⾬の⼊り明け(気象庁HP) https://www.data.jma.go.jp/fcd/yoho/baiu/kako_baiu09.html
Japan.R 2018 台⾵トラック(気象庁HP) http://www.data.jma.go.jp/fcd/yoho/typhoon/position_table/index.html
Japan.R 2018 上空データ(ワイオミング⼤学) http://weather.uwyo.edu/upperair/sounding.html
Japan.R 2018 不安定降⽔ 暖
Japan.R 2018
Japan.R 2018 午後の不安定降⽔有無 @PM ※1mm以上の降⽔を不安定降⽔有とした
Japan.R 2018 メカニズムの想像(不安定)
Japan.R 2018 メカニズムの想像(不安定) 寒 暖+湿潤
Japan.R 2018 メカニズムの想像(不安定) K −index = (T 850 −T 500
)+T d850 −(T 700 −T d700 )
Japan.R 2018 よく晴れた暑い⽇を抽出
Japan.R 2018 ✖︎ @AM よく晴れた暑い⽇を抽出
Japan.R 2018 よく晴れた暑い⽇を抽出 横浜 東京 熊⾕ 最⾼気温 ≧30℃
Japan.R 2018 低気圧や前線の影響を 受ける⽇を除外
Japan.R 2018 台⾵の影響がある⽇を除外 台⾵中⼼が 130~150°E 30~40°N の⽇を除外
Japan.R 2018 梅⾬期間を除外
Japan.R 2018 天気図を⾒ずに事例を 除外するのは難しい・・・
Japan.R 2018 モデル式 外れ値の 確率 ⼤気の 安定度 y =α ⋅
1 2 +(1−α)⋅logistic(β 0 + β 1 KI) ランダム⽣成 ロジスティック 関数で⽣成 降⽔の 有無
Japan.R 2018 MCMC MCMC MCMC MCMC MCMC MCMC MCMC MCMC
Japan.R 2018 推定結果 ※収束は確認済み ü 外れ値にロバストな推定値が得られた ü 外れ値は2割含む
Japan.R 2018
Japan.R 2018 Reference n 松浦 健太郎(2016)『StanとRでベイズ統計モ デリング』⽯⽥ 基広監修,共⽴出版 n John
K. Kruschke(2017)『ベイズ統計モデ リング: R,JAGS, Stanによるチュートリアル 原 著第⼆版』前⽥ 和寛 ・⼩杉 考司監訳,共⽴出版 n Nomura, S., and T. Takemi, 2011: Environmental stability for afternoon rain events in the Kanto plain in summer. SOLA, 7, 9-12.
Japan.R 2018 Enjoy!