ぎょーむ日誌 2015-03-(11-20)

2015 年 03 月 11 日 (水)

• 0600 起床． 朝飯・弁当の準備． 朝飯． コーヒー．
• 0840 自宅発． 昨日は雨から雪にかわり， 今日も雪． 0850 保育所着． 0920 研究室着．
• 来週木曜日の発表の準備がすすまない ……
  • RT @kenmetrics : 先週Gelmanがブログに書いていた並列化するstan()関数の実装 andrewgelman.com/2015/03/03/one… を手元で簡単に試してみたところ、たしかに通常のrstanの関数よりは速いけど、僕のMIRTモデルでは既存のdoParallelでの並列化の方がもっと速かったです。   09:53:46
  • 二週連続江戸講釈のあと，来週木曜日と再来週水曜日も出張講釈…準備が進捗しない．現実とーひ的に雑用あれこれ処分．再来週はGLMMがお題なんだが…尤度方程式の中で積分という説明，ウケが悪い．階層ベイズモデルのほうが説明が簡単ではないかな…   16:38:50
  • 
  • GLMM の尤度方程式の積分が不人気な理由（根拠ふりー憶測） - 積分とか好きではなかったし，高校卒業以降やったことがない - GLMで「何でも正規分布」が否定されて混乱してるのに，また正規分布が召喚されて「確率分布を混ぜて新しいのを作る」という邪悪な魔術みたいなことやってるから   16:49:29
  • RT @berobero11 : @KuboBook そうなんですよね、なのでまずはその簡単な階層ベイズでどんな分布が生成されそうか、グラフどりぶんな説明が分かりやすいかと思ってこうなりました。 heartruptcy.blog.fc2.com/blog-entry-154… まだ分かりにくい…。数式や階層ベイズとの関連を追加すべきなんでしょうね。   16:58:04
  • @berobero11 そうですね… GLMM の図示はなかなか難しいです．階層ベイズモデルだけ，なら…確率分布のまぜあわせという説明を回避できるかなー…という気もするのですが，回避するのもよくないような（…とうだうだしてみる）   17:26:31
  • いろいろな二項分布を並べてみました…という図を作ってみたが，何の説明にもなっていない…   17:28:05
  • 
  • まだ，こちらのほうがマシかも…という気もしたけどそれは錯覚で，「いろいろな二項分布あります」は誰でもわかるんだけど，正規分布の確率密度で重みをつけて足し合わせるというところが難関みたいで．   17:32:13
  • 
  • じつは再来週の GLMM について検討している場合ではなく，来週鹿児島で話さなければならない「どうして時系列データを直線回帰しちゃダメなの?」講釈の準備をしなければならないのだが…   17:34:28
  • 
  • 来週木曜日20時の生態学会・データ解析自由集会 hosho.ees.hokudai.ac.jp/~kubo/ce/EcoSj… ，久保の発表準備はできてないのに，懇親会の会場だけは確定しました．鹿児島大ちかくの「しながわ」 tabelog.com/kagoshima/A460… です．明日，参加予約フォームつくります!   22:14:04
• [今日の運動]
  • ラジオ体操: 0 回
• [今日の食卓]
  • 朝 (0640): 米 0.5 合． 豆腐の味噌汁． キャベツのサラダ． ブロッコリー． アスパラガス．
  • 昼 (1200): 研究室． 米 0.7 合． 豆腐． ブロッコリー． アスパラガス． ミニトマト．
  • 晩 (1900): えーと…?

2015 年 03 月 12 日 (木)

