ぎょーむ日誌 2011-01-(01-10)

2011 年 01 月 01 日 (土)

• 0800 起床． 札幌ではおだやかな天気の元旦． 朝飯． コーヒー．
• のんびり怠業．
• iTunes で映画レンタル， 「インセプション」． 脳くらっきんぐものは攻殻機動隊とか谷甲州 (!) とかでなれしたしんでいるので， 「あーよくあるやつだなぁ」 と思ってたんだけど， ノーラン監督の作癖というべきなのか 「ややこしい設定だけど，あとでつっこまれないように 事前に十分設定をツメて，しかもそれをハナシの展開に使おう」 ぽりしー + 画像技術者たちの現実捏造に気合が入っていて楽しめました． あえてコマかい指摘をすれば， この世界の「現実」の階層において夢を共有する device だけが突出したてくのろぢーと申しますか， ほとんど オーパーツ ぢみたしろもので， そのためこの「現実」のふわふわさが強められていたような (とすると，ラストシーンとも関係するかも) …… 製作側もその点に関連して何か演出的な意図があるのか， 計算機のたぐいは登場せづ， 電子機器・通信機器もびみょーに古いかんじのものばかり出していて， ですね．
• 昼飯．
• うろうろ．
• 晩飯の準備． 晩飯．
• ちょっとだけ仕事をススめる ……
• [今日の運動]
  • 朝のラジオ体操． うんどう休養日．
• [今日の食卓]
  • 朝 (0900): もち． 正月料理あれこれ．
  • 昼 (1330): もちはすまし汁の雑煮． 正月料理あれこれ．
  • 晩 (1900): ハクサイ・ダイコン・エノキダケ・豆腐・豚肉のキムチ鍋． 卵雑炊．

2011 年 01 月 02 日 (日)

• 0830 起床． ちょっとねぼう． 朝飯の準備． 朝飯． コーヒー．
• 1000 自宅発， AEON 桑園に買いものに． 天気曇天， 気温零度ぐらい． 今日の AEON はすかすか． あちこちに多種多様な福袋おいてあるけど， 買わずに逃げる． 1110 帰宅． ひろい店内をうろうろしてたので時間かかってしまった．
• 昼飯の準備． 昼飯．
• うーむ， 仕事がススまない日々 ……
• 1615 日没直後に自宅発 今日は ミニ大通公園 ではなく， JR 桑園方面に北上してから 北大構内に入って獣医学部-低温研前を走る． ちょっと雪． 1700 帰宅． 体重 69.8 kg．
• リンゴジャム作り．
• 晩飯の準備． 晩飯．
• [今日の運動]
  • 朝のラジオ体操． 雪上走 45 分間．
• [今日の食卓]
  • 朝 (0900): もち． 正月料理あれこれ．
  • 昼 (1150): ナメコ・ベイコンスパゲッティー．
  • 晩 (1850): 米 0.7 合． 手巻き寿司 …… 自分で準備したのは初めてだったんだけど， 意外とおてがるだった． ミツバ・豆腐のすまし汁．

2011 年 01 月 03 日 (月)

• 0700 起床． 朝飯の準備． 朝飯． コーヒー．
• 怠業．
• 1145 自宅発 晴． 今日も JR 桑園を経由して北大構内に侵入， 低温研前までいってひきかえすコース． 1225 帰宅．
• 昼飯の準備． 昼飯．
• 1400 自宅発． 晴． 1510 研究室着．
• 1/7 の発表準備， なげやりになると少し進捗した …… 最初から再利用 & きりはりに徹していれば， 苦しまずにすんだわけだ．
• 新春とびこみデータ解析こんさる． 大学院生， 例によって困難な地上部・地下部問題 + 形態のハナシ． 区画内個体は独立ではない， とゆーのをデータを前に指摘されてみればなるほどと思った．
• 新春とびこみ院生部屋ネット配線工事． うう．
• 1910 研究室発． 晩飯． 2140 帰宅．
• 自宅至近距離の温泉．
• [今日の運動]
  • 朝のラジオ体操． 雪上走 40 分間．
