ぎょーむ日誌 2012-06-(01-10)

2012 年 06 月 01 日 (金)

• 0620 起床． ラジオ体操． 弁当・朝飯の準備． 朝飯． コーヒー． 0750 自宅発． 晴． 0820 研究室着．
• 統計モデリング入門 の謹呈本， 昨日夕刻から本日にかけて， 日本のあちこちに到達したもよう． はい， 岩波書店で発送を「忘れられていた」あの謹呈本です．
• みなさんへのお返事かき． けっこう， Amazon とかで発注しているひとが多くて …… ひたすらおわび．
• 北大祭関係のメイルもめんどうなことに．
• 昼飯．
• また別の統計モデリング本でも書きませんか (分担執筆) という提案をいただいたんだけど， まあこれはお断りすることに．

  「これはぜひ!」と反射的にお返事を書こうとしたのです
  が，少し正気にもどりまして，やはり今回はやめておいた
  ほうがよいだろうと結論しました．
  
  私は執筆速度がきわめて遅く，今回の教科書かきには正味
  二年間以上を費やしてしまい，そのあいだ本業の研究はほ
  とんど停滞しておりました．このために，あちこちにいろ
  いろと学内・学外で問題が生じつつある状況です．しばら
  くは「解説かき」の仕事は断念し，研究仕事のほうに専念
  しようと考えております．
  
  せっかくの機会を与えていただきながら，もうしわけあり
  ません．
  

• なぜかしら， 甲山さんの Mac であれこれ設定作業． これは私が web server を勝手に北大くらうどに移転したことに由来するのだが …… Mac の scp クライアントってどうなってるの?
  • 「そんなの fugu だろ」 と思って universal 版をいれてみたのだけど (その Mac は intel 機だったので)， なぜかしらエラーをだして起動せづ (64 bit 機だから??) ……
  • FileZilla はよけーな機能ばかりで， ふつーに使えない
  • ともかく選択の余地がなさすぎる ……
  みなさん， Linux 使いみたいにターミナルから scp とかやってるのか?
• けっこうばてました．
• その後は Twitter にふりまわされてみたり ……
• 1700 研究室発． 晩飯． 2000 帰宅．
• [今日の運動]
  • 朝のラジオ体操．
• [今日の食卓]
  • 朝 (0700): シリアル． ヨーグルト．
  • 昼 (1240): 研究室． 米 0.7 合． 塩ザケ． ゴボ天． トマト． ミズナ．
  • 晩 (1830): なぜかしら S 2 W 8 のクネルで晩飯．

2012 年 06 月 02 日 (土)

• 0830 起床． よくねた． ラジオ体操． 朝飯． コーヒー．
• 怠業しつつ， 雑用メイルかき． とどこおっていたぎょーむ日誌かき．
• 正午をすぎても怠業つづく． シリアのアサド一族にすこしくわしくなってしまった． 昼飯のパンを買いにでる． 知事公館 の庭でおそい昼飯．
• また怠業．
• 1800 自宅発， ミニ大通公園走． 1840 帰宅．
• 晩飯の準備． 晩飯．
• [今日の運動]
  • 朝のラジオ体操． ミニ大通公園 走 40 分間．
• [今日の食卓]
  • 朝 (0900): シリアル． ヨーグルト．
  • 昼 (1520): 知事公館でモンクールのパン．
  • 晩 (1930): 米 0.5 合． キャベツ・ニンジン・タマネギ・ソーセイジのポトフ． アスパラガスのバター炒め．

2012 年 06 月 03 日 (日)

• 0740 起床． ラジオ体操． 朝飯． コーヒー．
• うーむ， 怠業がつづく ……
• 昼飯．
• 午後も怠業．
• 晩飯の準備． 晩飯．
• お， アマゾントップセラー という tweet になぜか掲載された． 瞬間的に， 分野内で一位になったみたいだ． これまたよくわからぬ順位づけだけど．

