2012年11月11日日曜日

tapplyとby

今年最後の軟禁。午後からネットが使用できない、ということなので、午前中に更新。

ヒノキ林データ整理の続き。

各調査地点において3地点ずつ全天写真を撮影し、樹冠開空度を測定している。


データセットはこんな感じ。nameが調査地点名、openが開空度の数値データ。


間伐直後とか、広葉樹とかあるので、当然、開空度に違いがある。箱ヒゲ図にするとこんな感じ。


で、各調査地点ごとの平均値も欲しいなと。キット良い方法があるだろうと調べたら、最適な動画を発見。

tapply(open,name,mean)、とすればOK。meanは平均値の意味で、これを変更するだけで、色々と求められる。by(open,name,mean)でも良いみたい。