2014年7月5日

曇り時々雨　梅雨らしい天気ですが、気温がかなり低い

これは、今年の米の収穫量が落ちるかもしれないですね・・・。雨も全体的に少ないし。

朝、少しだけ子どもとテレビ電話で話をしてから、職場へ。明日から出張のため、片付けられるところは少しでも片付けようと。

昨日からしたことは、gamを使って（package " mgcv "）、生物種毎のハビタットを推定するための計算。これは、まあまあ妥当と思われる結果が出てきました。

　　　model <- gam(y~s(x), family=binomial, data=data)
のような形。ここでは、二項分布を仮定。参考は、Zuur, A.F. A begginer's guide to generalized additive models with R.

データをbinomialにするところも、Rでコード化（大したことはない）しているので、データが追加されても、嫌がらずさくさく動くはず。他の膨大なデータにも応用可能なので、これからあちこちで使いまわす予定。

秋くらいに投稿を考えているので、その準備をしておきます。

で、16時過ぎまで論文の元原稿をあちこち修正しながら、少しずつ形を整えていっています。書いていて、分からないところは本で調べたりして・・・。まあ、論文を書くのなら、当たり前の作業なのですが、ここらへんが結構大変なのです。

朝：トースト
昼：幕の内弁当
夜：昨日と同じ、豚肉とピーマンの炒めもの

勉強：土曜ゼミ