一連の連続した時点での日射強度グラフの下の面積を見つけようとしています。基本的に、過去 24 時間の 7 日間の 1 時間ごとの日射量の積分 (過去 24 時間の移動合計) が必要です (土壌の温度は、過去 24 時間の日射量の結果であると思われます)。
ここに私が使用しているコードがありますが、機能しますが、統合の時間枠を簡単に変更できるように自動化したいと思います (12、18、24、36 時間を試してください)。関係があるかどうかを確認するために、時間ごとの温度データに対してプロットする可能性がある値)
ここで: Rg - 10 分間の測定での日射量 num - データフレームのエントリ番号 AUC_xxx - 過去 24 時間の総日射量
y<-as.numeric(xx$Rg[xx$num["2015-09-13 14:10"]:xx$num["2015-09-14 14:00"]])
x<-c(1:length(y))
id <- order(x)
AUC_s14_14 <- sum(diff(x[id])*rollmean(y[id],2))
私はrollapplyで試しましたが、私は再び立ち往生しています:
rollapply(xx$Rg[xx$num["2015-09-13 00:10"]:xx$num["2015-09-14 00:00"]], width = 144, by = 6, FUN = **"INTEGRAL"**, na.rm = TRUE, align = "left")
助けてくれてありがとう ！