・久しぶりに試験地に出勤。連絡調整だけで午前中は終わりとなってしまい、ううむ。ついにスギ花粉が飛び始めたらしい。が、ニュースを見るまでもなく、当方の鼻は敏感に反応しており、自分自身が検出キットになりそうである。既にマスクは装着しているのだが、これから2-3ヶ月は厳しいシーズンである。
・標高別造林試験地のデータ解析再び。標高と標高の2乗の両方を説明変数に入れると、二次関数での最適化ができる。stepAICで調べると、何種かは一次回帰の方が当てはまりがよいものの、ほとんどは二次関数の方がよくフィットしている。多少おかしな樹種もあるが、割合とよく推定できている。
・次に、樹高を測定した30個体のデータを用いて、DBHと樹高の関係式を求める。さらに、その推定式を使って全個体の樹高を推定し、最終的なvolumeを算定するためのコードを書く。相変わらず美しくないコードだが、何とかR上で樹種ごと・サイトごとの合計Volume(D2Hの総計)を計算できるようになった。なるほどこちらの方が、林学的には分かりやすい指標だ。
・GIS管理を行うため、試験地の区分けに3人でくるりと見回る。カラスが集団でカエデ林付近に集まっており、それなりにおどろおどろしい雰囲気をかもし出している。それにしても、せっかく樹林があるので、鳥の種類が分かるとずいぶんと楽しそうである。個人的に、水辺の鳥と山野の鳥という薄いハンドブックを購入したので、まずは初心者向きといわれる水辺の鳥を片手に遊水地辺りを歩いてみよう。
・標高別造林試験地のデータ解析再び。標高と標高の2乗の両方を説明変数に入れると、二次関数での最適化ができる。stepAICで調べると、何種かは一次回帰の方が当てはまりがよいものの、ほとんどは二次関数の方がよくフィットしている。多少おかしな樹種もあるが、割合とよく推定できている。
・次に、樹高を測定した30個体のデータを用いて、DBHと樹高の関係式を求める。さらに、その推定式を使って全個体の樹高を推定し、最終的なvolumeを算定するためのコードを書く。相変わらず美しくないコードだが、何とかR上で樹種ごと・サイトごとの合計Volume(D2Hの総計)を計算できるようになった。なるほどこちらの方が、林学的には分かりやすい指標だ。
・GIS管理を行うため、試験地の区分けに3人でくるりと見回る。カラスが集団でカエデ林付近に集まっており、それなりにおどろおどろしい雰囲気をかもし出している。それにしても、せっかく樹林があるので、鳥の種類が分かるとずいぶんと楽しそうである。個人的に、水辺の鳥と山野の鳥という薄いハンドブックを購入したので、まずは初心者向きといわれる水辺の鳥を片手に遊水地辺りを歩いてみよう。