ぎょーむ日誌 2014-10-(11-20)

2014 年 10 月 11 日 (土)

• とらだっこで背筋をすこし痛めたので， ここ数日， 夜中に何度か目がさめるのだが …… 2 時ごろおきたら， 何だかおどろくよーな依頼メイルがきていたので， ちょっと寝れなくなった． あれこれ検索してみたり ……
  • RT @berobero11 : 行動計量学会のStanの講師の岡田謙介先生「伝えるための心理統計」の著者でしたか。ベイズMDSのBUGSコードとか資料のページにあったので後でチェックします。www3.psy.senshu-u.ac.jp/~ken/misc.html   02:22:52
  • RT @izumeeeeeeeeeel : 行動計量学は統計の授業ではなく、ただマニュアル通りにパソコン操作をする授業だしな。   02:23:16
  • RT @millionsmile : 日本行動計量学会なるところで発表する予定。3月16日。チラシができてるようだ。[pdf] bit.ly/12Vh5tu   02:26:54
  • RT @kan_yukiko : #TokyoWebmining 第22回 にて「広告クリエイティブ最適化の実際」発表させていただきました。先日、行動計量学会にて山川と発表した内容とほぼ同じです。貴重な機会を賜り、ありがとうございました！　 www.slideshare.net/kan_yukiko/22-…   02:31:29
