Top

個別メモ

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

【第1421回】 疲れが残っている（2009年9月28日）

• 土曜は娘の中学の運動会に行って，帰ってきてからは地元の神社のお祭りで神輿を担いで練り歩き，何度か木遣りをやって力を使い果たして帰宅した。日曜は家事とか，午後は神輿などの片付けとか。本当は直会があったのだが，家の用事があるため失礼した。
• 月曜は5:45起床。具沢山の野菜スープを作って朝食。来月から長野市はゴミ出し制度が少し変わるからなのか，今朝は大量のゴミが出されていた。小雨交じりの天気だが自転車で長野駅に出て，往路あさま512号。年のせいか，一日経ってから神輿を担いでいた左肩に痛みが走るのが悲しい。
• K先生からボーディングパスの半券とパスポートのコピーが届き，ソロモン諸島調査がらみの他の書類仕事もだいたい終わった。来月の熱帯医学会出張関係の書類も終わった。これでデータ入力の続きができる（本当は本の原稿を書かねばならないという話もあるので，どちらを優先すべきか難しいところ）。
• 3週間ぶりのR実践活用勉強会だったが，競合リスクの解析のところで，舟尾さんのテキスト『[R流!]イメージで理解する統計処理入門』カットシステム，ISBN 978-4-87783-229-2（Amazon | bk1 | e-hon）に不備があった。競合リスクを考慮した解析の必要性の説明はわかりやすく素晴らしいのだが，本に掲載されているコードとデータでは，その次に書かれている結果は出てこないのだ。つまり，RのsurvivalライブラリでのSurv()関数で，右側打ち切りの場合の打ち切りフラグは0/1か1/2のどちらかしか許されていないので，現在の最新版では，2が最初に出現すると1/2で2がイベント発生と解釈されてしまうようなのだ（実は0/1/2の3値がある変数であっても）。2.6とかの古いバージョンだとエラーが出て実行できなくなってしまう場合もあり，読者が混乱するのは避けられないだろう。表8.3のデータを使って，図8.5と図8.6を描かせるには，

  MYDATA$CENSOR2<-ifelse(MYDATA$CENSOR==1,1,0)

  MYDATA$CENSOR3<-ifelse(MYDATA$CENSOR==2,1,0)

  として，打ち切りフラグを0/1に変換してやる必要がある。それでも図8.6は，喫煙者のみ表8.5で肺がん発生のはずだったのが心疾患発生に変わってしまっているので，本に掲載されている図とは違ってしまうが，論理的一貫性を得るためには図8.6の方を変えるべきと思う。その後のp.186のcmprskライブラリの使い方はわかりやすいが（競合リスクのCENSORのコードを2にしてCENSORを3値の変数にしてしまったのは，このライブラリのcuminc()関数とかcrr()関数を使うためだったと思うが，そのせいでsurvfit()関数を使うときとの整合性がなくなってしまった），欲をいえば，作図の完成度をもう少し上げてもいいんじゃないかと思うし，数学的な解説も割愛せずにつけて欲しかったところだ（ここを詳しく書いた日本語の本はほとんどないので）。
• 今回のソロモン諸島調査でいよいよ壊れてしまったHIDEO WAKAMATSUデザインの3WAYバッグをリプレースするために，今夜，アマゾンプライムでリュック型のPCザックが届くはずなのだが，まだ来ない。いつもだいたい19:00頃届くはずなんだが。交通事情でも悪いんだろうか。今日は，それが届いたら帰ろうと思っているのだが。
• 結局20:00頃届いて荷物を詰め込み，帰途に就いた。しかしバスから両毛線，新幹線の乗り継ぎが悪く，復路あさま553号になってしまった。荻原浩『四度目の氷河期』新潮文庫を読んだが，一風変わったビルドゥングスロマンであり，スポーツ小説であり，ラブストーリーでもあり，旅小説でもあり，不思議な読後感だった。学問語りも悪くないし。うん，いい小説だ。10年後くらいのエピローグが欲しかったような気もするし，また一方では，そういう蛇足はない方がいい終わり方であるような気もするが。

▼前【1420】（たぶん不在中に届いた仕事が溜まっているだろう（2009年9月25日） ）　▲次【1422】（今日も雨（2009年9月29日） ）　●Top

△Read/Write COMMENTS

Notice to cite or link here | [TOP PAGE]