データを分析し、より高い値が同じ間隔で繰り返されるインスタンスを見つける必要があります。例:
スペクトル関数を使用しようとしていますが、奇妙な結果が得られます。
これが私のサンプルデータです。1時間ごとに外れ値を挿入しました。
library(dplyr)
library(lubridate)
library(ggplot2)
set.seed(900)
data1 <-
data.frame(
datetime = seq.POSIXt(as.POSIXct("2020-12-26 10:00:00"), as.POSIXct("2020-12-26 10:00:00") + 15*50001, "15 sec"),
Value = sample(1:10, 50002, replace = T),
Instance = "A"
)
data1.1 <- data.frame(
datetime= seq.POSIXt(as.POSIXct("2020-12-26 10:00:00"), as.POSIXct("2020-12-26 10:00:00") + 15*50001, "hour"),
Value = sample(10:100, 209, replace = T),
Instance = "A"
)
data1 <- rbind(data1, data1.1) %>% group_by(datetime, Instance) %>% summarise(Value = max(Value)) %>% ungroup()
ggplot(data1, aes(x=datetime, y=Value, color = Instance)) +
geom_point()
spect <- spectrum(data1$Value, log="no", spans=c(5,5), plot=FALSE)
delta <- 1/4
specx <- spect$freq/delta
specy <- 2*spect$spec
plot(specx, specy, xlab="Period (minutes)", ylab="Spectral Density", type="l")
スペクトル密度が 60 分を示すスペクトル グラフが表示されることを期待していました。しかし、それは私が得たものです:
間隔を繰り返すピークを見つける方法 (私の例では 60 分)?