ぎょーむ日誌 2014-01-(21-31)

2014 年 01 月 21 日 (火)

• 夜中の 2 時ごろ， とら泣き． またまた手足がかゆいようなので， ワセリンを塗る．
• 0550 起床． コーヒー． 朝飯・弁当の準備． とらめし． 朝飯 …… 洗濯もの干しとかいろいろ作業があって， 0900 自宅発． 0910 保育所着．

[北大生協クラーク会館店りにゅーある]

ちょっと寄り道して， のぞいてみた． 半年以上， 工事してたんじゃないかな ……

[2F 書籍部]

建物は新しくなったけど， 1F と 2F の売り場面積そのものは以前とはあまりかわらないような (cf. フロア概要 PDF file) …… まあ， 「多目的スペース」とか追加されたようなので， そういう点が強化されたんでしょう．

• 1010 研究室着．
• 1030 より 研究室セミナー， 今日は野村さんの修論概要の発表みたいなかんじ ……
• 議論してたらけっこう長くなった． 昼飯．
• あれこれ作業．
• Hari さんのデータ解析こんさる …… 北大の天塩研究林のデータなのだが， まあ， 時系列構造とかあって， なんともめんどくさいもので． いろいろ議論して， randomization で信頼区間を評価できるんじゃなかろうかというほうこうで …ｚ
• それについての R coding ……
• 1730 研究室発． 保育所でとら回収． 1810 帰宅． とら風呂． とらめし． 晩飯．
• [今日の運動]
  • とらちゃんえくささいず．
• [今日の食卓]
  • 朝 (0630): パン．
  • 昼 (1230): 研究室． 米 0.7 合． フリーズドライスープ． リンゴ．
  • 晩 (1910): 米 0.9 合． カブの味噌汁． コマツナ・ブナシメジ・豚肉の炒めもの．

2014 年 01 月 22 日 (水)

• 0530 とら泣きで起床． 朝飯・弁当の準備． とらめし． 朝飯． コーヒー． 今日も洗濯ものにとりくんでから， 0830 自宅発． 雪． 0840 保育所着．
• 先週金曜日， 保育所の新年交流会のコマまわし担当メンバーのひとりだったので， とら送迎で保育所内をうろうろしてると， ぜんぜん知らない園児たち (ただしボーイズ限定) から「あのコマ，大きかったね，すげー」 「あれ回すのむずかしい?」 「おれも来年はまわすんだ (注: 最年長組の子供だけがコマまわししていいらしい?)」 …… などなど． 現時点では， 「コマまわしおじさん」と認知されてるよーで．
• 0920 研究室着． おや， いきなり晴れたな．
• 昨日の R coding のつづき …… けっこうややこしい．
• 昼飯．
• なんとかできた ……

• とりあえず， 報告メイルかき ……
• あとから考えると， この信頼区間の構成方法には問題がありそうだとわかってきた …… 「時系列そのもの」をリサンプリングしてから「時系列の平均」 をとるのはいいんだろうけど， 各時刻で CI を構成するのはダメではないかな …… たぶん， ある「時系列な区間」を設定すると， (たとえば) 95% の時系列はその区間内からまったくはみ出さないけれど， 残り 5 % は少なくとも 1 時刻ではみだす， といった区間構成が必要なんじゃなかろーか．
• このあたりの問題は理論とアルゴリズムの両方が難しいかも． たとえば， この本……

  統数研 web site で統計学本 PDF 版を無料配布している EBSA http://t.co/aDQREBJ2PJ たとえば 2012 年出版 5600 円の「21世紀の統計科学< Vol. III >数理・計算の統計科学」 http://t.co/b676SlifiJ

  — 久保拓弥 (@KuboBook) 2014, 1月 23

  これの下平さんの 章 を見ると， 時系列みたいな「構造のあるデータにおけるブートストラップは まだ発展段階にあり，いまのところデータに応じて十分な検討が必要である」 といったことも書かれてるしね．