• [今日の食卓]
  • 朝 (0900): チーズのせ焼バゲット． ベイコンエッグ． ヨーグルト． リンゴジャム．
  • 昼 (1300): うどん， ミツバ・ネギ・油揚の卵とじのせ．
  • 晩 (1940): なぜかしら Avion でおいしいものを食べてしまった ……

2011 年 01 月 04 日 (火)

• 0720 起床． 朝飯の準備． 朝飯． コーヒー． 0855 自宅発． 曇． 0920 研究室着．
• いきなり院生のとびこみデータ解析こんさる． 植物個体の空間分布の統計モデリング， ただし集中分布してるので「空地」が多い …… 難問． WinBUGS の car.normal() とかでガウス確率場を作って …… という手口， じつはゼロの多いデータではあまりうまくいかなかったりする． 少なくとも収束は悪くなる．
• 今回の場合は， じつはどの場所にどれぐらい (個体をあらわす)「点」 がちらばるかとゆーのは主問題ではなく …… とりあえずの問題としては， あたえられた点の配置のもとで， これらを二色 (二種類の植物) にわけなさい， という問題にしてしまえばよいような気もする． しかし， これまたゼロが多いと (species ごとの集中分布) MCMC に時間がかかる問題になるのかもしれない．
• 北大 の ArcGIS ライセンス， 研究室 LAN ルータのしたにある PC の場合は， ルータの IP アドレスとその下で使うマシンの台数を申告せよ， か．
• うー， アタマがなかなかきりかわらん ……
• 昼飯．
• 7 日のハナシ， ひたすらきりはり． いやはや．

[地球統計学とクリギング法]

間瀬さんの新著を買った． 空間統計本， いままでの間瀬さんの本をさらに丁寧にした内容． R の library(gstat) や library(geoR) の動作を理解したいときに役たちそう． まあ， この本が必要な生態学研究者はあまりいないかも． 出版社の編集はテぬきぎみで， 図の文字とか小さすぎたりして読みづらい．

[主要目次]
 1章　地球統計学とクリギング法
 2章　確率場と定常性
 3章　共分散関数とバリオグラム関数
 4章　共分散関数とバリオグラム関数の例
 5章　単純・通常クリギング法
 6章　ブロッククリギング法
 7章　普遍型クリギング法
 8章　多変量確率場とクロス共分散・クロス・バリオグラム関数
 9章　多変量クリギング法
10章　交差検証法によるモデル当てはめの良さの検証
11章　モデルの当てはめ
12章　確率場のシミュレーションと条件付きシミュレーション
13章　クリギング法による解析の実例
14章　パッケージgeoR を用いた地球統計学的解析の実際

• library(geoRglm) についての説明はない． help ファイルをつらつらと見くらべてみると (geoR と geoRglm)， 内容がけっこうちがっていて， たとえば library(geoR) には MCMC な関数とかない． 間瀬さん本をよく見てないんでよくわからないのだけど， ガウス確率場 + エラーも正規分布となっている場合には， 試行錯誤な推定計算とか必要ないのか?
• 現実とーひのため， R の HTML help まわりを改善．
• きりはり発表じゅんびはあとヒトいきでおわらづ …… 1810 研究室発． 雪． 買いもの． 1840 帰宅． 晩飯の準備． 晩飯．
• ばてたので， また近所温泉．
• [今日の運動]
  • 朝のラジオ体操． うんどう休養日．
• [今日の食卓]
  • 朝 (0810): チーズのせ焼バゲット． ゴボウ・ニンジンのコンソメスープ．
  • 昼 (1300): 研究室． 米 0.5 合． スープ． リンゴ．
  • 晩 (1930): 米 0.5 合． ゴボウ・ニンジンのコンソメスープ． ダイコン・ニンジンのなます． チンゲンサイ・ピーマン・モヤシ・シイタケ・豚肉の ショウガ・ニンニク・豆板醤・オイスターソース炒め．

2011 年 01 月 05 日 (水)