• 0600 起床． 朝飯・弁当の準備． 朝飯． コーヒー．
• 0810 自宅発． 0820 保育所着． 0850 研究室着．
• ラジオ体操．
  • 私のマト外れな質問に丁寧な返信いただいたおかげでハミルトニアンMCの理解が少し進んだ…かも．M-Hとちがって HMC には提案分布のたぐいは無く，陰的解法 leapfrog で次候補を決定論的に生成→単純なメトロポリス法，だったのか…   12:03:27
  • 
  • . @ibaibabaibai あ，すみません，またわかりにくい画像 tweet をやってしまいましたが，メトロポリスな受容確率の画像は stan-reference-2.6.0.pdf github.com/stan-dev/stan/… の p.440 です．   12:25:23
  • @ibaibabaibai 正直に言うと「単純なメトロポリス」の引用として画像を出したのですが，結果的に random momentum vector の記述も含まれていました．じつは勾配計算のランダムネスの重要性はわかっていなかったので…「だいじょうぶ」ではありませんでしたね…   12:38:44
  • @ibaibabaibai うーむ，Hamilton MC における乱数の役割とその軽重を理解できていませんでした…このあたり手がかりに勉強やりなおします．   13:40:41
  • 私の見当はずれ質問への返信メイルで HMC がちょっとわかった気がしたんだけど，twitter 上のやりとりで，じつは「わかってない」と判明…とりあえず jmlr.org/papers/v15/hof… とか読んでみようかな．   13:59:43
  • 
  • RT @wonosatoru : "「日本人の気質」と言われるものは、その時々の社会的状況に過ぎないんです" RT @HuffPostJapan 外国人歴史家が体験した3.11 huff.to/1Gtbv6z   14:15:21
  • いろいろ勉強しなければ…その前に来週の生態学会自由集会の発表準備が…さらにその前に自由集会の懇親会 hosho.ees.hokudai.ac.jp/~kubo/ce/EsjKo… 「おんらいん予約」の作業…メイルアドレスの Perl 正規表現（いいかげん版）は m/^[-\w\.]+@[-\w\.]+$/ か?   15:37:36
  • …ということで，ここで連絡してもヒトが集まるかわからんけど…来週木曜日の生態学会「データ解析」自由集会の懇親会の申し込みは hosho.ees.hokudai.ac.jp/~kubo/ce/EsjKo… からお願いいたします． 会場は鹿児島大ちかくの郷土料理「しながわ」 tabelog.com/kagoshima/A460…   15:46:38
  • あの日、Twitterのくじらが出なかったもう1つの理由 www.atmarkit.co.jp/ait/articles/1… "「Twitterにアクセスできなくなったら大変な二次的被害を被る人が…」" "そのとき、段ボールに詰まった未開封のサーバーがコンピュータールームにあることを思い出した"   15:59:52
  • 自分でもすっかり忘れていたんだけど Twitter を始めたのは 2011-03-11 の朝だった．まさかその日に大地震が起こるとは知らず，そのときに開催されていた生態学会札幌大会の「ねっと伝言板」としてスタートしたのだけど… twilog.org/KuboBook/date-…   16:04:31
  • 北大の情報基盤センターから連絡: 現在、学外からDoS攻撃と思われる大量のアクセスにより、 学外との通信で遅延が発生しております。 復旧まで今しばらくお待ちください。   16:29:31
  • RT @JSPP_qa : [Q1098] 質問：ヨーロッパの植生について jspp.org/hiroba/q_and_a… #jspp   18:17:38
  • @ibaibabaibai そういう原稿ありましたか…明日、さがしてみます。   22:13:24
  • RT @ibaibabaibai : STANの勉強まだしてないけど，思想は分かってきた．アルゴリズムありきではなくて，統計モデルを「多数の条件つき確率の集まり」と考えるのではなくて「変数全部の同時分布をギブス分布で書く→好きなlog prob（統計物理ならエネルギー）を書いてね」と理解するということなんだな！すげえ   22:13:46
  • RT @ibaibabaibai : 「変数全部の同時分布をギブス分布で書く→好きなlog probを書いてね」を実現するには，ランダムウォーク型のメトロポリス法を使えばスグできる．でも，それは行ったりきたりするしステップ幅に敏感．じゃあ，ハミルトニアンＭＣ，ステップ幅自動のNo U-turn型なら最強，って順番だな   22:13:59
  • RT @ibaibabaibai : しかし条件つき確率でモノを考えるのやめようよ．有向グラフのモデルも，同時確率をたまたまそう表現したと考えるべきだ．ＣＡＲモデル自体は問題ないけど，同時確率をギブス分布で書くべきで，久保本やゲルファント本のようなfull conditionalから導入するのは良くないよ．   22:14:13
• [今日の運動]
  • ラジオ体操: 0 回
• [今日の食卓]
  • また記録のだつらく

2015 年 03 月 13 日 (金)

• twilog しか記録がのこっていない ……
  • RT @ibaibabaibai : Log-Cauchy distribution　 en.wikipedia.org/wiki/Log-Cauch… 　やりすぎだろ，これ．   09:27:26
  • R package development のちーとしーと．たいへんすばらしい． twitter.com/StatGarrett/st…   09:41:15
  • 
  • 学会会員の年会費って科研費とかから支払っていいはず，だよね…? www.jsps.go.jp/j-first/qa.html   09:59:18
  • 
  • あ，まちがえた．これ科研費 FAQ ではなく最先端研究開発支援プログラム FAQ だ．しかし，どちらもぽりしーは同じかも?   10:01:47
  • @u_ribo うーむ，R2WBwrapper のようないいかげんな関数セットをわざわざ package 化する意味はほとんどないような……   10:07:11
  • 科研費 FAQ www.jsps.go.jp/j-grantsinaid/… では「学会」は「学会の大会や集会」といった意味，年会費をぶんどる学術団体の意味では使われてない．まあ，たいていの学術団体は会員でないと「大会や集会」で発表はできない…大会の旅費その他はOKなので年会費もOKなんだろう…   10:16:02
  • @TJO_datasci ありがとうございます．はい，たぶん OK なんですよ．問題は，お金を管理している事務担当のかたとの交渉かも．   10:18:56
  • 村上さんのコトバに反省させられる… "あまり他人の意見をあてにしていると、やがて痛い目にあうことになります。嗅覚を磨くことはとても大事です。野球場に行ったらスコアボードではなく、試合そのものをじっくり見ましょう。それと同じことです。" www.welluneednt.com/entry/2015/03/…   11:02:39
  • まだまだ雪ふりつづく札幌。外気温 +1℃ぐらい。   11:17:10
  • 
  • RT @nejiko : 男児が、概念としての電車や自動車だけでなく、実際の路線名や車種にこだわり、かつ、みているだけで嬉しい、という感じで線路脇に何十分も立っていられる、というのは割とよくあることだけど、うちの娘にはそういう傾向は見られない。一つのジャンルに対する偏執というのは今のところ感じられない   11:30:51
  • 北大ネット HINES の「学外から rsync 拒否」問題，担当のかたから連絡いただいた（現在調査中）: HINES機器のIPSではsshポートへの過度なアクセスについては 閾値を設けて一定時間通信を遮断します ※sshポートへのbruteforce attackなどの対策の為   12:20:28
  • １９世紀科学史あれこれ togetter.com/li/794391 "ガソリンはかつて産業廃棄物だった" おもしろい．しかし "モーターがそのまま発電機になることが偶然発見された"は「えっ」と… Wikipedia ながめると逆で，ダイナモに電流ながしてモーター発見とか…   12:38:28
  • 学外→北大への通信の保安体制 www.hines.hokudai.ac.jp/restriction/ IPS (不正侵入防止システム) なるものが稼働していて，たとえば ssh であれば「10 秒以内に 3 回の login 失敗」が記録されれば，そのアドレスからの通信は一時間ほど遮断される…とのこと．   14:17:22
  • どうも北大側に問題があるというより，自分で使っている ssh 経由の rsync に何か問題があるような気がしてきた．ある時点で rsync か ssh を update して，そのあとから挙動がおかしくなっている，とか……   14:18:52
  • @VARDIGA ありがとうございます．うーむ id_dsa には一個しか鍵が登録されていないのですが…しかし，たしかにこのあたりが怪しいので，調査してみたいと思います．   15:36:25
  • 今日はひたすらいろいろな雑用が続く．夕方はネット雑用．login user を調べる unix コマンド last とかひさしぶりに使ったよ．last -f /var/log/wtimp.1 とか．   17:22:32