• 0600 起床． ねぶそく感あり …… ねむい． コーヒー． 朝飯・弁当の準備． 朝飯． 「よくねた，よくねた」 と言いつつとら坊がおきてきたので， とらめし．
• とらの機嫌がよくなったり悪くなったりの変動がはげしくて， ちょっとしんどい …… 0825 自宅発． 晴． 風が強くて寒い． 0835 保育所着． とら， 土曜保育． 0910 研究室着．
• だめだ， ねむい …… せっかく土曜日も研究室にきたのに．
• ふちょーぎみなので (このあたり帰宅後・とら坊ねかしつけ後に書いてるわけだが)， また twilog はりでごまかすことに．
  • 地震情報 11時36分頃発生 最大震度：4 - 日本気象協会 tenki.jp #jishin @tenkijp www.tenki.jp/bousai/earthqu… 札幌はたぶん震度 2 ぐらいなんだろうけど (現時点でまだ情報なし)，A 棟 8F は震度 3 ぐらいのゆれだった   11:41:13
  • 結局，札幌市は震度 2 だったとのこと www.tenki.jp/bousai/earthqu… ← これはまだ更新されてないけど．   11:57:40
  • JAGS 並列化はそう簡単ではないとわかったけど時間切れ未解決放置な金曜日 - ぎょーむ日誌更新 hosho.ees.hokudai.ac.jp/~kubo/log/2014… #KuboLog   12:21:32
  • RT @teramonagi : 状態空間モデルまとめ(PDF)。シミュレーションスムーザの記述あり。　RT State Space Models in R - Journal of Statistical Software bit.ly/1C3ppZS #メモランダム   13:10:03
  • だれも使わぬ library(lattice) のハコひげ図 bwplot() わざ: 自分でハコ・ひげの位置を決めたい場合は panel 関数内の panel.bwplot() 引数で stats = my.boxplot.stats というふうに自作関数を呼び出すよう指定する   15:07:47
  • やっぱり，かくのごとき関数内の関数内で関数名を指定…みたいなことなんかを要求するから lattice 作図は普及しないんだろうか? ggplot2 ならもっとラクに描けるのかな…と lattice 作図脳なワタシは苦闘 1.5 時間の過程をふりかえりつつお茶を呑むのであった……   15:11:46
  • . @MayumiPod 自分で panel 関数を書くユーザー個人を見たことないなあ……もちろん印刷物・web とかに掲載されてる lattice 図とかで，そういう「仕事」がなされた出力を見たことはあるけれど (えがき人しらず)．   15:17:24
  • 先ほどまでの lattice 作図苦闘はとあるデータ解析の最終結果を図示するものなんだけど，結果は現物あつかっている皆さんの期待に反したものに……科学はそうしたオチになりがちなんだけど，私の作業によってそれが明確になってしまって，なんだかもうしわけないきぶんに．   15:37:14
  • まあ，あほみたいなコーディングミスとかないか，きちんと点検しよう……   15:37:39
  • . @berobero11 ありがとうございます! 自分でも少し ggplot2 勉強しようと docs.ggplot2.org/current/geom_b… ながめてみたんですが，たぶん stat 引数指定するだけで終了しそう…こういう簡単な図なので… pic.twitter.com/W2SX6lywGJ   16:31:58
  • 
  • やはり library(lattice) は歴史的経緯をひきずって複雑なのかな? library(ggplot2) のどくとくな "+" 演算子記法は病的とゆー印象だったんだけど，ひとつひとつの作図要素の複雑化が回避される…といったご利益があるのかなー? よくわからんけど…   16:36:55
  • @berobero11 まあ，科学の世界だとそこまで酷薄ではないかも……とデータ解析ねつぞうのはびこる生物学ぎょーかいに身を沈めているわたくしがいうのもあれですが……   16:40:52
  • おっと，のんびり土曜ぎょーむもここまでか……買いものして，保育所から土曜保育なとらを連れて帰る，と……   16:55:50
  • RT @berobero11 : @KuboBook 一応描いてみました。いつもまさにそういったマニュアルのページにあるコードをコピペしたりして作っています。data frameをreshapeするところだけが慣れないときつかった覚えがあります。 pic.twitter.com/uACMI3LCA3   16:59:08
  • @berobero11 ありがとうございます! reshape2 の melt() の使い方までわかり，おとくな気分です．あとで，自分でもやってみよう……   16:59:53
  • とら、よーやくねた… ちょっとふらふらするので検温してみたら、当方の体温37.3℃まで上昇…うーむ、ちょっと厚着してみるかな…   20:42:25
  • RT @ibaibabaibai : ちなみに，ベイズ流の検定なんて知らないと思うけど，あなたの使ってるかもしれないBICとかMDLとかいうのは実質それですよ．リンドレーパラドックスって知ってる？   21:49:15
  • RT @ibaibabaibai : リンドレーパラドックスは「ベイズ統計と統計物理」でちょっと触れている．無情報に近いときにベイズ因子でモデル選択することへの疑問，あるいは頻度論的検定との不一致のこと．   21:49:53
  • RT @ibaibabaibai : 「ベイズ統計と統計物理」は高級とされる話題も実はいろいろ入っていて，リンドレーパラドクスの件も某駒○先生にも「いやああいう普通は知らないことが書いてあって・もちろん（僕らは）知ってるわけですけど・入門の本にはなかなか書いてないことが出てるのは良いと」みたいに言われた記憶がある．   21:49:59
  • RT @ibaibabaibai : 「戦後のベイズ統計が華やかなりし頃」に書かれたリンドレーの教科書だって，普通の検定をベイズ流に置きかえる勇気はなくてそこは微妙になってる．というか読者が困るからだろう．主観確率についてはもっと過激にサヴェジ風にすればよかったと序文で書いてるけど．   21:57:42
  • RT @ibaibabaibai : 本当に新しい統計の入門書ってのは根本的に従来のものを見直した上で，自分とは違う分野での実用とか現に社会で行われている通法とかの折り合いをどうつけるかまで考えないといけないから尋常なことでは無理．なのでまず統計科学のフロンティアを出し，久保本を出し，上からリニューアルしていく計画に   21:59:51
  • RT @ibaibabaibai : 久保本もglm本やmcmc本でなく「新しい統計の入門書」にしたかったので，ガウス分布の例入れろ，ふつうの回帰分析ももっと扱ってくれ，と言ったんだけど，もうひとつ反応は鈍かった．「既成」へのアンチの面があると，「既成」を知らない人にはそこは意味がないから「既成」を置き換えられない．   21:59:56
• [今日の運動]
  • ラジオ体操: 1 回 …… とらのおかげで
• [今日の食卓]
  • 朝 (0630): 米 0.5 合． ゴボウ・ニンジン・鶏肉の炒めもの．
  • 昼 (1200): 研究室． 米 0.7 合． 野菜スープ． 豆腐．
  • 晩 (1840): 米 0.7 合． ゴボウ・ニンジン・鶏肉の炒めもの．