• [今日の運動]
  • 朝のラジオ体操．
• [今日の食卓]
  • 朝 (0750): シリアル． ヨーグルト．
  • 昼 (1200): 米 0.5 合． キャベツ・ニンジン・タマネギ・ソーセイジのポトフ．
  • 晩 (1930): 米 0.5 合． シジミの味噌汁． ヒラアジ (すごく小さいやつ) の煮物． ホウレンソウ・タマネギ・卵の炒めもの．

2012 年 06 月 04 日 (月)

• 0620 起床． ラジオ体操． 弁当・朝飯の準備． 朝飯． コーヒー． 0740 自宅発． 0810 研究室着．
• 今週の 岩波ベストテン， 統計モデリング入門 はまた理科系・単行本で 2 位． 「宇宙の謎」はつよい ……
• 北大祭関係の雑用メイルかき． けっこう時間かかるな ……
• ちょっとうんどう． 1200 研究室発． 北大構内走． 1230 研究室もどる． 昼飯．
• 午後は研究室紹介ポスター (これまた北大祭がらみ) を自動生成するしくみづくりの改訂作業 ……
• 統計モデリング入門 ってホントに本屋に並んでるのかなあ …… などとぶつぶつ言ってたら， 写真を送っていただきました!

[東京駅の丸善]

村瀬さんから送っていただきました (元の画像ファイル)． ありがとうございます． こんなポップを作ってくれた店員さんもありがとう ……

• 1750 研究室発． 買いもの． 1820 帰宅． 晩飯の準備． 晩飯．
• そうそう今晩は部分月食だ．

• 夜もポスター生成プログラミング ……
• [今日の運動]
  • 朝のラジオ体操．
• [今日の食卓]
  • 朝 (0700): シリアル． ヨーグルト．
  • 昼 (1200): 研究室． 米 0.7 合． 鶏からあげ． タコからあげ． インゲンマメ． ミニトマト． レタス．
  • 晩 (1900): 米 0.5 合． ネギ・豚肉の味噌汁． カツオのたたきとタマネギサラダ． マイタケ・ブナシメジのバター醤油いため． キュウリのピクルス．

2012 年 06 月 05 日 (火)

• 0620 起床． ラジオ体操． 弁当・朝飯の準備． 朝飯． コーヒー． 0740 自宅発． 0810 研究室着．
• 北大の就業規則とやらが変更されたとかいうので， そのweb page を見にいったんだけど， 給与削減 0.5% となっている． 文科省の「もっと下げろ」圧力に抵抗してるのかな?
• なぜかしら， 大阪大学からこういう 133.1.8.28 - - [03/Jun/2012:11:06:44 +0900] "GET / HTTP/1.1" ... "osaka university[QSC-0985]" アクセスが数万回も来襲． IP アドレスはたしかに大阪大学のもの． 何かネット徘徊プログラムが暴走しているのか?
• セミナーの予習． うーむ ……
• 研究室セミナー， 今日は雨谷君の論文紹介なのだが …… 予習段階でわかっていたことなんだけど， 4 本ともヘンなものだった． ひとつめは図をいくつか並べただけでスカンディナヴィア半島の 植生と環境変動の将来についてあれこれ述べているし (しかもその内容が図示された情報とは異なる)， ふたつめは統計モデリングがおかしい (分布・リンク関数がヘン， random effects のいれかたがおかしい) だけでなく二種類ある応答変数のつながりを考えてないので意味不明な結論となってる， みっつめは GIS + 分類木 + ROC/AUC なる 「環境系お役所論文ここにあり!」の典型で 地図ばかりだして版面をハデにしてるくせにデータ解析では空間構造もへちまもなく 説明変数を過剰に複雑怪奇に組み合わせた分類木をこしらえて「あてはまりの良さ」 に自己満足しているだけ (お役所は無意味なはざーどマップごときものを， ひたすら機械的に生成するのが仕事なので， これで問題ないのだろう)， よっつめは奇妙な仮定による樹齢推定とかひとつの枝からとったサンプルを 「いずれも完全独立な反復」とみなして「ゆーいゆーい」とか …… これが， これが「生態学・環境科学のデータ解析」だ!
• 昼飯．
• 午後はまた北大祭関係の雑用 ……
• 研究室紹介 web 版 と連動する紙版ポスター生成のしくみのつくりなおす …… 本来は作業を省力化するためのものなのだが， PStricks 関連をぬきとるのに苦闘する …… まだ不完全ではあるが， かなり単純化された． ただし， 画像まわりに何かまだ不備があるらしく， 生成した PDF ファイルは 49M なる巨大なしろものに ……