• お， もうじゅぎょーの時間だ．
• 1630 から 3 人あいてにじゅぎょー …… まあ， へんな総合科目みたいなやつなんで， 受講者が少ないのはどうでもよいのだが …… どうでもよくないのは， 自分が半年ほど前にいろいろ作った講義資料 (投影資料) の「ぼろぼろ」さ …… 当時， かなりアタマがおかしくなっていたのか?
• 育児ばてとかで 「作業できる量が激減したなあ」 とはわかっていたんだけど， 「作業の質もすごく低下していた」 ということに， 今さらながら， よーやくにして気づいた …… まあ， 私の場合， 量・質の改善どちらも時間と体力が必要なんで， ここ一年間以上のアウトプットがずたずたなのは， どうしもようないのかもしれない ……
• 1720 研究室発． 買いもの． 保育所でとら回収． 1750 帰宅． 晩飯の準備． とら風呂． とらめし． 晩飯．
• 体重 69.2 kg， 少し上昇しているな …… えーと， ここ数年は 70-72 kg ぐらいの範囲をうろうろしてたんだけど， 「とらばて」でしばらく 68 kg 台に落ちていた， という状況です． とらだっこでひへーした背筋はまだ痛くて， 湿布を毎日はっております．
• ついでに「だっこ測定」したとらの体重は 11 kg 弱かな …… 8 月以降， ずっと 10 kg 台のまま． まあ， 重くなったらだっこがますますしんどいですが ……
• [今日の運動]
  • とらちゃんえくささいず．
• [今日の食卓]
  • 朝 (0630): 米 0.7 合．
  • 昼 (1130): 研究室． 米 0.7 合． ミニトマト．
  • 晩 (1900): 米 0.7 合．

2014 年 01 月 23 日 (木)

• 5 時ごろとらなき …… あれ， 数分間のだっこでわりとあっさり寝た?
• 0600 起床． すぐにとら起き． とらめし． 朝飯・弁当の準備． 朝飯． コーヒー． 0830 自宅発． 晴． 0840 保育所着． 0910 研究室着．
• 「時系列の信頼区間」みたいなことをもう少しぐだぐだと考えてみる …… 進展なし．
• 昼飯．

[石田さんから本を送っていただきました]

ありがとうございます． らのべふうデータ解析本 第 2 巻， この本では因子分析その他あれこれをタテ書きで解説! …… 今回もすごく売れそうだな ……

• 時系列本やブートストラップ本， さらにクラスタリング本などながめる ……
• 15 時から雨谷君の時系列データ解析こんさる …… まあ， library(repeated) の GLMM 推定関数ではどうこうできる問題ではないわけで．
• とりあえずこういう複雑な構造のデータ解析をやるときには 図を描きましょうよ …… ということで， あれこれ作図法と data.frame 操作わざを実演しつつ説明 ……
• そんなこんなでおそくなってしまった …… 1730 研究室発． 保育所でとら回収． 1810 帰宅． とら風呂． とらめし． 晩飯．
• [今日の運動]
  • とらちゃんえくささいず．
• [今日の食卓]
  • 朝 (0630): パン．
  • 昼 (1130): 研究室． 米 0.7 合． 鶏肉のグリル． ブロッコリー． ミニトマト． リンゴ．
  • 晩 (1910): 米 0.7 合． ハクサイ・昆布の味噌汁． ニンジン・ゴボウ・タマネギ・鶏肉のバルサミコ炒め．

2014 年 01 月 24 日 (金)

• 0600 起床． こちらが起きたのを察知するのか， すぐにとら泣き起き． 朝飯・弁当の準備． とらめし． 朝飯． コーヒー． とらだっこで背筋がいたむので湿布． 0825 自宅発． 0835 保育所着． 0910 研究室着． 今日は 2 週間ぶりぐらいに「脱」真冬日， となるはず．
• 昨日の院生データの検討あれこれ ……
• 昼飯． ついでに散髪． 前回は 2013-11-29 で 2 ヶ月ちかくマがあいてしまった．
• 今度は Hari さんの修論データ， 天塩の森焼きあれこれ時系列データ解析の方法についての説明など ……
• …… でオワってしまった． 1720 研究室発． 買いもの． 保育所でとら回収． 1750 帰宅． 晩飯の準備． とら風呂． とらめし． 晩飯．
• [今日の運動]
  • とらちゃんえくささいず．
