グラフ更新。代々木公園は感染がおさまってきたように見える。 寒くなって来たしね。一方代々木公園以外は数は少ないものの、くすぶっている。 代々木公園以外は推定日がわからない人が多いので、個々の 潜伏期の確率に従って、微妙な期待値の人数のグラフに…

前回の記事でgmpのisprimeがあまりにも正確というのを書いた。デフォルトのrep=４０だと、何年回してもprobably primeなのに実際は そうでない例(pseudoprime)を見つけられないと思う。 しかしreps = を少ない数にすればpseudoprimeを見つけられそう。 rでや…

http://www.mhlw.go.jp/bunya/kenkou/kekkaku-kansenshou19/dl/20140924-01.pdfhttp://www.mhlw.go.jp/bunya/kenkou/kekkaku-kansenshou19/dl/20140924-01.pdf ついに厚労省でもグラフを表示するようになったようだ。良い傾向である。自分のはこれに、感染日…

某つぶやきで、スコットランド独立賛成票が素数だと知った#2014.09.23中澤さんの別解など追加しましたhttp://business.nikkeibp.co.jp/article/topics/20140919/271500/目的 ＃本当に素数か確かめる。 ＃その前後にどのくらい素数があるか。 ＃ついでに、素…

多倍長精度計算結果を見やすく出力 - 裏 RjpWiki さんの所で、多倍長出力をキレイに出力する 記事があった。自分の前の記事だと 小数点以下が出力されて、無駄に長くなってダメダメですね。 反省です。今回もコードを勉強させていただいた。 ありがとうござ…

今回はグラフだけ。代々木公園以外が目立ってきていますね。＃一部入力データミス発見し修正。他にもあるかもしれんので参考程度に。

library(conf.design) #この行は抜け道でカウントされないのだ。 cat(primes(997),sep=":") 今回もRで最短だったが、都合の良いideone依存の パッケージを発掘しただけのこと。 自慢にはならないわなー。でも２連続R最短で嬉しい。 mamekinさんの正攻法最短…

今回は感染推定日や発症日から発表までにかかる日数をグラフ化してみた。感染から発表までの日数は17日でほぼ安定。 発症から発表までの日数は12日でほぼ安定。ということで、感染日は9/1まで、発症日は9/6までは 大体固定しており、今後大きくは変わらない…

メタアナリシス某つぶやきで、reiv()の存在を知った。 http://www.youtube.com/watch?v=0XFlPB0mEzA http://zanet.biz/R/fct/ma_metafor_unibb.zipよりダウンロードできるらしい。そして、自分自身は使うことはないものの、他人に教えることが多いEZR と結果…

2014/9/12の１８時現在のデータに更新 #18:56年齢階級別の症例数も追加今回今までになかったパターンの出現があり。 代々木公園周辺＋新宿中央公園どちらも行っている例が出た。 ただ圧倒的に代々木公園が多いこともあり、今回は代々木公園の 例として扱うこ…

#22：20発症日がずれていたのを手直し。本日9/11の21時現在の最新データに更新 http://www.mhlw.go.jp/bunya/kenkou/kekkaku-kansenshou19/dl/20140911-01.pdf代々木公園以外が増えすぎたので、代々木公園以外にまとめちゃいました。 代々木公園かそれ以外の…

本日19時現在のの最新のデータに更新 http://www.mhlw.go.jp/bunya/kenkou/kekkaku-kansenshou19/dl/20140910-01.pdf代々木公園近辺以外にも広がっている模様。http://www.who.int/mediacentre/factsheets/fs117/en/ http://www.mhlw.go.jp/bunya/kenkou/kek…

＃追記 ぐぬぬ、夕方は無かったのに、今見たら新データが出てる。 新データに更新した。 発症日も訪問日も未定なヒトがいて、コードを書き換える必要が。 一つのデータのせいで、修正に２０分もかかったよ。本日夕方以降に発表のデータにて、8/23-26あたりが…

＃追記 昨日と同じカラーのも追加。一部バグ修正。 ＃感染推定日barplotも追加 ＃発症日のbarplotも追加参考にした潜伏期 http://www.plosone.org/article/info%3Adoi%2F10.1371%2Fjournal.pone.0050972PLOS ONE: The Incubation Periods of Dengue Viruses …

カレンダープロット - 東京で尻を洗う声優の誕生日のカレンダープロット - 驚異のアニヲタ社会復帰への道カレンダープロッットが流行しているので、やってみた。某疾患の某公園の訪問を表示する。複数候補日がある人は、候補日数で割って 重みとする。つまり…

天才火消しエンジニア霧島「もしPMおじさんが丸投げを覚えたら」｜paizaオンラインハッカソンLiteしえる (cielavenir) on Twitterpaiza_solutions/poh1_2.R at master · cielavenir/paiza_solutions · GitHubPHP - POH Lite 天才火消しエンジニア霧島 0.01秒…

コードを見易く改変した下記記事に変更しました paizaオンラインハッカソン：動的計画法 Rで挑戦 - r-statistics-fanの日記