2014 年 10 月 12 日 (日)

• 0600 起床． 今日の要約: とらと動物園に行ったんだけど， 途中から風邪症状がでてきて， 夕方以降は倒れていた．
• ばてぎみなので， また芸もなく Twilog はりつけ ……
  • RT @ibaibabaibai : @berobero11 「操作全体」での詳細つり合いをちゃんとチェックすること．自信の無いときはチューニングは予備のrunでやって本番は固定，というのが基本です．   14:45:22
  • RT @berobero11 : ブログ更新しました。でも読まなくていいです。自分のケジメのために更新してるんです。→[WinBUGS] 傾向スコアを使ったモデル 続き heartruptcy.blog.fc2.com/blog-entry-161…   14:50:32
  • ちょっとふちょーぎみだけど細胞生物学なデータ解析最終出力作図にとりくむ土曜日 - ぎょーむ日誌更新 hosho.ees.hokudai.ac.jp/~kubo/log/2014… #KuboLog   15:22:12
  • 午前はとらとバス移動&円山動物園、ホッキョクグマ・ララの水泳なるレアイヴェントを観測、帰路の中央バス「散策バス」に乗るあたりからふらふらしてきた。現在体温37.6℃、ふとんにもぐりこんでおります。とら「げんきになれよー」とのげきれい… pic.twitter.com/P7HcY21TWZ   16:23:22
  • 
  • RT @ibaibabaibai : 頻度論的な信頼区間を検定を含むように厳密に定義するのは意外と難しくて初心者向きではない．でもとりあえずなんちゃって誤差棒でいいのでは．ましてベイズの信頼区間であれば事後分布がわかれば即定義できるわけだし，とりあえずエラーバー主義もありでは．   17:09:55
  • RT @ibaibabaibai : 重回帰って１）作った回帰式全体で予測する　２）回帰式の係数に意味を持たせる，ってあって，その２つでかなり違うよね．もう少し内容的に言い換えると１’）もともとと似た状況での予測　２’）全然違う状況や自分が介入した場合についての予測，となるかもしれない．後者が因果推論ということになる   17:14:30
  • RT @berobero11 : 統数研みたいな研究所に勤めて、スパコン借りつつ、色々な企業や人からデータ解析コンサルを受けて解析して、たまに分からんときには近所の研究室に聞きに行くような生活が理想。   17:26:37
• [今日の運動]
  • ラジオ体操: 0 回
• [今日の食卓]
  • 朝 (0640): トウモロコシのパンケイキ．
  • 昼 (1200): ベイコン・コマツナのパスタ．
  • 晩 (???): えーと ……?

2014 年 10 月 13 日 (月)