• [今日の運動]
  • ラジオ体操: 0 回
• [今日の食卓]
  • 記録だつらく

2015 年 03 月 14 日 (土)

• また twilog だけ ……
  • ggplot2 の彩色ヒストグラム blog.datacamp.com/make-histogram…   06:44:23
  • 
  • 時系列モデルパラメーターのJAGS推定，いまさらじたばたしてるんだけど…「一本の時系列」みて，各年の決定論的な増分の推定って意外と難しい．乱歩にひきづられる．青が真の増加，赤が推定の中央値．独立な6反復の例．   15:14:31
  • 
  • 同じ設定でさらに 6 反復．ランダムウォークの標準偏差 0.1 で毎年の増分 0.02，これはわりと推定しやすい設定だと思うんだが…うまくいくかどうかはきまぐれ乱歩次第? 2パラメーター間の関係と推定の良さについての研究とかありそう．   15:20:47
  • 
  • うす赤な領域は95%ベイズ予測区間といったところで…時系列データという「点の数はたくさんあるけど各点独立ではないので情報は少ない」性質を反映して推定誤差大．独立と仮定する直線 glm()なら「p = 9.0×10^(-9) でゆーい! いえーい，温暖化傾向だ!」とかなるわけだが…   15:31:54
  • なンかおかしいな?…と思ってたんだけど，同じ図を二つもアップロードしていた．アタマがかなりおかしいのかもしれない…   15:46:23
• [今日の運動]
  • ラジオ体操: 0 回
• [今日の食卓]
  • きろくだつらく

2015 年 03 月 15 日 (日)

