[[過去の日録]] * 2008年 9 月 18 日 (木) 曇 [#bdd29fa2] - 06:10 起床 -- 台風が接近中 -- 07:30 研究室着 - 午前 -- 文長にポアソン分布を仮定して,Bayesで予測を行ってみる. ポアソンの共役事前分布はガンマ分布である(Gelman et al.: Bayesian Data Analysis, p.52, 2nd Edition). ---ポアソン分布 #mimetex(p(y|\theta) = \frac{\theta^y \, e^{-\theta}}{y!}) --- 尤度 #mimetex(\bf{y} = y_i, \dots , y_n) #mimetex(p(\bf{y}|\theta) = l(\theta|\bf{y}) = \Pi \frac{1}{y_i!} \theta^{y_i} \, e^{-\theta}) #mimetex( \propto \theta^{t(y)} \, e^{-n \theta}) 指数属で表現すると #mimetex( \propto e^{-n\theta}\, e^{t(y) log \theta} ) 文系学生相手にするので,指数法則を指摘する. #mimetex(\theta^y = e^{log(\theta^y)} = e^{y log \theta}) 指数属とは関数が次の形で表されること(Lee: Bayesian Statistics, pp.60). #mimetex(p(\bf{x}| \theta) = g(\bf{x}) \, h(\theta) \, e^{t(\bf{x}) \, \phi(\theta)}) 上の式では&mimetex(\frac{1}{y!});が消えているので注意 ---自然共役事前分布 #mimetex(p(\theta) \propto (e^{-\theta})^n e^{v log \theta} ) --- 事後分布 #mimetex(\theta|\bf{y} \sim Gamma(\alpha + n \bar{y}, \beta + n)) -- 医学関係ではポアソン分布を次のように書いている. #mimetex(y_i \sim Poisson(x_i \theta)) &mimetex(y_i); と &mimetex(x_i); が観測数で,&mimetex(\theta); が未知のパラメータだが,&mimetex(x); を exposure,&mimetex(\theta); をrate と表現しているので,分野違いの人間は,何度みても慣れないというか,馴染めない. * 2008年 9 月 17 日 (水) 曇 [#hd82accb] - 06:10 起床 -- 07:30 研究室着 - 午前 -- [[R-2.7.2の改良点:http://cran.md.tsukuba.ac.jp/src/base/NEWS]]に spirinf 関数は結果をUTF-8で返すことができるようになったとある. -- Bayesの計算に合間に,机周りなど掃除. - 午後. -- K氏N氏来訪.スピーカーBronze Br-1,2 やヘッドフォンSENNHEISER HD595,アンプCR-D2の話など. * 2008年 9 月 16 日 (火) 雨 -> 晴 [#f12f2a39] - 06:10 起床 -- 07:30 研究室着 -- じめじめしてかなわん. - 午前 -- 昨日までに校正した二つの原稿の梱包を終えた.午前中に送ろう. -- 二つ校正終わったので,また別の仕事に入る. - 午後.郵便局に行き,二つ校正原稿を送付して,帰りに蓬来でチャーハン小520と餃子350. -- Software for Data Analysis: Programming with R が届いているので,眺める.これはS社から翻訳が出ることが決まっている. -- R に[[tcl/tkオプションを付けてコンパイル,インストール>R_Tcl/TK]]したけど,今のところ,tcl/tkインターフェイスが起動するのは install.packagesを実行する時ぐらい.サイトを選択するのにわざわざ新規ウィンドウを出してもらってもなぁ. -- [[長さが100のベクトル生成]]をR備忘録に. * 2008年 9 月 15 日 (月) 雨 [#w4a1391e] - 06:30 起床 -- 07:50 研究室着 - 午前 -- S社第2校をもう一度確認.修正しても(そのままでも)良い点をいくつか見付けた. - 午後.麺. -- M 社初校の最終確認作業に入る. -- 17:15.一通り校正終了.えらい疲れた... -- 初校を返送する用意もしたし,明日午前中に送ろう. * 2008年 9 月 14 日 (日) 晴 [#u2375f11] - 06:30 起床 -- 朝方,自衛隊のヘリ2機が,西の方へ飛んでいった. -- 長いこと愛でていたサボテンの一つがへたってしまった.手術して,様子を見る. - S社の第二稿の確認を終わらせる. * 2008年 9 月 13 日 (土) 曇 [#f66d65ea] - 06:30 起床 -- 07:50 研究室着 --充電 - 午前 -- 「統計の基礎」の執筆 --[[xaxt_yaxt_パラメータと正規分布の図]]をRの備忘録に追記 --[[摂氏の表し方]]をTeXの備忘録に追記 -午後.麺 -- ひたすら校正作業 --[[RjpWiki:index.php]]に,[[Rの基礎とプログラミング技法:http://www.springer.jp/japan/math/j71218.html]]は原書より売れているという記事があった.そうなんですか. * 2008年 9 月 12 日 (金) 晴 [#od2c55c4] - 06:10 起床 -- 07:20 研究室着 - 午前 -- 今日も,ひたすらM社の初稿校正.S社の方は,明日集中的に行おう. -- 09:50 とりあえず,ざっと目を通した.今からもう一度読み直し,さらに『統計の基礎』の章の執筆に入る. - 午後.麺 -- うーーん,疲れた... -- 久しぶり晩酌 * 2008年 9 月 11 日 (木) 晴 [#k4137779] - 06:10 起床 -- 07:20 研究室着 - 午前 -- 初稿の校正. -- 電話があって,会議を開催中だという.いってみると出席者が全然足らない状況.それから40分ほど待って,ようやく定数に達して会議.会議そのものは10分程度で終了. - 午後.麺 -- 校正.うーん,編集部では細かいところまで,かなり手を入れてくれてるんだけど,その修正が印刷所の作業にはうまく反映されていないようで,もとの文意から微妙にずれてしまっているところが多い. --- ひたすら校正.さすがに疲弊してきた... -- 帰り際BookOffによって,文庫本を仕入れる. -- 帰宅すると,S 社からの第2稿が届いていた. * 2008年 9 月 10 日 (水) 曇 [#wfa2f880] - 07:10 起床 -- 09:10 研究室着 - 午前 -- 初稿の校正に入る.かなり時間がかかりそう. - 午後.麺 -- 「統計の基礎」を付録として追加することにした. -- factor の日本語にいつも迷う.これまで「要因」としてきたが,先月校正の終わった本では「因子」とした.でも,いま校正している本では,やっぱり「要因」のママにしておこうか.我ながら一貫性がないなぁ. -- [[RMeCab]]にかなり関数を追加しているので,それにあわせた追記が多くなる. -- 先週の学会出張の際,移動の間,気楽なもの読もうと思って,宮部みゆきや山本周五郎の文庫をBook Off でまとめ買いしたのだが,はまってしまって,帰宅すると,食事や入浴しながらも読み進めてしまう... * 2008年 9 月 9 日 (火) 晴 [#a3ea3d0f] - 07:10 起床 -- 自宅で,他の学会の原稿書き * 2008年 9 月 8 日 (月) 晴 [#wc26c248] - 07:10 起床 -- お昼ごろ,いったん研究室に赴く. -- M 社の初校が届いていた. * 2008年 9 月 7 日 (日) 晴 [#ia2a7cef] - 07:10 起床 -- 一日静養 * 2008年 9 月 6 日 (土) 晴 [#n854362d] - 07:10 起床 -- 夕方の飛行機なので,とりあえず新宿に出る.新宿は,自分のテリトリーだと思っていたのだが,目的の建物が見付からず(自分の子供の頃は,そんな建物は無かった),角の交番で尋ねる.この交番はちっとも変わっていない.道を教わって三越方面に戻ると,恐ろしく細いビルがあって,その上階が目的の店であった. -- 羽田の6階ラウンジで軽食を取って,搭乗。空港での着陸が、海側からではなく、陸側からの進入であった。珍しいが、何の都合かな?18時過ぎに自宅着. * 2008年 9 月 5 日 (金) 晴 [#oc005dd6] - 07:10 起床 -- 行動計量学会で発表するため,成蹊大学へ.吉祥寺は,本当に久しぶりである.中学高校の頃はしょっちゅう遊びに来ていた.3本立ての名画劇場が複数あったからだが,もう,取り壊されているんだろうな.駅前に昔の面影はない.周辺を散歩してみたかったが時間の余裕がない. -- バスで成蹊大学に向かう. -- D大のJ先生に久しぶりお会いする.今日,出てこられて,夜には戻られるそうだ. -- 発表を終える. -- 他の方々の発表で,認知意味論や生成文法の仮定を計量的に実証するという内容の話があった.個人的には,認知意味論や生成文法を捨てて確率言語論に入っているので,計量的に実証しようとしたら,逆に反証されてしまった,なんて心配はないのかなぁ,などと思ったりもした. -- 夜はまた中華料理. * 2008年 9 月 4 日 (木) 晴 [#fe913278] - 07:10 起床 -- 12時の飛行機で東京へ.そのまま目黒の白金台に宿泊。ベランダでビールでいっぱい。 -- 夜は中華料理を取り寄せたが、食べ過ぎた。 * 2008年 9 月 3 日 (水) 曇 [#d5034c91] - 07:10 起床 -- 09:10 研究室着 - 午前 -- 週末から行動計量学会だが,別の学会の論文要旨の作成に入る. - 午後.麺 -- とりあえず提出用の論文要旨の作成は完了した -- [[RMeCab]]のソースをだいぶ整理したので,WindowsとMacのそれぞれで一通り処理を行ってみた.直感的には,前よりも処理時間が短縮されたように思うのだが. -- S出版に本の要約を送付 -- 明日から東京出張なので,細々とした準備に入る. * 2008年 9 月 2 日 (火) 晴 [#gaac09aa] - 06:10 起床 - 午前 -- 自宅で原稿作成 -- 11:30 研究室着 - 午後.麺 -- [[RMeCab]] と [[RCaBoCha]]の整備 * 2008年 9 月 1 日 (月) 晴 [#rd10c4d2] - 06:10 起床 -- 自宅でスライド作成 -- ビールなど飲んで寝る,..