【求助】用R进行时间序列谱分析如何显著性检验

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

对一个千年的时间序列进行谱分析，假设为timeseries[1:1000]
现在通过R语言中的spec.pgram包进行谱分析，参数不多,具体如下图：

在R中我使用
ap=spec.pgram(timeseries[1:1000],detrend=T,plot=T)
进行绘图。

但是显著性检验的那条曲线如何添进去？这个软件包有提供自带的显著性检验吗？求有经验的大神指教。
下图是我在读文献时找到的，我也希望有这样的效果：


该图中黑线没话说就是频谱分析，红蓝线分别是95%置信度和99%置信度，我想把这两条线加上，请问如何分析。

提前谢过~

附加问题：


上图中横轴处理的非常漂亮，重点强带0.001-0.1区间的信息，请问这种该如何画？？？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏5 回帖

关键词：时间序列 谱分析 timeseries Detrend Series 如何

谱分析全用NCL实现了，坐标轴和显著性检验。放弃用R做谱分析了。
NCL谱分析的官网教程：
http://www.ncl.ucar.edu/Applications/spec.shtml

再加两个plot就可以了
你看这个例子
http://www.inside-r.org/r-doc/stats/spec.pgram