对一个千年的时间序列进行谱分析,假设为timeseries[1:1000]
现在通过R语言中的spec.pgram包进行谱分析,参数不多,具体如下图:
在R中我使用
ap=spec.pgram(timeseries[1:1000],detrend=T,plot=T)
进行绘图。
但是显著性检验的那条曲线如何添进去?这个软件包有提供自带的显著性检验吗?求有经验的大神指教。
下图是我在读文献时找到的,我也希望有这样的效果:
该图中黑线没话说就是频谱分析,红蓝线分别是95%置信度和99%置信度,我想把这两条线加上,请问如何分析。
提前谢过~
附加问题:
上图中横轴处理的非常漂亮,重点强带0.001-0.1区间的信息,请问这种该如何画???