Rでフーリエ解析

vsolj 10750 で紹介されているコードを走らせてみました。
後半の

　(sum(d$V2*ss))^2 + (sum(d$V2*cc))^2

の部分は普通にコピぺしてくるとエラーが出るので水色文字のようにしました。

テストで使ったのは VSOLJ の T Mon のデータ。
　
　MONT 200912212258 68 Onr
　　　　　　・
　　　　　　・
　　　　　　・
　MONT 201104252051 58 Heo

の範囲を使わせて頂きました。


グラフはこんな感じです。

取り出せた周期は以下の通り↓

　pgm$period[which.max(pgm$power)]
　[1] 26.96

GCVS (27.024649 d) よりちょっと短いですが (^^;