• 1750 研究室発． 買いもの． 1820 帰宅． 晩飯の準備． 晩飯．
• [今日の運動]
  • 朝のラジオ体操．
• [今日の食卓]
  • 朝 (0700): シリアル． ヨーグルト．
  • 昼 (1200): 研究室． 米 0.7 合． 鶏手羽先の赤ワイン・八角・ショウガ煮 (わけわからん)． タコからあげ． インゲンマメ． ミニトマト． レタス．
  • 晩 (2000): 米 0.5 合． レタス・油揚げのスープ． ズッキーニ・ニンジン・アスパラガス・トマト・エノキダケのラタトウィユ． カマスの干物．

2012 年 06 月 06 日 (水)

• 0620 起床． ラジオ体操． 弁当・朝飯の準備． 朝飯． コーヒー． 0750 自宅発． 晴． 0810 研究室着．
• 昨晩いただいた， おもしろい tweet ……

• また北大くらうどの設定にはまる． ntpdate ができないというだけなんだが …… port 123/udp がひらかん!
• その試行錯誤の過程でみつけた， なんとも脱力な北大おんらいんすとれいじ …… そのじったいは， これまた意味不明な北大 application server なので， このちんけな interface は無視して， rsync とかで自動バックアップ体制をつくればよさそう．

• 今日も研究室・北大祭ざつよう ……
• 昼飯． いきなりつよい雨．
• 論文改訂相談．
• おお， ようやく統計モデリング本 の出版契約書が岩波からきた …… えっ? 印税 8% なの? 10% だと思っていたんだが ……

[けいやくしょ]

まあ， アパートの部屋の契約書みたいなもので． しかしなあ …… R グラフィックス のときは， 作業にかかる前に共立出版から送られてきたような気もするのだが …… さてさて． なんでも「統計モデリング入門」 が電子化されれば印税は 12% だそーです． はたして， 岩波書店が理系本を電子化するかどうか， という点はともかく．

• 再来週の出張の手配をまったくやっていなかった …… ぢたばた． 統数研ののろわれな申請が ……
• 1730 研究室発． 買いもの． 1800 帰宅． 晩飯の準備． 晩飯．
• [今日の運動]
  • 朝のラジオ体操．
• [今日の食卓]
  • 朝 (0700): シリアル． ヨーグルト．
  • 昼 (1210): 研究室． 米 0.7 合． サバの塩焼き． いかねぎボール． ミニトマト． レタス．
  • 晩 (2030): 米 0.5 合． ネギ・モヤシ・豚肉の味噌汁． サバの塩焼き． ズッキーニ・ニンジン・アスパラガス・トマト・エノキダケのラタトウィユ． コマツナ・エノキダケのゴマあえ．

2012 年 06 月 07 日 (木)

• 0730 起床 ……… うわー! ねぼうした! 弁当なしの日でよかった! 0750 自宅発． 0820 研究室着． ぐったり朝飯 & コーヒー．
• やっぱり寝る時刻が重要なのかも …… まあ， 弁当なしの日ということで油断してましたが．
• Amazon のらんきんぐって移動平均か何かにもとづいているんだろうけど， 「マンガでわかる統計学」が強すぎ ……