• 風邪． ねたきりでした． twitter だけやってた， みたいな …… また芸もなく Twilog はりつけ ……
  • 熱があまりさがらないので、家の中ひとりで寝たきり状態です…少しはマシになってきたかな？まあ、今日はのんびり休もう…   12:22:41
  • @TJO_datasci おみまいありがとうございます。かぜには気をつけましょう…   12:27:53
  • @TJO_datasci なが患いでしたね…さっさと回復するべく今日はひたすら寝ます。   12:34:19
  • RT @mas__yamazaki : マレーシア航空１７便の撃墜事件で、オランダのティマーマンス外相と同国の検察当局は、乗客１人の遺体が酸素吸入マスクを装着していた状態で発見されたことを明らかに（CNN） bit.ly/1vVLIQg 「マスクを着ける時間があったことを意味し、即死しなかったことを示唆」   13:18:28
  • RT @ibaibabaibai : 「信頼区間」というのは「データを入れると信頼区間が出てくるルール」なんです．アルゴリズムとか関数といってもいい．推定値ってのはみんなそうじゃないか，っていわれそうだけど，信頼区間については特にそれを強調する価値があると思う．   14:47:19
  • RT @ibaibabaibai : 点推定値のまわりに「誤差棒」をつけるのなら，その推定値をモデルに入れて，何回もデータを発生させたと考えればよい．実際にそうやっても良いし，解析的な近似も可能かもしれない．これが私が「なんちゃって誤差棒」と呼んでいるもので普通はこれで十分なことも多い．   14:47:42
  • RT @ibaibabaibai : 本当は点推定値にも誤差があるので，「真の値」を入れたモデルからデータを発生させたい．でもわからない．大げさにいうと，自分の誤差を自分で評価できるのかという自己言及めいた問題になる．   14:49:48
  • 「信頼区間」というのは「データを入れると信頼区間が出てくるルール」 …という baibai さん指摘、生態学分野の glm() 使っている人はほぼ理解してない。むしろ事後分布みたいな脳内モデルが受容的で、そのセンに沿って説明すると直感にあうみたい…私の説明のしかたが悪いのかな？   15:07:50
  • RT @ibaibabaibai : 信頼区間についてはベイズ信頼区間（確信区間）のほうがわかりやすいと思う．なんちゃって誤差棒ですませばそれほど難しくはないように思うが，それでも多少は心理的抵抗がありそう．そもそも「調査」「抜き取り」とかでないのにサンプリング誤差とは何かというより根本的な疑問も出そうだし．   15:31:39
  • RT @ibaibabaibai : 検定を重視する分野だと有意性と誤差のサイズ・エフェクトの大きさが同時に表示できるのはありがたいだろうと思う．「なんちゃって誤差棒」だとそうはならない．   15:35:29
  • 生態学分野の場合、信頼区間生成プログラムを与えてもエフェクトサイズなんぞ気にかけてくれるかどうか…？95%区間内にゼロあるや否や、統計学は丁半バクチと心得、単純明快が信条のゆーい差決戦主義えころぢすとあなどれませんね…「ゆーいじゃない？もっとデータ取ろう」という往生際の悪さもあり   15:58:44
  • RT @ibaibabaibai : ただ「仮説検定など日常の論理ではない」となるとどうでしょう．　「この空の様子は普通ではありえない」という人は「いま見ている空が，私が「普通の空」と呼ぶ帰無仮説から生成された確率は大変小さい．よってこれは正常な空ではない」と言ってるんじゃないかな．   16:04:44
  • RT @ibaibabaibai : 確信区間でなく信用区間なのかな．覚えられないので「ベイズ信頼区間」でいいような気がしてきた．   16:33:12
  • それなりにこみいっていて日常感覚になじまない Neyman-Pearson な検定のわくぐみが生態学分野で普及したのは不思議…たぶん「p < 0.05 ならあんたのハナシは正しいと主張可能」なる誤解が集団内に固定し、p だけ見てればいいという「安心感が共有」されているんだろうね…   16:33:37
  • エフェクトサイズの検討が不人気な理由は、解釈にあたって使用している統計モデルの理解が必要とされるからかな？…p 値のみ重視するゆーい差決戦主義なら「キいてる」「キいてない（誤解釈）」だけで議論終了。かんたん。   17:28:43
  • . @TJO_datasci そうですね…主語は必要なので「○○はキいてる」と表現されてはいます。ただし、ごくまれには、解析者の意図と推定値の正負の符合が逆であることに気づいてないことも…   17:54:37
  • マイスターとら師匠の #今日のレゴ 作品 、作品名は「カバとヘリコプター」、みぎの緑・青のブロックがカバ、その他はヘリコプターだそうです。 pic.twitter.com/C4N76BRHPs   19:58:34
  • 
  • RT @TJO_datasci : 「彼女ができるデータマイニング」をお願いします — 彼女ぐらいデータマイニングなしで作れよ ask.fm/a/bc590qhj   20:03:23
  • これ、かなりおもしろいハナシでした。データベイスあるところデータ解析あり⁈ 私がオンラインデートを攻略した方法 www.ted.com/talks/amy_webb…   20:06:57
  • 生態学ゆーい差決戦主義史観ざつだん…10-20年ぐらい前(かしら?)，順位統計量をつかった検定いわゆるノンパラメトリック検定がはやり，これはある意味では「(人間に解釈可能な)統計モデルなし，検定だけ」だったので，「p<0.05 だけを確認すればヨシ」というナガレが加速したのかも?   21:40:02
  • 分散分析はいまだに生態学でよく使われるデータ解析，これは線形・正規分布モデルに特殊化した尤度比検定みたいなものとみなせる．ソフトウェアの出力は分散分析表なるもので，伝統的には「効果の大きさ」表示なし，分散比とF分布で算出されたpが示されるのみ．これまた統計モデル理解が不要の世界?   21:49:41
  • . @ibaibabaibai うーむ，粕谷さんはかなり早い時期から順位統計量検定の「万能ではない」論文を出したり，何しろ柔軟な勉強家なので一般化線形モデルの本 www.kyoritsu-pub.co.jp/bookdetail/978… も出しておられて…「のんぱら信仰」は安きに流れる解析者たちの暴走ですかね…   22:00:32
  • @ibaibabaibai 改訂されてないみたいですね．私の勝手な憶測ですが，たぶん改訂する必要なしと考えておられるのかも……ぴんく本目次 www.kinokuniya.co.jp/f/dsg-01-97848… をみても「のんぱらあおり本」というかんじでもないわけですし……   22:11:53
  • @ibaibabaibai 私も目次をひさしぶりにみて「あれっ? こんな構成だったっけ?」と意外に思いました．あした研究室にいったら，ぴんく本をながめてみることにします……   22:17:50
  • とりあえず，おかげさまで風邪気味状態を脱出しつつあり，ざつだん tweet をあれこれかけるぐらいになりました．明日も６時おきなので，さっさと寝ることにしよう……ふちょーぎみだと眠りも浅くなりがちだし．   22:19:55
  • @ibaibabaibai ……ということで，「順位統計量の確率分布」を使った次第に衰退しつつある（しかし絶滅しない）手法ではないかと思います．ならびかえや bootstrap を使った手法は発展していくと思いますが．   22:27:25
  • 歴史的ぢゃーごん補筆: かつてある時代において生態学をふくむいくつかの分野で以下のよーな用語が使われていた……ぱらめとりっく統計: 正規分布を使った検定のこと，のんぱらとりっく統計: 正規分布を使用しないで順位統計量を使った検定のこと．「統計」と「検定」はほぼ同義．   22:34:17
  • RT @ibaibabaibai : @KuboBook 考えてみると，ならびかえやサロゲーションは帰無仮説の作り方で，順位は検定統計量なので，論理的には独立でした．尤も，順位検定が好きな人というのは自分が「なにをランダム・意味なし」とするかをデザインしようとはあまり思わない傾向かな，という気はします．   22:36:37
  • 統計連合大会のチュートリアル「マッチングと統計解析」聞いてかなり勉強になったんだけど，「傾向スコアを使ってマッチングさせると p 値が0.06 から 0.04 になりました」みたいなオチの説明がわりとあって，「p 値好きな学問分野」は広大なのかもしれない，とゆー気分になりました．   22:48:17