• 0655 起床． 朝飯の準備． 朝飯． コーヒー． 0850 自宅発． 雪． 明日の出張， 鉄道・航空で無事に移動できるのだろうか． 0920 研究室着．
• うだうだと能率わるく出張の準備 …… ほんらいは岩波本のちぇっくだったのだけど， この発表準備がだらだらしていたので， 明日はあまり検討すべき項目が ……
• 生態学会雑用がしつこく ……
• 早めに昼飯くって， 散髪にでる． なンでこのタイミングで散髪なんかするの， とゆーかんぢだけど， 今月はごちゃごちゃしているので今やっとくのがよいと判断したのだが …… 雪のふるなかクラーク会館まで行ってみたら， 満員御礼状態だったのですごすごと撤退してきた． 平日の正午， 意外とこんでるねえ．
• 雪がどんどんふってきた． 明朝は早めの JR にのらないと空港に到達できないか?
• 研究室ねっと雑用とか …… いろいろわりこみが ……
• 7 日のいいかげんな 30 分とーく (いいかげんではよろしくないのだけど， まあチカラ不足で) の準備， ひととーりはできた． 時刻はすでに 16 時すぎ．

[寒波襲来]

明日からすごーく冷えこみます， とゆー予報． まあ， サムいのはいいんだけど， 寒波の襲来のせいで雪がどかどかふるんじゃないかと不安． とくに 7 日 (金) は 2235 新千歳空港着， 大雪で交通機関がストップしたら空港びばーく という楽しい事態になりかねん．

• 明日の夕方， 伊庭さんとの面談があるので， ひさびさの 岩波統計モデル本 作業．
• しかしあまりススまぬうちに撤退． 1720 研究室発． 雪． やはり今月はいろいろばたばたしているので， 気分転換もかねていまのうちに散髪しておく． 1730 昼間に撤退した 理髪店 に再度突入． 今回はがらがらですぐにやってくれた．
• 閉店時間寸前だったためか 50 分でしあげてくれた． 1820 同発． 買いもの． 1840 帰宅． 晩飯の準備． 晩飯．
• 晩飯後に岩波本作業をつづけようとしたんだけど …… どーにもこーにもテがうごかない． 自宅内雑用ごそごそ ……
• 明日の朝飯の準備とか．
• 岩波本作業， 少しだけつづける．
• [今日の運動]
  • 朝のラジオ体操． うんどう休養日．
• [今日の食卓]
  • 朝 (0750): 米 0.5 合． ネギ・ナメコ・豆腐の味噌汁． おろしナガイモ．
  • 昼 (1240): 研究室． 食パン． リンゴ．
  • 晩 (1920): 米 0.6 合． タマネギ・トマト・タコ・ホタテのカレー． ヨーグルト． リンゴジャム．

2011 年 01 月 06 日 (木)

• 0640 起床． 朝飯の準備． 朝飯． コーヒー． 0720 自宅発． 雪がどしどしふっている 早めにでることに
• 0735 JR 桑園． 0742 同発， のつもりだったけど 出発ホームまちがえた． 札幌方面いきは二番と四番に分離してるのか． 0746 同発．
• 0750 札幌着． うーむ …… この時間は快速エアポートは 15 分間隔の運転ではないのか． 20 分ほど待つことに． 寒いなあ， と思ってたら 0800 快速エアポートが入構． 車内で座れた． 0810 同発．
• 新千歳空港ふきんは晴． よしよし． おっと， 線路上の雪のかたまりのためしばらく停車． 0847 新千歳空港着．

[B-747-400]

翼面に湯噴射して氷結解除している． 0920 ANA 54 に登場． 席 15H． 満席．

• 0950 離陸． 1110 羽田着陸． 15 時すぎに立川の統数研につけばよいので， 空港内でのんびり無線LAN． 1149 京急に乗って同発． 東京， 晴れてて温暖．
• 品川で JR に乗り換えて 1217 同発． 新宿で中央線特快に乗り換えて 1237 同発． 1302 立川着．

[風にむかって]

とゆー名前の銅像らしいです． 立川飛行機 にちなんで， なのかな? JR 立川駅とモノレイル北駅の中間にある．

[今日の昼飯]

駅北のカフェで昼飯． 豚肉． 990 円．

[多摩モノレイル]

時間があるので， てけてけと 15 分ほど歩く． このモノレイル下は公園みたいになってて， 歩くと快適．

[統数研]

