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

 马上下载



1.png
2.png
3.png
4.png
二维码

扫码加我 拉你入群

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

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

全部回复
2014-6-23 15:20:29
自己试着给了一些参数,但有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数据给的不对?
二维码

扫码加我 拉你入群

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

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

2014-6-23 15:52:30
自己看手册,发现需要把时间序列通过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)
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群