2014年9月14日

晴れ

三連休の2日目。朝、少しだけのんびりと過ごした。8時過ぎに職場へ。下の階に職場の知り合いが引っ越ししてきた。

今日も0データを扱うモデルについてコードを書き加える。とりあえず、同僚も違うデータセットで試してみて、うまく進めているようだ。0と1を最初にbionominalモデルで切り分けて、各説明変数に対して予測値を返し、年効果を計算してlogit変換する。その次にGLMで1以上の部分について年効果を出し、それぞれを掛け合わせるというもの。それと、zero-inflated modelについて・・・。試してみたら、普通のnegative binomialも01切り分けモデルもzero-inflatedモデルもさほど結論は変わらない。これは、私と同僚がそれぞれ別のデータセットで行った場合でも変わらなかった。

同じデータセットで、別の種で確かめるために、あれこれコードを書き換えていたら・・・・思わずセーブして、書き換えのところで保存されてしまった。仕方ない、明日復旧させよう。別の種だと0を考慮したらかなり印象が異なってしまった。

データセットを少し見直す。生データをそのまま使うと、時間がかかり、しかもほとんどが0なので、やはり収束が悪い。したがって、環境条件を少し絞って、生データをあまり殺さずに、作り変えることに。

その結果、2種のうちの1種はほぼこれまで通りだが、もう1種もなんとなく見えてきた。単純にしただけのモデルでは見えなかったものが見えてきた感じ。明日も少し頑張ってみよう。きっと、あと少しで真実が見えてくるはず。

朝：トースト
昼：カレーライス
夜：カツオの刺し身

勉強：中学生は日曜ゼミ、小学生もちょっと休ませる。明日からまた再開。

Sep-6 2014-32