14 時ごろついたので， まずは 1F ロビーで岩波本のみなおし作業． 15 時に伊庭さんの研究室にいったのだけど， 室内から (先客との) 白熱した議論が聞こえてきたので， さらに 20 分ほどそのあたりの廊下のテイブル (統数研はこういうのがあちこちにあってよい!) で原稿みなおし．

• で， 岩波統計モデル本 そうだんなのだが …… 私のほうはぎょーむ日誌につらつらと記録されてるありさまで， 作業が進捗しておらず， しかも今回はたまたま伊庭さんもご多忙でチェックしてるヒマもなし， とゆーことで …… まあ， ふたりで一時間ちょっと雑談していただけでした．
• で， そのあと所内で 研究会 の案内はりのてつだい．
• 立川駅ちかくで晩飯．
• 1910 本日の 宿 着 …… ここは一泊 5500 円ぐらいの安いところなんだけど， 同じ建物内にこんびにがあっていろいろ便利で …… と思ったらツブれて撤退していた． うーむ， 立川のこんびに競争はキビしい ……
• 生態学会雑用その他雑用あれこれ ……

発表編成部会の皆さん

久保です．また仕事のお願いでもうしわけありませんが，
1/10 ぐらいまでに大会プログラムの口頭発表ペイジを
点検してもらえないでしょうか．

大会プログラムは

http://***

の「口頭発表一覧」からダウンロードしてください．
やっていただきたいことは……

・http://*** に示している
  それそれの分担の会場についてチェックしてください
・発表編成の web page に掲載されている情報が正しいものと考えてください
   - こちらで生成したデータをもとに手作業で一覧が作られたのですが，その
     転記の過程で生じたエラーをさがすことが目的です

   - 講演番号・時刻・題名が正しいか，発表者がぜんぜん別人になっていないか，
     をみてください
   - 登録者自身の入力ミスと考えられる不備は放置してください

・とくに問題なければ，連絡していただかなくてもかまいません
・何か問題ありましたら久保までメイルでお知らせください

これまでに発見された問題は物質循環 D1-12, D1-13 で，
・講演番号がまちがっていた
・同じ講演が二回くりかえされていた
といったまちがいで，すでに修正をお願いしています．

お手数をおかけしますが，よろしくお願いいたします．

• ぐったり．
• そして来週の水産統計モデリング講習会の準備が ……
• [今日の運動]
  • 朝のラジオ体操． うんどう休養日 …… 立川歩行?
• [今日の食卓]
  • 朝 (0700): 米 0.5 合． ダイコン・ネギ・豚肉の味噌汁．
  • 昼 (1330): 立川駅北のカフェのおすすめ定食．
  • 晩 (1800): 伊庭さんにごちそうになりました． ありがとうございました．

2011 年 01 月 07 日 (金)