• [今日の食卓]
  • 朝 (0630): 米 0.7 合． ハクサイ・豆腐の味噌汁． 納豆．
  • 昼 (1230): 北大生協．
  • 晩 (1900): 米 0.7 合． ハクサイ・豆腐の味噌汁． ホッケ． ホウレンソウのゴマあえ． ナガイモすりおろし．

2014 年 01 月 25 日 (土)

• 0550 起床． とらもすぐにおきてきたので， とらめし． 朝飯． コーヒー．

[とら散髪]

前髪がちょっとヘンかな …… ということで， これから補正するところ …… この作業は簡単だったけど， その後に散らばった髪の毛の掃除とかに時間をくって ……

• 朝はとらと遊んでいたら， だんだん眠たげなかんじになってきたので， 10 時すぎに寝かしつけ．
• ちょっとのんびり．
• かんたん昼飯 …… 食べ終わったら， とら泣き． 昼寝からめざめたもよう．
• とらめし …… よくたべるなあ．
• 14 時すぎ， とらおんぶして耳鼻科にいく． 一ヶ月以上マがあいてしまったわけだが …… いまのところ中耳炎再発している様子はないんだけど， なんだか不安な気分になったので．
• 土曜日の診察終了 20 分前という絶妙な時刻だったせいか， まち時間ほぼゼロで診察してもらえた …… 右耳は問題なし， 左耳は少し危なくなっている， ということで抗生物質を出してもらうことに．
• 買いものなどして， 16 時ごろ帰宅． コーヒー．
• 晩飯の準備， あれこれ．
• とらめし． 晩飯．
• 寝かしつけに 40 分ほど苦闘 …… 慢性的な筋肉痛 ……
• [今日の運動]
  • とらちゃんえくささいず．
• [今日の食卓]
  • 朝 (0630): 「どんぐり」の角食． チーズ． ベイコン．
  • 昼 (1150): かけラーメン．
  • 晩 (1910): 米 0.7 合． ネギ・豆腐の味噌汁． キャベツ・ソーセイジの蒸し煮． サトイモ・ニンジン・鶏肉の煮物．

2014 年 01 月 26 日 (日)

• 0550 起床． 朝飯の準備． とらめし． 朝飯． コーヒー．
• 今日もとら坊あいてにいろいろやっているうちに終了 ……

[とら画伯]

外は吹雪なので， 家の中にとじこもってばかり …… 初めて， とらに鉛筆をもたせてみたら， このように絵を書いてくれました． 指がちょっとぼろぼろしてるのは， あとぴーです (モクタール塗で改善中)．

• [今日の運動]
  • とらちゃんえくささいず．
• [今日の食卓]
  • 朝 (0630): 「どんぐり」の角食．
  • 昼 (1130): らでぃっしゅぼーやのパン． キャベツ・ソーセイジの蒸し煮．
  • 晩 (1920): 根菜蒸し煮のトマトソース鍋． チーズ雑炊．

2014 年 01 月 27 日 (月)

• 3 時すぎ， 「かゆいよー」 …… とら泣き． 15 分ほどだっこうろうろしたら寝てくれた ……
• 0550 起床． 朝飯・弁当の準備． とらめし． 朝飯． コーヒー． 0825 自宅発． 晴． 昨日からまた真冬日が再開している． 0835 保育所着． ひよこ組のみなさん， 園庭の除雪機が働いている様子を窓にはりついて観察している …… 0910 研究室着．
• …… えーい， 脱力しているばあいではない ……
• とはいえ， 進捗よろしくないまま昼飯．
• かたづけられるものから， かたづける ……
• ハイマツデータの統計モデリングにさきだち， いろいろと図を作らねばならない …… ところが， 院生の皆さん， なかなかこういうことをやらない or できないんだよね …… とりあえず， R coding 例もつけとくから， 参考にしてください．

panel.cor <- function(x, y, prefix="", method = "kendall"){
   r <- cor(x, y, method = method)
   txt <- sprintf("%.2f", r)
   txt <- paste(prefix, txt, sep="")
   par(usr = c(0, 1, 0, 1))
   text(0.5, 0.5, txt, cex = 1 + abs(r) * 4, col = "#0000ff")
}
panel.xy <- function(x, y) points(x, y, col = "#ff4000")
ds <- d[d$mountain == "AK",]
d.med <- data.frame(ss.med = tapply(ds$ss, ds$ye, median))
d.med <- cbind(d.med, unique(ds[,5:9]))
pairs(d.med, lower.panel = panel.cor, upper.panel = panel.xy)