• お礼状に対するお礼状かきとか ……
• 鍋嶋さんデータに少しとりくんでみる ……

• なぜかしら北大正門前の「駿」で昼飯． ついでに， N 9 W 3 (北大正門でて東にまっすぐ) にある小映画館 蠍座 の「蠍座通信」もとりにいく． なんと， 17 年間固定していた料金を， やむにやまれず値上げとか． たとえば一本だけ見るのは 800 円だったのが 900 円に …… Twitter でもいろいろ言及されている．
• 統計モデリング入門 直販サイトのプログラミング …… なぜかしら， すごく時間がかかっている．
• 1750 研究室発． 買いもの． 1820 帰宅． 晩飯の準備． 晩飯．
• [今日の運動]
  • 朝のラジオ体操．
• [今日の食卓]
  • 朝 (0700): シリアル． ヨーグルト．
  • 昼 (1200): 北大正門前の「駿」の昼定食． なまのビンチョウマグロをノリでまき， さらに卵を焼きながらその外にまきつけた 「ピカタふう」なる不思議な料理をたべた ……
  
  • 晩 (2000): 米 0.5 合． ネギ・ダイコン・シイタケの味噌汁． ゆで豚肉とミズナのサラダ． 卵焼き． キュウリのピクルス．

2012 年 06 月 08 日 (金)

• 0620 起床． ラジオ体操． 弁当・朝飯の準備． 朝飯． コーヒー． 0750 自宅発． 0820 研究室着．
• なぜかしら， なかなか進捗しない 統計モデリング入門 のちょくはんサイト準備 ……

[今日の店頭写真]

丸善丸の内本店 に陳列されている写真を送っていただきました． ありがとうございます．

• 昼飯．
• 午後はまず，雑用 ……
• オジロワシ繁殖の統計モデリングこんさる．
  1. データが少ないわりには， GIS な説明変数が過剰というありがちな状況 ……
  2. glm() + stepAIC() などやると， どう見ても無意味な (まったく解釈のできない) 説明変数を結合したようなモデルが AIC 最良となる (呪われな GIS 的説明変数!)
  3. とりあえず， データのとりかたをよく検討し， 「地域差」という効果があることがわかったので glmmML() を使うことに
  4. しかし現時点の最新の glmmML() …… これはまたしても以前のように stepAIC() できなくなってしまっている!
  5. しょーがないので 2003 年ごろ につくっていた自作の stepAIC.glmmML() 関数を復活させることに …… これ は少し修正しただけで使えた! (後記: この関数はかなりヘンなもので， stepAIC.glmmML("y ~ x1 + x2 + x3", cluster = d$id, data = d, ...) といった引数指定をします． 数式が formula ではなく文字列!!)
  6. しかし呪われな GIS 的説明変数の個数はなかなか減らない
  7. 呪われな GIS 的説明変数どうしの正負の相関を pairs() で図示して確認
  8. 論文の method section における「これこれこういう理由で」 という部分の正当化を考え， 「少数個の説明変数を使ったモデルたち」 という範囲で AIC によるモデル選択， これにて一件落着
  というかんじで， とりあえず何とかなった． しかし， また後日， これまた呪われぎみになりがちな「気象要因」データを 「どう捨てるか」について研究しなければならぬことに．
• 16 時 30 分ごろから 1.5 時間ほど， 北大祭の展示雑用 …… ぐったり．
• 1900 研究室発． 買いもの． 1930 帰宅． 晩飯の準備． 晩飯．
• [今日の運動]
  • 朝のラジオ体操．
• [今日の食卓]
  • 朝 (0700): シリアル． ヨーグルト．
  • 昼 (1220): 研究室． 米 0.7 合． 興部ハーブソーセイジ． アオノリ天． キュウリのピクルス． ミニトマト． ミズナ．
  • 晩 (2000): 米 0.7 合． 冷凍保存しといたチキン野菜カレー．

