- へろへろ JAGS 実験めも．cores = 4 の library(doMC) で n.chain = 3 の rjags を並列化計算させると，並列化しない場合に 90 秒を要する計算が 50 秒になった．もっと長い計算なら，たぶん 2 倍ぐらいの高速化になるかな? 08:51:10
- 2 倍になるのは WinBUGS でのサンプリング並列化と同じぐらい．JAGS の場合，各分業の結果をたばねて mcmc.list として出力してるんだけど，当方の細工がしょぼくて print(), plot() はできるけど gelman.diag() できない現状… 08:55:24
- 大会会場についても空き部屋でへろへろ MCMCglmm 実験のつづき…doMC の 4 core 並列化すると計算時間が15秒から9秒に短縮…というか、JAGSより6倍も速い⁉︎ さすがはGLMM専用のMCMCさんぷらーだな… 10:21:45
- RT @monotropastrum : @KuboBook むかしrjags並列化したときのコードですが、ご参考になりますでしょうか。 ito-hi.blog.so-net.ne.jp/2012-09-23 10:25:25
- @monotropastrum ありがとうございます！ 10:26:04
- せっかく学会会場にいるんだから少しは勉強するかということで、午前最後の一時間は群集生態学というか水産資源解析の最先端の統計モデリングの三講演を聴いた…temporal turnover、非対称ロジスティク回帰、MaxEnt の拡張計画… 12:25:45
- 昼飯はまた喫煙区画のあるカフェ。タバコ吸わないんだけど昼飯時間でもひろびろすかすか空間になってる。現時点ではケムリ発生源ゼロ。のんびりできるなー pic.twitter.com/FDHgdUBksB 12:29:32
-
- 午後の最初の20分間は統数研の江口さんの受賞講演を聞いてみることに。要旨によると1980年代から10年ごとぐらいに区切った研究変遷のハナシのようなので。キーワードは情報幾何、漸近有効性、選択性バイアス、漸近一致性の崩壊、機械学習、高次元データ解析、L1正則化の意外な発展……？ 13:12:45
- 今のところ日本統計学会が出す各種の賞あれこれの授与式が続いております。高校の先生たちは挨拶が長い… 13:15:58
- 受賞せれもにー、まだまだつづく…江口さんの受賞講演きいてる時間ないかも？ 13:20:49
- RT @ibaibabaibai : 単に逃亡しているように思われる久保さん． 13:23:50
- せれもにーがおわらないので、受賞講演は聞かないで、深谷君の群集動態モデルの改良の講演を聴きに行くことにする。 13:29:30
- 群集動態の統計モデリングのハナシを聞いて，ただちに会場から撤退．宿の部屋にひきこもって準備作業のつづき．いつまでたっても終わらん……というかぜんせん進んでいない…… 14:32:36
- 深谷君の潮間帯群集動態（フジツボとか固着性なやつら）の観測誤差のハナシをちゃんと聞くのは今回が初めてだったので，なかなかおもしろかった．「観測誤差によるバイアス」とは前回とはちょっとずれた点を見て「おや変化したな」と記録して群集動態の速度を実際よりも速いと推定してしまうこと． 14:41:27
- 今回の深谷モデルの改良点は，誤観測したときの「ホントはこれでしょ?」補完の改良．これはまでは広い範囲の平均値を使っていたけれど，今回は正規分布カーネルで観測点近傍から代替種を探すという方式にしたとのこと…「しましま」岩礁なので上下と左右の相関が異なるので二次元正規分布なカーネル． 14:56:29
- 深谷モデルはまあ正常進化・改訂なわけだが，すごいと思ったのはこのめんどうな計算を（よくわからない工夫によって）JAGS だけでこなしているところ．距離行列は事前に計算すればいいんだろうけど，局所群集構成は各ステップごとにいちいち評価する必要があるわけで…また教えてもらいますかね… 15:01:05
- ということで，自分の準備作業にもでるのであった…… 15:01:32
- #みどりぼん 第11章の空間統計学な問題は WinBUGS の car.normal() であつかっていたわけだが，「脱」ゐんばぐすということで，R の libaray(CARBayes) で代替できないかと…うーむ，出力されるオブジェクトがちょっとヘンぽくってやる気がうせる… 15:32:14
- あ，そういえば今日は祝日なのに東大正門前の「ルオー」 www.google.co.jp/search?q=%E5%9… が開店してたな，ちょっと入ってみようかな……と逃避な気分に……えーい，この部屋あついので，本郷６丁目住宅地どまんなかだけど窓を大きくあけて冷却してみる．外は曇天．気温26℃ぐらい? 15:36:31
- RT @ibaibabaibai : . @KuboBook ルオーのカレー美味でお勧めです．デフォルトで紅ショウガがつかなくなったんですが，頼むと来ます．紅ショウガ必須！ 15:46:23
- @ibaibabaibai ありがとうございます．晩ごはん第一候補は「ルオー」のカレー，紅ショウガ添え，ということで． 15:47:32
- . @ibaibabaibai はい，この空間統計モデルものまでぬいてしまうと，ハナシがどうしようもなくすかすかになりそうなので，一夜漬け勉強をつづけたいと思います． 15:50:44
- 講演も聞かずに作業に時間をわりふるものの進捗よろしくないしんどい日曜日 - ぎょーむ日誌更新 hosho.ees.hokudai.ac.jp/~kubo/log/2014… #KuboLog 16:10:38
- RT @ibaibabaibai : 統数研のYouTubeアカウントに伊庭のレプリカ交換MCMC講義の録画を追加しました．2012年度のニコ生講義をベースに新たに収録したオリジナル版です．局所的極大のある問題で効果的な手法を初歩から解説 www.ism.ac.jp/news/2014/anno… 17:59:36
- 東大正門前「ルオー」行ってみたら，懸念していたとーりとすでに本日は終了していた．祝日は土曜日と同じで 17 時までなのかな? 明日の昼ごろ再訪してみよう．晩ごはんは「ルオー」カレーではなく，Family Mart の納豆巻きと豆腐に変更 (なんとなく) …… 18:32:16
- @leeswijzer あいかわらず準備が進捗してなくておいこまれているので，フォレスト本郷の部屋にもどって作業などしております……明日の昼こそは「ルオー」のカレーとコーヒーを! 18:42:36
- PukiWiki の代わりに Gollum を導入してみた - あらびき日記 d.hatena.ne.jp/a_bicky/201409… なかなかおもしろそうだけど Ruby ですか……うーむ．ゴラムなる名前は「指輪物語」から ja.wikipedia.org/wiki/%E3%82%B4… 18:46:15
- @TJO_datasci まあ，今晩は（明日の晩も?）「楽しく夕食」といった気分ではないので…… 18:47:52
- #みどりぼん 「統計モデリング入門」の図 11.4 などで，以前から作図プログラムのまちがい (中央値を足し算など) を指摘していただいたのを思い出して，ようやくその作図 R コードを修正したみたのだが…あまり変わらないかな? pic.twitter.com/Y0GXWEa82B 20:44:24
-
- RT @ibaibabaibai : . @KuboBook 訂正謝罪会見をw 20:46:16
- みなさま…… #みどりぼん のまちがいは永久に不滅です……! hosho.ees.hokudai.ac.jp/~kubo/ce/Iwana… 20:47:31
- 修正することで，もうちょっと見てくれのよい図になってくれー……とホントは期待していたのですが，そういう問題ではなかったのでしょうか．ああ，さっきの訂正謝罪会見ぢゃないけど，まだまちがいとか残存してるのかな…… 20:50:50
- まあ，一次の conditional auto-regressive (CAR) モデルなので，実際のところ，こんなものなのかも……と以前は思考停止して，作図プログラミングのまちがいに気づかなかったわけですが…… 20:53:46