• [今日の運動]
  • ラジオ体操: 1 回 …… たしか， とらにあやつられて一回やったような?
• [今日の食卓]
  • 朝 (???): たしか， くってない
  • 昼 (???): シリアルすこしずつ， 時間不定
  • 晩 (1830): 米 0.7 合． コマツナ・豆腐の味噌汁． 鶏肉・サツマイモのショウガ炒め． ニンジンサラダ …… おかげさまで， だいぶ回復しました

2014 年 10 月 14 日 (火)

• 風邪ぎみ状態つづく …… ときどき， 目がさめて寝られなくなる断続的な睡眠．
• 0600 起床． 朝飯・弁当の準備． 朝飯． コーヒー．
• 0810 自宅発． 雨． 0820 保育所着． 0850 研究室着．
• 研究室の自己紹介セミナー．
• 新しく入った研究生・院生 5 人， 全員外国からきた人たち & つゆざきさん研究室所属， そのうち 4 人は女性， そのうち 3 人は中国あちこちからきた人たち， そのうち 2 人は日本語うまい ……
• 昼飯． いろいろなヒトに風邪を心配される． 声がかなりおかしいようで．
• 研究室のメンバー表など作りなおしてみて， 院生・研究生の過半数はすでに外国からきた人たち， というかんじ．
• ぢりぢりとぎょーむをすすめる ……
• 1730 研究室発． 買いもの． 1800 帰宅． 晩飯． とら風呂， 拒否． とら寝かしつけ …… ぐったり．
• [今日の運動]
  • ラジオ体操: 0 回