2012 年 06 月 09 日 (土)

• 0820 起床． よくねた． ラジオ体操． 朝飯． コーヒー． 0930 自宅発． 0950 研究室着．
• 北大祭の環境科学院展示の現場で， あれこれと …… これで私のぎょーむはだいたい終わり．
• 気がゆるんだので， 北大クラーク会館で散髪， 前回は 4/25 なので， 少しマがあいてしまった． いつもどーり 2625 円． さっぱりした．
• ついでに北大生協会館の書籍部によってみたら， 統計モデリング入門 が一冊だけおかれてるのを， よーやくにして目撃できた …… まあ， たくさん売れたあとの「補充待ち状態」だと信じたい． あとで twitter (@KuboBook) にその写真をあげてみたところ， さっそく伊庭さんが何やら思いつかれたようで ……

• 環境科学院にもどって， 午後の サイエンストーク を依頼してましまった小泉さんの陣中見舞い． 考えてみれば， 大学院生に依頼すればよかった． そうすれば， その院生には一日ぶんの謝金 (7000 円ぐらい?) が支払われるから．
• 1340 研究室発． ぱん吉でパンを購入． 1410 帰宅． 昼飯．
• ちょっと怠業．
• 怠業， つづく ……
• 1740 自宅発， ミニ大通公園走． たいへんすばらしい完全曇天． しかも涼しい． 1820 帰宅． 体重 73.8 kg． うんどう不足ぶとりだ．
• 晩飯の準備． 晩飯．
• いつものごとく仕事のすすまぬ土曜日．
• [今日の運動]
  • 朝のラジオ体操． ミニ大通公園 走 40 分間．
• [今日の食卓]
  • 朝 (0900): 米 0.5 合． 昨晩の冷凍カレーの残り．
  • 昼 (1420): ぱん吉のサンドイッチとパン．
  • 晩 (1920): 米 0.5 合． ネギ・油揚げの味噌汁． ダイコン・ニンジン・シイタケ・鶏肉の煮物． ヒメタラ． キュウリのピクルス． ミニトマト．

2012 年 06 月 10 日 (日)

• 0730 起床． よくねたけど， 最近なにやらストレスフルなのか， 寝てるあいだに歯ぎしりやっていて， これはものすごい筋力を使っているので， 朝になると首・肩が筋肉痛になっている …… ラジオ体操． 朝飯． 紅茶．
• 統計モデリング入門 の直販サイトづくり， ぢりぢり …… 注文うけつけの部分はもうできているんだけど， そのあとのフェイズのところでツマっている． 私が入金を確認し， 岩波の出荷担当た「送信しました」ボタンを押して終了するナガれなのだが ……

• おっと雲が日光をさえぎっているうちに走っておくか …… 1120 自宅発． ほぼ完全な曇天． 今日は宮の森の坂道をのぼったりおりたり …… 1220 帰宅．
• ミニ大通公園の東端にある倫敦館で昼飯．
• 買いものして帰宅． また家をでて， 大学へ． 北大祭展示の撤収． 私がやったことは， 次の展示である open campus (8 月) にそなえて各コースのポスターを強引に回収してまわっただけ …… 生態系生物学コース，つまり低温研には逃げられたが．
• 研究室にもどって， ちょっとめんどうなデータ解析こんさるメイルの読解． こりゃあ難物だなぁ．
• 1755 研究室発． 1810 帰宅． 晩飯の準備． 晩飯．
• [今日の運動]
  • 朝のラジオ体操． 宮の森坂道走 60 分間．
• [今日の食卓]
  • 朝 (0800): ぱん吉のパン． チーズ．
  • 昼 (1340): 倫敦館のポークソテーライス．
  • 晩 (2000): 米 0.5 合． ネギ・シジミの味噌汁． 豚肉・タマネギハンバーグ． 蒸し野菜 (アスパラガス・ジャガイモ・タマネギ)． レタス．