R 函数 spec.pgram() 的均匀间隔 1/freq

Question

目标是绘制两个时间序列之间的相干性（即相对于频率的相关系数）。如何让 x 轴上的 1/freq（即周期）均匀分布？

t <- 0:99
ts1 <- ts(2*cos((2*pi)/24*t))
ts2 <- ts(2*cos((2*pi)/48*t))
ts12 <- ts.intersect(ts1, ts2)

Coh <- spec.pgram(ts12, spans=3)
plot(Coh$freq, Coh$coh, type='l')
plot(1/Coh$freq, Coh$coh, type='l')  # how to get 1/freq to be evenly-space?

我试图修改函数 spec.pgram() 但没有成功。更具体地说，我替换了以下行：

freq <- seq.int(from = xfreq/N, by = xfreq/N, length.out = Nspec)

与：

freq.tmp <- seq.int(from = xfreq/N, by = xfreq/N, length.out = Nspec)
freq <- rev(1/seq(from=1/max(freq.tmp), to=1/min(freq.tmp), length.out=Nspec))  

还有其他人的运气更好吗？谢谢

Answer 1

您的意思是您只想用句点而不是频率重新标记 x 轴吗？这将保持 x 轴上值的间距，但代价是 x 值的非线性缩放。例如（使用ggplot2）：

library(ggplot2)

dat = as.data.frame(Coh[c("freq","spec","coh","phase")])

ggplot(dat, aes(freq, coh)) + 
  geom_line() +
  geom_point() +
  scale_x_continuous(breaks=dat$freq[seq(1,nrow(dat),3)], 
                     minor_breaks=dat$freq,
                     labels=round(1/dat$freq[seq(1,nrow(dat),3)],1)) +
  labs(x="Period")

您还可以将 x 值标签设置为整数周期：

breaks = c(1:10,15,25,50,100)

ggplot(dat, aes(freq, coh)) + 
  geom_line() +
  geom_point() +
  scale_x_continuous(breaks=1/breaks,
                     minor_breaks = 1/(breaks[-1] - 0.5 * diff(breaks)),
                     labels=breaks) +
  labs(x="Period")