• [今日の食卓]
  • 朝 (0630): 米 0.5 合． コマツナ・エノキダケの味噌汁．
  • 昼 (1200): 研究室． 米 0.7 合． モズクスープ．
  • 晩 (1850): 米 0.7 合． コマツナ・シイタケの味噌汁． コマツナ・タマゴ・豚肉のショウガ炒め． タラコ．

2014 年 10 月 15 日 (水)

• 0600 起床． 朝飯・弁当の準備． 朝飯． コーヒー．
• 0810 自宅発． 0820 保育所着． 0850 研究室着．
• 今日は R の 「ダイジョブ」 変換問題でもがくことに …… しかし library(stringi) を教えてもらえてよかった!
  • R のlibrary(Nippon) に数字英字の全角→半角変換関数 zen2han ある (カナは変換できない)．逆に半角カナの全角化は? > chartr("[ｱ-ﾝ]", "[ア-ン]", "ﾊﾙﾆﾚ") [1] "セナシニ" まだ工夫が必要．濁点・半濁点という難問もあり…   12:29:18
  • 呪われ半角カナの変換問題，ひょっとして前処理特攻隊のみなさんなら，ひょっとして解決策をすでに開発しておられるとか……?   12:33:50
  • さきほどのしくじり chartr() 解説． > chartr("[ｱ-ﾝ]", "[ア-ン]", "ｱｶｻﾀﾅ") [1] "アェガケザ" となる理由は半角・全角の文字コード表がこうなっているから． pic.twitter.com/pTwNaN6H06   12:46:39
  • 
  • @berobero11 たしかに R だけでやるのは，かなりたいへんそうな……   12:50:26
  • . @TJO_datasci @berobero11 ありがとうございます．ちなみに私が好きな Perl だと use Jcode; Jcode::convert(\$str, $ocode, $icode, "z"); ですよ! ……え? だれも Perl つかわない?!   12:57:58
  • 「R で半角→全角変換」について質問してきたヒトは，R での変換に見切りをつけてゑくせるの「半角→全角変換」 bit.ly/1DaSBAs を使うことにしたそうで……なんだろう，この敗北感は……   13:04:07
  • RT @kohske : @KuboBook こういうのも有ります。 > stringi::stri_trans_nfkd("ｱｶｻﾀﾅ") [1] "アカサタナ" windowsで使えるのかわかりません。   13:07:32
  • RT @kohske : @kubobook windowsでもﾀﾞｲｼﾞｮﾌﾞですね。   13:07:37
  • . @kohske ありがとうございます! library(stringi) cran.r-project.org/web/packages/s… 知りませんでした．今，ダウンロードしてコンパイルしてるんですけどかなり時間かかりますね…終わったら，いろいろいじってみます．あ，終わった…   13:09:59
  • RT @kohske : @kubobook しかし環境によっては濁点が独立１文字になってしまう場合もあるので、その点の対処は別途ひつようかもしれません。   13:11:20
  • RT @kohske : @kubobook この辺りの闇です。 d.hatena.ne.jp/miau/20110805/…   13:14:32
  • @kohske NFC/NFD という新しい闇を知ることになってしまいました……   13:16:38
  • ふむふむNFC/NFD

    > stri_trans_nfkd("ﾀﾞｲｼﾞｮﾌﾞ")
    [1] "タイショフ"
    > stri_trans_nfc("ﾀﾞｲｼﾞｮﾌﾞ")
    [1] "ﾀﾞｲｼﾞｮﾌﾞ"
    > stri_trans_nfkc("ﾀﾞｲｼﾞｮﾌﾞ")
    [1] "ダイジョブ"
    

      13:26:53
  • これは Ubuntu 14.04 (Linux) 上の R の library(stringi) の挙動の一例です．   13:27:56
  • @MayumiPod こちらではこんなかんじ

    > v <- c("ﾀﾞｲｼﾞｮﾌﾞ", "ダイジョーブ", "ﾀﾞｲｼﾞｮｰﾌﾞ"); stri_trans_nfkc(v)
    [1] "ダイジョブ"   "ダイジョーブ" "ダイジョーブ"
       
    13:35:23
    

  • RT @kohske : @KuboBook 裏でICUに投げてるので、JavaとかPerlとかPythonのICUバインディングと同じ結果になるとはおもいますが　、この辺りは正気の人間なら手を出してはいけない闇ですね。   14:07:47
  • @MayumiPod あらら，sapply() 使えばよかったのに……

    v <- c("ﾀﾞｲｼﾞｮﾌﾞ", "ダイジョーブ", "ﾀﾞｲｼﾞｮｰﾌﾞ") sapply(v, stri_trans_nfkc)
    

      14:16:09
  • . @kohske ICU en.wikipedia.org/wiki/Internati… というソフトウェアがあるんですね．調べてみると Perl の Encode も ICU の table を参照してるような記述が……現代におけるバベルの塔ですね．   14:25:12
  • @MayumiPod なるほど……R の version 依存ではないと思うんで，CVS ファイル化するときに文字コートを UTF-8 など unicode を指定すればいいのかな……じつは，このあたり，よくわかっていない．   14:26:44
