私が，Savitzky-Golayフィルタと出会ってから，もうすぐ１０年になります．私は，Savitzky-Golayフィルタをとても気に入っていましたが，残念ながら，さよならするときが来ました． 最近，Savitzky-Golayフィルタは理想のフィルタではないことを認識しました…

下の図のように，時系列解析では，パワースペクトルとか，コヒーレンスとか，DFAなどのスケーリング解析とかで，いくつかの結果のグラフをまとめて，平均的なグラフを求めたいことがあります．今回は，この処理をするRスクリプトを紹介します． x軸の各点で…