Top

個別メモ

Latest update on 2012年3月5日 (月) at 10:54:46.

【第725回】 一日遅れ（2007年2月1日）

• そういうわけで5:00から作業を継続している。ロジスティック回帰が収束しないので，悩んだ挙句に線型判別分析に逃げることにした。まあこれでも目的は達成できる。図示を多用したので，読みやすいものにはなったであろう。しかし，この報告書を真剣に読んでくれる人が何人くらいいるんだろうかと考えると，ちょっと空しかったりもする。
• 朝食もとらずに続けて9:30脱稿。さっそくメール添付で委員長に査読を依頼し，弁当を買いに行く。次は地域保健実習の一覧表を仕上げなくては。
• ……と思ったが，さっきメールで送ったpdfファイル中，図として挿入したメタファイルの半角英数字が文字化けしてしまっていることに（日本語は大丈夫だったのに）気づいて作り直し。OOo-2.1のバグか？　もしかしたら設定で回避できるのかもしれないが，試行錯誤するのが面倒なのでPDF reDirect v2に出力して済ませた。文字化けも起こらずイメージも崩れず，しかもOOoで作ったpdfよりサイズも小さくなった。念のためWord形式でも保存してみたが，案の定，ページ当たりの行数が変わってしまったり，表の体裁が崩れたりした。表の体裁が崩れた原因は，OOoではほとんど幅をとらなかった余分な半角スペースをWordはちゃんと幅をとってしまうためだったので，たぶんぼくが悪いんだが。表の元をRで出力させるのにsprintf("%5d",...)とかやったのがいけなくて，たんにsprintf("%d",...)なら余分なスペースが入ることはなかっただろう。
• 今回の仕事からメモしておくと，cor()をデータフレームに使うと相関係数行列ができるが，母相関係数がゼロであるという帰無仮説を検定するためのcor.test()は2つのベクトルしか引数にできないという点は，例えばデータフレームZを分析する場合なら，

  rZ <-cor(Z)

  と

  NZ <- nrow(Z)

  をしてから，

  tZ <- abs(rZ)*sqrt((NZ-2)/(1-rZ*rZ))

  でt値の行列を得て，最後に

  2*(1-pt(tZ,NZ-2))

  とすれば有意確率の行列が得られる。気づいてみれば当たり前だが，今までやったことがなかった。
• さて一覧表作業を始めよう。
• このサーバのCPUとマザーボードとメモリは去年の10月25日に買ったのだが，そのPC SUCCESSという店が昨日閉店したらしいのは残念。あれだけ安売り合戦をしていては小売業は厳しかろう。安い通販店で何でも買ってしまう自分がそんなこと言える立場じゃないが。
• 15:00を過ぎ，やや小腹が空いてきた。けれども，いま食べてしまうと眠くなりそうだし，まだ何とか耐えられそうなので耐えよう。
• 笹山掲示板の[5646]のChaosさんの書き込みで知ったが，厚労省が「新型インフルエンザ対策ガイドライン（フェーズ４以降）」のパブコメ募集をしている。コメントするかどうかは別として，公衆衛生部門の資料は読んでおかないとまずいな。
• 夕方，約１時間半の会議。その前に結局空腹に耐えかねてパンを買ってきて食べてしまった。一覧表はあと少し。
• 一覧表が完成して打ち出してチェックをしていたら21:00。今日は終電の1本前には帰りたい。
• というわけで，21:20に研究室を出て，復路あさま551+553号。

▼前【724】（いろいろ締め切り（2007年1月31日） ）　▲次【726】（積み残し（2007年2月2日） ）　●Top

△Read/Write COMMENTS

Notice to cite or link here | [TOP PAGE]