R中实现谱分析的相关参数设置问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

在R中有一个叫做multitaper的程序包里面有spec.mtm，不明觉厉地感到它的功能很强大。
但鉴于自己统计基础比较薄弱，谱分析中的各个参数并不知道是什么意思，所以贴出来让有经验的大神们指导一二。
我用到的数据是只有1列的时间序列，长度为1000行，代表着1000年，我想基于spec.mtm，直接搞定谱分析，跪求参数的一般设置。
multitaper.pdf (117.68 KB)

2014-6-23 14:47:58 上传

multitaper程序包说明

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：参数设置 谱分析 是什么意思 multi tape 程序 统计

自己试着给了一些参数，但有warning:
Warning message:
In spec.mtm(Ypyear1, nw = 4, k = 8, dtUnits = "year", sineSmoothFact = 0.2,  :
  Time series is not a ts object and deltat is not set. Frequency array and axes may be incorrect.
这是不是意味着我的time series数据给的不对？

自己看手册，发现需要把时间序列通过ts指令变成R认识的时间序列（我是1000年的数据）：
Ypyear1=array(scan('Ypyear_1.dat'),c(1000))
Yp1.ts<-ts(data=Ypyear1,start=c(850,1),deltat=1)
spec.mtm(Yp1.ts,nw=4.0,k=8,dtUnits="year",sineSmoothFact=0.2,jkCIProb=.95,deltat=1)