• 0600 起床． 朝飯の準備． 朝飯． コーヒー．
• 1010 自宅発． 1030 研究室着．
• 今日のぢたばた ……
  • パフォーマンスを向上させたいなら、カフェインを減らそう www.lifehacker.jp/2015/03/150314… "被験者のカフェイン使用をコントロールすることで…カフェインに関係するパフォーマンス向上はカフェインの禁断症状無しでは存在しないことを発見しました" 依存な私にはつらい結果…   12:08:13
  • 「ななめ直線まわりの乱歩」データへの状態空間モデルJAGSあてはめ，ニセ誤差はヤメて，ふつーに「システムノイズ」「観測ノイズ」を推定させたほうがマシな推定結果になっているような気がする．観測ノイズとやらは10倍ぐらい過大推定だけど…   15:18:53
  • 
  • しかし，この「傾き」の推定，ただ単に始点と終点を結ぶ直線の傾きを推定しているだけ? なぜか95%予測区間内に「ホントのかたむき」(青破線) をとらえてくれてはいるのだが…いっぽうで，この区間内に「傾きゼロ」が入る確率は25%ぐらいかな．   15:24:37
  • 今度は「逆」，システムノイズは0.01で観測ノイズとした場合．時間相関がすごく小さく，つまり「ふつーの直線あてはめ」モデルと対応する場合．先ほどと同じ状態空間モデルのJAGSあてはめ…まあ，妥当でしょう．システムノイズ過大推定だけど…   15:32:13
  • 
  • 「あたりまえ」すぎるから誰もいちいち指摘してないのかもしれないけれど，二種類のノイズを同時に推定する状態空間モデルは， ・自己回帰な乱歩 ・直線+各時点独立な正規乱数 (直線あてはめ) ふたつのモデルのどちらにも対応できる，というのが利点のひとつなのか?   15:39:19
  • 状態空間モデルの二種類のノイズ，システムノイズと観測ノイズの事後分布サンプルに負の相関がありそう．図はそれぞれの対数の散布図，対数をとると二変量正規分布っぽく見える（要・心眼），3000点つかって線形相関係数を計算してみると-0.17   15:49:51
  • 
  • RT @ibaibabaibai : @KuboBook それはそういうものです．データのぎざぎざを説明するのに両者に配分しているわけなので．   16:10:19
  • RT @ibaibabaibai : @KuboBook もし曲線のかたちのMAP推定値だけに注目するなら，両雑音の比が推定された曲線の形状をきめる．ところが，いまの理屈から周辺尤度がほぼ平らになるのは両雑音の積（？）が一定の曲線上．その曲線の上の微妙な凹凸で両雑音の比→曲線の形状がきまる．   16:10:22
  • . @ibaibabaibai まあ，library(dlm) で十分なのかもしれませんが，生態学会ではどうせ個体数とかオス・メス比とかの変動あつかいたいという要望があるので dlm ではダメ．CRAN packageさがしをがんばるなら，JAGS でやろうよといった煽動の準備中…   16:19:06
  • @TJO_datasci @ibaibabaibai 煽動おおくして船やまにのぼる，ですよ．ことわざの意味はよく知りませんが，いろいろあおっているうちに，尤度の高いところに到達できる…といった教えなのではないでしょうか?   16:30:50
  • 毎年の増分ゼロのランダムウォークへのあてはめ，状態空間モデルは無難なかんじで．事後分布の中央値は乱歩にひきづられるけど増分 95%区間にゼロが入っている．これが正規分布な直線あてはめだと，とうぜんながらどれも「ゆーいな傾き」となる…   17:17:28
  • 
  • @ibaibabaibai KFAS，ど忘れしないようメモしといたのにど忘れしてました…ちょっと cran.r-project.org/web/packages/K… など調べてみたのですが，時系列モデルの知識ゼロの聞き手が対象なので，BUGSの単純なモデル指定のほうが説明がラクそう…という気がしてます   22:49:42
  • かつて library(sspir) という非ガウスな時系列モデルあてはめ R package が存在していたとのことだけど，いまやそれは削除されてしまっている cran.r-project.org/web/packages/s…   22:53:38
  • 非ガウスな時系列あてはめできるKFAS cran.r-project.org/web/packages/K… 私が調べる時間があまりないのが最大の問題なんだけど…GLMMの関数の使いかた・JAGSの使いかた，どちらを勉強するのがおトク？みたいなハナシかな…   23:10:54
  • 
  • @ibaibabaibai う，そういうワナがあるんですか…   23:24:28
  • RT @ibaibabaibai : @KuboBook 分散が変動するから「張力」みたいなものも変動しますよね．収束の様子も大きく変化します．たとえば，データが不十分だと定数（１階）や直線（２階）になりますけど，その途中で動かなくなるなど．推定がうまくいってるときは割と大丈夫ですが，「ゆーい」でないと不安定に．   23:40:32
  • . @ibaibabaibai なるほど…うまく推定できているかどうか，ちゃんと確認するようにします．ありがとうございます．   23:42:04
  • そしてたったいま気づいたのでだけど，夕方にやってたあれこれ実験，分散がに種類あって，真の値がすごく小さく設定した s[2] の推定がうまくいってなかった，とは…そうか「過大推定」とかになってた原因はこのあたりにあるのかも．   23:44:32
  • 
  • RT @ibaibabaibai : @KuboBook いやダメなとき，本当にダメなのか自信もてない，みたいな．自分じゃなくて院生さんの体験ですけど．MCMCの陥ってる状況が本能的に察知できないと混乱してしまいますね．   23:46:23
  • 真夜中にじたばたしてもしょうがないので，本日はここで撤退．ああいう簡単な例題でコケるとは，なかなか油断ならないけれど…まあ，もうちょっと試行錯誤して，二種類の分散がうまく推定できるかどうか，調べてみましょう．どうしてもダメなら KFAS の使いかたもだいたいわかってきたことだし…   23:48:53
• [今日の運動]
  • ラジオ体操: 0 回
• [今日の食卓]
  • 朝 (0720): フレンチトースト．
  • 昼 (1200): 米 0.7 合． 鶏肉カレー． ミニトマト．
  • 晩 (1900): ハクサイ・ネギ・豚肉の鍋．

2015 年 03 月 16 日 (月)

