ぎょーむ日誌 2005-07-21
2005 年 07 月 21 日 (木)
-
0700 起床．
朝飯．
コーヒー．
洗濯．
0925 自宅発．
晴．
0940 研究室着．
-
朝から
R
での疎行列計算について検討している．
これは
Matrix
使うのが順当なんだろうな
……
``Introduction to the Matrix package''
読んでみる．
ふーむ，
class は
Matrix()
したときに勝手に指定されんのかしらん?
rcond()
はたんに
LU 分解
するだけみたいだし．
-
とはいえ，
私が興味があるのは行列とヴェクトルの積なんで LU 分解は関係ない．
そこで，
280 × 280
の行列の積を 100 回くりかえし計算させる実験やってみる．
ふつーの
matrix だと 2.7 秒ぐらいかかる計算が，
Matrix だと 2.2 秒ぐらいで片づくようだ．
行列とヴェクトルの積 300 回計算の実験やってみると，
matrix で 1.02 秒ぐらい，
Matrix だと 0.48 秒ぐらい，
という結果になった．
まあ，
ためしに使ってみるかな．
-
混合モデルの計算，
なんとか Newton-Raphson 法でケリがつかないものか，
とぢたばたしてみる．
だめっぽい．
ここで期待対数尤度の最大化，
と問題をすりかえてしまうというテもなくはないんだが
……
いや，
個体ごとに個別の random effects がつくわけだから，
これまたダメかもしれんな．
-
2005年度 統計関連学会連合大会の
プログラム
できたよーだ．
ひやかしに行きたいところなんだけど，
それこそひやかしに行くだけなんですけぢゅーる空けてる状態．
参加申し込みは 8/26 までか．
-
昼飯．
-
アカマツ原稿 (シミュレイションのほう)
残りがぜんぶもどってくる．
ということで，
文章修正作業．
しかしながら，
アタマがなかなか切りかわらない．
-
ぢりぢり進捗させていく．
-
抑圧的な気分で作業してると，
どこかでそれが発散してしまうのかもしれない．
D 論公聴会を翌日にひかえた
伊藤さん
「よりぬきカツラくん」
コメント欄で quasipoisson 議論 (ひょっとしてこれって「荒し」?!)．
よく見ると
Venables & Ripley (2002)
からの引用 (p.210) の久保訳が
あまり正しくないような気がする
……
ともあれ quasi-likelihood については明日のぎょーむ日誌にでも．
-
原稿修正，
全部を終える前に気力がつきてしまった
……
2130 研究室発．
2155 帰宅．
晩飯．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0820):
パン．
ニンジン・タマネギ・セロリ・豆腐のカレー．
- 昼 (1350):
研究室お茶部屋．
パン．
ニンジン・タマネギ・セロリ・豆腐のカレー．
- 晩 (2230):
米麦 0.8 合．
モヤシ・マイタケ・ニンニク・イカの炒めもの．
乾燥野菜海藻スープ．