ぎょーむ日誌 2009-08-17

2009 年 08 月 17 日 (月)

• 0815 起床． ねむい． 0825 自宅発． 曇． 0845 研究室着． 朝飯． コーヒー．
• モニターの色温度を下げるとよい， というハナシを読んだので 6500 にしてみる． もと記事が見つからなくなったので， 検索 してみた ……
• 岩波統計モデリング本がらみで， R の summary(glm()) をちょっと調べているんだが …… 以前 にこのあたr で何かややこしーことやってるんでは， とか書いてたわけだが …… うーむ， 当方のよくわからぬかんちがい (もしくは実行環境のちがい， R コンソール上で実行したか R --vanilla < script.R で実行したか) による混乱だった． 単純に options(show.signif.stars = TRUE) が設定されいれば

  Coefficients:
              Estimate Std. Error z value Pr(>|z|)    
  (Intercept)  1.26311    0.36963   3.417 0.000633 ***
  x            0.08007    0.03704   2.162 0.030620 *  
  fT          -0.03200    0.07438  -0.430 0.667035    
  ---
  Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
  

  となり， options(show.signif.stars = FALSE) ならこうなる， というだけのハナシだった (有効ケタ数が異なるのは下では options(digits = 5) としているため)．

  Coefficients:
              Estimate Std. Error z value Pr(>|z|)
  (Intercept)   1.2631     0.3696    3.42  0.00063
  x             0.0801     0.0370    2.16  0.03062
  fT           -0.0320     0.0744   -0.43  0.66703
  

• このあたりの問題， ネット上で検索してみるとおもしろかったので， summary(glm()) の星 メモをつくってみた． R-help 2000 年の一連の議論の Ripley 御大の発言， まったく同感です．
• えーい， おもしろいけど時間とられる ……
• へー， RPD (Restart Postdoctoral Fellowship) とゆー制度があるんだ． Restart とかいう名前はヘンだけど， 良い制度ですね．
• R 作図わざ本， 共立出版の制作な編集のかたとやりとりしつつ， 最後の作業． 「表記のゆれ」 統一作業， 私が作った日本語版さくいんの自動生成システムの説明， などなど．
• …… てなかんぢで本日も終了してしまいましたよ． 八甲田湿原のヤブに突入できなかった． 1840 研究室発北大構内走． 1910 研究室もどる． 1920 研究室発． R 作図わざ本の校正刷を紙袋につめて， クロネコヤマトの着払い荷札をつけて …… 北大 13 条門前のセヴンイレヴンで， この重量 1.8 kg の紙束をぽいっと送信してしまう． これで一件落着， であればよいのだが …… 1940 帰宅． 晩飯．
• [今日の運動]
  • 北大構内走 30 分間． ストレッチ．
• [今日の食卓]
  • 朝 (0900): 研究室． バナナ．
  • 昼 (1300): 研究室． バナナ． 食パン．
  • 晩 (2030): なんとなくごはんやデンスケでサバを食べてしまった．