• 0600 起床． 朝飯． コーヒー．
• 0800 自宅発． 今日もとらはバスを見ることができてごきげん． 0815 保育所着． 0845 研究室着．
• せとぎわ計算つづく ……
  • RT @ibaibabaibai : @KuboBook KFASはまあ線形近似ですからね．   06:17:29
  • RT @h_okumura : 🔖 Rで描いたクリーンなグラフ A Compendium of Clean Graphs in R shinyapps.org/apps/RGraphCom…   08:50:27
  • 良い証明は我々を賢くする証明である―ユーリ・マニンへのインタビュー | taro-nishinoの日記 | スラッシュドット・ジャパン slashdot.jp/~taro-nishino/…   10:22:54
  • 時系列データの状態空間モデル，二種類のノイズの標準偏差のうち「とても小さいほう」の事後分布推定，昨晩は収束がどーのこーの騒いでたけど，ただ単に「遅い」だけ．ただし推定値はどうやっても過大推定，これはたぶん原理的に難しいのではないかな…   11:22:36
  • 
  • 北大くらうど上さーばに rsync+ssh するとIPS (不正アクセス監視門) にはじかれる問題 (いまだ原因不明)，とりあえず担当者が IPS 監視対象外と設定してくれた．自分で不正アクセスを見張れ，ということで．えーと，/var/log/btmp の表示は lastb か…   11:28:49
  • RT @NatGeoJapan : カタツムリの「恋の矢」が相手の寿命短縮、東北大　子孫を残してくれるパートナーを傷つけるのはなぜ？　（ナショジオニュース） nationalgeographic.jp/nng/article/20…   12:36:38
  • 40歳を過ぎたら：村上さんのところ www.welluneednt.com/entry/2015/03/… "僕に言わせれば、治療法はとてもはっきりしています。日常的に身体を動かすことです" 集中力低下対策の回答．ここ2.5年ほどジョギングとかしてない…急いで帰宅するため15分ほど走ることはあるけど…   12:52:45
  • 明日，新千歳離陸なせとぎわ状態空間モデル，JAGS並列化がまだ理解できてなくてぱられるできん…計算おわらない…ふと気づいたのだけど，三とおりの場合について実験すれば良いので鎖ではなくジョブ三個を並列に実行させれば明朝には終了するはず!   17:34:11
  • 
  • 進行状況を確認してみると，場合わけジョブ並列化のおかげで，あと9時間ぐらいで終了するはず…まあ，何かミスがあれば全部やりなおしになってしまうのだけど．   17:38:28
  • RT @ibaibabaibai : まず．increment_log_prob(normal_log()で２次元CARとか実装する場合だが，これは結構「偶然合ってる」みたいな部分あり．状態変数xについては何の問題もないんだが，問題は正規化定数に入ってる事前分布の分散τ．   17:49:14
  • RT @ibaibabaibai : increment_log_prob(normal_log()にτが入ってるんだが，本当は正規化した条件つき分布を２回掛けるんじゃない．全部加えてから正規化するのがマルコフ場の正しいやり方．いまの場合は -2\log \tau^2になって偶然にも？正解．   17:49:18
  • RT @ibaibabaibai : 多変量正規でない一般のマルコフ場については，xについては任意のlog probを与えて計算できるんだけど，そこに入ってるハイパーパラメータを同時推定するのは，事前分布の正規化定数が明示的にかけないので無理．   17:49:21
  • RT @ibaibabaibai : ・・というだけならまだ簡単なんだが，実は２次元以上の非正規の場合には，ハイパーパラメータをベイズ推定すること自体がまずい可能性があるのが深淵なところ．「できないからダメ」って言っておけば簡単なんだけど，そうすると変分ベイズでやろうという人とか出てきて歴史的にも袋小路だった．   17:49:30
  • RT @ibaibabaibai : 北川本の非ガウスのような１次元は計算ができればかなり大丈夫．次元があがって非ガウスになると変わってくるというのが，いわゆる相転移とか複雑性とかいうこと．　＞次元以上の非正規の場合には，ハイパーパラメータをベイズ推定すること自体がまずい可能性がある   17:50:22
• 1810 研究室発． 買いもの． 1830 帰宅． 晩飯の準備． とら風呂． 晩飯．
• [今日の運動]
  • ラジオ体操: 0 回
• [今日の食卓]
  • 朝 (0630): 米 0.5 合． いろいろはいった味噌汁．
  • 昼 (1200): 研究室． 米 0.7 合． ミニトマト． 豆腐ハンバーグ．
  • 晩 (1900): 米 0.7 合． いろいろはいった味噌汁．

2015 年 03 月 17 日 (火)

• 0 時すぎの作業つづく ……
  • 3000 乱歩に対するあてはめ実験，3並列化ジョブのおかげで6 時間10分ほどで終了．3000 個の bugs オブジェクトから結果要約を…そうか，各 post.bugs の sims.matrix にすでに「まとめ」が格納されてて…   01:34:44
  • 
  • いまさら気づいたんだけど，CRAN の各 package のソースコードって， github.com に格納されているのかな? 例: print.bugs.R github.com/cran/R2WinBUGS…   01:43:17
  • せとぎわ真夜中作図．傾向のない (delta=0) 乱歩データ，直線回帰すると，「ゆーいな傾き」を発見する確率はほぼ100%．では状態空間モデルを使えば完璧かというと…6%ぐらいの確率で delta 事後95%区間にゼロが入らない…   04:09:45
  • 
  • 一階のランダムウォーク，おそるべし…ですよ．おやすみなさい．   04:10:13
  • . @ibaibabaibai はい，これが20%とか0.2%ではなくて安心はしたのですが…[自由自在なランダムウォーク] + [ランダム性の無い・少ない過程] で生成されたデータから何か推定すると，実現した乱歩にかなり依存してしまんだなあ，といった感想がにじんだ表現になりました…   09:27:40
  • 前回は飛行機に乗り遅れたので、今日は早めに札幌発。   10:55:12
  • 12:30 新千歳離陸の予定、羽田で乗り換えて鹿児島へ。今日は時間的よゆーあったので空港蕎麦、一灯園のもり。発表準備はススんでないけど。   12:13:36
  • 
  • 羽田でせとぎわ作業のつづき。このあと、ソラシドエア発搭乗、鹿児島へ。   15:12:28
  • 
  • 鹿児島空港着。鹿児島中央駅まで40分ぐらいとのこと。1250円。   17:19:44
  • 
  • 宿に到着、鹿児島中央駅のすぐ横、窓から九州新幹線の発着が見物できる…ちかごろはいよいよ鉄道熱が高まるとら坊が知ったら羨望するかも。バスもたくさん見物できるし。   18:14:02
  • 
  • RT @ibaibabaibai : [自由自在なランダムウォーク] + [ランダム性の無い・少ない過程] 　直線を出すなら２階差分のトレンドモデル，普通の乱歩は１階差分のレベルモデル，だと思うのだがその辺どうなってるんだろ．直線成分のあるのに１階差分を適用すると直線部分もシステムノイズの分散にはいるので過大推定の筈   18:27:28
  • 宿が駅に隣接してるので，晩飯はお気楽に駅デパ地下のカウンター式とんかつ屋へ．豚肉にはうるさい土地がらのようで，ソースつけずに食べると揚げ豚肉の良さがわかったような．   19:44:36
  • 
  • 駅デパ，何でもある…ちょっとふらふら…無印で文房具ちょっと買ったり，ビックカメラでレーザーポインターの単５電池を補充したり．鹿児島らしさはまるでないけど便利だな．道ゆく人たちのハナシことばに耳をかたむけてみたんだけど，薩摩イントネイションではないみたいなかんじで，ちょっと残念…   19:49:35
  • “ @masuiyuko : 息子、今週末生態学会デビュー。 www.esj.ne.jp/meeting/abst/6… #esj62 #生態学会 ” 変形菌の研究 変形体の「自他」を見分ける力…中学1年生でも研究発表できる生態学会大会、この点はすばらしいね。   20:00:53
  • RT @jack1972frost : 中１の変形菌研究家、増井真那さん。中２進級直前に生態学会で発表の快挙です！→ＲＴ　研究テーマの「変形体の自他を見分ける力」は、敷衍すれば、他者や対立する存在とどう向き合い得るのか、という根本的な問題への一視点にもなるのではないか、と期待しています。日本変形菌学の未来は明るいぞ！！   20:17:22
  • 札幌の今日の気温変動は 3.8-11.6℃だったんだけど，鹿児島は 13.1-23.3℃．暑い．宿の中はさらに室温たかい．Tシャツいっちょうでノンアルコールビールのみながら，せとぎわ作業をつづけております…発表は明後日の19時．   20:30:49
  • 
  • spBayes for Large Univariate and Multivariate Point-Referenced Spatio-Temporal Data Models www.jstatsoft.org/v63/i13 空間-時系列ベイズ+MCMCなspBayes改良版   20:57:56
  …… だめだ， もうねむい ……
• [今日の運動]
  • ラジオ体操: 0 回
• [今日の食卓]
  • 朝 (0720): 米 0.5 合． 豆腐の味噌汁． サケ． もずく酢・キュウリ．
  • 昼 (1150): 新千歳空港の蕎麦． 720 円．
  • 晩 (1900): 鹿児島中央駅の地下 1F のとんかつ屋 とんかつ かつ寿 アミュプラザ店 のカツカレー 980 円．

2015 年 03 月 18 日 (水)

• 鹿児島滞在 2 日目． 0600 起床． ぐーぐー寝たけど，ばててるかんじが …… 朝飯． コーヒー．
• 9 時すぎまで宿の部屋でせとぎわ作業 …… その後， 鹿児島大に移動． 1115 から生態学会の IT 委員会があるので．
  • RT @ibaibabaibai : 「モンテカルロ法の基本は大数の強法則でー」みたいに言っていても，デモの図見ると，そりゃ弱法則じゃん，みたいな．無限を含むような数学の命題をデモってどういう意味だが難しい訳ではあるけどね．いまの場合はたとえばコルモゴルフ不等式のデモみたいにするのが一案かな．   06:58:00
  • 鹿児島中央駅から郡元方面行きの市電に乗って5分ぐらいで唐湊（とそ）、料金固定で170円。通用門から鹿児島大に入れる。まっすぐ歩いて行くと共通棟一号館。   10:22:35
  • 共通教育棟の北側に学習交流プラザというのがあって、いろいろ作業とかできそう。2階というか2.5階の学習ラウンジは電源たくさん開放されてる。無線LANはsoftbank とか。FREESPOTはここでは電波弱いので移動する必要あり。   10:27:37
  • 
  • ということで鹿児島大生のいんふらを借用して、せとぎわ作業のつづき…   10:30:29
  • 
  • 鹿児島大の学習交流プラザ続報。学習ラウンジではないスペースでもたくさんの電源が開放されてます‼︎ 1F 食堂ですら、窓際の席では電源使用可能。すごいな…   10:42:03
  • 
• 学習交流プラザ 1F の食堂で大澤君と昼飯． まさか大澤君と育児について話しあう， といった「現実ばなれしている」としか言いようのない状況が実現してしまった ……
  • RT @dokinbi : 今回は、キンビのマメ知識をご紹介。玄関を抜けると、高さ12メートルの吹き抜けホール。正面を飾るのが、色とりどりの大理石を、横縞状に張り合わせた大壁面です。当館の設計者・太田実さんいわく、「カラー写真で夜景を横流しに撮ったイメージ... fb.me/2LcChNh3o   12:57:38
  • RT @teramonagi : 「えー、時系列本って言えばHamiltonもいいですけど、Brockwell&Davis本とか小学校で読むじゃないですかー？」ってゆるふわ数学勢に言われたのでちょっと見てみたら、2章のタイトルが「ヒルベルト空間」なってたので、もう何も言わずにそっ閉じ。   13:05:36
  • RT @teramonagi : @turuge いきなり完備な計量ベクトル空間を持ち出して、関数空間で議論する必要、ないと思いたいんです！！！   13:06:16
  • せとぎわ作業、オワリが見えない…気分転換にざつだん投稿。生態学会 #esj62 参加者の服装の圧倒的多数派は非スーツ。そして、相手がどなた様であろーが「先生」よばわりしない場合がほとんど。いつまでもこの文化は存続してほしいものです…   15:24:08
  • 
  • 人文系その他で、よく知りもしない相手をつらまえては、本人にことわりもなく「先生」よばわりしてみる文化にはなじめない…あんたの「先生」になった覚えないよ。しかし、本朝学術ぎょーかい全体では「ヒトをみたら先生」かるちゃーが多数派なのかも。   15:34:11
  • @berobero11 「御侍史」は若いお医者さんには違和感あるそうですよ… www.qlife.jp/square/hospita…   17:29:26
  • RT @NatGeo : Mysterious bright spots on dwarf planet Ceres may be water ice: on.natgeo.com/1GTpR0a   17:32:26
  • 宿にもどってせとぎわ準備のつづき．鹿児島大学構内を歩いていたら，ちょっと庭園っぽい場所をみつけた．あとで調べてみたら玉利池というところだそうで… www.google.co.jp/search?q=%E9%B…   17:39:37
  • 
  • 鹿児島大学ちかくの唐湊駅から 170 円市電に乗って鹿児島中央駅に移動．おきがるでいいな…   17:42:15
  • 
  • RT @Takashirouzu : 生態学会にDJIファントム2がいた．GIS・GPS技術と合わせてフィールド研究に活用しようという試み．   17:51:35
  • RT @sotetsuk : 線虫が尿嗅いでガン検知出来るんなら機械学習で出来るでしょw とか言ったら線虫舐めすぎという御言葉を頂きました…すみません…(´･_･`)   17:56:37
• その後もせとぎわ作業， つづく ……
  • @u_ribo ほんとに鹿児島大の学習交流プラザはいいですよね…   19:52:02
  • 時系列なデータと，直線回帰の GLM から生成されうるデータは性質がちがうんだ，といったことを説明する図…をいまだにちこちこ作っている．こんなことやってるから時間が足りなくなるのだが，やらないと，たぶん何も理解してもらえないだろうし…   19:53:55
  • 
  • えーと，明日夜の #esj62 自由集会後の懇親会 hosho.ees.hokudai.ac.jp/~kubo/ce/EsjKo… ，もうすぐ締め切ります…早めのもうしこみをお願いします．人数を店に連絡しないといけないんで． …さて，晩飯を食いたい気分でもないけど，何かまた駅デパ地下とかで豚肉でも食べようかな…   20:00:46
  • 宿でて鹿児島中央駅の新幹線改札口西側の（この付近では唯一の）公衆電話から明日の懇親会の店 tabelog.com/kagoshima/A460… に連絡して，すみません人数確定3人だけです，はい，お刺身注文しますといったハナシなどを．明日はそれほど混んでないとのこと．   20:56:57
  • 晩飯は駅デパ地下のそば茶屋 吹上庵で「黒豚つけそば」，人気の店のようで昨日は行列で入れなかった．薬味はわさびではなく柚子ごしょう．黒酢ドリンクがおまけでついてきた．このへんの和食系の店は卓上にカブのつけもの容器を置くのがデフォルト？   21:03:41
  • 
  • 関係ないけど，鹿児島市内はこんびにファミリーマートをやたらと見かけるような気がする…実際に多いらしい www.google.co.jp/search?q=%E9%B… エリアフランチャイズ方式？なんだかよくわからんな…   21:07:11
  • @kimegukurumi ありがとうございます！蕎麦好きというか，なぜか蕎麦についてだけは「いろいろな店で試してみたい」欲がすこしあるようです．出張先の外食→蕎麦，となるもうひとつの理由は「『何を食べようかな』と考えるのがイヤ」…というめんどくさがりな側面も無視できません…   21:11:03
• [今日の運動]
  • ラジオ体操: 0 回
• [今日の食卓]
  • まあ， うえのようなかんじで ……

2015 年 03 月 19 日 (木)

• 鹿児島出張 3 日目． 0600 起床． 朝飯． コーヒー．
• 今日の 18 時から 時間・時系列データの自由集会 ……
  • @ibaibabaibai そういう誤解があるとは知りませんでした…むしろ生態学分野の人たちのデータ解析みてると，ほぼ例外なく，Yが時系列だろうがなかろうが，glm(Y~t) の出力に含まれる p だの AIC だのをそのまま「ふつー」に解釈すればいいのだ…が多数派ですね．   01:05:49
  • @ibaibabaibai そうですね…構想としては「さあ，ここで時間的自己相関のありなしを見て，このデータに『むいた』モデルを作ろう」…とするのではなく，まず状態空間モデルをあてはめてみて二種類のノイズを同時に推定しようという方向性だったのですが…さて，どうなることやら，です．   01:14:16
  • @ibaibabaibai ここで想定している聴衆の考えている AIC 導出とは R の glm(Y~t) の出力に AIC が出力に含まれていれば，それは「正しい」AIC なのです．AIC そのものについても「あてはまりの良さ，ゼロに近いほどよい」ぐらいの人もけっこういますよ．   01:27:25
  • 鹿児島は今日も雨。宿の部屋でせとぎわ準備つづく…10時になったら大学の学習交流ラウンジ開くので、そちらに移動しよう。残り9時間…せとぎわすぎる…   09:15:26
  • 
  • RT @nowslope11 : 生態学会は、ジャージだったり登山靴だったり、熊除けの鈴付きリュックで来場する人もいるくらいだから、学会=スーツというイメージが全くない。   10:12:10
  • 短期時系列データを用いて非線形力学の因果関係を検出| Scientific Reports | Nature Publishing Group: www.natureasia.com/ja-jp/srep/abs… へー、どんな方法なんだろ…相原さんも共著   10:32:57
  • 
  • 生態学会大会 #esj62 会場である鹿児島大、今日も快適な学習交流プラザでせとぎわ作業進行中。カゴ大にこういう部外者か作業できる場があってよかった…北大だと図書館のラウンジとかかな？最近は北大生協も開放電源を設置してるのかな？   11:23:10
  • 
  • 生態学会 #esj62 のデータ解析自由集会，今回もたくさんのかたに来ていただき，ありがとうございました．粕谷さんは生存解析，私は時系列の状態空間モデル入門（しゃべりだす直前まで準備）のハナシをしました．けっこう需要あるかんじでした．   23:27:29
  • 
  • @ibaibabaibai クラシカルな生存解析の概要を聞きながら，現代の生存解析は R の library(survival) とはまた違ったものになっているのでは…という気がしました．Bayesian Survival Analysis といった本はかなり前に出てますね．   23:43:05
  • #esj62 自由集会，準備の時間もしゃべる時間もなかったので時系列 Y ~ 時系列 X な回帰については，みなさんの前でR実演やっただけで終了．「時間相関のある時系列データは見かけよりも情報が少ない」原理をお持ち帰りいただいて，と…   23:49:54
  • 
• まだ準備がおわってないのだが， 1730 になったので A 会場 に移動． 18 時まで会場の準備など．
• 18 時 時間・時系列データの自由集会 開始． みなさん， 200 人ぐらい集まってれたかな ……
• なんとか， 終了 …… みなさん，ありがとうございました．
• • 生態学会 #esj62 自由集会後の懇親会は鹿児島大学から徒歩 10 分ほどの「しながわ」 tabelog.com/kagoshima/A460… 薩摩・奄美の料理．キビナゴ天ぷら・豚なんこつ煮・トビウオつけあげ・鶏飯など．とても良いお店でした．   00:20:21
  • 
  • せとぎわな発表準備の日々が続いたので，ばてぎみ．にもかかわらず…明日もあさ早くに鹿児島空港に移動しないといけないんだけど，ネットとかながめてうだうだしております…いやー，やっぱりしんどいからさっさとねよう…   01:00:11
  [今日の運動]
  • ラジオ体操: 0 回
• [今日の食卓]
  • なにがなんだか ……

2015 年 03 月 20 日 (金)

• 鹿児島出張 4 日目． 0550 起床． ねむい． しんどい． 朝飯． すぐに宿を出て， 空港バスのりばへ． 40 分ちょいで到着． 0810 離陸．
  • せとぎわ準備と自由集会の本番・懇親会だけで終わった鹿児島滞在。早起きして B787 で北海道に戻ります。その前に、黒酢ドリンクを一杯ひっかけて、と…   07:52:12
  • 
  • 快速エアポートに乗って新千歳空港からの出発まち。鹿児島との気温差が10℃以上ある寒い世界に戻ってきた…お、とらの好きな旭川行きスーパーカムイな快速エアポートが向かいに入ってきた。   13:03:23
  • 
  • RT @Takashirouzu : 阿部永先生「たくさん動物を殺したけど，全部標本にしたよ」 #esj62   15:01:50
  • RT @hatebu : 安全なウェブサイトの作り方改訂第7版の変更点と変わらない点 | 徳丸浩の日記 (59 users) bit.ly/18V95SR 2件のコメント b.hatena.ne.jp/entry/blog.tok…   15:41:10
• 14 時ごろ札幌駅着， ちょっと本屋などをうろうろして， 15 時ごろ研室着． どーにもこーにも， もえつきうだうだ …… 1800 帰宅． 晩飯． とら風呂． 寝かしつけしてたら， こちらもあっというまに寝おち．
• [今日の運動]
  • ラジオ体操: 0 回
• [今日の食卓]
  • 朝 (0630): 宿の朝飯． 今日は鶏飯をいただいた． 鶏スープ飯．
  • 昼 : 羽田空港でカレーくったり， 機内でベイグルかじったり ……
  • 晩 (1900): 米 0.7 合． エノキダケの味噌汁． アジの塩焼き．