• 1720 研究室発． 買いもの． 1740 帰宅． 晩飯．
• とらねかしつけ． とらはエルゴでおんぶうろうろすると， わりと「すかっ」と寝てくれて助かるのだけど， 本人もおんぶするほうも「エルゴに依存」な状態が継続している …… ちょっとこのあたり， 不安．
• 保育所のお昼寝では保育士さんたちが何か魔術を使って， 10 人以上の二歳児を手早く寝かしつけているらしいんだけど (とらは保育所でも難物らしいのですが) …… ということで， 保育士さんたちほどうまくはできないだろうけど， 「おんぶ以外のなにか」 で寝かしつけることは可能， なはず ……
• 今晩の挑戦， 創作どきゅめんたりー童話「とらと動物園」， とらが円山動物園内を巡回しながら， いろいろなどうぶつたちと会話を楽しむハナシ …… 「ねないー! ねないー! おんぶー!!」 と絶叫していたのに， ハナシ始めたとたん沈黙， 動物園をひととーり巡回しおえたところで， 「もいっぺん，もいっぺん……」 とつぶやきながら寝おちてしまいました …… めでたし， めでたし．
• しかし， これおんぶ以上にくたびれるな ……
• [今日の運動]
  • ラジオ体操: 0 回
• [今日の食卓]
  • 朝 (0630): 米 0.5 合．
  • 昼 (1200): 研究室． 米 0.7 合． サバ水煮缶詰．
  • 晩 (1830): 何やらイタリアンそうざい， あれこれ ……

2014 年 10 月 16 日 (木)

• なんとなく 0520 起床． ちょっと朝ぐだぐだ …… 朝飯・弁当の準備． 朝飯． コーヒー．
• 0820 自宅発． 曇天． 0830 保育所着． いきなり， あひる組の女の子から「だめー!」連発のだめ出しをくらう …… とらの教育について， なにか問題を指摘されたのだろうか? よくわからぬまま， 0900 研究室着．
• ちょっとまた論文修正作業のうだうだ ……
• 昼飯．
• また別の修正作業をぢりぢり ……
• 1720 研究室発． 買いもの． 1750 帰宅． 晩飯の準備． 晩飯． とら風呂． とら寝かしつけ． また「どきゅめんたりー創作どうわ」 …… 「むかしむかし，あるところにとらちゃんがいました， とらは家の近くのバス亭から JR バスに乗って ……」
• [今日の運動]
  • ラジオ体操: 0 回
• [今日の食卓]
  • 朝 (0630): 米 0.5 合．
  • 昼 (1200): 研究室． 米 0.7 合． カレー豚肉． ネギ・シイタケの炒めもの． ミニトマト．
  • 晩 (1830): 米 0.7 合． ざんぎ．

