vsolj 10750 で紹介されているコードを走らせてみました。
後半の
(sum(d$V2*ss))^2 + (sum(d$V2*cc))^2
の部分は普通にコピぺしてくるとエラーが出るので水色文字のようにしました。
テストで使ったのは VSOLJ の T Mon のデータ。
MONT 200912212258 68 Onr
・
・
・
MONT 201104252051 58 Heo
の範囲を使わせて頂きました。
![](https://blogimg.goo.ne.jp/user_image/29/34/9fd729be901a6101f48b95dc4c165073.png)
グラフはこんな感じです。
取り出せた周期は以下の通り↓
pgm$period[which.max(pgm$power)]
[1] 26.96
GCVS (27.024649 d) よりちょっと短いですが (^^;
後半の
(sum(d$V2*ss))^2 + (sum(d$V2*cc))^2
の部分は普通にコピぺしてくるとエラーが出るので水色文字のようにしました。
テストで使ったのは VSOLJ の T Mon のデータ。
MONT 200912212258 68 Onr
・
・
・
MONT 201104252051 58 Heo
の範囲を使わせて頂きました。
![](https://blogimg.goo.ne.jp/user_image/29/34/9fd729be901a6101f48b95dc4c165073.png)
グラフはこんな感じです。
取り出せた周期は以下の通り↓
pgm$period[which.max(pgm$power)]
[1] 26.96
GCVS (27.024649 d) よりちょっと短いですが (^^;