• 1710 研究室発． 1740 帰宅． とらめし． 晩飯．
• [今日の運動]
  • とらちゃんえくささいず．
• [今日の食卓]
  • 朝 (0650): 米 0.7 合． ダイコン・油揚の味噌汁．
  • 昼 (1230): 研究室． 米 0.7 合． タマネギ天． サトイモ・ニンジンの煮物． ブロッコリー． ミニトマト．
  • 晩 (1900): 米 0.9 合． 根菜・鶏肉のスープ．

2014 年 01 月 28 日 (火)

• 0610 起床 …… ちょっと寝ぼう， アタマが混乱していて「今日はにちようびだ」と 朝飯・弁当の準備． 体内時計が 6 時すぎおき習慣に adjust されつつあるとら坊がおきてきたので， とらめし． 朝飯． あれこればたばた …… 0820 自宅発． 雪． 0830 保育所着． 0910 研究室着．
• ハイマツデータ作図のつづき．

• 昼飯．
• ここまでの作図わざや data.frame 操作わざを説明しつつ， 雨谷君とハイマツモデリング相談 …… この統計モデルを BUGS 言語で記述するためには， かなりの工夫が必要とされるような． いちばんやっかいなのは， 時系列データなのに， 開始時刻・終了時刻がサンプルごとにばらばらであること ……
• 1720 研究室発． 保育所でとら回収．

[キツネ?]

なぜか保育所の玄関にキツネの氷像がおいてありました ……

• 1800 帰宅． とらめし． 晩飯．
• [今日の運動]
  • とらちゃんえくささいず．
• [今日の食卓]
  • 朝 (0630): 米 0.7 合． ダイコン・油揚の味噌汁．
  • 昼 (1230): 研究室． 米 0.7 合． タマネギ天． ブロッコリー． ミニトマト． リンゴ．
  • 晩 (1910): 米 0.7 合． マッシュルーム・鶏肉のトマトスープ． 紅茶豚． ハクサイ・ベイコンの炒めもの．

2014 年 01 月 29 日 (水)

• 3 時すぎ， とら泣きおき． だっこうろうろ …… 30 分ほどで寝てくれた．
• 0600 起床． 朝飯・弁当の準備． とらめし． 朝飯． コーヒー． 0815 自宅発． 晴． 気温はかなり低い． 0825 保育所着． 0900 研究室着．
• R の package とか， ちょっと調べる ……

  library(glmulti) なる package があり http://t.co/ekS3v82DRw glm() 的な関数あれこれの線形予測子を取捨選別して「可能なモデル」一式を作ってAICなどを評価するみたい…自分では使いたいと思わないけど，院生に聞かれたとき用にメモ…

  — 久保拓弥 (@KuboBook) 2014, 1月 29

• データ解析こんさるメイル書き． GLM 的なモデルの説明変数が「方位」の場合はどうするか? …… その場で考えついた策としては， 方位を南北方向と東西方向に分割するとゆーもの …… 推定されたふたつの係数から「重要な方向」がわかるでしょう． この方式の欠点は，南北と東西に相関が生じること， そして「重要な方向」なるものがひとつに限定されてしまう， とか ……
• 抜本的な解決策は，いつものことながらベイズモデル化で， 角度を離散化 (1 度きざみとかに) して，「統計モデリング入門」 第 11 章の一次元空間モデルをリング状にしたものを考えて， 重要な方向みたいなものとかを発見すること ……
• 昼飯． なぜか， 卯和 でおいしいごはんを ……

[うーむ，まさかのビッグウォーズ再開?]

昼飯後に北大東の古本屋の前で見かけたポスター …… 展示・イヴェントの web site もできている …… これ， ずいぶんカネがかかっていそうだな． ナゾなかんじ． 道立文学館 は， なんだかじみなかんじで．

• Dell Latitude E6500 という古い PC で動作させている Ubuntu 13.10 の動きががくがくぎみのような気がするので， swap の設定とかをいろいろ試行錯誤してみたんだけど …… 改善されないなあ．
• 1730 研究室発． 保育所でとら回収． 1750 帰宅． とらめし． 晩飯．
• [今日の運動]
  • とらちゃんえくささいず．