• 東京出張二日目． 0640 起床． 朝飯．
• 生態学会雑用．
• 今日の発表内容のみなおし作業 ……
• 0905 宿発． 曇． 0930 統数研着．
• 10 時より 研究会 はじまる． 午前の部は 「説明変数に誤差のある回帰をめぐる話題」 なるもので最初は伊庭さん． 説明変数 x の誤差のおよぼす影響は， 線形ではなく非線形の問題で重要になるんでは， といったご指摘でカーネル回帰 + MCMC の工夫とか．
• で， 次が私で 2010 年 3 月の自由集会 の自分のハナシ (回帰はヤメて統計モデルをかえましょう) などでごまかす …… じつはこれが失敗で， 粕谷さん・伊東さんがされた直線回帰の例をくわしく説明すべきだった (いいわけ: 伊庭さん「線形の場合はそんなにおもしろくない」とのご指摘だったので) ……
• どうもいろいろな意味でハズしたハナシでしたな． しかしたまにはこういう居心地のワルい発表をすると， 自分が何モノなのか少し理解がススむような気がする．
• で， 午前の部のメインである の「セミパラメトリック統計数理の情報幾何」 (「セミパラメトリック」なのは各 x に対応する真の値がパラメーターになるから)． 直線回帰 (しかも傾きを推定するだけのやつ) でもバイアスのない推定はそれなりに難しい， といったことを数学的に示すもの． つまり 2010 年 3 月の自由集会 で伊東さんが実例で示されたあれこれで， 当時は「何かヘンだな?」と思っていたんだけど， じつはあのように「バイアスがなかなか解消されない」 というのは正しい結果なのだとわかった．
• 単純なスコア関数 (対数尤度をパラメーターで微分したもの) を使うと必ずバイアスが生じるという指摘があり， これを解消するためには nuisance parameter の軸に直交するような スコア関数を使って最尤推定しなさい，ということに ……
• で， これをやるために推定関数法というやりかたでスコア関数を作るわけだが (このあたりからよくわからなくなるわけだが) …… 尤度がわりと簡単なかたちでないと， 直交 (projected) スコア関数が explicit に導出できない， といったハナシになったりしてですね．
• ともかく， 切片ナシの傾き推定の場合， ある場合には TLS (直交最小二乗法の尤度) がよく， またある場合には (全 y の和) / (全 x の和) といった推定量がよく， 一般にはその中間に …… といった結果の例が示されたりですね． x について複数の観測が得られている場合にも (この状況は伊東さんの実験にもあり)， ひとつひとつの x に対する観測回数が少ない場合は ダメといったご指摘なども (n - 1 で割る不偏分散にすればよい?)．
• 甘利さんのハナシをきくのは今回が初めてなのだが， 内容だけでなく話しかたもたいへんよいものでした． ああいう軽妙なかんじで内容のあるハナシができるようになればよいのだけど ……
• 昼飯．
• 午後は「仮想データとその周辺」 というまた別の研究会です， という伊庭さんの紹介からはじまる． 関連きーわーど: 逆像問題，サロゲイション，データ秘匿， ……
• 最初は津田さん． graph kernel と pre-image (これの訳が逆像?) 問題， 例題は drug discovery である構造の化合物の特徴空間に写像する問題の逆問題， 特徴空間内のある位置にくるような化合物の構造を決める問題．
• 次は坂野さん． 手書き文字認識， GA で切りはりして， さらに ソフトウェアに質問を作らせる (そのソフトウェアが自分の能力を高めるための)， といった手法．
• 休憩．
• 星野さんによるデータ秘匿の review． お役所とかがもっているデータ， 公開したいけれど個人情報が特定されないようにしたい， といった状況でどういう統計学的な工夫が必要か， といったハナシ．
• 伊庭さんのレアイベントをサンプリングする MCMC 工夫のハナシ． 私はいままでレプリカ交換法とマルチカノニカルの区別がついて いなかった， と判明した ……
• 最後は平田さんで， 「ならびかえ」を駆使して「帰無仮説的な」 時系列データを生成するわざの解説． いったんフーリエ変換したりとか …… しかしこれとて「ならびかえ」に援用するためなのだが．
• ということで， 1745 ごろ終了． たいへん勉強になる， しかも元気のでる研究会でした．
• さて， 本日中に札幌に帰れるか? 1755 統数研発． 1805 多摩モノレイルで高松発． 1808 立川着． 500 円よぶんにはらって 1815 特急あずさに乗ってみたけど …… 1844 新宿． そんなに速くなかった? 1846 同発． 1906 品川． 京急にのりかえて 1907 同発． 1920 蒲田着． 空港いき急行にのりかえて 1928 同発． 1942 羽田．
• ちょっとぼーっとしてたら 20 時にいきなり あちこちの食べもの屋が閉店してしまった． 売店でベイグルサンドイッチを買って， 2020 晩飯．
• 2050 ANA4727 B-737-500 (AIR DO の機体) に搭乗． ほぼ満席． これが一番おそい便なんだけど， こっちにしといてたすかった． とゆーのも， ひとつ前のやつは機材がとどいてなくて， これよりも遅くなってしまったよーで． 2110 離陸． 2230 新千歳空港着． しかし空港内が大雪で混乱していて， なかなか降機できない． 空港内を走ったのだけど …… しかし雪のため最終の快速エアポートの到着も遅れていて 2250 発のはずが 2259 発． 吹雪． のろのろ運転． 2355 札幌． 56 分もかかった …… いつもより 20 分ぐらい遅い．
• このあたりで 1/7 (金) 終了， 桑園いき普通列車も遅れまくり． ホーム寒い． 2406 JR 札幌駅発． 2408 桑園着． 2420 帰宅．