2014 年 10 月 17 日 (金)

• 0600 起床． 朝飯・弁当の準備． 朝飯． コーヒー．
• 0805 自宅発． 雨． 0810 保育所着． 0850 研究室着．
• 今日もいまいちなかんじだな …… 風邪の症状はかなりおさまってきたんだけど ……
• おや， 科研費申請関連の連絡が ……
• そちらはひとくぎりつけて， 動物生態の小泉研 の院生・山崎君のヤツメウナギ産卵行動データのハナシを拝聴 (あぽなしなところが動物生態) …… 興味ぶかい， 一種の時系列データではあるんだが …… さて， どうなることやら．
• 雨谷君の ハイマツモデルの試行錯誤． 15 時から雨谷君に現状説明．
• うーむ， 強引に「各年・各枝独立のばらつき」 の大きさを指定してもダメらしい， ということがわかってきた …… 状態空間モデル， いつものことながら， 複数種類の「ばらつき」をうまく推定させるのが難しい．
• 1720 研究室発． 保育所でとら回収． 買いもの． 1800 帰宅． 晩飯． 今晩も話芸だけでとらを寝かしつけてみる ……
• [今日の運動]
  • ラジオ体操: 0 回
• [今日の食卓]
  • 朝 (0630): 米 0.5 合． ざんぎ．
  • 昼 (1200): 研究室． 米 0.7 合． 豆腐．
  • 晩 (1830): キムチなべ．

2014 年 10 月 18 日 (土)

• ちょっと風邪っぽい状態つづいているけれど， 午前はとら坊と円山動物園． 気温が低くなっているせいか， オランウータンなど熱帯の動物たちは， 室内あそびをしていた． トラとか寒い地方の動物たちは， 「やっとすずしくなった」 というかんじで元気に歩きまわっている．
• [今日の運動]
  • ラジオ体操: 0 回
• [今日の食卓]
  • 記録・記憶なし …… そうだ， ひるめしは円山「ざぼうあん」でうどんとか食べたな ……

2014 年 10 月 19 日 (日)

• 今日は 藻岩山 (標高 531 m) に行ってきました． 2011 年に新設されたミニケイブルカー「もーりす」に乗って山頂へ．

• 山頂からはなんと! 大雪山系の山々 …… 旭岳 (直線距離 142 km ぐらい) やトムラウシが見えました． 手前の芦別岳とかも． よくわからなかったけど日高山系の山も見えるとか．

• [今日の運動]
  • ラジオ体操: 0 回
• [今日の食卓]
  • えーと …… 昼ごはんは中腹駅で Burg のサンドイッチたべたなー ……

2014 年 10 月 20 日 (月)

• 0540 起床． 朝飯・弁当の準備． 朝飯． コーヒー．
• 0800 自宅発． 曇天． 0815 保育所着． 0840 研究室着．
• メイルかきあれこれ ……
• 昼飯．
• 細胞生物問題， どーしたものかと手のうごきが止まっているところ …… 一種類の「生起確率」だけなら， 群間比較はオッズ比をみるだけなので簡単なのだが， 今回比較しないといけないのは， この生起確率を三つもかけあわせた確率なので …… 単純にオッズ比に帰着させて， という策が使えない． うーむ ……
• 進展しないまま， 1720 研究室発． 買いもの． 1750 帰宅． 晩飯． ラジオ体操．
• いつものごとく， 「とらねかしつけ」 の創作どきゅめんたりー童話「むかしむかしあるところに，とらちゃんがいました……」． 寝るまで 20 分以上かかる場合が多いような気がするんだけど， まあおんぶうろうろで寝かしつけるよりは双方にとってよいだろう， ということで．
• [今日の運動]
  • ラジオ体操: 1 回
• [今日の食卓]
  • 朝 (0630): 米 0.7 合． サトイモのスープ．
  • 昼 (1200): 研究室． 米 0.7 合． ヒメタラ． ダイコンの煮物．
  • 晩 (1900): 海鮮丼．