• [今日の食卓]
  • 朝 (0630): 米 0.7 合．
  • 昼 (1130): 研究室． 米 0.7 合． ミニトマト．
  • 晩 (1900): 米 0.7 合．

2014 年 01 月 30 日 (木)

• 2 時すぎ， とらおき． しばらくだっこうろうろしたら， 寝てくれた．
• 4 時すぎ， とらおき． しばらくだっこうろうろしたら， 寝てくれた …… 何がとら安眠をさまたげているのか， ぜんぜんわからないので不安． かゆいのは， だいぶなおったみたいなんだけど ……
• 0550 起床． 朝飯・弁当の準備． とらおきの前にすばやく朝飯． その後， とらおき …… 今朝は牛乳がなくなっていたので， 豆乳をあげたら， すごくふきげんになって泣きつづける …… ミカンあげたらおさまった． とらめし． コーヒー． 0815 自宅発． 今日は「脱」 真冬日であたたかくなりそう …… 0825 保育所着． とら坊もふくめ， ひよこ組のみなさんごきげん． 0900 研究室着．
• どーにもこーにも， ちょうしが ……
• 昼飯．
• どーにもこーにも， ちょうしが ……
• 1720 研究室発． 買いもの． 保育所でとら回収， 書類提出あれこれ． 1800 帰宅． とらめし． 晩飯．
• 22 時すぎ， とら泣きおき． ドコノコノキノコ&コノキノコドコノ で再ねかしつけ ……
• [今日の運動]
  • とらちゃんえくささいず．
• [今日の食卓]
  • 朝 (0630): 米 0.7 合． カブ・エノキダケの味噌汁．
  • 昼 (1230): 研究室． 米 0.7 合． サケ． ブロッコリー． 加賀あげ． ミニトマト． ミカン．
  • 晩 (1900): 米 0.7 合． カブ・エノキダケの味噌汁． カスベの煮付け．

2014 年 01 月 31 日 (金)

• 0550 起床． 朝飯・弁当の準備． 早おきとら坊もおきてきたので， とらめし． 朝飯． コーヒー．
• さて， ここで …… とら坊先生がいきなりあれ (公式隠匿符丁は (う)) をふんばりはじめ (しかもおしりが痛いので泣く) …… 30 分以上がんばって， なんとか終了しました． いやはや， たいへんだね (今回は特に長かったわけですが) ……
• 0830 自宅発． 昨日は気温 +8 度とかひじょーしきな高温だったけど， 昨日夕方から単調に気温は低下しつづけており， 現在は -3 度ぐらい． 0840 保育所着． 0920 研究室着．
• 論文原稿よみぎょーむ．
• 昼飯．
• D 論予備審．
• 個体群生態学はとくに難しいハナシはないんだけど (多少の想像力は必要でしょうけど)， 個々のモデルがどのように・どれいぐらい「現実ばなれ」しているか， その距離感を正確に知っておかないとね ……
• 森林の毎木調査データとゆーのは， まあ簡単な構造なんだけど …… たとえば R なんかで自由自在にあつかうのは， 慣れていないとけっこう難しいようですな． 垂直構造・平面構造・時間構造・群集構造， これらをいろいろな断面ですぱすぱ切ったりできないと， 何がなんだかワケわかんないでしょう．
• 1730 研究室発． 買いもの． 保育所でとら回収． 1810 帰宅． とらめし． 晩飯の準備． 晩飯． とら風呂．
• iTunes …… というか iPhone 上の Music player 上で「どーしても消せない曲」 が出たときの 対処法あれこれ …… 壊れた曲ファイルを iCloud から上がきして， それをフリック消去するとゆーわざでうまくいった．
• [今日の運動]
  • とらちゃんえくささいず．
• [今日の食卓]
  • 朝 (0630): 米 0.7 合． ナメコ・カブ・油揚の味噌汁．
  • 昼 (1230): 研究室． 米 0.7 合． アサリチャウダー． 魚肉ソーセイジ．
  • 晩 (1900): 米 0.7 合． 紅茶ブタ・タマネギソース． ナメコ・カブ・油揚の味噌汁．