[桑園からの帰路]

こんなかんじで積雪してた．

• [今日の運動]
  • うんどう休養日， つづく．
• [今日の食卓]
  • 朝 (0730): ハムサンドイッチ． バナナ．
  • 昼 (1230): こんびににぎり飯．
  • 晩 (2030): 羽田空港のベイグルサンドイッチ．

2011 年 01 月 08 日 (土)

• 0800 起床． 雪がどんどんふってる． 朝飯の準備． 朝飯． コーヒー．
• 生態学会な雑用．
• 怠業．
• 昼飯の準備． 昼飯の買いだし． 昼飯．
• 怠業．
• 1/12 (水) の水産資源データ解析研修会の準備． しんどい．
• 買いもの． 晩飯の準備． 晩飯．
• 怠業．
• [今日の運動]
  • うんどう休養日， つづく ……
• [今日の食卓]
  • 朝 (0900): パンケイキ． ヨーグルト． リンゴジャム．
  • 昼 (1310): モンクールのパン． ダイコン・ニンジン・セロリ・油揚のコンソメスープ．
  • 晩 (1920): 米 0.5 合． コマツナ・牡蠣の炒めもの． シュンギク・ネギ・鶏肉のサラダ．

2011 年 01 月 09 日 (日)

• 0820 起床． ねぼうした． 大雪． 朝飯の準備． 朝飯．
• 雪がどかどかふってるので， 外出する気になれづ． また 1/12 (水) の水産資源データ解析研修会の準備のつづき， 自宅にて ……
• 昼飯の準備． 昼飯． ちょっと日がさしてきた．
• 午前 2 時間の 講義的な部分はいつもの階層ベイズモデルハナシでごまかすことにして …… 午後 3 時間の「水産データ解析らいぶ」 は HTML ペイジなどを準備して， と ……

• 日没ごろ， ばててきたので気分転換に AEON 桑園まで買いものにいく． ここの カルディコーヒーファーム の豆は安い! とりあえず 200g 買ってみた． おいしければよいのだけど．
• 晩飯の準備． 晩飯．
• スルメイカとマガレイとケガニのデータ解析をつづける ……
• 生態学会大会雑用， しつこく．
• [今日の運動]
  • 朝のラジオ体操だけ． うんどう休養日の日々， つづく．
• [今日の食卓]
  • 朝 (0940): モンクールのパン． ダイコン・ニンジン・セロリ・油揚のコンソメスープ．
  • 昼 (1300): ネギ・コマツナ・鶏肉うどん．
  • 晩 (1930): 米 0.7 合． ピーマン・ブナシメジ・鶏肉のタイカレー． ナガイモのポン酢あえ．

2011 年 01 月 10 日 (月)

• 0650 起床． 朝飯の準備． 朝飯． コーヒー．
• ちょっとまた雑用あれこれ．
• そして水産資源データ解析研修会の準備のつづき． 今日が資料の提出しめきりか? いちおー祝日なんだが …… あとで確認しよう．
• 本日提出の場合にそなえて， 最低限の部分だけはつくりあげてしまって， と．
• とりあえず， 最小限の内容を web 上で見られるようにしておいて …… 研修会の幹事のヒトに 「情報はここに何でもあります」 といった時間かせぎ的なメイルを送る．
• 反応がないので， やっぱり明日だせばよかったのかな?
• 買いもの． 昼飯． 晩飯の準備．
• ぢりぢりと水産データ解析研修会の準備作業のつづき …… しかしすでにばてぎみ．
• 晩飯の準備． 晩飯．
• 自宅ちかくの温泉．
• その後も明後日の研修会の準備のつづき．

• [今日の運動]
  • 今日はらぢお体操すらやらなかった ……
• [今日の食卓]
  • 朝 (0800): 米 0.5 合． ダイコン・ニンジン・セロリ・油揚のコンソメスープ． ピーマン・卵の炒めもの．
  • 昼 (1230): うどん．
  • 晩 (1830): 米 0.5 合． 豆腐・ワカメの味噌汁． ホッケのみりん干し． サトイモ・ゴボウ・ニンジン・コンニャク・鶏肉の煮物．