ぎょーむ日誌 2005-01-11

2005 年 01 月 11 日 (火)

• 0700 起床． 朝飯． コーヒー． 0820 自宅発． 晴． 0830 研究室着．
• 0905 研究室発． 正門横の学術交流会館で 0920 よりセンター試験監督説明会． 1025 終了． 今回はキレてしまうせんせーもなく比較的しづかに終った． 1030 研究室もどる．
• 一昨日来日した滞在研究者・操国光 (Cao GuoXing) さんの PC 設定など． A703 にあった PC を A801 に移す． 前にこの機械をつかってた岳明さんがかなり中国語な設定にしてくれるんで， かなりラク． しかし A801 の電源状況は心配．
• 本日のぎょーむ進捗． 3 人の (研究室外) M2 のだす難問にふりまわされて， あまり進捗しなかった， というか ……
• 一人目． 霜月さん が書いてる「見えない集団」問題． この Slatkin 論文はともかく， けっきょく Beerli の提案する集団間相関を coalescent process で推定する， という問題に帰着していくのである． これは MCMC お勉強モードのときにダウンロードしたけど， 未読だったな …… といま見かえすと， さすがに以前よりは理解できる． しかしここでまた没頭するのは危険なのでいったんおく．
• 二人目． 苫小牧 の平松さん． R の glm() とか線形モデルのたぐいで， 順序因子 (ordered なやつら) 使ったときに推定値として得られるアレはいったい …… というもの．

  	Coefficients:
  		     Estimate Std. Error t value Pr(>|t|)
  	(Intercept)     0.702      0.143    4.92  3.5e-05
  	         x.L    1.384      0.402    3.45   0.0018
  	         x.Q    0.919      0.349    2.63   0.0136
  	         x.C    0.440      0.286    1.54   0.1348
  	         x^4    0.345      0.206    1.68   0.1050
  

  ああコレは簡単で， k 水準の順序な変数 x を整数化しちゃって k 次の多項式にほうりこんでいて， その多項式の係数が …… と具体例をあげて説明しようとして 自分の理解が間違っていることに気づかされた． うーむ， 院生たちにはいつも教えられる．
• で， かとーさんもまきこんで調べてみると， S や R の中で順序因子は多項式は多項式でも 選点直交多項式系 というのに変換され， しかも私のニガ手な対比行列にからんでくるとわかった． この計算じたいは R の poly() 関数でできるんだが orthogonal polynomial がわかってないんで ススまない． えーい， これも宿題だ．
• 蛇足ながら， 単に作図したいだけならかかる直交多項式は理解する必要なく， 推定関数 results <- glm(...) などから得られる期待値 results$fitted をてきとうに選んで使えばよい …… というようなごまかし R プログラムつきお返事メイルを送信．
• 三人目． 苫小牧の平尾君で， 世の中では simultaneous (SAR) and conditional (CAR) autoregressions のうち， SAR 推定値がとくに注釈もなく使われていて， R の package (spdep とか) なんかではそのまま尤度計算してるんぢゃないでしょーか， というご指摘． つまり pseudo likelihood かどうか， なんて誰も気にしてない， というものだが．
• 私の理解ではこのあたりの区別は重要で， X という事象が観測される尤度が 「観測された X そのもの」 に依存するならそれは pseudo likelihood だ． たぶん定義としてはこれでまちがっていないだろう …… このあたりの故事来歴に関しては SAR にしろ CAR にしろ pseudo likelihood にしろいづれもことごとく Besag の 30 年前ぐらいの研究から来ているものなんで， このあたりのもと論文にあたっていくほかあるまい．
• それはいったんおくとして， 応用的な問題からすると maximum pseudo likelihood estimate (MPLE) が現象の説明にどれぐらい使えるだろうか …… というあたりなんだよねえ． こんにちでは MCMC 法などつかえば MPLE にもとづいた乱数も (なかばむりやり) 得られちゃったりしたりするわけで． そして MPLE については「やっぱりへん」 というのと「けっこう『良い』場合もある」 の両方が示されていて， ですね．
• えーい …… ともかく， これも宿題だ．
• さらに， 甲山さんからぱいぷ樹木原稿かえってきて， ですね． すけじゅーるの観点からするとこれが一番重要， となるわけだが．
• さらに， 操さんの中国語 IME インストールしてくれ問題にかなり苦闘して， ですね．
• ばてたので途中撤退． 1850 研究室発． 1910 帰宅． 晩飯くってさっさと寝てしまった． えー， 一種の時差ぼけです． うう．
• [今日の運動]
  • うんどう休養日．
• [今日の食卓]
  • 朝 (0730): シリアル．
  • 昼 (1410): 研究室お茶部屋． 米麦 0.7 合． タマネギ・ニンジン・ジャガイモ・豆腐のカレー．
  • 晩 (2000): 米麦 0.7 合． タマネギ・ニンジン・ジャガイモ・豆腐のカレー．